Resonance Raman spectra of uracil based on Kramers–Kronig relations using time-dependent density functional calculations and multireference perturbation theory

Time-Dependent Resonant Inelastic X-ray Scattering of Pyrazine at the Nitrogen K-Edge: A Quantum Dynamics Approach

We calculate resonant inelastic X-ray scattering spectra of pyrazine at the nitrogen K-edge in the time domain including wavepacket dynamics in both the valence and core-excited state manifolds. Upon resonant excitation, we observe ultrafast non-adiabatic population transfer between core-excited states within the core-hole lifetime, leading to molecular symmetry distortions. Importantly, our time-domain approach inherently contains the ability to manipulate the dynamics of this process by detuning the excitation energy, which effectively shortens the scattering duration. We also explore the impact of pulsed incident X-ray radiation, which provides a foundation for state-of-the-art time-resolved experiments with coherent pulsed light sources.

Sum-over-state expressions including second-order Herzberg-Teller effects for the calculation of absorption and resonance Raman intensities

The sum-over-state expressions are derived to calculate the second-order Herzberg-Teller (HT) effects in absorption and resonance Raman spectroscopies. These effects depend on the second derivatives of the transition dipole moment with respect to the vibrational coordinates. The method is applied to the molecule of 1,3-butadiene using density functional theory calculations. It is found that the second-order HT effects are significant for both absorption and resonance Raman intensities, and that the calculated spectra are in good agreement with the experimental data. The second-order HT effects originate from diagonal elements of the second derivatives matrix, whereas non-diagonal elements have a negligible impact on the intensities of 1,3-butadiene.

All-experimental analysis of doubly resonant sum-frequency generation spectra: Application to aggregated rhodamine films

A new method is proposed to analyze Doubly Resonant infrared-visible Sum-Frequency Generation (DR-SFG) spectra. Based on the transform technique, this approach is free from assumptions about vibronic modes, energies, or line widths and accurately captures through the overlap spectral function all required aspects of the vibronic structure from simple experimental linear absorption spectra. Details and implementation of the method are provided along with three examples treating rhodamine thin films about one monolayer thick. The technique leads to a perfect agreement between experiment and simulations of the visible DR-SFG line shapes, even in the case of complex intermolecular interactions resulting from J-aggregated chromophores in heterogeneous films. For films with mixed H- and J-aggregates, separation of their responses shows that the J-aggregate DR-SFG response is dominant. Our analysis also accounts for the unexplained results published in the early times of DR-SFG experiments.

Role of Explicit Solvation in the Simulation of Resonance Raman Spectra within Short-Time Dynamics Approximation

In short-time dynamics approximation, relative resonance Raman (RR) intensity of a vibrational mode primarily depends on the magnitude of square of the excited-state gradient along the corresponding normal coordinate, ground-state normal mode eigenvector, and harmonic vibrational wavenumbers. In this study, through simulation of RR spectra of guanosine-5'-monophosphate (GMP) in two ππ* singlet excited states, we analyze how the explicitly hydrogen-bonded local solvation structure of the chromophore dictates intensities of the RR active modes in an unprecedented manner. We show that the accuracy of the structural model of solvated chromophore plays a decisive role in determining an optimal theoretical method for prediction of the Franck-Condon region of the ππ* excited states. 9-Methylguanine (9-meG) in complex with six water molecules (9-meG·6H₂O) is found out to be the most accurate one for describing GMP in two different bright electronic states. We find that explicit hydrogen-bonded water molecules strongly influence computed RR intensities of GMP by modulating both the ground-state normal mode vectors and the excited-state energy gradients. We find that simultaneous inclusion of six explicit waters to describe the solute-solvent interaction near all hydration sites is essential for reliable prediction of the features of RR spectra in L_b and B_b electronic states of GMP.

Ehrenfest+R dynamics. II. A semiclassical QED framework for Raman scattering

In Paper I [Chen et al., J. Chem. Phys. 150, 044102 (2019)], we introduced Ehrenfest+R dynamics for a two-level system and showed how spontaneous emission can be heuristically included such that, after averaging over an ensemble of Ehrenfest+R trajectories, one can recover both coherent and incoherent electromagnetic fields. In the present paper, we now show that Ehrenfest+R dynamics can also correctly describe Raman scattering, whose features are completely absent from standard Ehrenfest dynamics. Ehrenfest+R dynamics appear to be quantitatively accurate both for resonant and off-resonant Raman signals, as compared with Kramers-Heisenberg-Dirac theory.

Time-Dependent Formulation of Resonance Raman Optical Activity Spectroscopy

In this work, we extend the theoretical framework recently developed for the simulation of resonance Raman (RR) spectra of medium-to-large sized systems to its chiral counterpart, namely, resonance Raman optical activity (RROA). The theory is based on a time-dependent (TD) formulation, with the transition tensors obtained as half-Fourier transforms of the appropriate cross-correlation functions. The implementation has been kept as general as possible, supporting adiabatic and vertical models for the PES representation, both in Cartesian and internal coordinates, with the possible inclusion of Herzberg-Teller (HT) effects. Thanks to the integration of this TD-RROA procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects can be included in an effective way. The implementation is validated on one of the smallest chiral molecule (methyloxirane). Practical applications are illustrated with three medium-size organic molecules (naproxen-OCD₃, quinidine and 2-Br-hexahelicene), whose simulated spectra are compared to the corresponding experimental data.

The role of Herzberg-Teller effects on the resonance Raman spectrum of trans-porphycene investigated by time dependent density functional theory

The S₁ excited state properties as well as the associated absorption and resonance Raman (RR) spectra of trans-porphycene are investigated by means of time dependent density functional theory calculations. The relative magnitude of the Franck-Condon (FC) contribution and of the Herzberg-Teller (HT) effects is evaluated for both the absorption and RR intensities. The accuracy of the calculated spectra is assessed by employing different theoretical approximations and by comparing with experimental data. The obtained results show that Duschinsky effects lead to noticeable modifications in the absorption intensities but are nearly negligible in the RR spectrum. By contrast, the HT effects are stronger for the RR intensities compared to the absorption intensities, and these effects significantly improve the agreement with the experimental RR spectrum. Moreover, the HT effects produce different values of the RR depolarization ratios, which can be used to quantify the relative importance of the FC and HT contributions. Generally, it is found that the HT effects have a significant role on the RR spectrum of trans-porphycene and that their inclusion in the computational scheme is mandatory to accurately predict the RR intensities.

Ultrafast Nuclear Dynamics of Photoexcited Guanosine-5'-Monophosphate in Three Singlet States

We report measurement of resonance Raman (RR) spectra of guanosine-5'-monophosphate (GMP), a DNA nucleotide at excitation wavelengths throughout its ππ* absorption band (B_b) in the 210-230 nm range. From these data, we constructed wavelength-dependent Raman intensity excitation profiles (REPs) for all observed modes. These profiles and the absorption spectrum were then modeled using self-consistent simulations based on the time-dependent wave packet propagation formalism. We inferred the initial structural dynamics of GMP immediately after photoexcitation in terms of dimensionless displacements. The simulations also provide linewidth-broadening parameters that in turn report on the time scale of dynamics. We compared deduced structural changes in the purine ring upon photoabsorption into the B_b state with those deduced for the two lowest lying ππ* (L_a and L_b at 280 and 248 nm, respectively) excited states of GMP. We find that excitation to the L_b state lengthens C₆-N₁ and C₂═N₃ bonds, which lie along the formation coordinate of various oxidative adducts but B_b excitation does not. We also find that photoabsorption by the B_b state weakens the C₈-N₉ bond and thus might assist imidazole ring opening via cleavage of the same bond. Electronic excitation to different ππ* states of the guanine chromophore results in contrasting structural changes; although absorption by the L_a and L_b states induces expansion of pyrimidine and contraction of imidazole rings, excitation results in overall shrinkage of both the rings. Computed absolute changes in internal coordinates imply that photoexcitation to any of the three singlet states of GMP does not lead directly to the formation of a cation radical of guanine.

Theoretical Assessment of Excited State Gradients and Resonance Raman Intensities for the Azobenzene Molecule

The ground state geometries and vibrational frequencies as well as the excitation energies and excited state gradients of the S₁(nπ*) and S₂(ππ*) states of trans- and cis-azobenzene are investigated by several DFT methods, namely B3LYP, PBE, M06-2X, CAM-B3LYP, and ωB97X. Excited state properties and in particular gradients are also assessed using the wave function based methods EOM-CCSD and RASPT2/RASSCF. Comparison with experimental data shows that the B3LYP functional gives the most accurate results for the ground state geometry and vibrational frequencies. The analysis of the vertical excitation energies reveals that the RASPT2 approach provides the most accurate excitation energies with deviations of the order of 0.1 eV. Among the TDDFT methods, the CAM-B3LYP functional shows the best performance on the excitation energies. By assessing the excited state gradients with respect to the reference RASPT2 data, the most accurate gradients are obtained with B3LYP, whereas other functionals as well as the EOM-CCSD and RASSCF calculations give less consistent results. Overall, despite the tendency of B3LYP to underestimate the excitation energies, this functional provides the most balanced description of both ground and excited state properties for both isomers of azobenzene in the Franck-Condon region.

Multiconfigurational Effects in Theoretical Resonance Raman Spectra

We analyze resonance Raman spectra of the nucleobase uracil in the short‐time approximation calculated with multiconfigurational methods. We discuss the importance of static electron correlation by means of density‐matrix renormalization group self‐consistent field (DMRG‐SCF) calculations. Our DMRG‐SCF results reveal that a minimal active orbital space that leads to a qualitatively correct description of the resonance Raman spectrum of uracil should encompass parts of the σ/σ* bonding/anti‐bonding orbitals of the pyrimidine ring. We trace these findings back to the considerable entanglement between the σ/σ* bonding/anti‐bonding as well as valence π/π* orbitals in the excited‐state electronic structure of uracil, which indicates non‐negligible non‐dynamical correlation effects that are less pronounced in the electronic ground state.

Vibronic-structure tracking: a shortcut for vibrationally resolved UV/Vis-spectra calculations

The vibrational coarse structure and the band shapes of electronic absorption spectra are often dominated by just a few molecular vibrations. By contrast, the simulation of the vibronic structure even in the simplest theoretical models usually requires the calculation of the entire set of normal modes of vibration. Here, we exploit the idea of the mode-tracking protocol [M. Reiher and J. Neugebauer, J. Chem. Phys. 118, 1634 (2003)] in order to directly target and selectively calculate those normal modes which have the largest effect on the vibronic band shape for a certain electronic excitation. This is achieved by defining a criterion for the importance of a normal mode to the vibrational progressions in the absorption band within the so-called "independent mode, displaced harmonic oscillator" (IMDHO) model. We use this approach for a vibronic-structure investigation for several small test molecules as well as for a comparison of the vibronic absorption spectra of a truncated chlorophyll a model and the full chlorophyll a molecule. We show that the method allows to go beyond the often-used strategy to simulate absorption spectra based on broadened vertical excitation peaks with just a minimum of computational effort, which in case of chlorophyll a corresponds to about 10% of the cost for a full simulation within the IMDHO approach.

A robust and effective time-independent route to the calculation of Resonance Raman spectra of large molecules in condensed phases with the inclusion of Duschinsky, Herzberg-Teller, anharmonic, and environmental effects

We present an effective time-independent implementation to model vibrational resonance Raman (RR) spectra of medium-large molecular systems with the inclusion of Franck-Condon (FC) and Herzberg-Teller (HT) effects and a full account of the possible differences between the harmonic potential energy surfaces of the ground and resonant electronic states. Thanks to a number of algorithmic improvements and very effective parallelization, the full computations of fundamentals, overtones, and combination bands can be routinely performed for large systems possibly involving more than two electronic states. In order to improve the accuracy of the results, an effective inclusion of the leading anharmonic effects is also possible, together with environmental contributions under different solvation regimes. Reduced-dimensionality approaches can further enlarge the range of applications of this new tool. Applications to imidazole, pyrene, and chlorophyll a1 in solution are reported, as well as comparisons with available experimental data.

Simplified Sum-Over-States Approach for Predicting Resonance Raman Spectra. Application to Nucleic Acid Bases

Resonance Raman spectra provide a valuable probe into molecular excited-state structures and properties. Moreover, resonance enhancement is of importance for the chemical contribution to surface-enhanced Raman scattering. In this work, we introduce a simplified sum-over-states scheme for computing Raman spectra and Raman excitation profiles. The proposed sum-over-states approach uses derivatives of electronic excitation energies and transition dipole moments, which can be efficiently computed from time-dependent density functional theory. We analyze and interpret the resonance Raman spectra and Raman excitation profiles of nucleic acid bases using the present approach. Contributions of individual excited states under strictly resonant and nonresonant conditions are investigated, and smooth interpolation between both limiting cases is obtained.

Selective calculation of high-intensity vibrations in molecular resonance Raman spectra

We present an intensity-driven approach for the selective calculation of vibrational modes in molecular resonance Raman spectra. The method exploits the ideas of the mode-tracking algorithm [M. Reiher and J. Neugebauer, J. Chem. Phys. 118, 1634 (2003)] for the calculation of preselected molecular vibrations and of Heller's gradient approximation [Heller et al., J. Phys. Chem. 86, 1822 (1982)] for the estimation of resonance Raman intensities. The gradient approximation allows us to construct a basis vector for the subspace iteration carried out in the mode-tracking calculation, which corresponds to an artificial collective motion of the molecule that contains the entire intensity in the resonance Raman spectrum. Subsequently, the algorithm generates new basis vectors from which normal mode approximations are obtained. It is then possible to provide estimates for (i) the accuracy of the normal mode approximations and (ii) the intensity of these modes in the final resonance Raman spectrum. This approach is tested for the examples of uracil and a structural motif from the E colicin binding immunity protein Im7, in which a few aromatic amino acids dominate the resonance Raman spectrum at wavelengths larger than 240 nm.

Analysis and prediction of absorption band shapes, fluorescence band shapes, resonance Raman intensities, and excitation profiles using the time-dependent theory of electronic spectroscopy

A general method for the simulation of absorption (ABS) and fluorescence band shapes, resonance-Raman (rR) spectra, and excitation profiles based on the time-dependent theory of Heller is discussed. The following improvements to Heller's theory have been made: (a) derivation of new recurrence relations for the time-dependent wave packet overlap in the case of frequency changes between the ground and electronically excited states, (b) a new series expansion that gives insight into the nature of Savin's preresonance approximation, (c) incorporation of inhomogeneous broadening effects into the formalism at no additional computational cost, and (d) derivation of a new and simple short-time dynamics based equation for the Stokes shift that remains valid in the case of partially resolved vibrational structure. Our implementation of the time-dependent theory for the fitting of experimental spectra and the simulation of model spectra as well as the quantum mechanical calculation of the model parameters is discussed. The implementation covers all electronic structure approaches which are able to deliver ground- and excited-state energies and transition dipole moments. The technique becomes highly efficient if analytic gradients for the excited-state surface are available. In this case, the computational cost for the simultaneous prediction of ABS, fluorescence, and rR spectra is equal to that of a single excited-state geometry optimization step while the limitations of the short-time dynamics approximation are completely avoided. As a test case we discuss the well-known case of the strongly allowed 1 (1)A(g) --> 1 (1)B(u) transition in 1,3,5 trans-hexatriene in detail using method ranging from simple single-reference treatments to elaborate multireference electronic structure approaches. At the highest computational level, the computed spectra show the best agreement that has so far been obtained with quantum chemical methods for this problem.