From NWChem
Viewed 1296 times, With a total of 3 Posts

Clicked A Few Times
Threads 3
Posts 5


6:19:52 PM PDT  Wed, Nov 2nd 2016 

Is it possible to calculate the AO overlap matrix between two different geometries? Something like eq. 11 in :
[http://pubs.acs.org/doi/abs/10.1021/acs.jctc.5b01148
Thanks,
Craig




Forum Regular
Threads 1
Posts 189


6:40:58 AM PDT  Thu, Nov 3rd 2016 

This is possible in a couple different ways. The easiest way (i.e. the way that doesn't require modifying the source code) is to concatenate your two geometries into one geometry and then run a calculation on the combined geometry. For example, if you were interested in the AO overlap matrix between H2 with a bond length of 1.0 and a bond length of 1.1, you could set up the combined geometry as
set geom:dont_verify t
geometry noautosym nocenter
H 0.0 0.0 0.5
H 0.0 0.0 0.5
H 0.0 0.0 0.5
H 0.0 0.0 0.6
end
Important things to note, you will probably need to tell the code not to check if you are using a reasonable geometry since I imagine you will have overlapping (or nearly overlapping) atoms. This is accomplished with the "set geom:dont_verify t" directive. Also it is probably a good idea, if not required, to turn off symmetry and centering of the system in order to make sure you are getting the overlap that you want.
Then you just need to add
print "ao overlap"
to your DFT/SCF input block to get the code to print the overlap matrix to the output file. The ordering of the elements in the overlap matches the ordering of the atoms in the geometry input. So the printed overlap matrix will have a structure like
AA AB
BA BB
where AA corresponds to the AO overlap matrix for the first geometry, BB is for the second geometry, and AB/BA is for the overlap of the first geometry with the second geometry. The dimensions of each block are the number of basis functions for each geometry. For example, if we used the STO3G basis with the H2 example above (2 basis functions per molecule), S(1:2,1:2) would correspond to the AO overlap matrix for the H2 molecule with bond length 1.0; S(3:4,3:4) would correspond to the AO overlap matrix for the H2 molecule with bond length 1.1; and S(1:2,3:4)=S(3:4,1:2) would be the AO overlap matrix you are interested in.




Clicked A Few Times
Threads 3
Posts 5





Just Got Here
Threads 1
Posts 2


2:43:47 PM PDT  Thu, Oct 31st 2019 

Ahhh, yes! What a brilliant soln. Thanks as always Sean, and shoutsout to Craig for raising the question back before this place became a ghosttown. Going to save me a nice chunk of time...
Lingerfelt



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