SUMMARY: Potential parameter fitting
I received a few attempts to summarise the responses, so here is a
brief summary.
I asked what programs were avialable for potential parameter fitting
to energy hypersurfaces, such as ab-initio energies. The main
responses were
1. NBP by the late Prof D. E. Williams,
http://www.louisville.edu/~dewill01/nbp.html
This seems to be about the best suited program for such tasks.
It used to be obtainable directly from Prof Williams, but none of
my respondents had any idea how to obtain it now, as he recently
passed away. I hope this does not become lost to the community.
2. A number of respondents pointed me to GULP by Julian Gale
http://www.ch.ic.ac.uk/gale/Research/gulp.html
which is slightly embarrassing since I thought I knew that program
fairly well :) It can indeed be used to fit to an energy hypersurface
(which I erroneously thought it wasn't well suited for).
GULP uses a 3D periodic model, which isn't what I need for this
particular calculation, (A simple isolated moecular cluster), but
of course it's not to hard to achieve this by judicious use of a
large supercell and interaction cutoffs and turning off the
long-ranged part of coulomb interactions.
3. Though not what I was looking for Andrew Horsfield
(horsfield()at()fecit.co.uk)suggested that The OXON tight binding
package also fits (and implements) Tersoff type potentials.
4. "Kjeld Rasmussen studied with Lifson and has modified that august
program over the years to be a parameter-fitting tool. It can fit to
structural parameters and vibrational frequencies but i do not believe
to a hypersurface. Kjeld readily gives the program away. Some of the
background to his program is given at:
[1]http://struktur.kemi.dtu.dk/cff/cffhome.html"
5. "Arnold Hagler had a similar program (I believe called
"Probe") that
also came from Lifson's program and it was modified to fit
hypersurface data but I think it is now a proprietary commercial
product for MSI consortium members."
6. Marcelo Zaldini has been "developing a methodology, called SJBR,
to do the parametrization of the intermolecular terms (coulomb +
lennard_jones), based upon ab initio calculations, in an automatic
fashion." which he hopes to release at some future time, when it
is finished.
7. Other co-respondents recommended general least squares fitting
methods such as the least-squares module in Scientific Python,
available from
http://dirac.cnrs-orleans.fr/programs/scientific.html,
or even gnuplot! Although of course the built-in fit is only
1-dimensional and therefore not suitable, the gnuplot FAQ contains
links to a number of other general purpose leasy-squares programs,
such as Fudgit.
None of these really fulfil the spirit of my request which was for a
program which handles the bookkeeping and set-up for molecular
models.
Thanks to
Konrad Hinsen <hinsen()at()cnrs-orleans.fr>,
Manuel Melle-Franco <manu()at()ciamserv.ciam.unibo.it>,
Irena Efremenko <chrirena()at()techunix.technion.ac.il>,
Andrew Rohl <andrew()at()power.curtin.edu.au>,
Brian J. Teppen" <teppen()at()msu.edu>,
Marcelo Zaldini Hernandes <zaldini()at()npd.ufpe.br>,
Daquan Gao <dgao()at()chem.iupui.edu>
Andrew Horsfield <horsfield()at()fecit.co.uk>
for the suggestions.
Keith Refson
--
Dr Keith Refson, "Paradigm is a word too often used by those who
would
Dept of Earth Sciences like to have a new idea but cannot think of
one."
Parks Road, -- Mervyn King, Deputy Governor, Bank of England
Oxford OX1 3PR, UK
Keith.Refson()at() Tel: 01865 272026
earth.ox.ac.uk Fax: 01865 272072