If you have problems during the execution of MRCC, please attach the output with an adequate description of your case as well as the followings:
  • the way mrcc was invoked
  • the way build.mrcc was invoked
  • the output of build.mrcc
  • compiler version (for example: ifort -V, gfortran -v)
  • blas/lapack versions
  • as well as gcc and glibc versions

This information really helps us during troubleshooting :)

Higher-order coupled cluster calculations

  • bakowies
  • Topic Author
  • Offline
  • New Member
  • New Member
More
5 years 6 months ago #629 by bakowies
Dear Mihaly,
dear staff,

checking some higher order coupled cluster calculations, which I
had done some time ago with older versions of MRCC, I noted
seemingly strange behavior:

The CCSDT(Q) is larger in absolute terms than the CCSDT[Q] energy,
and both differ from CCSDT. The same is seen for CCSDTQP(6) and
CCSDTQP[6] relative to CCSDTQP and this is what I do expect.

CCSDTQ(P) and CCSDTQ[P] energies, on the other hand, are the
same to machine precision, and both somewhat larger in absolute
terms than CCSDTQ.

I observed the described behavior for a range of very small
molecules.

Now I reran calculations for CH4 (cc-pVDZ) using a recent
version of MRCC (2017, statically linked) and reproduced this
behavior. See attached file.

Is there any reason to expect E(CCSDTQ(P)) = E(CCSDTQ[P])? After
all CH4 has 8 valence electrons to correlate and so we still have
some way to go to FCI.

Thanks a lot,
best regards,

Dirk

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

  • kallay
  • Offline
  • Administrator
  • Administrator
  • Mihaly Kallay
More
5 years 6 months ago #630 by kallay
Replied by kallay on topic Higher-order coupled cluster calculations
Dear Dirk,
Please note that CCSDTQ(P) and CCSDTQ[P] are identical by definition, see the discussion before Eq. 35 in JCP 123, 214105 (2005).

Best regards,
Mihaly Kallay

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

Time to create page: 0.038 seconds
Powered by Kunena Forum