Computation of the S1 ← S0 Vibronic Absorption Spectrum of Formaldehyde by Variational Gaussian Wavepacket and Semiclassical IVR Methods

A time averaged semiclassical approach to IR spectroscopy

We propose a new semiclassical approach to the calculation of molecular IR spectra. The method employs the time averaging technique of Kaledin and Miller upon symmetrization of the quantum dipole-dipole autocorrelation function. Spectra at high and low temperatures are investigated. In the first case, we are able to point out the possible presence of hot bands in the molecular absorption line shape. In the second case, we are able to reproduce accurate IR spectra as demonstrated by a calculation of the IR spectrum of the water molecule, which is within 4% of the exact intensity. Our time averaged IR spectra can be directly compared to time averaged semiclassical power spectra as shown in an application to the CO2 molecule, which points out the differences between IR and power spectra and demonstrates that our new approach can identify active IR transitions correctly. Overall, the method features excellent accuracy in calculating absorption intensities and provides estimates for the frequencies of vibrations in agreement with the corresponding power spectra. In perspective, this work opens up the possibility to interface the new method with the semiclassical techniques developed for power spectra, such as the divide-and-conquer one, to get accurate IR spectra of complex and high-dimensional molecular systems.

Finite-temperature vibronic spectra from the split-operator coherence thermofield dynamics

We present a numerically exact approach for evaluating vibrationally resolved electronic spectra at finite temperatures using the coherence thermofield dynamics. In this method, which avoids implementing an algorithm for solving the von Neumann equation for coherence, the thermal vibrational ensemble is first mapped to a pure-state wavepacket in an augmented space, and this wavepacket is then propagated by solving the standard, zero-temperature Schrödinger equation with the split-operator Fourier method. We show that the finite-temperature spectra obtained with the coherence thermofield dynamics in a Morse potential agree exactly with those computed by Boltzmann-averaging the spectra of individual vibrational levels. Because the split-operator thermofield dynamics on a full tensor-product grid is restricted to low-dimensional systems, we briefly discuss how the accessible dimensionality can be increased by various techniques developed for the zero-temperature split-operator Fourier method.

High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics

Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but they still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multidimensional, nonseparable coupled Morse potential.

DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.

Homogeneous Dephasing in Photosynthetic Bacterial Reaction Centers: Time Correlation Function Approach

A new tractable linear electronic transition dipole moment time correlation function (ETDMTCF) that accurately accounts for electronic dephasing, asymmetry, and width of 1-phonon profile, which the zero-phonon line (ZPL) contributes to it, in Rhodopseudomonas viridis bacterial reaction center is derived. This time correlation function proves to be superior to other frequency-domain expressions in case of strong electron-phonon coupling (which is often the case in bacterial RCs and pigment-protein complexes), many vibrational modes involved, and high temperature, whereby more vibronic and electronic (sequence) transitions would arise. The Fourier transform of this ETDMTCF leads to asymmetric multiphonon profiles composed of Lorentzian distribution and Gaussian distribution on the high- and low-energy sides, respectively, whereby the overtone widths fold themselves with that of the one-phonon profile. This ETDMTCF also features expedient computation in large systems using asymmetric phonon profiles to account correctly for dephasing and pigment-protein interaction (electron-phonon coupling). The derived ETDMTCF allows computing all nonlinear optical signals in both time and frequency domains, through the nonlinear dipole moment time correlation functions (as guided by nonlinear optical response theory) in line with the eight Liouville space pathways. The linear transition dipole moment time correlation function is of a central value as the nonlinear transition dipole moment time correlation function is expressed in terms of the linear transition dipole moment time correlation function, derived herein. One of the great advantages of presenting this ETDMTCF is its applicability to nonlinear transition dipole moment time correlation functions in line with the eight Liouville space pathways needed in computing nonlinear signals. As such, there is more to the utility and applicability of the presented ETDMTCF besides computational expediency and efficiency. Results show good agreement with the reported literature. The intimate connection between a one-phonon profile and the corresponding bath spectral density in photosynthetic complexes is discussed.

Isotope Effects on the Electronic Spectra of Ammonia from Ab Initio Semiclassical Dynamics

Despite its simplicity, the single-trajectory thawed Gaussian approximation has proven useful for calculating the vibrationally resolved electronic spectra of molecules with weakly anharmonic potential energy surfaces. Here, we show that the thawed Gaussian approximation can capture surprisingly well even more subtle observables, such as the isotope effects in the absorption spectra, and we demonstrate it on the four isotopologues of ammonia (NH3, NDH2, ND2H, and ND3). The differences in their computed spectra are due to the differences in the semiclassical trajectories followed by the four isotopologues, and the isotope effects─narrowing of the transition band and reduction of the peak spacing─are accurately described by this semiclassical method. In contrast, the adiabatic harmonic model shows a double progression instead of the single progression seen in the experimental spectra. The vertical harmonic model correctly shows only a single progression but fails to describe the anharmonic peak spacing. Analysis of the normal-mode activation upon excitation provides insight into the elusiveness of the symmetric stretching progression in the spectra.

Mixed Quantum-Classical Liouville Equation Treatment of Electronic Spectroscopy of Condensed Systems: Harmonic and Anharmonic Electron-Phonon Coupling

This Review integrates the use of electronic optical response function theory and the mixed quantum-classical (MQC) Liouville equation (MQCLE), thereby leading to electronic spectroscopy in MQC media. It further sheds light on the applicability, utility, and efficiency of the mixed quantum-classical dynamics (MQCD) formalism, which starts off with the MQCLE, in probing spectroscopy and dynamics of condensed systems, whereby quantum mechanics and classical mechanics are combined systematically. The author has been exploring and implementing MQCD to investigate electron-phonon coupling effects on electronic dephasing in harmonic and anharmonic systems by calculating linear and nonlinear optical transition analytically and numerically dipole moment time correlation functions in an MQC environment, thereby presenting an in depth spectral profile analysis and their shape and symmetry. The distinctive capability of the MQC time correlation functions is that ergodicity and stationarity properties are inherently satisfied as part of the mixed quantum-classical dynamics (MQCD) framework, unlike classical correlation functions. While some research groups have applied MQCLE to calculate vibrational spectra to study hydrogen-bonded complexes in a MQC environment and other groups calculated Optical response function to probe electron transfer dynamics using the basis mapping technique, the approach, purpose, rigor, applications, and path to the end results reported herein are different. Finally, the same framework is employed to study dissipative systems in the MQC limit, whereby the zero-phonon line adopts the correct width and eliminates its asymmetry. While the full quantum mechanical model, like the multimode Brownian oscillator (MBO) model, yields the correct width and inaccurate shape in the low-temperature limit, the MQCD formalism seems to produce an accurate zero-phonon profile. Nonlinear optical signals are also reviewed in MQC media to show the applicability and utility of this approach. The vibronic optical response functions developed here will account for geometry change, frequency change, and anharmonicity upon electronic excitation to accurately probe electronic dephasing, electron-phonon coupling, shape, and symmetry of profiles and present differences and similarities to the MBO model on pure electronic dephasing. Frequency change and anharmonicity are vitally crucial for accurately assessing electron-phonon coupling upon electronic excitation. This is an additional unique result obtained by the author to further demonstrate the applicability and utility of this approach over other approximation schemes in probing electronic dephasing including that of the MBO model.

Computational Design of Molecular Probes for Electronic Preresonance Raman Scattering Microscopy

Recently developed electronic preresonance stimulated Raman scattering (epr-SRS) microscopy, in which the Raman signal of a dye is significantly boosted by setting the incident laser frequency near the electronic excitation energy, has pushed the sensitivity of SRS microscopy close to that offered by confocal fluorescence microscopy. Prominently, the maintained narrow line-width of epr-SRS also offers high multiplexity that breaks the "color barrier" in optical microscopy. However, detailed understanding of the fundamental mechanism in these epr-SRS dyes still remains elusive. Here, we combine experiments with theoretical modeling to investigate the structure-function relationship, aiming to facilitate the design of new probes and expanding epr-SRS palettes. Our ab initio approach employing the displaced harmonic oscillator (DHO) model provides a consistent agreement between simulated and experimental SRS intensities of various triple-bond bearing epr-SRS probes with distinct scaffolds. We further review two popular approximate expressions for epr-SRS, namely the short-time and Albrecht A-term equations, and compare them to the DHO model. Overall, the theory allows us to illustrate how the observed intensity differences between molecular scaffolds stem from the coupling strength between the electronic excitation and the targeted vibrational mode, leading to a general design strategy for highly sensitive next-generation vibrational imaging probes.

Applicability of the Thawed Gaussian Wavepacket Dynamics to the Calculation of Vibronic Spectra of Molecules with Double-Well Potential Energy Surfaces

Simulating vibrationally resolved electronic spectra of anharmonic systems, especially those involving double-well potential energy surfaces, often requires expensive quantum dynamics methods. Here, we explore the applicability and limitations of the recently proposed single-Hessian thawed Gaussian approximation for the simulation of spectra of systems with double-well potentials, including 1,2,4,5-tetrafluorobenzene, ammonia, phosphine, and arsine. This semiclassical wavepacket approach is shown to be more robust and to provide more accurate spectra than the conventional harmonic approximation. Specifically, we identify two cases in which the Gaussian wavepacket method is especially useful due to the breakdown of the harmonic approximation: (i) when the nuclear wavepacket is initially at the top of the potential barrier but delocalized over both wells, e.g., along a low-frequency mode, and (ii) when the wavepacket has enough energy to classically go over the low potential energy barrier connecting the two wells. The method is efficient and requires only a single classical ab initio molecular dynamics trajectory, in addition to the data required to compute the harmonic spectra. We also present an improved algorithm for computing the wavepacket autocorrelation function, which guarantees that the evaluated correlation function is continuous for arbitrary size of the time step.

Excited state distribution function for probing Herzberg-Teller vibronic coupling using linear optical response theory: Application to glassy pheophytin a

The goal of the present work is to develop an excited-state distribution function that can be used to calculate electronic transition dipole moment time correlation functions at a considerably low computational cost. An additional merit of the distribution function is its capability to probe the Hertzberg-Teller vibronic coupling effect in terms of the previously reported Condon correlation functions in the literature without having to start from the equilibrium density operator to probe spectral non-Condon effects, thereby exploring Hertzberg-Teller vibronic coupling by building on the Condon regime. It is easily extendable to anharmonic systems. Model calculations are reported to show the high degree of accuracy and computational efficiency of the presented approach. Application to a photosynthetic system such as pheophytin a in triethylamine is provided.

Finite-Temperature, Anharmonicity, and Duschinsky Effects on the Two-Dimensional Electronic Spectra from Ab Initio Thermo-Field Gaussian Wavepacket Dynamics

Accurate description of finite-temperature vibrational dynamics is indispensable in the computation of two-dimensional electronic spectra. Such simulations are often based on the density matrix evolution, statistical averaging of initial vibrational states, or approximate classical or semiclassical limits. While many practical approaches exist, they are often of limited accuracy and difficult to interpret. Here, we use the concept of thermo-field dynamics to derive an exact finite-temperature expression that lends itself to an intuitive wavepacket-based interpretation. Furthermore, an efficient method for computing finite-temperature two-dimensional spectra is obtained by combining the exact thermo-field dynamics approach with the thawed Gaussian approximation for the wavepacket dynamics, which is exact for any displaced, distorted, and Duschinsky-rotated harmonic potential but also accounts partially for anharmonicity effects in general potentials. Using this new method, we directly relate a symmetry breaking of the two-dimensional signal to the deviation from the conventional Brownian oscillator picture.

Born-Oppenheimer approximation in optical cavities: from success to breakdown

The coupling of a molecule and a cavity induces nonadiabaticity in the molecule which makes the description of its dynamics complicated. For polyatomic molecules, reduced-dimensional models and the use of the Born-Oppenheimer approximation (BOA) may remedy the situation. It is demonstrated that contrary to expectation, BOA may even fail in a one-dimensional model and is generally expected to fail in two- or more-dimensional models due to the appearance of conical intersections induced by the cavity.

On-the-fly ab initio semiclassical evaluation of vibronic spectra at finite temperature

To compute and analyze vibrationally resolved electronic spectra at zero temperature, we have recently implemented the on-the-fly ab initio extended thawed Gaussian approximation [A. Patoz et al., J. Phys. Chem. Lett. 9, 2367 (2018)], which accounts for anharmonicity, mode-mode coupling, and Herzberg-Teller effects. Here, we generalize this method in order to evaluate spectra at non-zero temperature. In line with thermo-field dynamics, we transform the von Neumann evolution of the coherence component of the density matrix to the Schrödinger evolution of a wavefunction in an augmented space with twice as many degrees of freedom. Due to the efficiency of the extended thawed Gaussian approximation, this increase in the number of coordinates results in nearly no additional computational cost. More specifically, compared to the original, zero-temperature approach, the finite-temperature method requires no additional ab initio electronic structure calculations. At the same time, the new approach allows for a clear distinction among finite-temperature, anharmonicity, and Herzberg-Teller effects on spectra. We show, on a model Morse system, the advantages of the finite-temperature thawed Gaussian approximation over the commonly used global harmonic methods and apply it to evaluate the symmetry-forbidden absorption spectrum of benzene, where all of the aforementioned effects contribute.

Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene

Azulene is a prototypical molecule with an anomalous fluorescence from the second excited electronic state, thus violating Kasha's rule, and with an emission spectrum that cannot be understood within the Condon approximation. To better understand the photophysics and spectroscopy of azulene and other nonconventional molecules, we developed a systematic, general, and efficient computational approach combining the semiclassical dynamics of nuclei with ab initio electronic structure. First, to analyze the nonadiabatic effects, we complement the standard population dynamics by a rigorous measure of adiabaticity, estimated with the multiple-surface dephasing representation. Second, we propose a new semiclassical method for simulating non-Condon spectra, which combines the extended thawed Gaussian approximation with the efficient single-Hessian approach. S₁ ← S₀ and S₂ ← S₀ absorption and S₂ → S₀ emission spectra of azulene, recorded in a new set of experiments, agree very well with our calculations. We find that accuracy of the evaluated spectra requires the treatment of anharmonicity, Herzberg-Teller, and mode-mixing effects.

Reduced rovibrational coupling Cartesian dynamics for semiclassical calculations: Application to the spectrum of the Zundel cation

We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from the work of Huang et al. [J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows for the physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid initially exciting lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of ∼29 cm^-1 with respect to available multiconfiguration time-dependent Hartree calculations and MAE ∼ 14 cm^-1 when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster systems.

Single-Hessian thawed Gaussian approximation

To alleviate the computational cost associated with on-the-fly ab initio semiclassical calculations of molecular spectra, we propose the single-Hessian thawed Gaussian approximation in which the Hessian of the potential energy at all points along an anharmonic classical trajectory is approximated by a constant matrix. The spectra obtained with this approximation are compared with the exact quantum spectra of a one-dimensional Morse potential and with the experimental spectra of ammonia and quinquethiophene. In all cases, the single-Hessian version performs almost as well as the much more expensive on-the-fly ab initio thawed Gaussian approximation and significantly better than the global harmonic schemes. Remarkably, unlike the thawed Gaussian approximation, the proposed method conserves energy exactly, despite the time dependence of the corresponding effective Hamiltonian, and, in addition, can be mapped to a higher-dimensional time-independent classical Hamiltonian system. We also provide a detailed comparison with several related approximations used for accelerating prefactor calculations in semiclassical simulations.