Capturing anharmonic effects in single vibronic level fluorescence spectra using local harmonic Hagedorn wavepacket dynamics
Zhang, Visegrádi, VanÃÄek
Hagedorn wavepacket dynamics yields exact single vibronic level (SVL) fluorescence spectra in global harmonic models. To partially describe the effects of anharmonicity, important in the spectra of real molecules, we describe a combination of the Hagedorn wavepacket approach to SVL spectroscopy with the local harmonic approximation. In a proof-of-principle study [Phys. Rev. A 111, L010801 (2025)], we successfully demonstrated the utility of this method by computing the SVL spectra of difluorocarbene, a floppy molecule with moderately anharmonic potential. Here, we describe the theory in detail and analyse the method more thoroughly. To assess the accuracy of the method independently of electronic structure errors, we use a two-dimensional Morse-type potential for which exact quantum benchmarks are available, and show that the local harmonic approach yields more accurate results than global harmonic approximations, especially for the emission spectra from higher initial vibrational levels. Next, we compare the global and local harmonic SVL spectra of anthracene, where the more expensive local harmonic corrections turn out to be less important as long as the correct global harmonic model is used. We also present additional local harmonic results for difluorocarbene, where treating anharmonicity is essential for accurate evaluation of the spectra. Yet, we also show that the structure of the difluorocarbene spectra can be explained qualitatively (but not quantitatively) with a reduced-dimensional harmonic model, for which the spectral intensities can be evaluated analytically.
academic
Catturare effetti anarmonici negli spettri di fluorescenza a singolo livello vibrazionale utilizzando la dinamica del pacchetto d'onda di Hagedorn armonica locale
Questo articolo propone un metodo che combina la dinamica del pacchetto d'onda di Hagedorn con l'approssimazione armonica locale per calcolare gli effetti anarmonici negli spettri di fluorescenza a singolo livello vibrazionale (SVL). Sebbene la dinamica del pacchetto d'onda di Hagedorn fornisca spettri SVL esatti nel modello armonico globale, gli effetti anarmonici delle molecole reali hanno un'influenza significativa sugli spettri. Gli autori verificano l'accuratezza del metodo mediante test di riferimento quantistico esatto su potenziali di tipo Morse bidimensionali e lo applicano al calcolo spettrale di molecole di antracene e difluorocarbene, dimostrando l'efficacia del metodo nel trattamento degli effetti anarmonici.
Problema Centrale: Nella spettroscopia di fluorescenza a singolo livello vibrazionale, è necessario propagare pacchetti d'onda non gaussiani, mentre i metodi tradizionali di dinamica gaussiana, sebbene esatti in potenziali armonici, non riescono a gestire gli effetti anarmonici delle molecole reali.
Importanza del Problema:
La spettroscopia SVL è uno strumento importante per lo studio del rilassamento intramolecolare, la caratterizzazione della struttura vibrazionale e l'identificazione conformazionale
Il calcolo accurato degli spettri elettronici a risoluzione vibrazionale di molecole poliatomatiche richiede generalmente la considerazione dell'anarmonicità della superficie di energia potenziale molecolare
Per molecole flessibili con vibrazioni anarmoniche di grande ampiezza, il trattamento quantistico è spesso necessario
Limitazioni dei Metodi Esistenti:
I metodi tradizionali indipendenti dal tempo per il calcolo dei fattori di Franck-Condon comportano un carico computazionale enorme nei sistemi ad alta dimensionalità
L'approssimazione armonica globale non riesce a catturare gli effetti anarmonici
I metodi di traiettoria semiclassica esistenti, sebbene possano gestire l'anarmonicità, hanno difficoltà nel trattare stati vibrazionali eccitati arbitrari
Motivazione della Ricerca: Sviluppare un metodo efficiente in grado di calcolare gli spettri SVL di livelli vibrazionali arbitrari da una singola traiettoria gaussiana e includere parzialmente gli effetti anarmonici.
Metodo Teorico: Propone il metodo della dinamica del pacchetto d'onda di Hagedorn armonico locale, combinando il metodo del pacchetto d'onda di Hagedorn con l'approssimazione armonica locale
Vantaggi Algoritmici: Realizza il calcolo degli spettri SVL per tutti i livelli vibrazionali iniziali da una singola traiettoria gaussiana
Verifica dell'Accuratezza: Verifica l'accuratezza del metodo mediante test di riferimento quantistico esatto su potenziali Morse bidimensionali accoppiati
Applicazioni Pratiche: Applica con successo il metodo ai calcoli ab initio di molecole di antracene e difluorocarbene, dimostrando l'utilità pratica del metodo
Analisi Teorica: Fornisce una spiegazione analitica della scissione dell'inviluppo spettrale del difluorocarbene
Invarianza dei Coefficienti: Sotto l'approssimazione armonica locale, i coefficienti cK nel pacchetto d'onda di Hagedorn rimangono indipendenti dal tempo
Calcolo a Singola Traiettoria: Gli spettri SVL per tutti i livelli vibrazionali possono essere ottenuti da una singola traiettoria gaussiana
Equazione di Schrödinger Non Lineare: La funzione di Hagedorn è una soluzione esatta dell'equazione di Schrödinger non lineare dipendente dal tempo
Integrali di Sovrapposizione Ricorsivi: Utilizza espressioni algebriche ricorsive esatte per calcolare gli integrali di sovrapposizione delle funzioni di Hagedorn
Emissione dello Stato Fondamentale (10²0): Eccellente accordo tra il modello armonico locale e armonico adiabatico con i risultati quantistici esatti
Emissione dello Stato Eccitato: Il metodo armonico locale mostra prestazioni migliori nella regione a bassa frequenza (corrispondente a transizioni di livelli vibrazionali elevati)
Accuratezza della Posizione dei Picchi: Il modello armonico verticale mostra le peggiori prestazioni nella posizione dei picchi
Mediante un modello analitico di oscillatore armonico spostato unidimensionale, si verifica il fenomeno della scissione dell'inviluppo spettrale del difluorocarbene:
Fattore di Huang-Rhys: S = 6.15
L'intensità di transizione segue polinomi di Laguerre associati modulati da una distribuzione di Poisson
Il metodo di Hagedorn dipendente dal tempo è completamente coerente con la formula analitica