Coherent control from quantum commitment probabilities

We introduce a general definition of a quantum committor in order to clarify reaction mechanisms and facilitate control in processes where coherent effects are important. With a quantum committor, we generalize the notion of a transition state to quantum superpositions and quantify the effect of interference on the progress of the reaction. The formalism is applicable to any linear quantum master equation supporting metastability for which absorbing boundary conditions designating the reactant and product states can be applied. We use this formalism to determine the dependence of the quantum transition state on coherences in a polaritonic system and optimize the initialization state of a conical intersection model to control reactive outcomes, achieving yields of the desired state approaching 100%. In addition to providing a practical tool, the quantum committor provides a conceptual framework for understanding reactions in cases when classical intuitions fail.

Managing temperature in open quantum systems strongly coupled with structured environments

In non-perturbative non-Markovian open quantum systems, reaching either low temperatures with the hierarchical equations of motion (HEOM) or high temperatures with the Thermalized Time Evolving Density Operator with Orthogonal Polynomials Algorithm (T-TEDOPA) formalism in Hilbert space remains challenging. We compare different ways of modeling the environment. Sampling the Fourier transform of the bath correlation function, also called temperature dependent spectral density, proves to be very effective. T-TEDOPA [Tamascelli et al., Phys. Rev. Lett. 123, 090402 (2019)] uses a linear chain of oscillators with positive and negative frequencies, while HEOM is based on the complex poles of an optimized rational decomposition of the temperature dependent spectral density [Xu et al., Phys. Rev. Lett. 129, 230601 (2022)]. Resorting to the poles of the temperature independent spectral density and of the Bose function separately is an alternative when the problem due to the huge number of Bose poles at low temperatures is circumvented. Two examples illustrate the effectiveness of the HEOM and T-TEDOPA approaches: a benchmark pure dephasing case and a two-bath model simulating the dynamics of excited electronic states coupled through a conical intersection. We show the efficiency of T-TEDOPA to simulate dynamics at a finite temperature by using either continuous spectral densities or only all the intramolecular oscillators of a linear vibronic model calibrated from ab initio data of a phenylene ethynylene dimer.

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.

Nonadiabatic transition paths from quantum jump trajectories

We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. We detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.

Steady State Photoisomerization Quantum Yield of Model Rhodopsin: Insights from Wavepacket Dynamics?

We simulate the nonequilibrium steady state cis-trans photoisomerization of retinal chromophore in rhodopsin on the basis of a two-state, two-mode model coupled to a thermal environment. By analyzing the systematic trends within an inhomogeneously broadened ensemble of systems, we find that the steady state reaction quantum yield (QY) correlates strongly with the excess energy above the crossing point of the system, in agreement with the prediction of the short-time dynamical wavepacket picture. However, the nontrivial dependence of the QY on the system-environment interaction indicates that a pure dynamical picture is insufficient and that environment-induced partial internal energy redistribution takes place before the reaction concludes. These results imply that a proper treatment of the photoisomerization reaction, particularly its high QY, must account for the redistribution and dissipation of energy beyond the dynamical wavepacket motion that is typically employed in the literature and that is appropriate only in the transient regime.

The Energy Transfer Yield between Carotenoids and Chlorophylls in Peridinin Chlorophyll a Protein Is Robust against Mutations

The energy transfer (ET) from carotenoids (Cars) to chlorophylls (Chls) in photosynthetic complexes occurs with almost unitary efficiency thanks to the synergistic action of multiple finely tuned channels whose photophysics and dynamics are not fully elucidated yet. We investigated the energy flow from the Car peridinin (Per) to Chl a in the peridinin chlorophyll a protein (PCP) from marine algae Amphidinium carterae by using two-dimensional electronic spectroscopy (2DES) with a 10 fs temporal resolution. Recently debated hypotheses regarding the S₂-to-S₁ relaxation of the Car via a conical intersection and the involvement of possible intermediate states in the ET were examined. The comparison with an N89L mutant carrying the Per donor in a lower-polarity environment helped us unveil relevant details on the mechanisms through which excitation was transferred: the ET yield was conserved even when a mutation perturbed the optimization of the system thanks to the coexistence of multiple channels exploited during the process.

The fate of molecular excited states: modeling donor-acceptor dyes

We investigate the relaxation of a coherently excited molecule in the Redfield approximation. The molecular model, parametrized to describe donor-acceptor dyes that represent a large family of molecules of interest...

Funneling dynamics in a phenylacetylene trimer: Coherent excitation of donor excitonic states and their superposition

Funneling dynamics in conjugated dendrimers has raised great interest in the context of artificial light-harvesting processes. Photoinduced relaxation has been explored by time-resolved spectroscopy and simulations, mainly by semiclassical approaches or referring to open quantum systems methods, within the Redfield approximation. Here, we take the benefit of an ab initio investigation of a phenylacetylene trimer, and in the spirit of a divide-and-conquer approach, we focus on the early dynamics of the hierarchy of interactions. We build a simplified but realistic model by retaining only bright electronic states and selecting the vibrational domain expected to play the dominant role for timescales shorter than 500 fs. We specifically analyze the role of the in-plane high-frequency skeletal vibrational modes involving the triple bonds. Open quantum system non-adiabatic dynamics involving conical intersections is conducted by separating the electronic subsystem from the high-frequency tuning and coupling vibrational baths. This partition is implemented within a robust non-perturbative and non-Markovian method, here the hierarchical equations of motion. We will more precisely analyze the coherent preparation of donor states or of their superposition by short laser pulses with different polarizations. In particular, we extend the π-pulse strategy for the creation of a superposition to a V-type system. We study the relaxation induced by the high-frequency vibrational collective modes and the transitory dissymmetry, which results from the creation of a superposition of electronic donor states.

Efficient Steady-State Solver for the Hierarchical Equations of Motion Approach: Formulation and Application to Charge Transport through Nanosystems

An iterative approach is introduced, which allows the efficient solution of the hierarchical equations of motion (HEOM) for the steady-state of open quantum systems. The approach combines the method of matrix equations with an efficient preconditioning technique to reduce the numerical effort of solving the HEOM. Illustrative applications to simulate nonequilibrium charge transport in single-molecule junctions demonstrate the performance of the method.

Exploring the Many-Body Dynamics Near a Conical Intersection with Trapped Rydberg Ions

Conical intersections between electronic potential energy surfaces are paradigmatic for the study of nonadiabatic processes in the excited states of large molecules. However, since the corresponding dynamics occurs on a femtosecond timescale, their investigation remains challenging and requires ultrafast spectroscopy techniques. We demonstrate that trapped Rydberg ions are a platform to engineer conical intersections and to simulate their ensuing dynamics on larger length scales and timescales of the order of nanometers and microseconds, respectively; all this in a highly controllable system. Here, the shape of the potential energy surfaces and the position of the conical intersection can be tuned thanks to the interplay between the high polarizability and the strong dipolar exchange interactions of Rydberg ions. We study how the presence of a conical intersection affects both the nuclear and electronic dynamics demonstrating, in particular, how it results in the inhibition of the nuclear motion. These effects can be monitored in real time via a direct spectroscopic measurement of the electronic populations in a state-of-the-art experimental setup.

Extreme Parametric Sensitivity in the Steady-State Photoisomerization of Two-Dimensional Model Rhodopsin

We computationally studied the photoisomerization reaction of the retinal chromophore in rhodopsin using a two-state two-mode model coupled to thermal baths. Reaction quantum yields at the steady state (10 ps and beyond) were found to be considerably different than their transient values, suggesting a weak correlation between transient and steady-state dynamics in these systems. Significantly, the steady-state quantum yield was highly sensitive to minute changes in system parameters, while transient dynamics was nearly unaffected. Correlation of such sensitivity with standard level spacing statistics of the nonadiabatic vibronic system suggests a possible origin in quantum chaos. The significance of this observation of quantum yield parametric sensitivity in biological models of vision has profound conceptual and fundamental implications.

Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives

We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

Photoisomerization-coupled electron transfer

Photochromic molecular structures constitute a unique platform for constructing molecular switches, sensors, and memory devices. One of their most promising applications is as light-switchable electron acceptor or donor units. Here, we investigate a previously unexplored process that we postulate may occur in such systems: an ultrafast electron transfer triggered by a simultaneous photoisomerization of the donor or the acceptor moiety. We propose a theoretical model for this phenomenon and, with the aid of density functional theory calculations, apply it to the case of a dihydropyrene-type photochromic molecular donor. By considering the wavepacket dynamics and the photoisomerization yield, we show that the two processes involved, electron transfer and photoisomerization, are in general inseparable and need to be treated in a unified manner. We finish by discussing how the efficiency of photoisomerization-coupled electron transfer can be controlled experimentally.

Cooperative Conical Intersection Dynamics of Two Pyrazine Molecules in an Optical Cavity

Hybrid light-matter states in optical cavities, known as polaritons, offer a novel means of manipulating and controlling photochemical processes. We investigate the cooperative cavity photochemistry of two pyrazine molecules undergoing conical intersection dynamics and interacting with a single cavity photon mode by exact quantum dynamics. When the cavity mode is coupled to the electronic transition between the ground and excited states, we find an enhanced polaritonic splitting and collective dark states. These features dominate the cooperative polariton dynamics and can be observed in the transient absorption spectrum.

Manipulating nonadiabatic conical intersection dynamics by optical cavities

Optical cavities hold great promise to manipulate and control the photochemistry of molecules. We demonstrate how molecular photochemical processes can be manipulated by strong light-matter coupling. For a molecule with an inherent conical intersection, optical cavities can induce significant changes in the nonadiabatic dynamics by either splitting the pristine conical intersections into two novel polaritonic conical intersections or by creating light-induced avoided crossings in the polaritonic surfaces. This is demonstrated by exact real-time quantum dynamics simulations of a three-state two-mode model of pyrazine strongly coupled to a single cavity photon mode. We further explore the effects of external environments through dissipative polaritonic dynamics computed using the hierarchical equation of motion method. We find that cavity-controlled photochemistry can be immune to external environments. We also demonstrate that the polariton-induced changes in the dynamics can be monitored by transient absorption spectroscopy.

Mapping of Wave Packet Dynamics at Conical Intersections by Time- and Frequency-Resolved Fluorescence Spectroscopy: A Computational Study

Monitoring of wave packet dynamics at conical intersections by time- and frequency-resolved fluorescence spectroscopy has been investigated theoretically for a three-state two-mode model of a conical intersection coupled to a dissipative environment. The ideal and the actually measurable time- and frequency-gated fluorescence spectra are accurately and efficiently simulated by combining the hierarchy equations-of-motion method for dissipative quantum dynamics with the methodology of the equation-of-motion phase-matching approach for the calculation of spectroscopic signals. It is shown that time- and frequency-resolved fluorescence spectra reveal essential aspects of the wave packet dynamics at conical intersections and the effects of environment-induced dissipation. The results of the present work indicate that fluorescence up-conversion spectroscopy with femtosecond time resolution is an efficient tool for the characterization of ultrafast dynamics at conical intersections.

Orientational relaxation of a quantum linear rotor in a dissipative environment: Simulations with the hierarchical equations-of-motion method

We study the effect of a dissipative environment on the orientational relaxation of a three-dimensional quantum linear rotor. We provide a derivation of the Hamiltonian of a linear rotor coupled to a harmonic bath from first principles, confirming earlier conjectures. The dynamics generated by this Hamiltonian is investigated by the hierarchical equations-of-motion method assuming a Drude spectral density of the bath. We perform numerically accurate simulations and analyze the behavior of orientational correlation functions and the rotational structures of infrared absorption and Raman scattering spectra. We explore the features of orientational correlation functions and their spectra for a wide range of system-bath couplings, bath memory times, and temperatures. We discuss the signatures of the orientational relaxation in the underdamped regime, the strongly damped regime, and the librational regime. We show that the behavior of orientational correlation functions and their spectra can conveniently be analyzed in terms of three characteristic times, which are explicitly expressed in terms of the parameters of the Hamiltonian.

An efficient approach to the quantum dynamics and rates of processes induced by natural incoherent light

In many important cases, the rate of excitation of a system embedded in an environment is significantly smaller than the internal system relaxation rates. An important example is that of light-induced processes under natural conditions, in which the system is excited by weak, incoherent (e.g., solar) radiation. Simulating the dynamics on the time scale of the excitation source can thus be computationally intractable. Here we describe a method for obtaining the dynamics of quantum systems without directly solving the master equation. We present an algorithm for the numerical implementation of this method and, as an example, use it to reconstruct the internal conversion dynamics of pyrazine excited by sunlight. Significantly, this approach also allows us to assess the role of quantum coherence on biological time scales, which is a topic of ongoing interest.