Molecular theory and simulation of coherence transfer in metal carbonyls and its signature on multidimensional infrared spectra

Effects of the Intramolecular Group and Solvent on Vibrational Coupling Modes and Strengths of Fermi Resonances in Aryl Azides: A DFT Study of 4-Azidotoluene and 4-Azido-N-phenylmaleimide

The high transition dipole strength of the azide asymmetric stretch makes aryl azides good candidates as vibrational probes (VPs). However, aryl azides have complex absorption profiles due to Fermi resonances (FRs). Understanding the origin and the vibrational modes involved in FRs of aryl azides is critically important toward developing them as VPs for studies of protein structures and structural changes in response to their surroundings. As such, we studied vibrational couplings in 4-azidotoluene and 4-azido-N-phenylmaleimide in two solvents, N,N-dimethylacetamide and tetrahydrofuran, to explore the origin and the effects of intramolecular group and solvent on the FRs of aryl azides using density functional theory (DFT) calculations with the B3LYP functional and seven basis sets, 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), and 6-311++G(df,pd). Two combination bands consisting of the azide symmetric stretch and another mode form strong FRs with the azide asymmetric stretch for both molecules. The FR profile was altered by replacing the methyl group with maleimide. Solvents change the relative peak position and intensity more significantly for 4-azido-N-phenylmaleimide, which makes it a more sensitive VP. Furthermore, the DFT results indicate that a comparison among the results from different basis sets can be used as a means to predict more reliable vibrational spectra.

Tensor-Train Thermo-Field Memory Kernels for Generalized Quantum Master Equations

The generalized quantum master equation (GQME) approach provides a rigorous framework for deriving the exact equation of motion for any subset of electronic reduced density matrix elements (e.g., the diagonal elements). In the context of electronic dynamics, the memory kernel and inhomogeneous term of the GQME introduce the implicit coupling to nuclear motion and dynamics of electronic density matrix elements that are projected out (e.g., the off-diagonal elements), allowing for efficient quantum dynamics simulations. Here, we focus on benchmark quantum simulations of electronic dynamics in a spin-boson model system described by various types of GQMEs. Exact memory kernels and inhomogeneous terms are obtained from short-time quantum-mechanically exact tensor-train thermo-field dynamics (TT-TFD) simulations and are compared with those obtained from an approximate linearized semiclassical method, allowing for assessment of the accuracy of these approximate memory kernels and inhomogeneous terms. Moreover, we have analyzed the computational cost of the full and reduced-dimensionality GQMEs. The scaling of the computational cost is dependent on several factors, sometimes with opposite scaling trends. The TT-TFD memory kernels can provide insights on the main sources of inaccuracies of GQME approaches when combined with approximate input methods and pave the road for the development of quantum circuits that implement GQMEs on digital quantum computers.

Electronic Absorption Spectra from Off-Diagonal Quantum Master Equations

Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrodinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrodinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrodinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained in this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion (HEOM) method.

On simulating the dynamics of electronic populations and coherences via quantum master equations based on treating off-diagonal electronic coupling terms as a small perturbation

Quantum master equations provide a general framework for describing the dynamics of electronic observables within a complex molecular system. One particular family of such equations is based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. In this paper, we show how different choices of projection operators, as well as whether one starts out with the time-convolution or the time-convolutionless forms of the generalized quantum master equation, give rise to four different types of such off-diagonal quantum master equations (OD-QMEs), namely, time-convolution and time-convolutionless versions of a Pauli-type OD-QME for only the electronic populations and an OD-QME for the full electronic density matrix (including both electronic populations and coherences). The fact that those OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. To address this, we also extend a procedure for extracting Schrödinger picture electronic coherences from interaction picture populations recently introduced by Trushechkin in the context of time-convolutionless Pauli-type OD-QME to the other three types of OD-QMEs. The performance of the aforementioned four types of OD-QMEs is explored in the context of the Garg-Onuchic-Ambegaokar benchmark model for charge transfer in the condensed phase across a relatively wide parameter range. The results show that time-convolution OD-QMEs can be significantly more accurate than their time-convolutionless counterparts, particularly in the case of Pauli-type OD-QMEs, and that rather accurate Schrödinger picture coherences can be obtained from interaction picture electronic inputs.

Ultrafast vibrational dynamics of the tyrosine ring mode and its application to enkephalin insertion into phospholipid membranes as probed by two-dimensional infrared spectroscopy

Enkephalins are small opioid peptides whose binding conformations are catalyzed by phospholipid membranes. Binding to opioid receptors is determined by the orientation of tyrosine and phenylalanine side chains. In this work, we investigate the effects of different charged phospholipid headgroups on the insertion of the tyrosine side chain into a lipid bilayer using a combination of 2D IR spectroscopy, anharmonic DFT calculations, and third order response function modeling. The insertion is probed by using the ∼1515 cm^-1 tyrosine ring breathing mode, which we found exhibits rich vibrational dynamics on the picosecond timescale. These dynamics include rapid intramolecular vibrational energy redistribution (IVR), where some of the energy ends up in a dark state that shows up as an anharmonically shifted combination band. The waiting-time dependent 2D IR spectra also show an unusual line shape distortion that affects the extraction of the frequency-frequency correlation function (FFCF), which is the dynamic observable of interest that reflects the tyrosine side chain's insertion into the lipid bilayer. We proposed three models to account for this distortion: a hot-state exchange model, a local environment dependent IVR model, and a coherence transfer model. A qualitative analysis of these models suggests that the local environment dependent IVR rate best explains the line shape distortion, while the coherence transfer model best reproduced the effects on the FFCF. Even with these complex dynamics, we found that the tyrosine ring mode's FFCF is qualitatively correlated with the degree of insertion expected from the different phospholipid headgroups.

Observation of Ultrafast Coherence Transfer and Degenerate States with Polarization-Controlled Two-Dimensional Electronic Spectroscopy

Optical spectroscopy is a powerful tool to interrogate quantum states of matter. We present simulation results for the cross-polarized two-dimensional electronic spectra of the light-harvesting system LH2 of purple bacteria. We identify a spectral feature on the diagonal, which we assign to ultrafast coherence transfer between degenerate states. The implication for the interpretation of previous experiments on different systems and the potential use of this feature are discussed. In particular, we foresee that this kind of feature will be useful for identifying mixed degenerate states and for identifying the origin of symmetry breaking disorder in systems like LH2. Furthermore, this may help identify both vibrational and electronic states in biological systems such as proteins and solid-state materials such as hybrid perovskites.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Conducting Excitation and Emission Spectra in the IR Regime: Frequency-Domain Time-Resolved Vibrational Four Wave Mixing Spectroscopy

We report new features of recently developed ultrafast coherent multidimensional spectroscopy (CMDS), an optical analogue to multidimensional NMR. By using both frequency- and time-domain nonlinear four wave mixing methods, CMDS is able to directly observe coherence transfer (CT), the coherent quantum mechanical analogue of population relaxation. Using a mixture of acetonitrile and magnesium perchlorate (1.0 M) as a model system, we demonstrated that this one color- and population-involving CT process makes CMDS capable of measuring samples with features that mimic excitation and emission spectral measurements in fluorescence spectroscopy. With the new capabilities, one might develop CMDS into a versatile vibrational tool for revealing the role of coherence as a design element in realizing a function. Furthermore, CT-based vibrational resonance energy transfer (VRET) methods may be developed for label-free biosensing and imaging, such as those demonstrated by fluorescence resonance energy transfer (FRET).

Multispectral multidimensional spectrometer spanning the ultraviolet to the mid-infrared

Multidimensional spectroscopy is the optical analog to nuclear magnetic resonance, probing dynamical processes with ultrafast time resolution. At optical frequencies, the technical challenges of multidimensional spectroscopy have hindered its progress until recently, where advances in laser sources and pulse-shaping have removed many obstacles to its implementation. Multidimensional spectroscopy in the visible and infrared (IR) regimes has already enabled respective advances in our understanding of photosynthesis and the structural rearrangements of liquid water. A frontier of ultrafast spectroscopy is to extend and combine multidimensional techniques and frequency ranges, which have been largely restricted to operating in the distinct visible or IR regimes. By employing two independent amplifiers seeded by a single oscillator, it is straightforward to span a wide range of time scales (femtoseconds to seconds), all of which are often relevant to the most important energy conversion and catalysis problems in chemistry, physics, and materials science. Complex condensed phase systems have optical transitions spanning the ultraviolet (UV) to the IR and exhibit dynamics relevant to function on time scales of femtoseconds to seconds and beyond. We describe the development of the Multispectral Multidimensional Nonlinear Spectrometer (MMDS) to enable studies of dynamical processes in atomic, molecular, and material systems spanning femtoseconds to seconds, from the UV to the IR regimes. The MMDS employs pulse-shaping methods to provide an easy-to-use instrument with an unprecedented spectral range that enables unique combination spectroscopies. We demonstrate the multispectral capabilities of the MMDS on several model systems.

2D THz-THz-Raman Photon-Echo Spectroscopy of Molecular Vibrations in Liquid Bromoform

Fundamental properties of molecular liquids are governed by long-range interactions that most prominently manifest at terahertz (THz) frequencies. Here we report the detection of nonlinear THz photon-echo (rephasing) signals in liquid bromoform using THz-THz-Raman spectroscopy. Together, the many observed signatures span frequencies from 0.5 to 8.5 THz and result from couplings between thermally populated ladders of vibrational states. The strongest peaks in the spectrum are found to be multiquantum dipole and 1-quantum polarizability transitions and may arise from nonlinearities in the intramolecular dipole moment surface driven by intermolecular interactions.

Two-Dimensional Vibrational Spectroscopy of a Dissipative System with the Optimized Mean-Trajectory Approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions representing radiation-matter interactions. Here we apply this method to an anharmonic chromophore coupled to a harmonic bath. A forward-backward trajectory implementation of the OMT method is described that addresses the numerical challenges of applying the OMT to large systems with disparate frequency scales. The OMT is shown to well reproduce line shapes and waiting time dynamics in the pure dephasing limit of weak coupling to an off-resonant bath. The OMT is also shown to describe a case where energy transfer is the predominant source of line broadening.

Vibrational coherence and energy transfer in two-dimensional spectra with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions that represent radiation-matter interactions. Here, we extend the OMT to include additional vibrational coherence and energy transfer processes. This generalized approximation is applied to a pair of anharmonic chromophores coupled to a bath. The resulting 2D spectra are shown to reflect coherence transfer between normal modes.

Direct observation of the coherent nuclear response after the absorption of a photon

How molecules convert light energy to perform a specific transformation is a fundamental question in photophysics. Ultrafast spectroscopy reveals the kinetics associated with electronic energy flow, but little is known about how absorbed photon energy drives nuclear motion. Here we used ultrabroadband transient absorption spectroscopy to monitor coherent vibrational energy flow after photoexcitation of the retinal chromophore. In the proton pump bacteriorhodopsin, we observed coherent activation of hydrogen-out-of-plane wagging and backbone torsional modes that were replaced by unreactive coordinates in the solution environment, concomitant with a deactivation of the reactive relaxation pathway.

Two-dimensional spectroscopy of coupled vibrations with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical representation of the nonlinear vibrational response function used to compute multidimensional infrared spectra. In this method, response functions are calculated from a sequence of classical trajectories linked by discontinuities representing the effects of radiation-matter interactions, thus providing an approximation to quantum dynamics using classical inputs. This approach was previously formulated and assessed numerically for a single anharmonic degree of freedom. Our previous work is generalized here in two respects. First, the derivation of the OMT is extended to any number of coupled anharmonic vibrations by determining semiclassical approximations for pairs of double-sided Feynman diagrams. Second, an efficient numerical procedure is developed for calculating two-dimensional infrared spectra of coupled anharmonic vibrations in the OMT approximation. The OMT approximation is shown to reproduce the fundamental features of the quantum response function including both coherence and population dynamics.

Quantum process tomography quantifies coherence transfer dynamics in vibrational exciton

Quantum coherence has been a subject of great interest in many scientific disciplines. However, detailed characterization of the quantum coherence in molecular systems, especially its transfer and relaxation mechanisms, still remains a major challenge. The difficulties arise in part because the spectroscopic signatures of the coherence transfer are typically overwhelmed by other excitation-relaxation processes. We use quantum process tomography (QPT) via two-dimensional infrared spectroscopy to quantify the rate of the elusive coherence transfer between two vibrational exciton states. QPT retrieves the dynamics of the dissipative quantum system directly from the experimental observables. It thus serves as an experimental alternative to theoretical models of the system-bath interaction and can be used to validate these theories. Our results for coupled carbonyl groups of a diketone molecule in chloroform, used as a benchmark system, reveal the nonsecular nature of the interaction between the exciton and the Markovian bath and open the door for the systematic studies of the dissipative quantum systems dynamics in detail.

Determining equilibrium fluctuations using temperature-dependent 2D-IR

We demonstrate the capability of temperature-dependent 2D-IR to characterize sources of vibrational population transfer. In a model system of iron diene tricarbonyl "piano stool" complexes, this approach reveals symmetry breaking associated with equilibrium fluctuations and differentiates these from fluxional rearrangement. Tricarbonyl(1,3-butadiene)iron and tricarbonyl(1,5-cyclooctadiene)iron are shown to undergo intramolecular vibrational redistribution (IVR) coupled to the wagging motion of their carbonyl ligands. In the case of both molecules, these equilibrium fluctuations are distinguished from chemical exchange behaviors by their temperature dependence and arguments of molecular symmetry.

Quantum Beats and Coherence Decay in Degenerate States Split by Solvation

Coherent dynamics of degenerate quantum states symmetry-broken on the femtosecond timescale is found to exhibit the phenomenon of quantum beats. Frequency-resolved and polarization-selective heterodyned transient grating spectroscopy enabled us to retrieve the oscillation pattern characteristic of the beating in systems undergoing ultrafast dynamical processes. This methodology applies to the general phenomena of coherence dynamics which is important in any ultrafast multidimensional spectroscopy. A particular application to the vibrational spectroscopy of coherence in the degenerate normal modes of the tricyanomethanide anion solvated in water is explored in this study. The relaxation of the cross-polarization transient grating anisotropy is shown to reflect the loss of the vibrational coherence, which is caused by ultrafast dynamics of the water solvation shell.

Photoinduced Charge Transfer from Titania to Surface Doping Site

We evaluate a theoretical model in which Ru is substituting for Ti at the (100) surface of anatase TiO2. Charge transfer from the photo-excited TiO2 substrate to the catalytic site triggers the photo-catalytic event (such as water oxidation or reduction half-reaction). We perform ab-initio computational modeling of the charge transfer dynamics on the interface of TiO2 nanorod and catalytic site. A slab of TiO2 represents a fragment of TiO2 nanorod in the anatase phase. Titanium to ruthenium replacement is performed in a way to match the symmetry of TiO2 substrate. One molecular layer of adsorbed water is taken into consideration to mimic the experimental conditions. It is found that these adsorbed water molecules saturate dangling surface bonds and drastically affect the electronic properties of systems investigated. The modeling is performed by reduced density matrix method in the basis of Kohn-Sham orbitals. A nano-catalyst modeled through replacement defect contributes energy levels near the bottom of the conduction band of TiO2 nano-structure. An exciton in the nano-rod is dissipating due to interaction with lattice vibrations, treated through non-adiabatic coupling. The electron relaxes to conduction band edge and then to the Ru cite with faster rate than hole relaxes to the Ru cite. These results are of the importance for an optimal design of nano-materials for photo-catalytic water splitting and solar energy harvesting.

An optimized semiclassical approximation for vibrational response functions

The observables of multidimensional infrared spectroscopy may be calculated from nonlinear vibrational response functions. Fully quantum dynamical calculations of vibrational response functions are generally impractical, while completely classical calculations are qualitatively incorrect at long times. These challenges motivate the development of semiclassical approximations to quantum mechanics, which use classical mechanical information to reconstruct quantum effects. The mean-trajectory (MT) approximation is a semiclassical approach to quantum vibrational response functions employing classical trajectories linked by deterministic transitions representing the effects of the radiation-matter interaction. Previous application of the MT approximation to the third-order response function R(3)(t3, t2, t1) demonstrated that the method quantitatively describes the coherence dynamics of the t3 and t1 evolution times, but is qualitatively incorrect for the waiting-time t2 period. Here we develop an optimized version of the MT approximation by elucidating the connection between this semiclassical approach and the double-sided Feynman diagrams (2FD) that represent the quantum response. Establishing the direct connection between 2FD and semiclassical paths motivates a systematic derivation of an optimized MT approximation (OMT). The OMT uses classical mechanical inputs to accurately reproduce quantum dynamics associated with all three propagation times of the third-order vibrational response function.