NQCDynamics.jl: A Julia Package for Nonadiabatic Quantum Classical Molecular Dynamics in the Condensed Phase

Abstract

Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in condensed phase. To facilitate effective development, code sharing and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package which provides a framework for established and emerging methods for performing semiclassical and mixed quantum-classical dynamics in condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semi-classical and mixed quantum-classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the Molecular Dynamics with Electronic Friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems.

19 pages, 8 figures; updated figures, added scaling data, added model E results

References in corpus (12)

• Efficient stochastic thermostatting of path integral molecular dynamics
• Generalized spin mapping for quantum-classical dynamics
• Ab-initio tensorial electronic friction for molecules on metal surfaces: nonadiabatic vibrational relaxation
• Mode specific electronic friction in dissociative chemisorption on metal surfaces: H${}_2$ on Ag(111)
• Unified Formulation of Phase Space Mapping Approaches for Nonadiabatic Quantum Dynamics
• Symmetry-Adapted High Dimensional Neural Network Representation of Electronic Friction Tensor of Adsorbates on Metals
• Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation
• Negative Zero-Point-Energy Parameter in the Meyer-Miller Mapping Model for Nonadiabatic Dynamics
• Commutator Matrix in Phase Space Mapping Models for Nonadiabatic Quantum Dynamics
• State Dependent Ring Polymer Molecular Dynamics for Investigating Excited Nonadiabatic Dynamics
• Electronic friction coefficients from the atom-in-jellium model for $Z=1-92$
• High-order geometric integrators for representation-free Ehrenfest dynamics