Coordination numbers as reaction coordinates in constrained molecular dynamics

Abstract

Reaction coordinates for chemical reactions in solution are often of a collective nature involving solvent degrees of freedom. A recent generalization of the blue-moon ensemble technique has enabled us to determine the free energy required to change these variables, using them as constraints in a molecular dynamics simulation. We discuss a first application to the free energy of coordination fluctuations of an argon atom in liquid argon.

References

1. A. Curioni, M. Sprik, W. Andreoni, J. Hutter, M. Parrinello and H. Schiffer, J. Am. Chem. Soc., 1997, 119, 7218.
2. E. J. Meijer and M. Sprik, J. Am. Chem. Soc., 1998, 120, 6345.
3. C. J. Cerjan and W. H. Miller, J. Chem. Phys., 1980, 75, 7800.
4. T. R. Walsh and D. J. Wales, J. Chem. Soc., Faraday Trans., 1996, 92, 2505.
5. J. P. Valleau, in Computer Simulation in Material Science, ed. M. Meyer and V. Pontikis, Kluwer, Dordrecht, 1991, p. 67.
6. E. A. Carter, G. Ciccotti, J. T. Hynes and R. Kapral, Chem. Phys. Lett., 1989, 156, 472.
7. M. J. Ruiz-Montero, D. Frenkel and J. J. Brey, Mol. Phys., 1997, 90, 925.
8. D. Frenkel and B. SmitUnderstanding Molecular Simulation, Academic Press, San Diego, CA 1996.
9. T. Mülders, P. Krüger, W. Swegat and J. Schlitter, J. Chem. Phys., 1996, 104, 4869.
10. M. Sprik and G. Ciccotti, to be published.
11. J. P. Ryckaert and G. Ciccotti, J. Chem. Phys., 1983, 78, 7368.