| 9:28:35 AM PST - Tue, Mar 6th 2012 |
|
Hi there,
I am new to NWchem and would like to learn how to use the electron transfer module.
I read the relevant section of the manual, but it was not clear to me how to go about doing what I want to do, namely calculate the transfer integral between two sections of a conjugated molecule.
If I have understood correctly, I need to generate charge-localised vectors for the products and reactant states.
To me, the most obvious way of doing this would be to chop the molecule in two, add capping hydrogrens and optimise their neutral and charged geometries, then rebuild the entire molecule using these fragments.
If I am lucky, when I run the SCF, the charge will localise onto the part which has the relaxed charged geometry. It might be a good idea to reoptimise the whole structure, as long as the charge does not become delocalised.
However, if this does not work, I think I will have to try to force the charge to localise by superimposing fragment molecular orbitals and perhaps not running many SCF cycles following superposition. My question is how to do this if I have added hydrogens to cap my fragments? Alternatively, how to do the fragment calculations without adding hydrogens to cap the fragments?
Many thanks in advance for your help,
Mike
|