Marcus Theory and The Condon Approximation Revisited I: E-SHAKE and Seam Sampling

Basic Information

• Paper ID: 2510.11810
• Title: Marcus Theory and The Condon Approximation Revisited I: E-SHAKE and Seam Sampling
• Authors: D. Vale Cofer-Shabica, Jennifer R. DeRosa, Joseph E. Subotnik
• Classification: physics.chem-ph
• Publication Date: October 13, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.11810v1

Abstract

Marcus theory is an essential tool in theoretical chemistry for predicting charge and energy transfer rates. For the naphthalene-bridged-biphenyl and naphthalene-bridged-benzophenone systems studied by Piotrowiak, Miller, and Closs, Marcus theory shows excellent agreement with experimental results for electron transfer and triplet energy transfer. However, discrepancies between theory and experiment persist. This paper revisits a critical inconsistency: the molecule C-13-ae (3,equatorial-naphthalene-cyclohexane-1,axial-benzophenone). To better understand the theory-experiment discrepancy, the authors introduce and employ a novel sampling scheme to sample the seam between two nonadiabatic electronic states (E-SHAKE). Through this method, they reveal the failure of the Condon approximation and the presence of conical intersections in the C-13-ae molecule, and predict isotope effects on triplet-triplet energy transfer rates.

Research Background and Motivation

Problem Background

1. Importance of Marcus Theory: Marcus theory is the core theoretical tool for predicting electron transfer (ET) and triplet energy transfer (TET) rates, expressed as: $k_{D→A} = \frac{2π}{\hbar}|H_{DA}|^2(4πλk_BT)^{-1/2}e^{-(λ+ΔG^0)^2/(4λk_BT)}$
2. Theoretical Assumptions: Marcus theory is based on several key assumptions:
   • Thermalized harmonic oscillator bath
   • Activated crossing limit
   • Applicability of the golden rule
   • Condon Approximation: Electronic coupling treated as a constant
3. Anomalies in Experimental Observations: In the Closs molecule series, the 'ee' (equatorial-equatorial) system exhibits perfect exponential decay relationships, while 'ae' (axial-equatorial) and 'ea' (equatorial-axial) systems display more complex behavior.

Research Motivation

1. Theory-Experiment Discrepancy: Particularly significant differences between experimental rates and theoretical predictions for the C-13-ae molecule
2. Validity of the Condon Approximation: Need to verify whether the Condon approximation remains applicable in "flexible" molecules
3. Molecular Rigidity vs. Flexibility: Explore how molecular flexibility affects energy transfer processes

Core Contributions

1. Development of the E-SHAKE Method: A novel electronic constrained molecular dynamics sampling technique capable of sampling seams between nonadiabatic electronic states
2. Revelation of Condon Approximation Failure: Discovery of dramatic variations in electronic coupling in the C-13-ae molecule through seam sampling
3. Identification of Conical Intersections: Recognition of energetically accessible conical intersection points in the C-13-ae molecule
4. Prediction of Isotope Effects: Prediction that deuteration will produce significant effects on triplet energy transfer rates
5. Provision of New Theoretical Insights: Explanation of why certain Closs molecules deviate from Marcus theory predictions

Methodology Details

Construction of Localized Nonadiabatic States

Nonadiabatic-Adiabatic Transformation

Construction of nonadiabatic states through rotation of adiabatic states: $|Ξ_i⟩ = \sum_{j=1}^{N_{states}} (U)_{ji} |Φ_j⟩$

Boys and BoysOV Methods

• Boys Criterion: Maximization of dipole moment differences $f_{Boys}(U) = \sum_{i,j=1}^{N_{states}} |⟨Ξ_i|\hat{μ}|Ξ_i⟩ - ⟨Ξ_j|\hat{μ}|Ξ_j⟩|^2$
• BoysOV Method: Separate localization of occupied and virtual orbitals $f_{BoysOV}(U) = \sum_{i,j=1}^{N_{states}} |⟨Ξ_i|\hat{μ}_{occ}|Ξ_i⟩ - ⟨Ξ_j|\hat{μ}_{occ}|Ξ_j⟩|^2 + |⟨Ξ_i|\hat{μ}_{virt}|Ξ_i⟩ - ⟨Ξ_j|\hat{μ}_{virt}|Ξ_j⟩|^2$

E-SHAKE Algorithm

Core Concept

E-SHAKE is the electronic version of the SHAKE/RATTLE algorithm, used to impose electronic constraints in molecular dynamics.

Constraint Conditions

• Position Constraint: $σ(t) = E_D(R(t)) - E_A(R(t)) = 0$
• Velocity Constraint: $\dot{σ}(t) = \dot{R}(t) · ∇σ|_{R(t)} = 0$

Algorithm Implementation

Using the constrained version of the velocity Verlet integrator: $R(t+h) = R(t) + hQ$$\dot{R}(t+h) = Q + \frac{1}{2}M^{-1} · [hF(t+h) + kG(t+h)]$

where Lagrange multipliers are determined numerically to satisfy constraint conditions.

Technical Innovations

1. Implementation of Electronic Constraints: First extension of the SHAKE algorithm to electronic degrees of freedom
2. Seam Space Sampling: Systematic exploration of the entire seam space of nonadiabatic state crossings
3. Real-Time Coupling Monitoring: Real-time monitoring of electronic coupling changes during dynamics
4. Generality: Applicable to arbitrary electronic constraints, not limited to energy difference constraints

Experimental Setup

Computational Details

• Electronic Structure Method: CIS/6-31G* basis set, using Q-Chem 6.0
• Constraint Precision: Energy difference constraint precision of 4×10^-5 a.u. (≈1 meV)
• Temperature Control: Random corrected Berendsen method with time constant τ = 10 fs
• Integration Time Step: Standard molecular dynamics time step

Molecular Systems

Study objects from the Closs molecule series:

• C-13-ee, C-13-ea, C-13-ae
• C-14-ee, C-14-ea, C-14-ae
• D-26, D-27 (decahydronaphthalene-bridged systems)

Analysis Methods

1. Seam Energy Analysis: Analysis of coupling variations in high-energy regions with the lowest T1 energy as zero point
2. Atomic Contribution Decomposition: Determination of atoms dominating coupling changes through gradient analysis
3. Attachment-Detachment Density Analysis: Analysis of changes in donor-acceptor wavefunction overlap

Experimental Results

Major Findings

1. Significant Differences in Electronic Coupling

• C-13-ae vs C-13-ea: Electronic coupling in C-13-ae is two orders of magnitude smaller than in C-13-ea
• Seam Sampling Results: Coupling in C-13-ae spans two orders of magnitude variation across seam space

2. Failure of the Condon Approximation

• C-13-ee and C-13-ea: Relatively stable coupling, consistent with Condon approximation
• C-13-ae: Coupling distributed around zero, severe violation of Condon approximation

3. Identification of Conical Intersections

When $H_{DA} → 0$ and $E_D → E_A$, the system transitions from nonadiabatic seam to conical intersection: $\lim_{H_{DA}→0, E_D→E_A=E} \begin{pmatrix} E_D & H_{DA} \\ H_{DA} & E_A \end{pmatrix} = \begin{pmatrix} E & 0 \\ 0 & E \end{pmatrix}$

4. Effects of Molecular Flexibility

• RMSD Analysis: RMSD of C-13-ea (0.13 Å) is twice that of C-13-ee (0.06 Å)
• Barrier Analysis: Similar ground state isomerization barriers for all systems (9-11 kcal/mol)

Isotope Effect Predictions

Atomic contribution analysis reveals that motion of equatorial hydrogen atoms (He3) dominates coupling changes. Prediction that replacement of this hydrogen atom with deuterium will produce significant isotope effects, particularly for the C-13-ae system.

Related Work

Development of Marcus Theory

1. Classical Marcus Theory: Based on harmonic oscillator model and Condon approximation
2. Efforts Beyond Condon Approximation:
   • Stuchebrukhov's two-mode model approach
   • Jang and Newton's non-Condon Hamiltonian
3. Modern Nonadiabatic Dynamics Methods: FSSH, ab initio multiple spawning, etc.

Related Seam Sampling Methods

1. Biased Potential Methods: Using artificial potentials to modify electronic Hamiltonian
2. Enhanced Sampling: Constrained extensions of nudged elastic band methods
3. Advantages of This Work: Generic and exhaustive exploration of seam space

Conclusions and Discussion

Main Conclusions

1. Anomalous Behavior of C-13-ae Molecule: Due to presence of conical intersections and failure of Condon approximation
2. Limitations of Marcus Theory Applicability: Inapplicable in flexible molecules and near conical intersections
3. Experimental Predictions: C-13-ae should exhibit significantly slower TET rates than C-13-ea

Limitations

1. Computational Method Constraints: Use of CIS method and moderate basis set may affect quantitative accuracy
2. Approximate Gradients: Gradients constructed using strict nonadiabatic approximation
3. Temperature Effects: Seam sampling temperature does not represent physical temperature

Future Directions

1. Extension to Solution Systems: Combining solvent effects to study seam dynamics
2. More Accurate Electronic Structure Methods: Using higher-level quantum chemistry methods
3. Experimental Verification: Verification of theoretical predictions through ultrafast spectroscopy
4. Isotope Effect Experiments: Verification of predicted deuteration effects

In-Depth Evaluation

Strengths

1. Methodological Innovation: E-SHAKE algorithm provides a new tool for studying nonadiabatic processes
2. Deep Theoretical Insights: Reveals microscopic mechanisms of Marcus theory failure
3. Concrete Experimental Predictions: Proposes verifiable isotope effect predictions
4. Computational Efficiency: Moderate computational cost compared to full quantum dynamics calculations

Limitations

1. Computational Accuracy: CIS method may lack sufficient precision for quantitative predictions
2. System Limitations: Primarily focused on specific Closs molecules; generalizability remains to be verified
3. Neglect of Solvent Effects: Actual experiments typically occur in solution
4. Statistical Sampling: Statistical sufficiency of seam sampling requires further verification

Impact

1. Theoretical Contribution: Provides new theoretical framework for understanding nonadiabatic processes
2. Methodological Value: E-SHAKE method widely applicable to other systems
3. Experimental Guidance: Provides new perspective for reinterpreting historical experimental data
4. Educational Value: Effectively demonstrates limitations of approximation methods in theoretical chemistry

Applicable Scenarios

1. Energy Transfer in Flexible Molecules: Particularly suitable for studying effects of conformational changes on electronic coupling
2. Conical Intersection Research: Provides new method for systematic exploration of conical intersection seams
3. Isotope Effect Prediction: Can be used to design and interpret isotope-labeled experiments
4. Theoretical Method Verification: Can serve as benchmark for validating other nonadiabatic dynamics methods

References

The paper cites multiple classical references in theoretical chemistry and molecular dynamics, including:

1. Nittan, A. Chemical dynamics in condensed phases (classic textbook)
2. Pioneering experimental work by Closs et al. on Marcus theory verification
3. Tully's surface hopping method
4. Development history of various nonadiabatic dynamics methods

Overall Assessment: This is a high-quality theoretical chemistry paper that deeply investigates the applicability boundaries of Marcus theory through development of new computational methods. The paper not only provides important theoretical insights but also offers concrete experimental predictions, exemplifying the research paradigm combining theory and experiment. The development of the E-SHAKE method opens new avenues for nonadiabatic dynamics research and possesses significant methodological value.