Dear Q-Chem developers (and users),
I’d like to calculate a certain numbers of EE singlets, EE/SF triplets, IP_BETA doublets and EA_ALPHA doublets, all from the singlet CCSD reference.
In addition to energies of the states, I’d like to get all (non-zero) alpha and beta state and transition densities of states/between pairs of states in fchk file. Here, I face two issues:
-
The CC_TRANS_PROP keyword gives either only ground-to-excited (for value 1) or excited-to-excited (for value 2) transition densities, and I need both. Is the only way to get both to repeat the same calculation twice, once with CC_TRANS_PROP = 1 and once with CC_TRANS_PROP = 2? If that’s the case, can Q-Chem dump T and R amplitudes to some file in the first call, and then read them and only do the different analysis in the second call? I don’t see that option in the manual.
-
I’d like to calculate triplets with MS = 1, so I assume SF_STATES is the option I need, rather than EE_TRIPLETS (which gives MS = 0). I want this because then I get alpha/beta T(MS=1) → T(MS=1) state/transition densities, which cannot be constructed by Wigner-Eckart theorem from the corresponding T(MS=0) → T(MS=0) densities (while vice versa can be done). However, if I do SF_STATES, then I don’t get S → T(MS=0) alpha/beta transition densities. Is the only way to get both T(MS=1) → T(MS=1) state/transition densities and S → T(MS=0) transition densities to specify both EE_TRIPLETS and SF_STATES?
Further I’d like to get all Dyson orbitals between (EOM-EE-)CCSD and EOM-IP/EA-CCSD states. Here, I face an issue that Q-Chem 6.2 simply freezes after SCF when I add CC_DO_DYSON to $rem, if STATE_ANALYSIS is also present. I assume this is a bug, and the only quick work-around I see is to do the calculation of all the states twice, once to get densities, and once to get Dyson orbitals? If that’s the case, dumping/reading the amplitudes would again be useful.
Thanks in advance for the help.
Best regards,
Tomislav Piteša