[molpro-user] EOM-CCSD second moments
Tatiana Korona
tania at tiger.chem.uw.edu.pl
Thu May 17 14:08:31 BST 2007
Dear Joseph,
This is not the problem of EOM, but a general issue with definitions of
moments calculated by Molpro and Dalton, for the conventions used in
Molpro see http://www.molpro.net/molpro-user/archive/all/msg01427.html.
Since you see differences for the ground state, you can test on, e.g. SCF
moments as well.
Best wishes,
Tatiana
Dr. Tatiana Korona http://tiger.chem.uw.edu.pl/staff/tania/index.html
Quantum Chemistry Laboratory
University of Warsaw
Pasteura 1, PL-02-093 Warsaw, POLAND
`The man who makes no mistakes does not usually make anything.'
Edward John Phelps (1822-1900)
On Fri, 18 May 2007, Joseph Lane wrote:
> I am having some difficulty understanding the one-electron properties
> calculated with EOM-CCSD in molpro as compared to the CCSD linear response in
> Dalton. I realise that the EOM-CCSD properties are calculated as expectation
> and the linear response as energy derivatives, hence there should be small
> numeric differences between the two.
>
> For ground state calculations, the dipole moments and quadrupole moments
> calculated with both methods are the same, however the second moments are
> different.
>
> Molpro:
>
> !CCSD_nonrel dipole moment : 0.00000000 0.00000000
> -0.64447126
>
> !CCSD(NOREL) expec <1.1|QMXX|1.1> 2.937794363335
> !CCSD(NOREL) expec <1.1|QMYY|1.1> -3.662513450822
> !CCSD(NOREL) expec <1.1|QMZZ|1.1> 0.724719087487
> !CCSD(NOREL) expec <1.1|QMRR|1.1> -53.392633597836
>
> !CCSD(NOREL) expec <1.1|XX|1.1> -15.839014957055
> !CCSD(NOREL) expec <1.1|YY|1.1> -20.239220166494
> !CCSD(NOREL) expec <1.1|ZZ|1.1> -17.314398474288
>
> Dalton:
>
> Total Molecular Dipole Moment
> au Debye
> z -0.64448139 -1.63810812
>
> Total Molecular quadrupole moment
>
> X Y Z
> Column 1 Column 2 Column 3
> 1 2.93725559 0.00000000 0.00000000
> 2 0.00000000 -3.66309263 0.00000000
> 3 0.00000000 0.00000000 0.72583703
>
>
> Electronic second moment of charge
>
> X Y Z
>
> 1 15.83901789 0.00000000 0.00000000
> 2 0.00000000 107.37034240 0.00000000
> 3 0.00000000 0.00000000 32.24924174
>
>
> For excited state calculations, the dipole moments are again the same with
> both codes but the quadrupole and second moments are completely different.
>
>
>
> Any help would be greatly appreciated.
>
>
>
>
>
>
>
> --
> Joseph Lane
> Department of Chemistry
> University of Otago
> P.O. Box 56
> Dunedin, New Zealand
> ph +64 3 479 5378
> email : jlane at alkali.otago.ac.nz
More information about the Molpro-user
mailing list