Atom-Atom partitioning of El. Energy
Hello world,
To my knowledge there is no option for GAMESS (US) to do atom-atom
partitionings of the electronic energy.
By this partitioning I mean that the total electronic energy, as can be
obtained using the classic (i.e. well known) RHF formula, is split up
into different atom-atoms terms (i.e. atom1-atom1, atom1-atom2, and
atom2-atom2 for a binuclear molecule.)
Because I thought it might be interesting for me to have a look at these
terms, I headed off to find out if I could not write a small program to
do this, or to change the FORTRAN code of gamess (US) as to fit my needs.
However, I am not an experienced programmer in FORTRAN (yet). So I ran
into trouble (I have studied the first part of my first book on FORTRAN
for the first time), so it is not easy for me change the GAMESS code.
Here is what I did...
I have always been programmig in Basic, or Pascal or C(++), or some even
more exotic languages, so as we have a RS6000 workstation I started
writing a C program.
I have changed the NPRINT=3 option in the GAMESS code and the rungms
executable, so that files are saved with the density, fock and hcore
matrices.
These files are read in by the program and changed as to yield some
temporary files that may be used in calculations.
So when I now make a calculation I specify NPRINT=3 and get the normal
output file plus a $JOB.at file in which I find the atom-atom parts.
Everything is written in a very simple C program, because I needed the
results asap, so I did not care for any mumbo-jumbo optimising, or
elegance (except for the elegance of using indents and such).
All this stuff has been done on the 1991 version of GAMESS (they wouldn't
let me change the new one.)
My question is now : does anybody have an example of a calculation (with
at least 10 **ten** basisfunctions) so that I can control my results.
I have confronted my findings with some theoretical facts that have to be
obeyed, and they do obey...( I off course need to know which Basis Set is
used...)
So if anybody could fax me a result, or if anybody has a reference to an
article in which the atom-atom splittimgs are proposed, I would be very
grateful if she/he could inform me on this matter...
P.S. when I have studied the whole FORTRAN book I will try to change the
1993 version in the right way (namely through the FORTRAN code of GAMESS(US)).
Any help is mighty welcome,
Regards,
|-----------------------------------------------------------------------|
| C-C Patrick Bultinck |
| / \ Dept. Physical & Inorganic Chemistry |
| C-O O-C Section Quantum Chemistry |
| | | University of Ghent |
| C-O O-C Krijgslaan 281 (S-3) |
| \ / 9000 Gent |
| C-C Belgium |
| Tel. Int'l code/32/9/264.44.44 |
| Macrocycles Fax. Int'l code/32/9/264.49.83 |
| Quantum Chemical E-mail : Patrick.Bultinck ^at^ rug.ac.be |
| Calculations |
|-----------------------------------------------------------------------|