Isomorphic classical molecular dynamics model for an excess electron in a supercritical fluid
J. Chem. Phys. 129, 194502 (2008); doi:10.1063/1.3013357
Published 17 November 2008
You are logged in to this journal.
Ring polymer molecular dynamics (RPMD) is used to directly simulate the dynamics of an excess electron in a supercritical fluid over a broad range of densities. The accuracy of the RPMD model is tested against numerically exact path integral statistics through the use of analytical continuation techniques. At low fluid densities, the RPMD model substantially underestimates the contribution of delocalized states to the dynamics of the excess electron. However, with increasing solvent density, the RPMD model improves, nearly satisfying analytical continuation constraints at densities approaching those of typical liquids. In the high-density regime, quantum dispersion substantially decreases the self-diffusion of the solvated electron. In this regime where the dynamics of the electron is strongly coupled to the dynamics of the atoms in the fluid, trajectories that can reveal diffusive motion of the electron are long in comparison to 
.
©2008 American Institute of Physics

.
©2008 American Institute of Physics
| History: | Received 4 August 2008; accepted 14 October 2008; published 17 November 2008 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/129/194502/1 |
REFERENCES (63)
-
I. R. Craig and D. E. Manolopoulos, J. Chem. Phys. 121, 3368 (2004). [MEDLINE]
-
B. J. Braams and D. E. Manolopoulos, J. Chem. Phys. 125, 124105 (2006). [MEDLINE]
-
J. Cao and G. A. Voth, J. Chem. Phys. 100, 5106 (1994).
-
S. Jang and G. A. Voth, J. Chem. Phys. 111, 2371 (1999).
-
A. D. McLachlan, Mol. Phys. 8, 39 (1964). [ISI]
-
H. D. Meyer and W. H. Miller, J. Chem. Phys. 72, 2272 (1980). [ChemPort]
-
J. C. Tully, J. Chem. Phys. 93, 1061 (1990). [ChemPort]
-
J. C. Tully, Faraday Discuss. 110, 407 (1998).
-
S. Nielsen, R. Kapral, and G. Ciccotti, J. Chem. Phys. 115, 5805 (2001).
-
P. V. Parandekar and J. C. Tully, J. Chem. Phys. 122, 094102 (2005). [MEDLINE]
-
J. R. Schmidt, P. V. Parandekar, and J. C. Tully, J. Chem. Phys. 129, 044104 (2008). [MEDLINE]
-
H. Wang, X. Sun, and W. H. Miller, J. Chem. Phys. 108, 9726 (1998).
-
X. Sun, H. Wang, and W. H. Miller, J. Chem. Phys. 109, 7064 (1998). [ISI] [ChemPort]
-
R. Hernandez and G. A. Voth, Chem. Phys. 233, 243 (1998).
-
J. Shao and N. Makri, J. Phys. Chem. A 103, 7753 (1999). [ISI] [ChemPort]
-
J. Shao and N. Makri, J. Phys. Chem. A 103, 9479 (1999). [ISI] [ChemPort]
-
Q. Shi and E. Geva, J. Phys. Chem. A 107, 9059 (2003).
-
J. A. Poulsen, G. Nyman, and P. J. Rossky, J. Chem. Phys. 119, 12179 (2003).
-
J. Liu and W. H. Miller, J. Chem. Phys. 127, 114506 (2007). [ISI] [MEDLINE]
-
J. O. Hirschfelder, C. F. Curtiss, and R. B. Bird, Molecular Theory of Gases and Liquids (Wiley, New York, 1954).
-
N. R. Kestner, J. Jortner, M. H. Cohen, and S. A. Rice, Phys. Rev. 140, A56 (1965).
-
D. F. Coker, B. J. Berne, and D. Thirumalai, J. Chem. Phys. 86, 5689 (1987). [ISI] [ChemPort]
-
D. F. Coker and B. J. Berne, J. Chem. Phys. 89, 2128 (1988). [ISI] [ChemPort]
-
B. Space and D. F. Coker, J. Chem. Phys. 94, 1976 (1991). [ISI] [ChemPort]
-
B. Space and D. F. Coker, J. Chem. Phys. 96, 652 (1992). [ISI] [ChemPort]
-
E. Gallicchio and B. J. Berne, J. Chem. Phys. 101, 9909 (1994). [ISI] [ChemPort]
-
E. Gallicchio and B. J. Berne, J. Chem. Phys. 105, 7064 (1996). [ISI] [ChemPort]
-
R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill, New York, 1965).
-
D. Chandler and P. G. Wolynes, J. Chem. Phys. 74, 4078 (1981).
-
M. Parrinello and A. Rahman, J. Chem. Phys. 80, 860 (1984).
-
B. De Raedt, M. Sprik, and M. L. Klein, J. Chem. Phys. 80, 5719 (1984). [ISI] [ChemPort]
-
D. Chandler, in Liquids, Freezing and Glass Transition, edited by D. Levesque, J. P. Hansen, and J. Zinn-Justin (Elsevier, New York, 1990).
-
T. F. Miller III and D. E. Manolopoulos, J. Chem. Phys. 122, 184503 (2005). [MEDLINE]
-
T. F. Miller III and D. E. Manolopoulos, J. Chem. Phys. 123, 154504 (2005). [MEDLINE]
-
I. R. Craig and D. E. Manolopoulos, J. Chem. Phys. 122, 084106 (2005).
-
I. R. Craig and D. E. Manolopoulos, J. Chem. Phys. 123, 034102 (2005).
-
M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, Oxford, 1987).
-
Centroid molecular dynamics has also been employed with the Einstein relationship to obtain the diffusion coefficient. See, for example, J. Lobaugh and G. A. Voth, J. Chem. Phys. 106, 2400 (1997) [ISI] [ChemPort]
-
R. Collepardo-Guevara, I. R. Craig, and D. E. Manolopoulos, J. Chem. Phys. 128, 144502 (2008). [MEDLINE]
-
S. Jang, Y. Pak, and G. A. Voth, J. Phys. Chem. A 103, 10289 (1999). [ISI]
-
E. Geva, Q. Shi, and G. A. Voth, J. Chem. Phys. 115, 9209 (2001).
-
G. Baym and N. D. Mermin, J. Math. Phys. 2, 232 (1961). [ISI]
-
J. E. Gubernatis, M. Jarrell, R. N. Silver, and D. S. Sivia, Phys. Rev. B 44, 6011 (1991). [MEDLINE]
-
J. Skilling and S. F. Gull, SIAM-AMS Proc. 14, 167 (1984).
-
B. Buck and V. A. Macaulay, Maximum Entropy in Action (Clarendon, Oxford, 1991).
-
M. Jarrell and J. E. Gubernatis, Phys. Rep. 269, 133 (1996).
-
S. Habershon, B. J. Braams, and D. E. Manolopoulos, J. Chem. Phys. 127, 174108 (2007). [MEDLINE]
-
G. Krilov and B. J. Berne, J. Chem. Phys. 111, 9147 (1999). [ISI] [ChemPort]
-
R. K. Bryan, Eur. Biophys. J. 18, 165 (1990). [Inspec] [ChemPort]
-
C. L. Lawson and R. J. Hanson, Solving Least Square Problems (Prentice-Hall, Englewood Cliffs, NJ, 1974).
-
K. Miller, SIAM J. Math. Anal. 1, 52 (1970).
-
G. J. Martyna, M. E. Tuckerman, D. J. Tobias, and M. L. Klein, Mol. Phys. 87, 1117 (1996).
-
D. Frenkel and B. Smit, Understanding Molecular Simulation, 2nd ed. (Academic, New York, 2002).
-
M. E. Tuckerman, B. J. Berne, and G. J. Martyna, J. Chem. Phys. 97, 1990 (1992).
-
P. W. Anderson, Phys. Rev. 109, 1492 (1958).
-
D. Chandler and K. Leung, Annu. Rev. Phys. Chem. 45, 557 (1994). [ISI] [ChemPort]
-
A. L. Nichols III and D. Chandler, J. Chem. Phys. 87, 6671 (1987). [ISI] [ChemPort]
-
G. A. Voth and T. D. Hone, J. Chem. Phys. 122, 057102 (2005). [ISI]
-
Estimates of the DHe obtained during the RPMD simulations of the electron (see Table I) and during the purely classical simulations of the solvated electron (not shown) agree to within 2% at each density. The small difference in these results is a finite system size effect; since the RPMD model correctly includes quantum dispersion of the electron, it is effectively larger than in the classical simulations. This small difference in the size of the electron changes the available space available to the solvent atoms and thus slightly effects the solvent self-diffusion coefficient.
-
B. Dünweg and K. Kremer, J. Chem. Phys. 99, 6983 (1993).
-
I. -C. Yeh and G. Hummer, J. Phys. Chem. B 108, 15873 (2004).
-
We find another minor system size effect at the lowest fluid density of
*=0.1. Given the very high mobility of the electron at this density, the velocity autocorrelation function of the helium atoms does not fully decay to zero on the timescale required for the electron to travel the sidelength of the periodic box. The electron can thus encounter helium atoms with which it is dynamically correlated. -
J. Lobaugh and G. A. Voth, J. Chem. Phys. 106, 2400 (1997). [ISI] [ChemPort]








