Hi NWChem users,
I am trying to do DFT-NEB calculation of C-N reaction. It needs ‘mult’ keyword since there are odd number of electrons in the system. ‘mult’ is the spin multiplicity of the system. I found there are two
ways to calculate the 'mult' values.
One is ‘number of alpha electrons (i.e., electron with the spin quantum number of 1/2) minus beta electrons (i.e., electron with the spin quantum number of -1/2) plus 1’ other is ‘2S+1’ where S is the 0.5*(number of unpaired electrons).
Using the both definition, I got mult = 6 for the C-N system. I am wondering whether I am doing any mistake. Here is the DFT-NEB input script.
start neb-c-n
basis
* library 6-311++G**
end
dft
xc b3lyp
maxiter 5000
cgmin
mult 6
end
- define the start geometry ####
geometry units angstroms print xyz noautosym
C 39.59843 39.09966 39.40598
N 39.26613 41.94704 40.23538
end
- define the end geometry ####
geometry endgeom units angstroms print xyz noautosym
C 39.20759 39.88065 39.43266
N 39.08069 41.28918 39.77152
end
neb
maxiter 1000
nbeads 20
stepsize 0.10
kbeads 0.1
end
task dft neb ignore
I will appreciate your comments regarding the mult value calculation and input script.
Thank you.
kind regards,
Sanjib
|