The excitation energy on EOM-CCSD QM/MM calculation?

From NWChem

You are viewing a single post from the thread title above
Jump to: navigation, search

Click here for full thread
Clicked A Few Times
Threads 12
Posts 43
Hi everyone,

I performed EOM-CCSD QM/MM calculation for H2O molecules as a test case. My question is how to get the excitation energy? NWChem does not print it explicitly.

The following is the input file:
start e4h2o

memory total 4500 Mb

charge -1

basis spherical
* library 6-31++g** 

thresh 1.0e-10
tol2e 1.0e-10

system e4h2o_opt
cutoff 1.0 qmmm 1.0
noshake solute

mm_charges expand all

freeze atomic core
tilesize 8
attilesize 60
2emet 14
diis 5
thresh 1.0d-5
nroots 2
eomsol 2
idiskx 1

task tce qmmm energy

The following is part of the output file, if you interested in the file, I can email it to you.
EOMCCSD root nr.                     1

Singles contributions
22a (alpha) --- 21a (alpha) 0.8957902018
25a (alpha) --- 21a (alpha) -0.1010208625
26a (alpha) --- 21a (alpha) 0.1495866218
27a (alpha) --- 21a (alpha) 0.2654077772
28a (alpha) --- 21a (alpha) -0.1891541423
34a (alpha) --- 21a (alpha) 0.1191698732

Doubles contributions
EOMCCSD root nr. 2

Singles contributions
22a (alpha) --- 21a (alpha) 0.1078353808
23a (alpha) --- 21a (alpha) 0.8192252691
24a (alpha) --- 21a (alpha) -0.2875936030
26a (alpha) --- 21a (alpha) -0.2268543499
27a (alpha) --- 21a (alpha) -0.3199715289
29a (alpha) --- 21a (alpha) 0.1373032935

Doubles contributions

Parallel integral file used     436 records with       0 large values

Bq nuclear interaction energy =   6.89703372905328     

QM/MM Energy

quantum energy                      -305.169313270 (-.801222E+06  kjoule/mol)
quantum energy adjusted -305.169313270 (-.801222E+06 kjoule/mol)
quantum energy internal -312.066346999 (-.819330E+06 kjoule/mol)
Bq-nuclear energy 6.897033729 (0.181082E+05 kjoule/mol)
Bq-electron energy 0.000000000 (0.000000E+00 kjoule/mol)
classical energy -4.752909506 (-.124788E+05 kjoule/mol)
total qmmm energy -309.922222776 (-.813701E+06 kjoule/mol)

Thanks in advance.

Who's here now Members 0 Guests 0 Bots/Crawler 0

AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC