Simulation of Nonlinear Electronic Spectroscopy in the Condensed Phase. In: Prigogine I, Rice SA, editors. Advances in Chemical Physics. Hoboken: John Wiley & Sons, Inc.; 1993. pp. 435-516. [DOI: 10.1002/9780470141427.ch6]

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.

Linear and nonlinear spectroscopy from quantum master equations

We investigate the accuracy of the second-order time-convolutionless (TCL2) quantum master equation for the calculation of linear and nonlinear spectroscopies of multichromophore systems. We show that even for systems with non-adiabatic coupling, the TCL2 master equation predicts linear absorption spectra that are accurate over an extremely broad range of parameters and well beyond what would be expected based on the perturbative nature of the approach; non-equilibrium population dynamics calculated with TCL2 for identical parameters are significantly less accurate. For third-order (two-dimensional) spectroscopy, the importance of population dynamics and the violation of the so-called quantum regression theorem degrade the accuracy of TCL2 dynamics. To correct these failures, we combine the TCL2 approach with a classical ensemble sampling of slow microscopic bath degrees of freedom, leading to an efficient hybrid quantum-classical scheme that displays excellent accuracy over a wide range of parameters. In the spectroscopic setting, the success of such a hybrid scheme can be understood through its separate treatment of homogeneous and inhomogeneous broadening. Importantly, the presented approach has the computational scaling of TCL2, with the modest addition of an embarrassingly parallel prefactor associated with ensemble sampling. The presented approach can be understood as a generalized inhomogeneous cumulant expansion technique, capable of treating multilevel systems with non-adiabatic dynamics.

A comparison between different semiclassical approximations for optical response functions in nonpolar liquid solutions

The temporal behavior of optical response functions (ORFs) reflects the quantum dynamics of an electronic superposition state, and as such lacks a well-defined classical limit. In this paper, we consider the importance of accounting for the quantum nature of the dynamics when calculating ORFs of different types. To this end, we calculated the ORFs associated with the linear absorption spectrum and the nonlinear two-pulse photon-echo experiment, via the following approaches: (1) the semiclassical forward-backward approach; (2) an approach based on linearizing the path-integral forward-backward action in terms of the difference between the forward and backward paths; (3) an approach based on ground state nuclear dynamics. The calculations were performed on a model that consists of a two-state chromophore solvated in a nonpolar liquid. The different methods were found to yield very similar results for the absorption spectrum and "diagonal" two-pulse photon echo (i.e., the homodyne-detected signal at time t=t(0) after the second pulse, where t(0) is the time interval between the two pulses). The different approximations yielded somewhat different results in the case of the time-integrated photon-echo signal. The reasons for the similarity between the predictions of different approximations are also discussed.

Optical linear response function with linear and diagonal quadratic electron-vibration coupling in mixed quantum-classical systems

Optical linear response function of linearly and quadratically coupled mixed quantum-classical condensed phase systems is derived. The linear response function is derived using Kapral's formalism of statistical mechanics in mixed quantum-classical systems. Our mixed quantum-classical linear dipole moment correlation function J(t) is compared with the full quantum J(t) [Y. J. Yan and S. Mukamel, J. Chem. Phys. 85, 5908 (1986)] in the high temperature limit. Model calculations and discussion of our results are presented. Various formulas of Franck-Condon factors for both linear and quadratic coupling are discussed.

On the mode-coupling theory of vibrational line broadening in near-critical fluids

Molecular-dynamics simulations of a neat atomic fluid, coupled with a simple model for vibrational frequency perturbations, are used to investigate vibrational line broadening near the liquid-gas critical point. All features of our simulations are in qualitative agreement with recent Raman experiments on nitrogen. We also use our simulation results to assess the validity of the mode-coupling theories that have been used to analyze experiment. We find that the theoretical results are not in good agreement with simulation, both for the temperature dependence of the linewidth, and for the frequency time-correlation functions. However, the mode-coupling prediction that critical line broadening is due to the diverging correlation time of the frequency fluctuations is shown to be correct.

Combined electronic structure/molecular dynamics approach for ultrafast infrared spectroscopy of dilute HOD in liquid H2O and D2O

We present a new approach that combines electronic structure methods and molecular dynamics simulations to investigate the infrared spectroscopy of condensed phase systems. This approach is applied to the OH stretch band of dilute HOD in liquid D2O and the OD stretch band of dilute HOD in liquid H2O for two commonly employed models of water, TIP4P and SPC/E. Ab initio OH and OD anharmonic transition frequencies are calculated for 100 HOD x (D2O)n and HOD x(H2O)n (n = 4-9) clusters randomly selected from liquid water simulations. A linear empirical relationship between the ab initio frequencies and the component of the electric field from the solvent along the bond of interest is developed. This relationship is used in a molecular dynamics simulation to compute frequency fluctuation time-correlation functions and infrared absorption line shapes. The normalized frequency fluctuation time-correlation functions are in good agreement with the results of previous theoretical approaches. Their long-time decay times are 0.5 ps for the TIP4P model and 0.9 ps for the SPC/E model, both of which appear to be somewhat too fast compared to recent experiments. The calculated line shapes are in good agreement with experiment, and improve upon the results of previous theoretical approaches. The methods presented are simple, and transferable to more complicated systems.

Multidimensional femtosecond correlation spectroscopies of electronic and vibrational excitations

Femtosecond visible and infrared analogues of multiple-pulse nuclear magnetic resonance techniques provide novel snapshot probes into the structure and electronic and vibrational dynamics of complex molecular assemblies such as photosynthetic antennae, proteins, and hydrogen-bonded liquids. A classical-oscillator description of these spectroscopies in terms of interacting quasiparticles (rather than transitions among global eigenstates) is developed and sets the stage for designing new pulse sequences and inverting the multidimensional signals to yield molecular structures. Considerable computational advantages and a clear physical insight into the origin of the response and the relevant coherence sizes are provided by a real-space analysis of the underlying coherence-transfer pathways in Liouville space.