Symmetry-dependent answers in open shell RELCCSD

21 views
Skip to first unread message

Alexander Oleynichenko

unread,
Jun 5, 2022, 4:51:14 AMJun 5
to dirac-users
Dear DIRAC users and developers,

I'm trying to calculate MP2 and CCSD energy of the open shell nitrogen monooxide (NO) molecule. Basically I'm following the manual (http://diracprogram.org/doc/release-22/tutorials/highspin_cc/O2.html).

It seems rather weird that I obtain different MP2 and CCSD correlation energies in different symmetries (output files are attached):

symmetry   MP2                     CCSD
C1:                -0.34256475      -0.33838921
Cs:                -0.34254696      -0.33828588
C2:                -0.34269326      -0.33838921
C2v:              -0.34269326      -0.33838921
Cinfv:            -0.34269326     -0.33838921

It seems like smth is wrong with the reconstructed Fock matrix and re-calculated diagonal Fock matrix elements. I've also printed transformed one-electron part of the Hamiltonian stored in the MRCONEE file and found that all matrix elements are real and coincide for Cs and Cinfv symmetries except the only one pair of spinors (number 19 and 25, the matrix elements between them is complex in Cs, but real parts are the same as in the Cinfv case).

Please, could you take a look, what can be wrong?

Best regards,
Alexander Oleynichenko

CC_NO-Cinfv.out
CC_NO-C2v.out
CC_NO-C2.out
CC_NO-C1.out
CC_NO-Cs.out

Ilias Miroslav, doc. RNDr., PhD.

unread,
Jun 5, 2022, 5:39:35 AMJun 5
to dirac...@googlegroups.com
Hello,


Few recalculated diagonal elements are different , as seen on the output  (C1:  -0.2016263919    0.1460295075; C2v:  -0.2016263919    0.1093294917 ) - these are the open shell pi-orbitals.


Best,

Miro



From: dirac...@googlegroups.com <dirac...@googlegroups.com> on behalf of Alexander Oleynichenko <alexvole...@gmail.com>
Sent: Sunday, June 5, 2022 10:51
To: dirac-users <dirac...@googlegroups.com>
Subject: [dirac-users] Symmetry-dependent answers in open shell RELCCSD
 
--
You received this message because you are subscribed to the Google Groups "dirac-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dirac-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dirac-users/04f07fa3-9efc-4e9d-92f5-21e88a36796an%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages