GPU-Accelerated Large-Scale Nonadiabatic Excited State Molecular Dynamics
- 13 November 2026
- journal article
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
Abstract
Nonadiabatic molecular dynamics is essential for simulating photoinduced processes, but its application to large molecular systems is often limited by the cost of repeatedly evaluating excited-state energies, forces, and nonadiabatic couplings. Here, we present a GPU-accelerated implementation of trajectory surface hopping within PYSEQM (https://github.com/lanl/PYSEQM/), a PyTorch-based semiempirical quantum chemistry framework. The module combines semiempirical model Hamiltonians with configuration-interaction singles excited states and introduces the nonadiabatic coupling vectors and time-derivative couplings required for surface-hopping dynamics. The implementation further includes treatments of trivial crossings and electronic decoherence, providing a complete framework for propagating coupled electronic and nuclear dynamics. To improve GPU utilization, independent trajectories are propagated concurrently in batches. Benchmark calculations show that this strategy substantially increases aggregate simulation throughput as the batch size is increased for small molecules. We demonstrate the applicability of the method to large molecular systems by simulating the excited-state population dynamics of methyl- and tert-butyl-substituted molecular nanostar dendrimers, which contain dense manifolds of electronically excited states. These calculations show that semiempirical surface-hopping dynamics can be extended to molecular systems that remain challenging for conventional ab initio approaches. In addition, the differentiable PyTorch formulation provides a data-science route for optimizing semiempirical Hamiltonian parameters against reference excited-state properties and, ultimately, dynamical observables. PYSEQM therefore provides a platform for high-throughput nonadiabatic molecular dynamics and for the development of data-driven semiempirical models for photophysics and photochemistry.All Related Versions
Funding Information
- Center for Integrated Nanotechnologies
- Los Alamos National Laboratory
- Basic Energy Sciences (FWP LANLE8B3)
- National Nuclear Security Administration (892333218NCA000001)
This publication has 72 references indexed in Scilit:
- Dynamics of Excitons in Conjugated Molecules and Organic Semiconductor SystemsChemical Reviews, 2022
- Surface hopping methods for nonadiabatic dynamics in extended systemsWIREs Computational Molecular Science, 2019
- Semiempirical molecular orbital models based on the neglect of diatomic differential overlap approximationInternational Journal of Quantum Chemistry, 2018
- Recent Advances and Perspectives on Nonadiabatic Mixed Quantum–Classical DynamicsChemical Reviews, 2018
- Ab Initio Nonadiabatic Quantum Molecular DynamicsChemical Reviews, 2018
- Ultrafast Elementary Photochemical Processes of Organic Molecules in Liquid SolutionChemical Reviews, 2016
- Implementation of surface hopping molecular dynamics using semiempirical methodsChemical Physics, 2008
- Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elementsJournal of Molecular Modeling, 2007
- Optimization of parameters for semiempirical methods I. MethodJournal of Computational Chemistry, 1989
- Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanesJournal of Computational Physics, 1977