The MaxFlux algorithm for calculating variationally optimized reaction paths for conformational transitions in many body systems at finite temperature
J. Chem. Phys. 107, 5000 (1997); doi:10.1063/1.474863
Issue Date: 1 October 1997
You are not logged in to this journal. Log in
An algorithm for the calculation of reaction paths between known reactant and product states in systems of low or high dimension is described. The optimal reaction path is defined as the path of maximum flux for a diffusive dynamics assuming isotropic friction. The resulting reaction path is temperature dependent and independent of the magnitude of the friction. Comparison is made with a number of algorithms designed for the calculation of minimum-energy reaction paths. Applications to two model potentials and an extended atom model of a dipeptide are presented. The applications demonstrate the ability of the algorithm to isolate a temperature-dependent optimal reaction path for a high dimensional molecular system. ©1997 American Institute of Physics.
| History: | Received 29 April 1997; accepted 26 June 1997 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/107/5000/1 |
KEYWORDS and PACS
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (28)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- C. J. Cerjan and W. H. Miller, J. Chem. Phys. 75, 2800 (1981);
- K. Müller and L. D. Brown, Theor. Chim. Acta (Berlin) 53, 75 (1979).
- K. Müller,
Angew. Chem. Int. Ed. Engl. 19, 1 (1980 ). - S. Fischer and M. Karplus, Chem. Phys. Lett. 96, 5272 (1992);
- A. Ulitsky and D. Shalloway, J. Chem. Phys. 106, 10099 (1997).
- R. Elber and M. Karplus,
Chem. Phys. Lett. 139, 375 (1987 ). - L. R. Pratt, J. Chem. Phys. 85, 5045 (1986).
- R. Czerminski and R. Elber,
Int. J. Quantum Chem. 24, 167 (1990 ). - R. Czerminski and R. Elber, J. Chem. Phys. 92, 5580 (1990).
- A. Ulitsky and R. Elber, J. Chem. Phys. 92, 1510 (1990).
- G. Mills, H. Jónsson, and G. K. Schenter,
Surf. Sci. 324, 305 (1995 );
O. S. Smart, - R. Elber, in Recent Developments in Theoretical Studies of Proteins, edited by R. Elber (World Scientific, Singapore, 1996).
- P. G. Wolynes, in Complex Systems, SFI Studies in the Sciences of Complexity, edited by D. L. Stein, (Addison-Wesley Longman, 1989).
- L. Onsager and S. Machlup,
Phys. Rev. 91, 1505, 1512 (1953 ). - R. Olender and R. Elber, J. Chem. Phys. 105, 9299 (1996).
- R. E. Gillian and K. R. Wilson, J. Chem. Phys. 97, 1757 (1992).
- J. E. Straub, B. J. Berne, and B. Roux, J. Chem. Phys. 93, 6804 (1990);
- D. G. Truhlar and B. C. Garrett,
Acc. Chem. Res. 13, 440 (1980 ). - C. W. Gardiner, Handbook of Stochastic Methods (Springer, New York, 1983).
- B. J. Berne, M. Borkovec, and J. E. Straub,
J. Phys. Chem. 92, 3711 (1988 ). - A. M. Berezhkovskii, L. M. Berezhkovskii, and V. Yu. Zitzerman,
Chem. Phys. 130, 55 (1989 ). - R. S. Larson,
Physica A 137, 295 (1986 ). - M. Berkowitz, J. D. Morgan, J. A. McCammon, and S. H. Northrup, J. Chem. Phys. 79, 5563 (1983).
- A more general definition of the path of least resistance was presented for the case of isotropic, spatially dependent friction
(r) (Ref. 23). It was defined as the path which minimizes the line integral =
![[integral]](http://scitation.aip.org/stockgif3/int.gif)
e
dl(r). The effective potential
(r) = U(r)+kBT ln[
(r)/
(rR)], where
(rR) is the friction at an arbitrarily chosen reference point. For the case of spatially independent friction
considered in this work (r) = U(r). In general, it is best to optimize
rather than
. However, while methods exist for the calculation of spatially dependent diffusion tensors (Refs. 17 and 28) the calculation of the friction along all coordinates in a multidimensional space remains computationally forbidding.
- P. Amara and J. E. Straub,
J. Phys. Chem. 99, 14 840 (1995 ). - K. Fukui, S. Kato, and H. Fujimoto,
J. Am. Chem. Soc. 97, 1 (1975 ). - B. R. Brooks, R. Bruccoleri, B. Olafson, D. States, S. Swaninathan, and M. Karplus,
J. Comput. Chem. 4, 187 (1983 ). - J. E. Straub, M. Borkovec, and B. J. Berne,
J. Phys. Chem. 91, 4995 (1987 ).
D. T. Nguyen and D. A. Case,
T. Lazaridis, D. J. Tobias, C. L. Brooks III, and M. E. Paulaitis, J. Chem. Phys. 95, 7612 (1991).



HCN(X 1
+)+H2(X 1




