[molpro-user] Molpro energy variables
Pfeifle, Mark Claudius Gerhard
mpfeifle at anl.gov
Fri Oct 6 20:18:10 CEST 2017
Cayo,
Have you added the " singles correction for state 1" to the cluster-corrected energies listed in the log file?
After that, they should be identical to the saved variables.
Best regards,
Mark
-----Original Message-----
Date: Thu, 05 Oct 2017 21:19:30 +0000
From: Cayo Emilio <cayoquimica at gmail.com>
To: "molpro-user at molpro.net" <molpro-user at molpro.net>
Subject: [molpro-user] Molpro energy variables
Content-Type: text/plain; charset="utf-8"
Dear MOLPRO users and developers,
When saving the MRCI-F12 energies into variables, the one with the Davidson correction is not the one given automatically in the output (neither the fixed, relaxed or the rotated correction).
For instance, consider the input:
basis=vqz-f12
nosym
r=2.33480912 bohr
geometry={C1,
N1,C1,r}
{multi;wf,13,1,1}
{mrci-f12;wf,13,1,1}
E0=energd0
E1=energd1
E2=energd2
E=energd
And the output is:
Reference energy -92.31349173 <(92)%203134-9173>
...
!CI-F12 STATE 1.1 Energy -92.548624174407
...
Cluster corrected energies -92.56542331 <(92)%205654-2331>
(Davidson,
fixed reference)
Cluster corrected energies -92.56526146 <(92)%205652-6146>
(Davidson,
relaxed reference)
...
SETTING E0 = -92.56531799 <(92)%205653-1799> AU
SETTING E1 = -92.56547984 <(92)%205654-7984> AU
SETTING E2 = -92.56547984 <(92)%205654-7984> AU
SETTING E = -92.56531799 <(92)%205653-1799> AU
Which energy is the variable energd0, energd1 and energd2 actually saving?
I am using the 2012.1 version.
Thanks in advance for any help,
Cayo Goncalves
More information about the Molpro-user
mailing list