Related Experiment Video
Updated: May 25, 2026

Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package
Published on: September 17, 2021
Nonadiabatic excited-state molecular dynamics: numerical tests of convergence and parameters
Tammie Nelson1, Sebastian Fernandez-Alberti, Vladimir Chernyak
1Theoretical Division, Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
Abstract:
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 compared. In addition to influencing the speed of the simulation, the chosen parameters combined with the frequently reduced number of trajectories can sometimes lead to unanticipated changes in the accuracy of the simulated dynamics. We have previously developed a nonadiabatic excited state molecular dynamics methodology employing Tully's fewest switches surface hopping algorithm. In this study, we seek to investigate the impact of the number of trajectories and the various parameters on the simulation of the photoinduced dynamics of distyrylbenzene (a small oligomer of polyphenylene vinylene) within our developed framework. Various user-defined parameters are analyzed: classical and quantum integration time steps, the value of the friction coefficient for Langevin dynamics, and the initial seed used for stochastic thermostat and hopping algorithms. Common approximations such as reduced number of nonadiabatic coupling terms and the classical path approximation are also investigated. Our analysis shows that, at least for the considered molecular system, a minimum of ~400 independent trajectories should be calculated in order to achieve statistical averaging necessary for convergence of the calculated relaxation timescales.
More Related Videos
Related Concept Videos
¹H NMR of Conformationally Flexible Molecules: Temporal Resolution
¹H NMR of Conformationally Flexible Molecules: Variable-Temperature NMR
Reaction Mechanisms: The Steady-State Approximation
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
Adiabatic Processes for an Ideal Gas
Real Gases: Effects of Intermolecular Forces and Molecular Volume Deriving Van der Waals Equation

