<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" id="owaParaStyle">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body fpstyle="1" ocsi="0">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear Molpro-Team,<br>

<br>

I need a hint on a rather simple problem. Currently I am trying to get into MOLPROs capabilities of calculating anharmonic vibrational spectrum (e.g. VSCF) of isotopologues with the PROGRAMS: SURF, POLY and PESTRANS.

<br>

<br>

The specific example is to obtain the spectrum of deuterated water "HDO", by transforming the PES of a previous calculation of non-deuterated water "H<sub>2</sub>O".

<br>

<br>

I use the following input for the calculation:<br>

<br>

<font size="1">-----------------------------------------------------------------------</font><br>

<font size="1">***,h2o vibrational test<br>

memory,64,m<br>

gprint,orbitals,basis<br>

<br>

geometry={<br>

O<br>

H1,O,B1<br>

H2,O,B1,H1,A1<br>

}<br>

B1=  0.95907103 ANG<br>

A1=  103.51462331 DEGREE<br>

<br>

orient,mass<br>

mass,iso<br>

<br>

basis=vtz<br>

<br>

hf<br>

mp2<br>

optg<br>

<br>

label1<br>

{hf<br>

start,atden<br>

}<br>

mp2<br>

<br>

{surf,start1D=label1<br>

disk,where='scr',extern='h2o_1D2D3D_mp2_di_vrc.pot'}<br>

<br>

poly,dipole=1,vrc=1<br>

<br>

{pestrans,vrc=1,plot=3<br>

intensity,dipole=1<br>

mass,2,2.01410178d0}<br>

<br>

poly,dipole=1,vrc=0<br>

<br>

vscf,type=poly,dipole=1,combi=2,ubound=10000<br>

vci,type=poly,vam=3,dipole=1,combi=2,ubound=10000<br>

<br>

put,irspec,irspec.dat<br>

</font><br>

<font size="1">-----------------------------------------------------------------------</font><br>

<br>

This yields the correct results, but as you can see, i need to transform the initial PES three times! I think there should be a more sensible way of doing this?<br>

<br>

<br>

</div>

</body>

</html>