CCL: intermolecular potential
Hello,
I have a very basic question -- please bear with me. I am still trying to
learn the basics, and unfortunately, there is no one around me who can
answer this question.
I plan to do a simulation of a solution, with one solute molecule
surrounded by solvent molecules. Before I do the simulation, I need to
obtain parameters for the intermolecular and intramolecular potential
models.
I plan to use the intramolecular potential models available in the
literature for both the solute and the solvent. For the intermolecular
potential, I will just use the Lennard-Jones and electrostatic interaction
terms, but I still need to obtain the parameters.
I would like to obtain those parameters by fitting the ab initio
interaction energies to this intermolecular potential. This is where I
have a question (finally). When I calculate the interaction energies by
ab initio methods, should I keep the geometries of the solute and the
solvent molecules fixed to the equilibrium geometries as specified by the
intramolecular potential models that I chose? I am inclined to do so, but
I wanted to get the opinions of the experts...
Thank you very much for your patience.
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
Naoko Akiya nakiya (+ at +) engin.umich.edu
University of Michigan
Dept. of Chemical Engineering phone (734) 764-7121
Ann Arbor, MI 48109-2136 fax (734) 763-0459
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=