How to compute overlap between two wavefunctions in NWChem

From NWChem

Viewed 2905 times, With a total of 3 Posts
Jump to: navigation, search

Just Got Here
Threads 1
Posts 2
Hi,

I am a new user to NWChem.
As an initial toy problem, I am trying to compute an overlap between <phi_1|phi_2>. Where phi_1 is a wavefunction of a diatomic molecule A-B relaxed without any constraint and phi_2 is a wavefunction of the same diatomic molecule relaxed with constrained DFT when B is in anionic state and A is in a cationic state [i.e., A(+)-B(-) ].
I was able to printout both the wavefunctions phi_1 and phi_2 in ASCII format using mov2asc utility. However, I am still confused about the format of the printed wavefunctions. Could someone please shed some light on how to compute <phi_1|phi_2>.

Any response will be much appreciated. Thanks.

Regards,
Ghanshyam

  • Bert Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
    Profile
    Send PM
Forum Vet
Threads 4
Posts 597
You would need more than just the orbitals. You would need the overlap matrix for the AO basis too.

There is a quick way to get what you need. Look at http://www.nwchem-sw.org/index.php/Release61:Electron_Transfer . Here the overlap between two wave functions is calculated. Hence, just feeding the two wave functions to ET should give you what you need in the output.

Just Got Here
Threads 1
Posts 2
Thanks Bert, that was very useful.

Regards,
Ghanshyam

Just Got Here
Threads 0
Posts 3
frozen-density embedding
Dear all,

Is the frozen-density embedding calculation method available in NWCHEM?
Regards


Forum >> NWChem's corner >> General Topics



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


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