The difference between molpro and mrcc

jamcas
Topic Author
Offline
New Member

4 years 7 months ago #879 by jamcas

The difference between molpro and mrcc was created by jamcas

I found a question that puzzled me, that is, for the same closded shell system, the triples (T) contribution for mrcc and molpro is totally different:.

molpro:
E(CCSD);-2980.779852732739
CCSD(T) correction:-0.324599006747
E(CCSD(T):-2981.10445174

mrcc:
E(CCSD);-2980.779847028194
CCSD(T) correction:-0.369297289025
E(CCSD(T):-2981.149144317219

My input file:

Code:

basis=aug-cc-pvdz
calc=ccsd(t)
scftype=rhf
dft=user
2
1.0 b88
1.0 pw91c
mem=60GB
dfbasiscor=none
dfbasis_scf=none
ccmaxit=50000
scfmaxit=20000
mult=1
charge=-1
scfdamp=0.7
scflshift=5
verbosity = 3

unit=angs
geom=xyz
7

C                -0.35087300    1.64374600    0.00250500
C                 1.01404000    1.70160700    0.00213000
C                 1.01401300   -1.70160600   -0.00211500
C                -0.35091000   -1.64372500   -0.00255200
V                 0.54132900    0.00190000   -1.06181900
V                 0.54127000   -0.00191700    1.06185300
V                -1.42858300    0.00001200   -0.00002500

I also tested N2 and I got the same result:

molpro:
E(CCSD);-109.278009380004
CCSD(T) correction:-0.015147862165
E(CCSD(T):-109.29315724

mrcc:
E(CCSD);-109.278009034846
CCSD(T) correction:-0.015291141772
E(CCSD(T):-109.293300176618

Please Log in or Create an account to join the conversation.

Nike
Offline
Premium Member

4 years 7 months ago #880 by Nike

Replied by Nike on topic The difference between molpro and mrcc

Dear jamcas,
1) One difference is in the integral transformation for frozen core systems, but in your closed-shell system I would hope the difference wouldn't be this big.
2) I believe MRCC does UCCSD whereas in MOLPRO you can specify RCCSD or UCCSD, so you might want to show us your input file (it won't be very big anyway).
3) Are your SCF (Hartree-Fock) energies the same?
With best wishes,
Nike

Please Log in or Create an account to join the conversation.

jamcas
Topic Author
Offline
New Member

4 years 7 months ago #881 by jamcas

Replied by jamcas on topic The difference between molpro and mrcc

2. This is my molpro input:

Code:

memory,300,m
file,2,r-m1-34n-avdz-ccsdt.wfu,old

gprint,basis,orbitals,civector;
gthresh,printci=0.05;

geometry={angstrom;
C                -0.35087300    1.64374600    0.00250500
C                 1.01404000    1.70160700    0.00213000
C                 1.01401300   -1.70160600   -0.00211500
C                -0.35091000   -1.64372500   -0.00255200
V                 0.54132900    0.00190000   -1.06181900
V                 0.54127000   -0.00191700    1.06185300
V                -1.42858300    0.00001200   -0.00002500
}

basis={
c=avdz
v=avdz
}

{rks,b,pw91c;
closed,47;
occ,47;
wf,94,1,0;
ACCU,14;
}

{rccsd(t);
shift,1,1;
core,31;
closed,47;
occ,47;
wf,94,1,0;
}

3.Yes, both the reference energy and the ccsd correlation energy are same.
Reference energy
molpro:-2979.279074961804
mrcc:-2979.279083044167

Please Log in or Create an account to join the conversation.

kallay
Offline
Administrator
Mihaly Kallay

4 years 7 months ago #882 by kallay

Replied by kallay on topic The difference between molpro and mrcc

Probably the reason is that the contribution of the occupied-virtual block of the Fock-operator to the (T) correction is also calculated in molpro. This contribution is zero for canonical HF. It is not zero for a KS reference but it is usually small.

Best regards,
Mihaly Kallay

Please Log in or Create an account to join the conversation.

jamcas
Topic Author
Offline
New Member

4 years 7 months ago #883 by jamcas

Replied by jamcas on topic The difference between molpro and mrcc

Thanks for your reply! How to get a open-shell singlet solution in mrcc? I know gaussian can performe stable check, but I don't seem to see similar keywords in mrcc.

Please Log in or Create an account to join the conversation.

kallay
Offline
Administrator
Mihaly Kallay

4 years 7 months ago #884 by kallay

Replied by kallay on topic The difference between molpro and mrcc

You can get an open-shell singlet state with MCSCF.

Best regards,
Mihaly Kallay

Please Log in or Create an account to join the conversation.

Time to create page: 0.042 seconds

The difference between molpro and mrcc

MRCC Login