Time-reversible ab initio molecular dynamics

Niklasson, Anders M. N.; Tymczak, C. J.; Challacombe, Matt
April 2007
Journal of Chemical Physics;4/14/2007, Vol. 126 Issue 14, p144103
Academic Journal
Time-reversible ab initio molecular dynamics based on a lossless multichannel decomposition for the integration of the electronic degrees of freedom [Phys. Rev. Lett. 97, 123001 (2006)] is explored. The authors present a lossless time-reversible density matrix molecular dynamics scheme. This approach often allows for stable Hartree-Fock simulations using only one single self-consistent field cycle per time step. They also present a generalization, introducing an additional “forcing” term, that in a special case includes a hybrid Lagrangian, i.e., Car-Parrinello-type, method, which can systematically be constrained to the Born-Oppenheimer potential energy surface by using an increasing number of self-consistency cycles in the nuclear force calculations. Furthermore, in analog to the reversible and symplectic leapfrog or velocity Verlet schemes, where not only the position but also the velocity is propagated, the authors propose a Verlet-type density velocity formalism for time-reversible Born-Oppenheimer molecular dynamics.


Related Articles

  • Semiclassical multistate Liouville dynamics in the adiabatic representation. Donoso, Arnaldo; Martens, Craig C. // Journal of Chemical Physics;3/1/2000, Vol. 112 Issue 9 

    In this paper, we describe implementation of the semiclassical Liouville method for simulating molecular dynamics on coupled electronic surfaces in the electronic adiabatic representation. We cast the formalism in terms of semiclassical motion on Born-Oppenheimer potential energy surfaces with...

  • Nonadiabatic excited-state molecular dynamics: Numerical tests of convergence and parameters. Nelson, Tammie; Fernandez-Alberti, Sebastian; Chernyak, Vladimir; Roitberg, Adrian E.; Tretiak, Sergei // Journal of Chemical Physics;2/4/2012, Vol. 136 Issue 5, p054108 

    Nonadiabatic molecular dynamics simulations, involving multiple Born-Oppenheimer potential energy surfaces, often require a large number of independent trajectories in order to achieve the desired convergence of the results, and simulation relies on different parameters that should be tested and...

  • Finite-temperature electronic simulations without the Born-Oppenheimer constraint. Mazzola, Guglielmo; Zen, Andrea; Sorella, Sandro // Journal of Chemical Physics;10/7/2012, Vol. 137 Issue 13, p134112 

    The adiabatic approximation, typically assumed when performing standard Born-Oppenheimer (BO) molecular dynamics, can become unreliable at finite temperature, and specifically when the temperature is larger than the electronic energy gap between the ground state and the low-lying excited states....

  • Response to “Comment on ‘Curvy-steps approach to constraint-free extended-Lagrangian ab initio molecular dynamics, using atom-centered basis functions: Convergence toward Born–Oppenheimer trajectories’ ” [J. Chem. Phys. 123, 027101 (2005)] Herbert, John M.; Head-Gordon, Martin // Journal of Chemical Physics;7/8/2005, Vol. 123 Issue 2, p027102 

    In response to the Comment by Iyengar et al. [(unpublished)], we reiterate the differences between our “curvy-steps” implementation of extended-Lagrangian molecular dynamics and their “atom-centered density matrix propagation” technique. We further clarify our conclusions...

  • Electronic decoherence time for non-Born-Oppenheimer trajectories. Jasper, Ahren W.; Truhlar, Donald G. // Journal of Chemical Physics;8/8/2005, Vol. 123 Issue 6, p064103 

    An expression is obtained for the electronic decoherence time of the reduced density electronic matrix in mixed quantum-classical molecular-dynamics simulations. The result is obtained by assuming that decoherence is dominated by the time dependence of the overlap of minimum-uncertainty packets...

  • Hydrated proton clusters: Ab initio molecular dynamics simulation and simulated annealing. Dongqing Wei; Salahub, D.R. // Journal of Chemical Physics;4/8/1997, Vol. 106 Issue 14, p6086 

    Describes an ab initio molecular dynamics simulation technique employing the Born-Oppenheimer approach in the framework of a Gaussian implementation of John-Sham density functional theory. Simulation results for H[sub 5]O[sub 2]; Optimization of structures of large hydrated proton clusters;...

  • Thermal behavior of a 13-molecule hydrogen cluster under pressure. Santamaria, Rubén; Soullard, Jacques; Jellinek, Julius // Journal of Chemical Physics;Mar2010, Vol. 132 Issue 12, p124505 

    The thermal behavior of a 13-molecule hydrogen cluster is studied as a function of pressure and temperature using a combination of trajectory and density functional theory simulations. The analysis is performed in terms of characteristic descriptors such as caloric curve, root-mean-square bond...

  • Communications: Evaluation of the nondiabaticity of quantum molecular dynamics with the dephasing representation of quantum fidelity. Zimmermann, Tomásˇ; Vanícˇek, Jirˇí // Journal of Chemical Physics;6/28/2010, Vol. 132 Issue 24, p241101 

    We propose an approximate method for evaluating the importance of non-Born–Oppenheimer effects on the quantum dynamics of nuclei. The method uses a generalization of the dephasing representation (DR) of quantum fidelity to several diabatic potential energy surfaces and its computational...

  • Symmetry-invariant reaction-path potentials. Ischtwan, Josef; Collins, Michael A. // Journal of Chemical Physics;6/1/1991, Vol. 94 Issue 11, p7084 

    A general method is presented for constructing a molecular potential-energy surface, based on a reaction path, which is invariant to operations of the complete nuclear permutation inversion group. The method is illustrated with a number of examples, including an approximate potential for NH+3+D2...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics