Assessing the performance of coupled-trajectory schemes on full-dimensional two-state linear vibronic coupling models
Peter Schürger1, Lea M Ibele1, David Lauvergnat1
1CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France.
The Journal of Chemical Physics
|March 14, 2025
Summary
Coupled-trajectory methods improve nonadiabatic molecular dynamics simulations. The combined coupled-trajectory Tully surface hopping (CCTTSH) algorithm resolves inconsistencies, but coupled-trajectory mixed quantum-classical (CTMQC) struggles with strong nonadiabaticity.
Area of Science:
- Computational chemistry
- Quantum dynamics
- Photochemistry
Background:
- Nonadiabatic molecular dynamics simulations are crucial for understanding photochemical processes.
- Coupled-trajectory methods offer a potential improvement over independent-trajectory approaches.
- Accurate simulation of excited-state dynamics requires robust theoretical frameworks.
Purpose of the Study:
- To evaluate the performance of coupled-trajectory methods for simulating the photodynamics of DMABN and fulvene.
- To compare the coupled-trajectory Tully surface hopping (CCTTSH) and coupled-trajectory mixed quantum-classical (CTMQC) algorithms against independent-trajectory methods.
- To identify strengths and weaknesses of these methods in regions of strong nonadiabaticity.
Main Methods:
- Simulations of 4-(dimethylamino)benzonitrile (DMABN) and fulvene photodynamics.
- Application of linear vibrational coupling models for electronic structure.
- Comparison of coupled-trajectory mixed quantum-classical (CTMQC) and combined coupled-trajectory Tully surface hopping (CCTTSH) algorithms.
- Analysis of electronic populations and nuclear properties in position and momentum space.
Main Results:
- The CCTTSH algorithm successfully addresses internal inconsistencies found in previous coupled-trajectory Tully surface hopping methods.
- DMABN simulations reveal a significant limitation of the CTMQC algorithm when trajectories linger in strongly nonadiabatic regions.
- Both DMABN and fulvene exhibit distinct behaviors under the tested simulation methods.
Conclusions:
- The CCTTSH algorithm represents a significant advancement for coupled-trajectory surface hopping methods.
- The CTMQC method shows limitations in simulating systems with prolonged interactions in regions of strong nonadiabaticity.
- Careful selection of simulation methods is essential for accurate photodynamic studies, particularly in complex systems.
More Related Videos
Related Concept Videos
Spin–Spin Coupling: Three-Bond Coupling (Vicinal Coupling)
1.0K
Vicinal or three-bond coupling is commonly observed between protons attached to adjacent carbons. Here, nuclear spin information is primarily transferred via electron spin interactions between adjacent C‑H bond orbitals. This generally favors the antiparallel arrangement of spins, so 3J values are usually positive.
The extent of coupling depends on the C‑C bond length, the two H‑C‑C angles, any electron-withdrawing substituents, and the dihedral angle between the...
The extent of coupling depends on the C‑C bond length, the two H‑C‑C angles, any electron-withdrawing substituents, and the dihedral angle between the...
1.0K
Spin–Spin Coupling: Two-Bond Coupling (Geminal Coupling)
936
Two NMR-active nuclei bonded to a central atom can be involved in geminal or two-bond coupling. Geminal coupling is commonly seen between diastereotopic protons in chiral molecules and unsymmetrical alkenes, among others.
The central atom need not be NMR-active because its electrons are affected by the electron polarization of the spin-active atoms. However, spin information is transmitted less effectively than in one-bond coupling, and 2J values are usually weaker than 1J values. The energy of...
The central atom need not be NMR-active because its electrons are affected by the electron polarization of the spin-active atoms. However, spin information is transmitted less effectively than in one-bond coupling, and 2J values are usually weaker than 1J values. The energy of...
936
Spin–Spin Coupling Constant: Overview
862
In bromoethane, the three methyl protons are coupled to the two methylene protons that are three bonds away. In accordance with the n+1 rule, the signal from the methyl protons is split into three peaks with 1:2:1 relative intensities. The methylene protons appear as a quartet, with the relative intensities of 1:3:3:1.
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must...
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must...
862
Spin–Spin Coupling: One-Bond Coupling
918
Coupling interactions are strongest between NMR-active nuclei bonded to each other, where spin information can be transmitted directly through the pair of bonding electrons. While nuclei polarize their electrons to the opposite spins, the bonding electron pair has opposite spins. Configurations with antiparallel nuclear spins are expected to be lower in energy. When coupling makes antiparallel states more favorable, J is considered to have a positive value. The one-bond coupling constant, 1J,...
918
¹H NMR: Long-Range Coupling
1.6K
The coupling interactions of nuclei across four or more bonds are usually weak, with J values less than 1 Hz. While these are usually not observed in spectra, the presence of multiple bonds along the coupling pathway can result in observable long-range coupling.
In alkenes, spin information is communicated via σ–π overlap, as seen in allylic (four-bond) and homoallylic (five-bond) couplings. These coupling interactions are stronger when the σ bond is parallel to the alkene...
In alkenes, spin information is communicated via σ–π overlap, as seen in allylic (four-bond) and homoallylic (five-bond) couplings. These coupling interactions are stronger when the σ bond is parallel to the alkene...
1.6K
Hybridization of Atomic Orbitals II
31.6K
sp3d and sp3d 2 Hybridization
31.6K


