Extended Lagrangian molecular dynamics on vibronic surfaces in the nuclear-electronic orbital framework
Joseph A Dickinson1,2, Mathew Chow1,2, Eno Paenurk2
1Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA.
The Journal of Chemical Physics
|June 15, 2026
Summary
We developed new methods for simulating proton transfer dynamics using the nuclear-electronic orbital (NEO) framework. These techniques accurately capture quantum effects, improving simulations of chemical reactions involving proton movement.
Area of Science:
- Quantum Chemistry
- Chemical Dynamics
- Computational Chemistry
Background:
- Proton transfer is crucial in many chemical processes.
- Simulating proton transfer requires accounting for nuclear quantum effects like tunneling and zero-point energy.
- Existing methods may not fully capture these quantum mechanical aspects.
Purpose of the Study:
- Introduce novel simulation methods within the nuclear-electronic orbital (NEO) framework for proton transfer dynamics.
- Incorporate nuclear quantum effects, such as quantum proton treatment, into molecular dynamics simulations.
- Enhance the efficiency and accuracy of simulating complex proton transfer systems.
Main Methods:
- Utilize nuclear-electronic orbital density functional theory (NEO-DFT) to treat protons quantum mechanically.
- Formulate a NEO extended Lagrangian molecular dynamics (NEO-ELMD) approach.
- Implement density matrix extrapolation and purification to accelerate the NEO self-consistent field procedure.
Main Results:
- Demonstrated the fidelity and efficiency of the NEO-ELMD method for intramolecular proton transfer in malonaldehyde.
- Successfully simulated nonequilibrium single and double proton transfer dynamics in larger benzimidazole-phenol systems.
- Validated the accelerated techniques for efficient proton transfer simulations.
Conclusions:
- The developed NEO-ELMD approach provides a robust foundation for simulating proton transfer dynamics.
- The accelerated techniques significantly improve the efficiency of NEO-DFT calculations.
- This work enables future methodologies for simulating proton transfer with nonadiabatic effects.
Related Concept Videos
Molecular Orbital Theory II
Molecular Orbital Energy Diagrams
Molecular Orbital Theory I
Overview of Molecular Orbital Theory
Hybridization of Atomic Orbitals II
sp3d and sp3d 2 Hybridization
Valence Bond Theory and Hybridized Orbitals
According to valence bond theory, a covalent bond results when: (1) an orbital on one atom overlaps an orbital on a second atom, and (2) the single electrons in each orbital combine to form an electron pair. The strength of a covalent bond depends on the extent of overlap of the orbitals involved. Maximum overlap is possible when the orbitals overlap on a direct line between the two nuclei.
A σ bond (single bond in a Lewis structure) is a covalent bond in which the electron density is...
A σ bond (single bond in a Lewis structure) is a covalent bond in which the electron density is...
Atomic Orbitals
An atomic orbital represents the three-dimensional regions in an atom where an electron has the highest probability to reside. The radial distribution function indicates the total probability of finding an electron within the thin shell at a distance r from the nucleus. The atomic orbitals have distinct shapes which are determined by l, the angular momentum quantum number. The orbitals are often drawn with a boundary surface, enclosing densest regions of the cloud.
Hybridization of Atomic Orbitals I
The mathematical expression known as the wave function, ψ, contains information about each orbital and the wavelike properties of electrons in an isolated atom. When atoms are bound together in a molecule, the wave functions combine to produce new mathematical descriptions that have different shapes. This process of combining the wave functions for atomic orbitals is called hybridization and is mathematically accomplished by the linear combination of atomic orbitals. The new orbitals that...


