Efficient real-space configuration-interaction method for the simulation of multielectron mixed quantum and classical nonadiabatic molecular dynamics in the condensed phase

Journal of Chemical Physics

Volumn 119, Issue 15, 2003, Pages 7672-7684

  Larsen, Ross E ^a   Schwartz, Benjamin J ^a  

^a UNIVERSITY OF CALIFORNIA   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

DEGREES OF FREEDOM (MECHANICS); ELECTRONS; MOLECULAR DYNAMICS; POTENTIAL ENERGY;

CONDENSED PHASES;

QUANTUM THEORY;

EID: 0242352597     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.1610438     Document Type: Article

References (73)

1. An overview of various mixed QM/CM methods can be found in Classical and quantum dynamics in condensed phase simulations: proceedings of the International School of Physics "Computer simulation of rare events and dynamics of classical and quantum condensed-phase systems": Euroconference on "technical advances in particle-based computational material sciences," edited by B. J. Berne, G. Cicotti, and D. F. Coker (World Scientific, Singapore, 1998).
2. O. V. Prezhdo, J. Chem. Phys. 117, 2995 (2002).
3. W. H. Miller, J. Phys. Chem. A 105, 2942 (2001).
4. E. J. Heller, J. Chem. Phys. 62, 1544 (1975); 75, 2923 (1981).
5. E. J. Heller, J. Chem. Phys. 62, 1544 (1975); 75, 2923 (1981).
6. G. D. Billing, Phys. Chem. Chem. Phys. 4, 2865 (2002); Chem. Phys. Lett. 343, 130 (2001); C. Coletti and G. D. Billing, ibid. 368, 289 (2003).
7. G. D. Billing, Phys. Chem. Chem. Phys. 4, 2865 (2002); Chem. Phys. Lett. 343, 130 (2001); C. Coletti and G. D. Billing, ibid. 368, 289 (2003).
8. G. D. Billing, Phys. Chem. Chem. Phys. 4, 2865 (2002); Chem. Phys. Lett. 343, 130 (2001); C. Coletti and G. D. Billing, ibid. 368, 289 (2003).
9. R. Car and M. Parrinello, Phys. Rev. Lett. 55, 2471 (1985).
10. M. Dahlbom, W. Beenken, V. Sundström, and T. Pullerits, Chem. Phys. Lett. 364, 556 (2002).
11. F. Webster, P. J. Rossky, and R. A. Friesner, Comput. Phys. Commun. 63, 494 (1991).
12. O. V. Prezhdo and P. J. Rossky, J. Phys. Chem. 100, 17094 (1996).
13. E. Neria, A. Nitzan, R. N. Barnett, and U. Landman, Phys. Rev. Lett. 67, 1011 (1991).
14. B. J. Schwartz and P. J. Rossky, J. Phys. Chem. 98, 4489 (1994); 99, 2953 (1995); J. Chem. Phys. 101, 6902 (1994); 101, 6917 (1994); Phys. Rev. Lett. 72, 3282 (1994); J. Mol. Liq. 65-66, 23 (1995); B. J. Schwartz, E. R. Bittner, O. V. Prezhdo, and P. J. Rossky, J. Chem. Phys. 104, 5942 (1996).
15. B. J. Schwartz and P. J. Rossky, J. Phys. Chem. 98, 4489 (1994); 99, 2953 (1995); J. Chem. Phys. 101, 6902 (1994); 101, 6917 (1994); Phys. Rev. Lett. 72, 3282 (1994); J. Mol. Liq. 65-66, 23 (1995); B. J. Schwartz, E. R. Bittner, O. V. Prezhdo, and P. J. Rossky, J. Chem. Phys. 104, 5942 (1996).
16. B. J. Schwartz and P. J. Rossky, J. Phys. Chem. 98, 4489 (1994); 99, 2953 (1995); J. Chem. Phys. 101, 6902 (1994); 101, 6917 (1994); Phys. Rev. Lett. 72, 3282 (1994); J. Mol. Liq. 65-66, 23 (1995); B. J. Schwartz, E. R. Bittner, O. V. Prezhdo, and P. J. Rossky, J. Chem. Phys. 104, 5942 (1996).
17. B. J. Schwartz and P. J. Rossky, J. Phys. Chem. 98, 4489 (1994); 99, 2953 (1995); J. Chem. Phys. 101, 6902 (1994); 101, 6917 (1994); Phys. Rev. Lett. 72, 3282 (1994); J. Mol. Liq. 65-66, 23 (1995); B. J. Schwartz, E. R. Bittner, O. V. Prezhdo, and P. J. Rossky, J. Chem. Phys. 104, 5942 (1996).
18. B. J. Schwartz and P. J. Rossky, J. Phys. Chem. 98, 4489 (1994); 99, 2953 (1995); J. Chem. Phys. 101, 6902 (1994); 101, 6917 (1994); Phys. Rev. Lett. 72, 3282 (1994); J. Mol. Liq. 65-66, 23 (1995); B. J. Schwartz, E. R. Bittner, O. V. Prezhdo, and P. J. Rossky, J. Chem. Phys. 104, 5942 (1996).
19. B. J. Schwartz and P. J. Rossky, J. Phys. Chem. 98, 4489 (1994); 99, 2953 (1995); J. Chem. Phys. 101, 6902 (1994); 101, 6917 (1994); Phys. Rev. Lett. 72, 3282 (1994); J. Mol. Liq. 65-66, 23 (1995); B. J. Schwartz, E. R. Bittner, O. V. Prezhdo, and P. J. Rossky, J. Chem. Phys. 104, 5942 (1996).
20. B. J. Schwartz and P. J. Rossky, J. Phys. Chem. 98, 4489 (1994); 99, 2953 (1995); J. Chem. Phys. 101, 6902 (1994); 101, 6917 (1994); Phys. Rev. Lett. 72, 3282 (1994); J. Mol. Liq. 65-66, 23 (1995); B. J. Schwartz, E. R. Bittner, O. V. Prezhdo, and P. J. Rossky, J. Chem. Phys. 104, 5942 (1996).
21. M. D. Hack, A. M. Wensmann, D. G. Truhlar, M. Ben-Nun, and T. J. Martinez, J. Chem. Phys. 115, 1172 (2001).
22. S. Hammes-Schiffer and J. C. Tully, J. Chem. Phys. 101, 4657 (1994).
23. W.-S. Sheu and P. J. Rossky, J. Phys. Chem. 100, 1295 (1996).
24. C. J. Smallwood, R. E. Larsen, W. B. Bosma, and B. J. Schwartz, J. Chem. Phys. (in press).
25. A. Staib and D. Borgis, J. Chem. Phys. 103, 2642 (1995).
26. W. Stier and O. V. Prezhdo, J. Phys. Chem. B 106, 8047 (2002).
27. One notable exception to the single-electron trend is a recent semiempirical CI calculation of photoreaction dynamics, by M. Persico, G. Granucci, S. Inglese, T. Laino, and A. Toniolo, J. Mol. Struct.: THEOCHEM 621, 119 (2003).
28. S. E. Bradforth and P. Jungwirth, J. Phys. Chem. A 106, 1286 (2002).
29. A. Szabo and N. S. Ostlund, Modem Quantum Chemistry: Introduction to Advanced Electronic Structure Theory (Dover, New York, 1989).
30. O. V. Prezhdo and P. J. Rossky, J. Chem. Phys. 107, 825 (1997).
31. J. C. Tully, J. Chem. Phys. 93, 1061 (1990).
32. K. F. Wong and P. J. Rossky, J. Phys. Chem. A 105, 2546 (2001); J. Chem. Phys. 116, 8418 (2002); 116, 8429 (2002).
33. K. F. Wong and P. J. Rossky, J. Phys. Chem. A 105, 2546 (2001); J. Chem. Phys. 116, 8418 (2002); 116, 8429 (2002).
34. K.F. Wong and P. J. Rossky, J. Phys. Chem. A 105, 2546 (2001); 116, 8418 (2002); 116, 8429 (2002).
35. J. Schnitker and P. J. Rossky, J. Chem. Phys. 86, 3462 (1987).
36. P. J. Rossky and J. Schnitker, J. Phys. Chem. 92, 4277 (1988).
37. J. Schnitker, K. Motakabbir, P. J. Rossky, and R. Friesner, Phys. Rev. Lett. 60, 456 (1988); D. Borgis and A. Staib, Chem. Phys. Lett. 230, 405 (1994) ; D. Borgis and S. Bratos, J. Mol. Struct. 437, 537 (1997); C. Nicolas, A, Boutin, B. Lévy, and D. Borgis, J. Chem. Phys. 118, 9689 (2003).
38. J. Schnitker, K. Motakabbir, P. J. Rossky, and R. Friesner, Phys. Rev. Lett. 60, 456 (1988); D. Borgis and A. Staib, Chem. Phys. Lett. 230, 405 (1994) ; D. Borgis and S. Bratos, J. Mol. Struct. 437, 537 (1997); C. Nicolas, A, Boutin, B. Lévy, and D. Borgis, J. Chem. Phys. 118, 9689 (2003).
39. J. Schnitker, K. Motakabbir, P. J. Rossky, and R. Friesner, Phys. Rev. Lett. 60, 456 (1988); D. Borgis and A. Staib, Chem. Phys. Lett. 230, 405 (1994) ; D. Borgis and S. Bratos, J. Mol. Struct. 437, 537 (1997); C. Nicolas, A, Boutin, B. Lévy, and D. Borgis, J. Chem. Phys. 118, 9689 (2003).
40. J. Schnitker, K. Motakabbir, P. J. Rossky, and R. Friesner, Phys. Rev. Lett. 60, 456 (1988); D. Borgis and A. Staib, Chem. Phys. Lett. 230, 405 (1994) ; D. Borgis and S. Bratos, J. Mol. Struct. 437, 537 (1997); C. Nicolas, A, Boutin, B. Lévy, and D. Borgis, J. Chem. Phys. 118, 9689 (2003).
41. An overview of the configuration interaction method can be found in Chap. 4 of Ref. 20
42. 12 cannot, of course, be represented as a direct product of single-electron operators: interactions between the electrons will not be diagonal in a one-electron basis.
43. Because the Coulomb interaction does not affect the spin of an electron, spin flips can only be induced by terms in the electron-solvent Hamiltonians (e.g., by magnetic fields induced by moving charges in the solvent) or by direct spin-spin interactions. In this work, we will neglect such spin interactions, but we note that the extension to such cases is straightfor-ward.
44. D. Finocchiaro, M. Pellegrini, and P. Bientinesi, J. Comput. Phys. 146, 707 (1998).
45. ij results from the symmetry of the cubic lattice.
46. -8, therefore, we can use our tabulated values to compute ø between cubes within 16 grid points of each other, approximating all other values of ø by 1/r.
47. A. Alavi, J. Chem. Phys. 113, 7735 (2000).
48. M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, London, 1992).
49. N. Basco, G. A. Kenney, and D. C. Walker, Chem. Commun. 917 (1969); N. Basco, G. A. Kenney-Wallace, S. K. Vidyarthi, and D. C. Walker, Can. J. Chem. 50, 2059 (1972); D. Meisel, G. Czapski, M. S. Matheson, and W. A. Mulac, Int. J. Radiat. Phys. Chem. 7, 233 (1975).
50. N. Basco, G. A. Kenney, and D. C. Walker, Chem. Commun. 917 (1969); N. Basco, G. A. Kenney-Wallace, S. K. Vidyarthi, and D. C. Walker, Can. J. Chem. 50, 2059 (1972); D. Meisel, G. Czapski, M. S. Matheson, and W. A. Mulac, Int. J. Radiat. Phys. Chem. 7, 233 (1975).
51. N. Basco, G. A. Kenney, and D. C. Walker, Chem. Commun. 917 (1969); N. Basco, G. A. Kenney-Wallace, S. K. Vidyarthi, and D. C. Walker, Can. J. Chem. 50, 2059 (1972); D. Meisel, G. Czapski, M. S. Matheson, and W. A. Mulac, Int. J. Radiat. Phys. Chem. 7, 233 (1975).
52. K. J. Schmidt and D. M. Bartels, Chem. Phys. 190, 145 (1995).
53. See D. F. Feng, K. Fueki, and L. Kevan, J. Chem. Phys. 58, 3281 (1978), and references therein.
54. H.-P. Kaukonen, R. N. Barnett, and U. Landman, J. Chem. Phys. 97, 1365 (1992).
55. Z. H. Deng, G. J. Martyna, and M. L. Klein, Phys. Rev. Lett. 68, 2496 (1992); 71, 267 (1993); J. Chem. Phys. 100, 7590 (1994); G. J. Martyna, Z. H. Deng, and M. L. Klein, ibid. 98, 555 (1993).
56. Z. H. Deng, G. J. Martyna, and M. L. Klein, Phys. Rev. Lett. 68, 2496 (1992); 71, 267 (1993); J. Chem. Phys. 100, 7590 (1994); G. J. Martyna, Z. H. Deng, and M. L. Klein, ibid. 98, 555 (1993).
57. Z. H. Deng, G. J. Martyna, and M. L. Klein, Phys. Rev. Lett. 68, 2496 (1992); 71, 267 (1993); J. Chem. Phys. 100, 7590 (1994); G. J. Martyna, Z. H. Deng, and M. L. Klein, ibid. 98, 555 (1993).
58. Z. H. Deng, G. J. Martyna, and M. L. Klein, Phys. Rev. Lett. 68, 2496 (1992); 71, 267 (1993); J. Chem. Phys. 100, 7590 (1994); G. J. Martyna, Z. H. Deng, and M. L. Klein, ibid. 98, 555 (1993).
59. K. Toukan and A. Rahman, Phys. Rev. B 31, 2643 (1985).
60. F. Webster, E. T. Wang, P. J. Rossky, and R. A. Friesner, J. Chem. Phys. 100, 4835 (1994).
61. J. J. Sakurai, Modern Quantum Mechanics, revised ed. (Addison-Wesley, Reading, MA, 1994), Chap. 2.
62. We define energy conservation by requiring any long-term drift in the total energy to be within the rms fluctuation over the course of a 2 ps run. The rms energy fluctuation is 0.08 eV for △t = 1.0 and 0.02 eV for △t = 0.5fs.
63. update= 4 fs criterion is actually overly generous.
64. R. E. Larsen and B. J. Schwartz (unpublished).
65. A simple uncertainty-principle argument shows that the kinetic energy reduction associated with going from a cavity 4.0 Å in diameter to one 4. 8 Å in diameter is less than ∼0.1 eV.
66. This change of shape from an "s-like" symmetry to a "p-like" symmetry is similar to what happens when a single hydrated electron is excited (see Refs. 9, 11, and 23), so it is not suprising that the time scale and mechanism of the Stokes shift are similar for the hydrated electron and the singlet dielectron.
67. The two cavities of the excited state dielectron cannot properly be thought of as two separated single electrons because the exchange energy of ∼0.5 eV indicates that each one-electron basis state has amplitude in both holes.
68. i). At first glance, the latter term appears quite similar to the Pulay force (Ref. 51) that corrects the Hellmann-Feynman force when the basis functions used to describe the quantum system vary with, e.g., nuclear postion. A difference between the adiabatic and nonadiabatic cases, however, is that this nonadiabatic "Pulay-type" term comes from changes in the adiabatic expansion coefficients rather than from changes in the basis. Unlike the Pulay force, therefore, the time dependence of the CI expansion coefficients is physically meaningful, and does not vanish in the limit that the basis set (here |nm〉±) is complete.
69. i〉) = 0, applying the product rule for derivatives and integrating by parts once.
70. P. Pulay, in Applications of Electronic Structure Theory, edited by H. F. Schaefer III (Plenum, New York, 1977); M. Di Ventra and S. T. Pantelides, Phys. Rev. B 61, 16207 (2000).
71. P. Pulay, in Applications of Electronic Structure Theory, edited by H. F. Schaefer III (Plenum, New York, 1977); M. Di Ventra and S. T. Pantelides, Phys. Rev. B 61, 16207 (2000).
72. One difference between the usual single-electron force and that for our CI case is that the CI expansion coefficients change in time. This change gives rise to a correction to the HF force, the Pulay force, Ref. 51, whenever a finite CI basis set is used (as will always be the case in practice). The specific form of the Pulay force for CI wave functions is straightfor-ward, though tedious, to derive using first-order perturbation theory and the approach of Di Ventra and Pantelides (Ref. 51). We have examined the Pulay force for the hydrated dielectron system discussed in Sec. IV and found it to be negligible in comparison to the HF forces themselves, so we do not bother to explicitly display it here.
73. C. C. Huang, G. S. Couch, E. F. Pettersen, and T. E. Ferrin, in Proceedings of the Pacific Symposium on Biocomputing, 1996, Vol. 1, p. 724. The CHIMERA code is freely available on the worldwide web at http://www.cgl.ucsf.edu/chimera