<div dir="ltr">Dear users,<br>I have a question about the theory behind the calculation of dipole moment for charged systems as implemented in molpro.<br>My system is CO++ for e.g. triplet sigma plus state the input is :<br>
<br> ***,CO++ dipole moment<br> memory,46,m<br> geometry={<br>  x,y<br>  a1,C<br>  a2,O,,,r(i)}<br> basis=av5z<br>  i=0; t=1.7<br> start:<br>   i=i+1; r(i)=t;<br>     int<br>     gexpec,dm<br>     rhf;occ,4,1,1;wf,12,1,0;<br>
     multi;occ,6,2,2;closed,0;<br>     config<br>     wf,12,1,2;lquant,0;<br>      maxiter,150,200;<br>    ci;wf,12,1,2;save,3024.2;                                                    !3sgpl<br>    e_ci2(i)=energy<br>    dipz(i)=dmz<br>
    dipy(i)=dmy<br>    dipx(i)=dmx<br>    t=t+0.2<br>    if (t.ge.10.0) t=t+1.00<br>   t0=r(i); t5=e_ci(i); t4=e_ci2(i);<br>   text,fi,$t0,$t5,$t4<br> if (t.le.20) goto,start:<br>   table r,dipz,dipy,dipx<br><br>So I did not put any constrains for dipole moment.<br>
Unfortunately I did not manage to find, which are the default<br>parameters for this calculation. So I would like to ask you if<br>someone could answer me these questions:<br><br>1- Do the dipole moments include both the nuclear and electronic<br>
contributions?<br><br>2- Have the calculations of d.m. been performed with respect to the nuclear center<br>of mass or the geometric center?<br><br>Or send me a reference where these are mentioned?<br>I supposed that  dipole moments include both the nuclear and electronic<br>
contributions and it is performed with respect to the nuclear center<br>of mass. But I would like to see it somewhere written.<br>I will be appreciate to your help.<br>Thank you very much,<br><br><br>Tereza Uhlikova (born Sedivcova)<br>
Physics dept. and INFM - University of Milano<br>- Via Celoria 16 - 20133 Milano (Italy)<br>  TEL.+39 02 50317379<br>and European Theoretical Spectroscopy Facility (ETSF) <a href="http://www.etsf.eu">http://www.etsf.eu</a><br>
<br><br><br></div>