Non-adiabatic coupling in EOM-CCSD

Questions
#1

I want to get non-adiabatic coupling between two excited states in EOM-CCSD, but I get an error saying
Q-Chem fatal error occurred in module /scratch/qcdevops/jenkins/workspace/build_qchem_linux_distrib/tags/qc601rc1/qchem/ccman2/qchem/ccman2_main.C, line 26:

Key ‘i_o1o1o1o2’ is not known.

Is NACV between EOM-CCSD excited states unavailable?

#2

Please copy–paste a sample input. NAC calculations are available with EOM-CCSD methods.

#3

$molecule
0 1
O -0.340078 0.171394 -0.000000
H 0.627988 0.226708 0.000000
H -0.611448 1.102297 -0.000000
$end

$rem
INPUT_BOHR FALSE
MAX_SCF_CYCLES 300
METHOD EOM-CCSD
BASIS cc-pVDZ
SYM_IGNORE TRUE
EE_SINGLETS 2
EE_TRIPLETS 4
PRINT_ORBITALS 10
MEM_TOTAL 150000
JOBTYPE SP
SCF_ALGORITHM DIIS
CC_SYMMETRY FALSE
STATE_ANALYSIS TRUE
CC_EOM_PROP TRUE
CC_TRANS_PROP TRUE
CC_CONVERGENCE 5
CALC_NAC true
$end

$trans_prop
state_list
ee_singlets 1 1
ee_singlets 1 2
ref
end_list
calc nac
$end

#4

CALC_NAC should be set to 2 for EOM-CC methods. Please refer to the Q-Chem Manual for more information (7.10.21 EOM-CC Optimization and Properties Job Control‣ 7.10 Coupled-Cluster Excited-State and Open-Shell Methods ‣ Chapter 7 Open-Shell and Excited-State Methods ‣ Q-Chem 6.2 User’s Manual).

1 Like