Many-Body Approaches for Simulating Coherent Nonlinear Spectroscopies of Electronic and Vibrational Excitons

Coherence in Chemistry: Foundations and Frontiers

Coherence refers to correlations in waves. Because matter has a wave-particle nature, it is unsurprising that coherence has deep connections with the most contemporary issues in chemistry research (e.g., energy harvesting, femtosecond spectroscopy, molecular qubits and more). But what does the word "coherence" really mean in the context of molecules and other quantum systems? We provide a review of key concepts, definitions, and methodologies, surrounding coherence phenomena in chemistry, and we describe how the terms "coherence" and "quantum coherence" refer to many different phenomena in chemistry. Moreover, we show how these notions are related to the concept of an interference pattern. Coherence phenomena are indeed complex, and ambiguous definitions may spawn confusion. By describing the many definitions and contexts for coherence in the molecular sciences, we aim to enhance understanding and communication in this broad and active area of chemistry.

The development and applications of multidimensional biomolecular spectroscopy illustrated by photosynthetic light harvesting

The parallel and synergistic developments of atomic resolution structural information, new spectroscopic methods, their underpinning formalism, and the application of sophisticated theoretical methods have led to a step function change in our understanding of photosynthetic light harvesting, the process by which photosynthetic organisms collect solar energy and supply it to their reaction centers to initiate the chemistry of photosynthesis. The new spectroscopic methods, in particular multidimensional spectroscopies, have enabled a transition from recording rates of processes to focusing on mechanism. We discuss two ultrafast spectroscopies - two-dimensional electronic spectroscopy and two-dimensional electronic-vibrational spectroscopy - and illustrate their development through the lens of photosynthetic light harvesting. Both spectroscopies provide enhanced spectral resolution and, in different ways, reveal pathways of energy flow and coherent oscillations which relate to the quantum mechanical mixing of, for example, electronic excitations (excitons) and nuclear motions. The new types of information present in these spectra provoked the application of sophisticated quantum dynamical theories to describe the temporal evolution of the spectra and provide new questions for experimental investigation. While multidimensional spectroscopies have applications in many other areas of science, we feel that the investigation of photosynthetic light harvesting has had the largest influence on the development of spectroscopic and theoretical methods for the study of quantum dynamics in biology, hence the focus of this review. We conclude with key questions for the next decade of this review.

Interconnection between Polarization-Detected and Population-Detected Signals: Theoretical Results and Ab Initio Simulations

Most of spectroscopic signals are specified by the nonlinear laser-induced polarization. In recent years, population-detection of signals becomes a trend in femtosecond spectroscopy. Polarization-detected (PD) and population-detected signals are fundamentally different, because they are determined by photoinduced processes acting on disparate time scales. In this work, we consider the fluorescence-detected (FD) N-wave-mixing (NWM) signal as a representative example of population-detected signals, derive a rigorous expression for this signal, and discuss its approximate variants suitable for numerical simulations. This leads us to the definition of the phenomenological FD (PFD) signal, which contains as a special case all definitions of FD signals available in the literature. Then we formulate and prove the population-polarization equivalence (PPE) theorem, which states that PFD NWM signals produced by (possibly strong) laser pulses can be evaluated as conventional PD signals in which the effective polarization is determined by the PFD transition dipole moment operator. We use the PPE theorem for the construction of the ab initio protocol for the simulation of PFD 4WM signals. As an example, we calculate electronic two-dimensional (2D) PFD spectra of the gas-phase pyrazine and compare them with the corresponding PD 2D spectra.

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.

Parameter estimation in ultrafast spectroscopy using probability theory

Ultrafast spectroscopy is a powerful technique that utilizes short pulses on the femtosecond time scale to generate and probe coherent responses in molecular systems. While the specific ultrafast methodologies vary, the most common data analysis tools rely on discrete Fourier transformation for recovering coherences that report on electronic or vibrational states and multi-exponential fitting for probing population dynamics, such as excited-state relaxation. These analysis tools are widely used due to their perceived reliability in estimating frequencies and decay rates. Here, we demonstrate that such "black box" methods for parameter estimation often lead to inaccurate results even in the absence of noise. To address this issue, we propose an alternative approach based on Bayes probability theory that simultaneously accounts for both population and coherence contributions to the signal. This Bayesian inference method offers accurate parameter estimations across a broad range of experimental conditions, including scenarios with high noise and data truncation. In contrast to traditional methods, Bayesian inference incorporates prior information about the measured signal and noise, leading to improved accuracy. Moreover, it provides estimator error bounds, enabling a systematic statistical framework for interpreting confidence in the results. By employing Bayesian inference, all parameters of a realistic model system may be accurately recovered, even in extremely challenging scenarios where Fourier and multi-exponential fitting methods fail. This approach offers a more reliable and comprehensive analysis tool for time-resolved coherent spectroscopy, enhancing our understanding of molecular systems and enabling a better interpretation of experimental data.

Ab initio simulation of peak evolutions and beating maps for electronic two-dimensional signals of a polyatomic chromophore

By employing the doorway-window (DW) on-the-fly simulation protocol, we performed ab initio simulations of peak evolutions and beating maps of electronic two-dimensional (2D) spectra of a polyatomic molecule in the gas phase. As the system under study, we chose pyrazine, which is a paradigmatic example of photodynamics dominated by conical intersections (CIs). From the technical perspective, we demonstrate that the DW protocol is a numerically efficient methodology suitable for simulations of 2D spectra for a wide range of excitation/detection frequencies and population times. From the information content perspective, we show that peak evolutions and beating maps not only reveal timescales of transitions through CIs but also pinpoint the most relevant coupling and tuning modes active at these CIs.

The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions

This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Absolute excited state molecular geometries revealed by resonance Raman signals

Ultrafast reactions activated by light absorption are governed by multidimensional excited-state (ES) potential energy surfaces (PESs), which describe how the molecular potential varies with the nuclear coordinates. ES PESs ad-hoc displaced with respect to the ground state can drive subtle structural rearrangements, accompanying molecular biological activity and regulating physical/chemical properties. Such displacements are encoded in the Franck-Condon overlap integrals, which in turn determine the resonant Raman response. Conventional spectroscopic approaches only access their absolute value, and hence cannot determine the sense of ES displacements. Here, we introduce a two-color broadband impulsive Raman experimental scheme, to directly measure complex Raman excitation profiles along desired normal modes. The key to achieve this task is in the signal linear dependence on the Frank-Condon overlaps, brought about by non-degenerate resonant probe and off-resonant pump pulses, which ultimately enables time-domain sensitivity to the phase of the stimulated vibrational coherences. Our results provide the tool to determine the magnitude and the sensed direction of ES displacements, unambiguously relating them to the ground state eigenvectors reference frame.

Ab Initio Quasiclassical Simulation of Femtosecond Time-Resolved Two-Dimensional Electronic Spectra of Pyrazine

Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.

The configuration effect on the exciton dynamics of zinc chlorin aggregates

Excitonic energy transfer among the zinc chlorin molecules is significant for the photovoltaic process because of their high sensitivities to harvesting sunlight. Zinc chlorin monomers and dimers can be synthesized experimentally, and they can form various self-assembled structures. Using the realistic parameters of zinc chlorin molecules, we assume that 20 molecules with J-, H- or J-H aggregation are arranged in a line and we investigate their dipole configuration effect on exciton dynamics. The expectation value approximation of operators is applied to derive the equations of motion of multi-exciton states. The temporal evolution of multi-exciton states is analyzed in the scheme of density matrix theory. Our simulations show that the inter-molecular coupling results in an exciton band and the wave-packet progressing excited by the resonant laser pulse exhibits attractive or repulsive behavior at the exciton level due to the dipole configuration effect. In the defined J-H coupling, the coherent wave-packet cannot overcome the configuration barrier to the no-excited part. The exciton dynamics revealed here might be helpful to better understand the energy transfer process in organic photovoltaic devices.

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.

Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.

The Relationship between the Spatial Arrangement of Pigments and Exciton Transition Moments in Photosynthetic Light-Harvesting Complexes

Considering bacteriochlorophyll molecules embedded in the protein matrix of the light-harvesting complexes of purple bacteria (known as LH2 and LH1-RC) as examples of systems of interacting pigment molecules, we investigated the relationship between the spatial arrangement of the pigments and their exciton transition moments. Based on the recently reported crystal structures of LH2 and LH1-RC and the outcomes of previous theoretical studies, as well as adopting the Frenkel exciton Hamiltonian for two-level molecules, we performed visualizations of the LH2 and LH1 exciton transition moments. To make the electron transition moments in the exciton representation invariant with respect to the position of the system in space, a system of pigments must be translated to the center of mass before starting the calculations. As a result, the visualization of the transition moments for LH2 provided the following pattern: two strong transitions were outside of LH2 and the other two were perpendicular and at the center of LH2. The antenna of LH1-RC was characterized as having the same location of the strongest moments in the center of the complex, exactly as in the B850 ring, which actually coincides with the RC. Considering LH2 and LH1 as supermolecules, each of which has excitation energies and corresponding transition moments, we propose that the outer transitions of LH2 can be important for inter-complex energy exchange, while the inner transitions keep the energy in the complex; moreover, in the case of LH1, the inner transitions increased the rate of antenna-to-RC energy transfer.

Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex

Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead?

A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.

Effects of high pulse intensity and chirp in two-dimensional electronic spectroscopy of an atomic vapor

The effects of high pulse intensity and chirp on two-dimensional electronic spectroscopy signals are experimentally investigated in the highly non-perturbative regime using atomic rubidium vapor as clean model system. Data analysis is performed based on higher-order Feynman diagrams and non-perturbative numerical simulations of the system response. It is shown that higher-order contributions may lead to a fundamental change of the static appearance and beating-maps of the 2D spectra and that chirped pulses enhance or suppress distinct higher-order pathways. We further give an estimate of the threshold intensity beyond which the high-intensity effects become visible for the system under consideration.

Distinguishing Energy- and Charge-Transfer Processes in Layered Perovskite Quantum Wells with Two-Dimensional Action Spectroscopies

Interest in photovoltaic devices based on layered perovskites is motivated by their tunable optoelectronic properties and stabilities in humid conditions. In these systems, quantum wells with different sizes are organized to direct energy and charge transport between electrodes; however, these relaxation mechanisms are difficult to distinguish based on conventional transient absorption techniques. Here, two-dimensional "action spectroscopies" are employed to separately target processes that lead to the production of photocurrent and energy loss due to fluorescence emission. These measurements show that energy transfer between quantum wells dominates the subnanosecond time scale, whereas electron transfer occurs at later times. Overall, this study suggests that while the intense exciton transitions promote light harvesting, much of the absorbed energy is lost by way of spontaneous emission. This limitation may be overcome with alternate layered perovskite systems that combine smaller exciton binding energies with large absorbance cross sections in the visible spectral range.

Origin of Unexpectedly Simple Oscillatory Responses in the Excited-State Dynamics of Disordered Molecular Aggregates

Unraveling the many facets of coherent and incoherent exciton motion in an ensemble of chromophores is an inherently complex quantum mechanical problem that has triggered a vivid debate on the role of quantum effects in molecular materials and biophysical systems. Here the dynamics of a statistical ensemble of molecular aggregates consisting of identical chromophores is investigated within a new theoretical framework. Taking account of intrinsic properties of the system, the Hamiltonian of the aggregate is partitioned into two mutually commuting vibrational and vibronic operators. This representation paves the way for an analysis that reveals the role of static disorder in ensembles of aggregates. Using analytical methods, it is demonstrated that after a critical time τ_D ≃ 2π/σ (σ being the dispersion of the disorder) any dynamic variable of the aggregate exhibits purely vibrational dynamics. This result is confirmed by exact numerical calculations of the time-dependent site populations of the aggregate. These findings may be useful for the interpretation of recent femtosecond spectroscopic experiments on molecular aggregates.

Extreme and Topological Dissipative Solitons with Structured Matter and Structured Light

Structuring of matter with nanoobjects allows one to generate soliton-like light bundles with extreme characteristics-temporal duration and spatial dimensions. On the other hand, structuring of light gives the possibility to form light bundles with complicated internal structure; their topology could be used for information coding similar to that in self-replicating RNA molecules carrying genetic code. Here we review the both variants of structuring. In the first variant, we consider a linear molecular chain and organic film interacting resonantly with laser radiation. Demonstrated are optical bistability, switching waves, and dissipative solitons, whose sizes for molecular J-aggregates can reach the nanometer range. We also discuss some theoretical approaches to take into account multi-particle interaction and correlations between molecules. In the second variant, light structuring in large-size laser medium with saturable amplification and absorption is achieved by preparation of the initial field distribution with a number of closed and unclosed vortex lines where the field vanishes. Various types of topological solitons, parameter domains of their stability, and transformation of the solitons with slow variation of the scheme parameters are presented.

Simulation of Femtosecond Phase-Locked Double-Pump Signals of Individual Light-Harvesting Complexes LH2

Recent phase-locked femtosecond double-pump experiments on individual light-harvesting complexes LH2 of purple bacteria at ambient temperature revealed undamped oscillatory responses on a time scale of at least 400 fs [ Hildner et al. Science 2013 , 340 , 1448 ]. Using an excitonic Hamiltonian for LH2 available in the literature, we simulate these signals numerically by a method that treats excitonic couplings and exciton-phonon couplings in a nonperturbative manner. The simulations provide novel insights into the origin of coherent dynamics in individual LH2 complexes.

Expanded Theory of H- and J-Molecular Aggregates: The Effects of Vibronic Coupling and Intermolecular Charge Transfer

The electronic excited states of molecular aggregates and their photophysical signatures have long fascinated spectroscopists and theoreticians alike since the advent of Frenkel exciton theory almost 90 years ago. The influence of molecular packing on basic optical probes like absorption and photoluminescence was originally worked out by Kasha for aggregates dominated by Coulombic intermolecular interactions, eventually leading to the classification of J- and H-aggregates. This review outlines advances made in understanding the relationship between aggregate structure and photophysics when vibronic coupling and intermolecular charge transfer are incorporated. An assortment of packing geometries is considered from the humble molecular dimer to more exotic structures including linear and bent aggregates, two-dimensional herringbone and "HJ" aggregates, and chiral aggregates. The interplay between long-range Coulomb coupling and short-range charge-transfer-mediated coupling strongly depends on the aggregate architecture leading to a wide array of photophysical behaviors.

UV-Light-Induced Vibrational Coherences: The Key to Understand Kasha Rule Violation in trans-Azobenzene

We combine sub-20 fs transient absorption spectroscopy with state-of-the-art computations to study the ultrafast photoinduced dynamics of trans-azobenzene (AB). We are able to resolve the lifetime of the ππ* state, whose decay within ca. 50 fs is correlated to the buildup of the nπ* population and to the emergence of coherences in the dynamics, to date unobserved. Nonlinear spectroscopy simulations call for the CNN in-plane bendings as the active modes in the subps photoinduced coherent dynamics out of the ππ* state. Radiative to kinetic energy transfer into these modes drives the system to a high-energy planar nπ*/ground state conical intersection, inaccessible upon direct excitation of the nπ* state, that triggers an ultrafast (0.45 ps) nonproductive decay of the nπ* state and is thus responsible for the observed Kasha rule violation in UV excited trans-AB. On the other hand, cis-AB is built only after intramolecular vibrational energy redistribution and population of the NN torsional mode.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

On the Simulation of Two-dimensional Electronic Spectroscopy of Indole-containing Peptides

A benchmark study of low-cost multiconfigurational CASSCF/CASPT2 schemes for computing the electronic structure of indole is presented. This facilitates the simulation of near-ultraviolet (UV) pump visible (VIS) probe (i.e. two-color) two-dimensional electronic spectra (2DES) of homo- and hetero-aggregates as well as for processing of multiple snapshots from molecular dynamics simulations. Fingerprint excited-state absorption signatures of indole are identified in a broad spectral window between 10 and 25 k cm^-1 . The 18-24 k cm^-1 spectral window which has no absorption of the monomer and noninteracting aggregates is ideally suited to embed charge-transfer signatures in stacked aggregates. The small peptide Trp-cage, containing a tryptophan and a tyrosine amino acids, having indole and phenol as side chains, respectively, serves to prove the concept. Clear charge-transfer signatures are found in the proposed spectral window for an interchromophore distance of 5 Å making near-UV pump VIS probe 2DES a suitable technique for resolving closely packed aggregates. We demonstrate that 2DES utilizing ultra-short pulses has the potential to resolve the nature of the spectroscopically resolved electronic states and that the line shapes of the excited-state absorption signals can be correlated to the polarity of the relevant states.

Quantum electron-vibrational dynamics at finite temperature: Thermo field dynamics approach

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on the thermo field dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. A comparison with the theoretically equivalent density matrix formulation shows the key numerical advantages of the present approach. The solution of thermo field dynamics equations with a novel technique for the propagation of tensor trains (matrix product states) is discussed. Numerical applications to model spin-boson systems show that the present approach is a promising tool for the description of quantum dynamics of complex molecular systems at finite temperature.

Photocell Optimization Using Dark State Protection

Conventional photocells suffer a fundamental efficiency threshold imposed by the principle of detailed balance, reflecting the fact that good absorbers must necessarily also be fast emitters. This limitation can be overcome by "parking" the energy of an absorbed photon in a dark state which neither absorbs nor emits light. Here we argue that suitable dark states occur naturally as a consequence of the dipole-dipole interaction between two proximal optical dipoles for a wide range of realistic molecular dimers. We develop an intuitive model of a photocell comprising two light-absorbing molecules coupled to an idealized reaction center, showing asymmetric dimers are capable of providing a significant enhancement of light-to-current conversion under ambient conditions. We conclude by describing a road map for identifying suitable molecular dimers for demonstrating this effect by screening a very large set of possible candidate molecules.

Nonlinear Optical Materials for the Smart Filtering of Optical Radiation

The control of luminous radiation has extremely important implications for modern and future technologies as well as in medicine. In this Review, we detail chemical structures and their relevant photophysical features for various groups of materials, including organic dyes such as metalloporphyrins and metallophthalocyanines (and derivatives), other common organic materials, mixed metal complexes and clusters, fullerenes, dendrimeric nanocomposites, polymeric materials (organic and/or inorganic), inorganic semiconductors, and other nanoscopic materials, utilized or potentially useful for the realization of devices able to filter in a smart way an external radiation. The concept of smart is referred to the characteristic of those materials that are capable to filter the radiation in a dynamic way without the need of an ancillary system for the activation of the required transmission change. In particular, this Review gives emphasis to the nonlinear optical properties of photoactive materials for the function of optical power limiting. All known mechanisms of optical limiting have been analyzed and discussed for the different types of materials.