From NWChem
You are viewing a single post from the thread title above
-
Bert Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
|
|
Forum Vet
Threads 4
Posts 597
|
|
2:52:46 PM PDT - Wed, May 30th 2012 |
|
You may want to try "task dft et" instead of "task scf et". As the orbitals are not optimized, changing the Hamiltonian used to calculate the ET could give weird results.
Bert
Quote:Adsd04 May 26th 10:27 pmHi,
I am newer with nwchem and i triying to run some imput test to coupled ET with C-DFT
similiar to some recent reports in the literature. The problems is that the VRP values obtained in my calculations are larger than the values reported for similar molecules in similar calculation conditions. I has been used different inputs models (see atached files) but the problem continue. That is why i turn to you waiting obtain any recommendation to continue with my work
My input is:
start Fc2
title "Fc2 cationico"
echo
charge 1
geometry
C 3.19286700 -1.79797600 -0.10018600
C 3.63730800 -1.21812200 1.12353800
C 4.51916500 -0.13996100 0.80164200
C 4.60983600 -0.05370300 -0.61878400
C 3.78473000 -1.07851400 -1.17622400
H 3.35407700 -1.53305700 2.11824900
H 3.62435200 -1.26239300 -2.22962700
H 5.01567100 0.50814100 1.51077800
H 5.18393300 0.67312200 -1.17679800
H 2.48778600 -2.61286900 -0.19593200
Fe 2.61500000 0.24507600 0.00078900
C 1.11014000 1.12649000 1.14889200
C 2.02985100 2.11921800 0.70707800
C 2.03480400 2.11659900 -0.71854000
C 1.11809800 1.12226500 -1.16307700
C 0.47961600 0.53873700 -0.00825000
H 0.89577800 0.88000500 2.17929600
H 2.64786700 2.73629600 1.34444400
H 2.65717900 2.73148100 -1.35379400
H 0.91222000 0.87051900 -2.19396500
C -0.47996300 -0.53955600 -0.00819600
C -1.11826400 -1.12330600 -1.16301800
C -2.03532600 -2.11728500 -0.71842100
C -2.03076500 -2.11946100 0.70718700
C -1.11093100 -1.12682800 1.14896200
H -0.91208900 -0.87189800 -2.19392800
H -0.89679500 -0.88006700 2.17934700
H -2.65765200 -2.73222100 -1.35367000
H -2.64909300 -2.73622000 1.34455900
Fe -2.61514600 -0.24527600 0.00033200
C -4.51666800 0.14317500 0.80545200
C -4.61174900 0.05344100 -0.61447400
C -3.78752500 1.07617800 -1.17696100
C -3.19176800 1.79783300 -0.10459700
C -3.63294400 1.22140400 1.12197600
H -5.01151900 -0.50277800 1.51769800
H -5.18807700 -0.67434300 -1.16892000
H -3.34652000 1.53867400 2.11502600
H -3.63008900 1.25715400 -2.23131400
H -2.48631300 2.61188700 -0.20457500
end
basis
H library 6-31g**
C library 6-31g**
Fe library lanl2dz_ecp
end
ecp
Fe library lanl2dz_ecp
end
dft
xc b3lyp
iterations 500
convergence nolevelshifting
odft
mult 2
cdft 1 20 charge 1.0
cdft 21 40 charge 0.0
vectors output Fc2dador2.mo
end
task dft
dft
xc b3lyp
iterations 500
convergence nolevelshifting
odft
mult 2
cdft 1 20 charge 0.0
cdft 21 40 charge 1.0
vectors output Fc2aceptor2.mo
end
tasK dft
et
vectors reactants Fc2dador2.mo
vectors products Fc2aceptor2.mo
end
task scf et
My output is :
-------------------
Electron Transfer Calculation
-----------------------------
MO vectors for reactants: Fc2dador2.mo
MO vectors for products : Fc2aceptor2.mo
Electronic energy of reactants H(RR) -3167.1137335703
Electronic energy of products H(PP) -3167.1137385630
Reactants/Products overlap S(RP) : 1.14D-01
Reactants/Products interaction energy:
-------------------------------------
One-electron contribution H1(RP) -652.3015883720
Beginning calculation of 2e contribution
Two-electron integral screening (tol2e) : 1.14D-08
Two-electron contribution H2(RP) 292.3372635365
Total interaction energy H(RP) -359.9643248355
Electron Transfer Coupling Energy |V(RP)| 0.8387068798
184074.857 cm-1
22.822388 eV
526.297 kcal/mol
Thank you very much
Best regards
Andrei
|
|
|
AWC's:
2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC