Capturing anharmonic effects in single vibronic level fluorescence spectra using local harmonic Hagedorn wavepacket dynamics

Basic Information

• Paper ID: 2408.11991
• Title: Capturing anharmonic effects in single vibronic level fluorescence spectra using local harmonic Hagedorn wavepacket dynamics
• Authors: Zhan Tong Zhang, Máté Visegrádi, Jiří J. L. Vaníček
• Classification: physics.chem-ph quant-ph
• Publication Date: October 14, 2025
• Paper Link: https://arxiv.org/abs/2408.11991

Abstract

This paper proposes a method combining Hagedorn wavepacket dynamics with local harmonic approximation to calculate anharmonic effects in single vibronic level (SVL) fluorescence spectra. While Hagedorn wavepacket dynamics yields exact SVL spectra in the global harmonic model, anharmonic effects in real molecules significantly influence the spectra. The authors validate the method's accuracy through exact quantum benchmarks on two-dimensional Morse-type potentials and apply it to spectral calculations of anthracene and difluorocarbene molecules, demonstrating the method's effectiveness in handling anharmonic effects.

Research Background and Motivation

1. Core Problem: In single vibronic level fluorescence spectroscopy, non-Gaussian wavepacket propagation is required. While traditional Gaussian wavepacket dynamics methods are exact under harmonic potentials, they cannot handle anharmonic effects in real molecules.
2. Problem Significance:
   • SVL fluorescence spectroscopy is an important tool for studying intramolecular relaxation, vibrational structure characterization, and conformational identification
   • Accurate calculation of vibrationally resolved electronic spectra of polyatomic molecules typically requires consideration of anharmonicity in molecular potential energy surfaces
   • For flexible molecules with large-amplitude anharmonic vibrations, quantum treatment is often necessary
3. Limitations of Existing Methods:
   • Traditional time-independent methods for calculating Franck-Condon factors are computationally prohibitive in high-dimensional systems
   • Global harmonic approximations cannot capture anharmonic effects
   • Existing semiclassical trajectory methods, while capable of handling anharmonicity, struggle with arbitrary vibrational excited states
4. Research Motivation: Develop an efficient method capable of calculating SVL spectra for arbitrary vibrational levels from a single Gaussian trajectory while partially incorporating anharmonic effects.

Core Contributions

1. Theoretical Method: Proposes local harmonic Hagedorn wavepacket dynamics, combining the Hagedorn wavepacket method with local harmonic approximation
2. Algorithmic Advantage: Enables calculation of SVL spectra for all initial vibrational levels from a single Gaussian trajectory
3. Accuracy Verification: Validates method accuracy through exact quantum benchmarks on two-dimensional coupled Morse potentials
4. Practical Application: Successfully applies the method to ab initio calculations on anthracene and difluorocarbene molecules, demonstrating practical utility
5. Theoretical Analysis: Provides analytical explanation for spectral envelope splitting in difluorocarbene

Methodology Details

Task Definition

Calculate the emission rate of single vibronic level fluorescence spectra: $\sigma_{em}(\omega) = \frac{4\omega^3}{3\pi\hbar c^3}|\mu_{ge}|^2 \text{Re}\int_0^{\infty} C(t) \exp[it(\omega - \omega_{e,K})] dt$

where the autocorrelation function is: $C(t) = \langle K| \exp(-iH_gt/\hbar)|K\rangle$

Model Architecture

Hagedorn Function Construction

The initial vibrational wavepacket is represented as a Hagedorn function: $\phi_K = (K!)^{-1/2}(A^{\dagger})^K\phi_0$

where the raising operator is defined as: $A^{\dagger} := \frac{i}{\sqrt{2\hbar}}(P_t^{\dagger} \cdot (\hat{q} - q_t) - Q_t^{\dagger} \cdot (\hat{p} - p_t))$

Local Harmonic Approximation

The real potential V is replaced by an effective quadratic potential: $V_{LHA}(q; q_t) := V(q_t) + V'(q_t) \cdot x + x^T \cdot V''(q_t) \cdot x/2$

Equations of Motion

Wavepacket parameters evolve according to: $\dot{q}_t = m^{-1} \cdot p_t, \quad \dot{p}_t = -V'(q_t)$$\dot{Q}_t = m^{-1} \cdot P_t, \quad \dot{P}_t = -V''(q_t) \cdot Q_t$

Technical Innovations

1. Invariant Coefficients: Under local harmonic approximation, coefficients $c_K$ in Hagedorn wavepackets remain time-independent
2. Single-Trajectory Calculation: SVL spectra for all vibrational levels can be obtained from a single Gaussian trajectory
3. Nonlinear Schrödinger Equation: Hagedorn functions are exact solutions of the nonlinear time-dependent Schrödinger equation
4. Recursive Overlap Integrals: Employs exact recursive algebraic expressions for calculating Hagedorn function overlaps

Experimental Setup

Test Systems

1. Two-Dimensional Coupled Morse Potential:
   • Parameters: $\omega_g = (0.0041, 0.005)$, $\chi = (0.005, 0.002)$
   • Coupling term: $d' = 0.08$, $\chi' = (0.001, 0.001)$
   • Used for exact benchmarking without electronic structure errors
2. Anthracene Molecule:
   • 66 vibrational degrees of freedom
   • Calculated at PBE0/def2-TZVP level
   • $^1B_{2u} \rightarrow ^1A_g$ transition
3. Difluorocarbene (CF₂):
   • Flexible molecule with 3 vibrational modes
   • Calculated at PBE0/aug-cc-pVTZ level
   • Bending mode dominates spectral features

Comparison Methods

• Vertical Harmonic Approximation (VHA)
• Adiabatic Harmonic Approximation (AHA)
• Exact Quantum Calculation (Morse potential only)

Computational Details

• Propagation time: 80000 a.u. (Morse potential), 80000 a.u. (molecular systems)
• Time step: 4 a.u. (Morse potential), 8 a.u. (molecular systems)
• Integrator: Second-order TVT geometric integrator

Experimental Results

Main Results

Two-Dimensional Morse Potential System

• Ground State Emission (10²0): Local harmonic and adiabatic harmonic models show excellent agreement with exact quantum results
• Excited State Emission: Local harmonic method performs better in low-frequency regions (corresponding to high vibrational level transitions)
• Peak Position Accuracy: Vertical harmonic model shows poorest performance in peak positions

Anthracene Molecule

• Local harmonic spectra highly similar to adiabatic harmonic approximation results
• Vertical harmonic spectra show significant differences
• Anharmonic effects are relatively small in anthracene

Difluorocarbene

• Local harmonic method shows excellent agreement with experimental results
• Global harmonic approximation deviates severely quantitatively
• Large displacement of bending mode dominates spectral features

Ablation Studies

Analytical verification through one-dimensional displaced harmonic oscillator model confirms difluorocarbene spectral envelope splitting phenomenon:

• Huang-Rhys factor: S = 6.15
• Transition intensities follow Poisson-modulated associated Laguerre polynomials
• Time-dependent Hagedorn method completely consistent with analytical formula

Experimental Findings

1. Accuracy Hierarchy: Local harmonic > Adiabatic harmonic > Vertical harmonic
2. Excited State Effects: High vibrational excited states are more sensitive to anharmonicity
3. Molecular Dependence: Importance of anharmonicity varies by molecule (CF₂ >> anthracene)
4. Artificial Peaks: Local harmonic method may produce non-physical negative intensity peaks

Related Work

Main Research Directions

1. Time-Independent Methods: Franck-Condon factor calculations based on variational principles or perturbation theory
2. Time-Dependent Methods: Spectral calculations based on wavepacket propagation
3. Semiclassical Methods: Trajectory methods such as thawed Gaussian approximation (TGA)

Advantages of This Work

• Handles arbitrary vibrational excited states (not limited to single excitations)
• Calculates all spectra from single trajectory
• Partially incorporates anharmonic effects
• Suitable for combination with ab initio dynamics

Conclusions and Discussion

Main Conclusions

1. Local harmonic Hagedorn method effectively captures anharmonic effects in SVL spectra
2. For moderately anharmonic systems, this method is significantly more accurate than global harmonic approximation
3. Method importance varies with molecular flexibility
4. Single-trajectory property makes it particularly suitable for expensive ab initio calculations

Limitations

1. Assumes excited state potential energy surface is harmonic
2. Accuracy decreases at high vibrational excitations
3. May produce non-physical peaks
4. Requires Hessian matrix calculations, which are computationally expensive

Future Directions

1. Project anharmonic excited states onto Hagedorn basis
2. Apply Hessian interpolation techniques to improve efficiency
3. Extend to other spectroscopic techniques requiring non-Gaussian wavepackets
4. Develop higher-accuracy methods for anharmonic treatment

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Solid mathematical foundation based on Hagedorn functions
2. Computational Efficiency: Single-trajectory calculation for all vibrational level spectra
3. Thorough Validation: Verified through exact quantum benchmarks and experimental comparisons
4. Strong Practicality: Suitable for combination with ab initio methods
5. Detailed Analysis: Provides physical interpretation of spectral features

Weaknesses

1. Applicable Range: Primarily suitable for moderately anharmonic systems
2. Artificial Peak Problem: Local harmonic approximation may produce non-physical results
3. Computational Cost: Requires Hessian calculation at each time step
4. Excited State Limitation: Assumption of harmonic excited states may be insufficiently accurate

Impact

1. Academic Value: Provides new theoretical framework for SVL spectral calculations
2. Practical Value: Particularly suitable for spectral prediction of flexible molecules
3. Methodological Contribution: Demonstrates potential of local harmonic approximation in quantum dynamics

Applicable Scenarios

• Vibrationally resolved fluorescence spectral calculations of flexible molecules
• Spectral predictions requiring consideration of anharmonic effects
• Quantum dynamics simulations of high-dimensional systems
• On-the-fly dynamics combined with ab initio electronic structure calculations

References

This paper cites over 70 relevant references, primarily including:

• Theoretical foundations of Hagedorn wavepacket dynamics
• Experimental and theoretical studies of SVL spectroscopy
• Quantum chemistry methods for handling anharmonic effects
• Classical literature on time-dependent spectroscopy