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.
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: