Shadow excited state molecular dynamics with the ΔSCF method
O Jonathan Fajen1,2,3, Oscar Grånäs4, Todd J Martínez1,2
1Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, USA.
The Journal of Chemical Physics
|March 5, 2026
Summary
We developed a new shadow excited-state molecular dynamics (ESMD) method using density functional theory. This approach improves computational efficiency, stability, and robustness for excited state dynamics simulations.
Area of Science:
- Computational Chemistry
- Quantum Chemistry
- Theoretical Chemistry
Background:
- Born-Oppenheimer molecular dynamics (MD) is crucial for simulating molecular behavior.
- Simulating excited-state molecular dynamics (ESMD) presents significant computational challenges.
- Existing methods often face limitations in accuracy, stability, or computational cost.
Purpose of the Study:
- To extend the shadow extended Lagrangian Born-Oppenheimer molecular dynamics (ESMD) method for excited state dynamics.
- To develop a more computationally efficient and robust approach for ESMD.
- To demonstrate the method's efficacy using self-consistent charge density functional tight binding (SCC-DFTB) theory.
Main Methods:
- Introduced a shadow ESMD approach replacing iterative solutions with a single-step solution of an approximate shadow excited-state potential.
- Linearized the energy functional about an approximate excited state density.
- Propagated the approximate excited-state (charge) density as an additional dynamical variable within an extended Lagrangian framework.
Main Results:
- The shadow ESMD method offers significant improvements in computational cost compared to direct ESMD.
- The new method demonstrates enhanced stability and robustness over its "exact" counterpart.
- The implementation in SCC-DFTB theory shows broad generalizability to other electronic structure methods.
Conclusions:
- The shadow ESMD method provides a computationally efficient, stable, and robust alternative for excited-state molecular dynamics.
- This approach is adaptable to various quantum chemistry methods, including ab initio techniques.
- The findings pave the way for more accurate and accessible simulations of excited-state phenomena.
Related Concept Videos
Molecular Geometry and Dipole Moments
13.6K
The VSEPR theory can be used to determine the electron pair geometries and molecular structures as follows:
13.6K
Bond Dissociation Energy and Activation Energy
8.6K
Bond energy is the energy required to break a bond homolytically. These values are usually expressed in units of kcal/mol or kJ/mol and are referred to as bond dissociation energies when given for specific bonds or average bond energies when indicated for a given type of bond over many compounds. Firstly, the bond dissociation energy for a single bond is weaker than that of a double bond, which in turn is weaker than that of a triple bond. Secondly, hydrogen forms relatively strong bonds with...
8.6K
π Electron Effects on Chemical Shift: Overview
1.5K
An applied magnetic field causes loosely bound π-electrons in organic molecules to circulate, producing a local or induced diamagnetic field over a large spatial volume. As the molecules tumble in solution, the field generated by π-electrons in spherical substituents results in a zero net field. However, the net field generated by π-electrons in non-spherical substituents is not zero. The effect of this induced field depends on the orientation of the molecule with respect to B0,...
1.5K
Atomic Force Microscopy
3.1K
Atomic force microscopy (AFM) is a type of scanning probe microscopy that can analyze topographic details of various specimens like ceramics, glass, polymers, and biological samples. AFM offers over 1000 times more resolution than the optical imaging system. Images generated from AFM are three-dimensional surface profiles, offering an advantage over the flat, two-dimensional images from other imaging techniques.
The AFM Probe
The probe is regarded as the heart of any AFM setup and comprises the...
The AFM Probe
The probe is regarded as the heart of any AFM setup and comprises the...
3.1K


