Quantum corrections in vibrational and electronic condensed phase spectroscopy: line shapes and echoes

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. The strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

Taming the third order cumulant approximation to linear optical spectroscopy

The second order cumulant method offers a promising pathway to predicting optical properties in condensed phase systems. It allows for the computation of linear absorption spectra from excitation energy fluctuations sampled along molecular dynamics (MD) trajectories, fully accounting for vibronic effects, direct solute-solvent interactions, and environmental polarization effects. However, the second order cumulant approximation only guarantees accurate line shapes for energy gap fluctuations obeying Gaussian statistics. A third order correction has recently been derived but often yields unphysical spectra or divergent line shapes for moderately non-Gaussian fluctuations due to the neglect of higher order terms in the cumulant expansion. In this work, we develop a corrected cumulant approach, where the collective effect of neglected higher order contributions is approximately accounted for through a dampening factor applied to the third order cumulant term. We show that this dampening factor can be expressed as a function of the skewness and kurtosis of energy gap fluctuations and can be parameterized from a large set of randomly sampled model Hamiltonians for which exact spectral line shapes are known. This approach is shown to systematically remove unphysical contributions in the form of negative absorbances from cumulant spectra in both model Hamiltonians and condensed phase systems sampled from MD and dramatically improves over the second order cumulant method in describing systems exhibiting Duschinsky mode mixing effects. We successfully apply the approach to the coumarin-153 dye in toluene, obtaining excellent agreement with experiment.

The Influence of Electronic Polarization on Nonlinear Optical Spectroscopy

The environment surrounding a chromophore can dramatically affect the energy absorption and relaxation process, as manifested in optical spectra. Simulations of nonlinear optical spectroscopy, such as two-dimensional electronic spectroscopy (2DES) and transient absorption (TA), will be influenced by the computational model of the environment. We here compare a fixed point charge molecular mechanics model and a quantum mechanical (QM) model of the environment in computed 2DES and TA spectra of Nile red in water and the chromophore of photoactive yellow protein (PYP) in water and protein environments. In addition to simulating these nonlinear optical spectra, we directly juxtapose the computed excitation energy correlation function to the dynamic Stokes shift function often used to analyze environment dynamics. Overall, we find that for the three systems studied here the mutual electronic polarization provided by the QM environment manifests in broader 2DES signals, as well as a larger reorganization energy and a larger static Stokes shift due to stronger coupling between the chromophore and the environment.

Nuclear Quantum Effects from the Analysis of Smoothed Trajectories: Pilot Study for Water

Nuclear quantum effects have significant contributions to thermodynamic quantities and structural properties; furthermore, very expensive methods are necessary for their accurate computation. In most calculations, these effects, for instance, zero-point energies, are simply neglected or only taken into account within the quantum harmonic oscillator approximation. Herein, we present a new method, Generalized Smoothed Trajectory Analysis, to determine nuclear quantum effects from molecular dynamics simulations. The broad applicability is demonstrated with the examples of a harmonic oscillator and different states of water. Ab initio molecular dynamics simulations have been performed for ideal gas up to the temperature of 5000 K. Classical molecular dynamics have been carried out for hexagonal ice, liquid water, and vapor at atmospheric pressure. With respect to the experimental heat capacity, our method outperforms previous calculations in the literature in a wide temperature range at lower computational cost than other alternatives. Dynamic and structural nuclear quantum effects of water are also discussed.

Dephasing and Decoherence in Vibrational and Electronic Line Shapes

Absorption and emission line shapes of vibrational and electronic transitions in liquids are broadened by interactions with the "bath" (in this case, the rotational and translational degrees of freedom of all the molecules in the liquid). If these degrees of freedom are treated classically, the broadening process is often known as dephasing. If, on the other hand, the bath degrees of freedom are instead treated quantum mechanically, there is additional broadening due to what is known in the chemical-physics literature as decoherence. The question addressed in this paper is the relative importance of decoherence (bath quantum effects) and dephasing. We present general developments of this subject for absorption and emission line shapes, discover several new relationships connecting classical and quantum treatments of the bath, and also consider the Stokes shift (difference in peak frequencies in absorption and emission). We next draw some general conclusions by considering a model system whose transition-frequency time-correlation function has only one bath time scale. We then consider a realistic system of the vibrational OH stretch transition of dilute HOD in liquid D₂O at room temperature. For this system, we conclude that bath quantum effects are not very important, except for the Stokes shift. More generally, we argue that this is the case for many vibrational and most electronic transitions in room-temperature liquids.

Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches

Simulating optical spectra in the condensed phase remains a challenge for theory due to the need to capture spectral signatures arising from anharmonicity and dynamical effects, such as vibronic progressions and asymmetry. As such, numerous simulation methods have been developed that invoke different approximations and vary in their ability to capture different physical regimes. Here, we use several models of chromophores in the condensed phase and ab initio molecular dynamics simulations to rigorously assess the applicability of methods to simulate optical absorption spectra. Specifically, we focus on the ensemble scheme, which can address anharmonic potential energy surfaces but relies on the applicability of extreme nuclear-electronic time scale separation; the Franck-Condon method, which includes dynamical effects but generally only at the harmonic level; and the recently introduced ensemble zero-temperature Franck-Condon approach, which straddles these limits. We also devote particular attention to the performance of methods derived from a cumulant expansion of the energy gap fluctuations and test the ability to approximate the requisite time correlation functions using classical dynamics with quantum correction factors. These results provide insights as to when these methods are applicable and able to capture the features of condensed phase spectra qualitatively and, in some cases, quantitatively across a range of regimes.

Thermal weights for semiclassical vibrational response functions

Semiclassical approximations to response functions can allow the calculation of linear and nonlinear spectroscopic observables from classical dynamics. Evaluating a canonical response function requires the related tasks of determining thermal weights for initial states and computing the dynamics of these states. A class of approximations for vibrational response functions employs classical trajectories at quantized values of action variables and represents the effects of the radiation-matter interaction by discontinuous transitions. Here, we evaluate choices for a thermal weight function which are consistent with this dynamical approximation. Weight functions associated with different semiclassical approximations are compared, and two forms are constructed which yield the correct linear response function for a harmonic potential at any temperature and are also correct for anharmonic potentials in the classical mechanical limit of high temperature. Approximations to the vibrational linear response function with quantized classical trajectories and proposed thermal weight functions are assessed for ensembles of one-dimensional anharmonic oscillators. This approach is shown to perform well for an anharmonic potential that is not locally harmonic over a temperature range encompassing the quantum limit of a two-level system and the limit of classical dynamics.

Can far-IR action spectroscopy combined with BOMD simulations be conformation selective?

The combination of conformation selective far-IR/UV double resonance spectroscopy with Born–Oppenheimer molecular dynamics (BOMD) simulations is presented here for the structural characterization of the Ac-Phe-Pro-NH₂ peptide in the far-infrared spectral domain, i.e. for radiation below 800 cm⁻¹.

2D IR spectra of cyanide in water investigated by molecular dynamics simulations

Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN(-) solvated in D2O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN(-) molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN(-) and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T1 times are sensitive to the van der Waals ranges on the CN(-) is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm(-1) vs. 14.9 cm(-1)) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements.

On the role of nonbonded interactions in vibrational energy relaxation of cyanide in water

The vibrationally excited cyanide ion (CN(-)) in H2O or D2O relaxes back to the ground state within several tens of picoseconds. Pump-probe infrared spectroscopy has determined relaxation times of T1 = 28 ± 7 and 71 ± 3 ps in H2O and D2O, respectively. Atomistic simulations of this process using nonequilibrium molecular dynamics simulations allow determination of whether it is possible at all to describe such a process, what level of accuracy in the force fields is required, and whether the information can be used to understand the molecular mechanisms underlying vibrational relaxation. It is found that, by using the best electrostatic models investigated, absolute relaxation times can be described rather more qualitatively (T1(H2O) = 19 ps and T1(D2O) = 34 ps) whereas the relative change in going from water to deuterated water is more quantitatively captured (factor of 2 vs 2.5 from experiment). However, moderate adjustment of the van der Waals ranges by less than 20% (for NVT) and 7.5% (for NVE), respectively, leads to almost quantitative agreement with experiment. Analysis of the energy redistribution establishes that the major pathway for CN(-) relaxation in H2O or D2O proceeds through coupling to the water-bending plus libration mode.

Bond connectivity measured via relaxation-assisted two-dimensional infrared spectroscopy

The relaxation-assisted two-dimensional infrared (RA 2DIR) method is a novel technique for probing structures of molecules, which relies on vibrational energy transport in molecules. In this article we demonstrate the ability of RA 2DIR to detect the bond connectivity patterns in molecules using two parameters, a characteristic intermode energy transport time (arrival time) and a cross-peak amplification coefficient. A correlation of the arrival time with the distance between the modes is demonstrated. An 18-fold amplification of the cross-peak amplitude for the modes separated by approximately 11 A is shown using RA 2DIR; larger cross-peak amplifications are expected for the modes separated by larger distances. The RA 2DIR method enhances the applicability of 2DIR spectroscopy by making practical the long-range measurements using a variety of structural reporters, including weak IR modes. The data presented demonstrate the analytical power of RA 2DIR which permits the speedy structural assessments of the bond connectivity patterns.

Frequency-frequency correlation functions and apodization in two-dimensional infrared vibrational echo spectroscopy: a new approach

Ultrafast two-dimensional infrared (2D-IR) vibrational echo spectroscopy can probe structural dynamics under thermal equilibrium conditions on time scales ranging from femtoseconds to approximately 100 ps and longer. One of the important uses of 2D-IR spectroscopy is to monitor the dynamical evolution of a molecular system by reporting the time dependent frequency fluctuations of an ensemble of vibrational probes. The vibrational frequency-frequency correlation function (FFCF) is the connection between the experimental observables and the microscopic molecular dynamics and is thus the central object of interest in studying dynamics with 2D-IR vibrational echo spectroscopy. A new observable is presented that greatly simplifies the extraction of the FFCF from experimental data. The observable is the inverse of the center line slope (CLS) of the 2D spectrum. The CLS is the inverse of the slope of the line that connects the maxima of the peaks of a series of cuts through the 2D spectrum that are parallel to the frequency axis associated with the first electric field-matter interaction. The CLS varies from a maximum of 1 to 0 as spectral diffusion proceeds. It is shown analytically to second order in time that the CLS is the T(w) (time between pulses 2 and 3) dependent part of the FFCF. The procedure to extract the FFCF from the CLS is described, and it is shown that the T(w) independent homogeneous contribution to the FFCF can also be recovered to yield the full FFCF. The method is demonstrated by extracting FFCFs from families of calculated 2D-IR spectra and the linear absorption spectra produced from known FFCFs. Sources and magnitudes of errors in the procedure are quantified, and it is shown that in most circumstances, they are negligible. It is also demonstrated that the CLS is essentially unaffected by Fourier filtering methods (apodization), which can significantly increase the efficiency of data acquisition and spectral resolution, when the apodization is applied along the axis used for obtaining the CLS and is symmetrical about tau=0. The CLS is also unchanged by finite pulse durations that broaden 2D spectra.

Theoretical characterization of carbon monoxide vibrational spectrum in sperm whale myoglobin distal pocket

In this article we use the perturbed matrix method and an extended molecular dynamics sampling of the carbon monoxide (CO) in the myoglobin distal pocket to characterize the CO vibrational spectrum and hence to relate its spectroscopic features with the atomic-molecular behavior. Results show the accuracy of the method employed and confirm the assignment of the spectroscopic B1 and B2 states proposed by Lim et al.

Quantal cumulant dynamics: General theory

The authors have derived coupled equations of motion of cumulants that consist of a symmetric-ordered product of the position and momentum fluctuation operators in one dimension. The key point is the utilization of a position shift operator acting on a potential operator, where the expectation value of the shift operator is evaluated using the cumulant expansion technique. In particular, the equations of motion of the second-order cumulant and the expectation values of the position and momentum operators are given. The resultant equations are expressed by those variables and a quantal potential that consists of an exponential function of the differential operators and the original potential. This procedure enables us to perform quantal (semiclassical) dynamics in one dimension. In contrast to a second-order quantized Hamilton dynamics by Prezhdo and Pereverzev which conserves the total energy only with an odd-order Taylor expansion of the potential [J. Chem. Phys. 116, 4450 (2002); 117, 2995 (2002)], the present quantal cumulant dynamics method exactly conserves the energy, even if a second-order approximation of the cumulants is adopted, because the present scheme does not truncate the given potential. The authors propose three schemes, (i) a truncation, (ii) a summation of derivatives, and (iii) a convolution method, for evaluating the quantal potentials for several types of potentials. The numerical results show that although the truncation method preserves the energy to some degree, the trajectory obtained gradually deviates from that of the summation scheme after 2000 steps. The phase space structure obtained by the truncation scheme is also different from that of the summation scheme in a strongly anharmonic region.

Fifth-order two-dimensional Raman spectroscopy of liquid water, crystalline ice Ih and amorphous ices: Sensitivity to anharmonic dynamics and local hydrogen bond network structure

The theoretical study of off-resonant fifth-order two-dimensional (2D)-Raman spectroscopy is made to analyze the intermolecular dynamics of liquid and solid water. The 2D-Raman spectroscopy is susceptible to the nonlinear anharmonic dynamics and local hydrogen bond structure in water. It is found that the distinct 2D-Raman response appears as the negative signal near the t(2) axis. The origin of this negative signal for t(2)<15 fs is from the nonlinear polarizability in the librational motions, whereas that for 30 fs<t(2)<150 fs is attributed to the anharmonic translational motions. It is found that the mechanical anharmonicity and nonlinear polarizability couplings among modes clearly can be observed as the sum- and difference-frequency peaks in the 2D-Raman spectrum (i.e., Fourier transforms of the response). The 2D-Raman spectroscopies of ice Ih and amorphous ices, i.e., low density, high density, and very high density amorphous ices, are also investigated. It is found that the 2D-Raman spectroscopy is very sensitive to the anisotropy of the structure of ice Ih. The strong hydrogen bond stretching band is seen in the 2D-Raman spectroscopy of the polarization directions parallel to the c axis, whereas the contributions of the librational motion can be also seen in the spectrum with the polarization directions parallel to the a axis. The 2D-Raman spectroscopy is also found to be also very sensitive to the differences in local hydrogen bond network structures in various amorphous phases.

Evaluation of quantum correlation functions from classical data: Anharmonic models

The previously introduced method of evaluating quantum mechanical time correlation functions using as input only classical simulation data is generalized and applied to two anharmonic model systems, as a further test. The quantum correction approach utilizes the relation between a general quantum correlation function and its classical analog. For the tested models, we obtain numerical results of nonlinear correlation functions with comparable accuracy to that of the centroid molecular dynamics method, although the present method is much simpler to implement and not limited to real valued quantum correlation functions.

CONNECTING CHEMICAL DYNAMICS IN GASES AND LIQUIDS

▪ Abstract  Modern ultrafast spectroscopic techniques provide new opportunities to study chemical reaction dynamics in liquids and hold the possibility of obtaining much of the same detailed information available in gases. Vibrational energy transfer studies are the most advanced of the investigations and demonstrate that it is possible to observe state-specific pathways of energy flow within a vibrationally excited molecule (intramolecular vibrational relaxation) and into the surrounding solvent molecules (intermolecular energy transfer). Energy transfer in liquids and gases share many common aspects, but the presence of the solvent also alters the relaxation in both obvious and subtle ways. Photodissociation is amenable to similarly detailed study in liquids, and there are informative new measurements. Bimolecular reactions have received the least attention in state-resolved measurements in liquids, but the means to carry them much further now exist. Studying photodissociation and bimolecular reaction of molecules prepared with initial vibrational excitation in liquids is a realistic, but challenging, goal.

IR spectra of N-methylacetamide in water predicted by combined quantum mechanical/molecular mechanical molecular dynamics simulations

We applied the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulation method in assessing IR spectra of N-methylacetamide and its deuterated form in aqueous solutions. The model peptide is treated at the Austin Model 1 (AM1) level and the induced dipole effects by the solvent are incorporated in fluctuating solute dipole moments, which are calculated using partial charges from Mulliken population analyses without resorting to any available high-level ab initio dipole moment data. Fourier transform of the solute dipole autocorrelation function produces in silico IR spectra, in which the relative peak intensities and bandwidths of major amide bands are quantitatively compatible with experimental results only when both geometric and electronic polarizations of the peptide by the solvent are dealt with at the same quantum-mechanical level. We cast light on the importance of addressing dynamic charge fluctuations of the solute in calculating IR spectra by comparing classical and QM/MM MD simulation results. We propose the adjustable scaling factors for each amide mode to be directly compared with experimental data.