Optimal control of ultrafastcis-transphotoisomerization of retinal in rhodopsin via a conical intersection

XFEL Beamline Optical Instrumentation for Ultrafast Science

Free electron lasers operating in the soft and hard X-ray regime provide capabilities for ultrafast science in many areas, including X-ray spectroscopy, diffractive imaging, solution and material scattering, and X-ray crystallography. Ultrafast time-resolved applications in the picosecond, femtosecond, and attosecond regimes are often possible using single-shot experimental configurations. Aside from X-ray pump and X-ray probe measurements, all other types of ultrafast experiments require the synchronized operation of pulsed laser excitation for resonant or nonresonant pumping. This Perspective focuses on the opportunities for the optical control of structural dynamics by applying techniques from nonlinear spectroscopy to ultrafast X-ray experiments. This typically requires the synthesis of two or more optical pulses with full control of pulse and interpulse parameters. To this end, full characterization of the femtosecond optical pulses is also highly desirable. It has recently been shown that two-color and two-pulse femtosecond excitation of fluorescent protein crystals allowed a Tannor-Rice coherent control experiment, performed under characterized conditions. Pulse shaping and the ability to synthesize multicolor and multipulse conditions are highly desirable and would enable XFEL facilities to offer capabilities for structural dynamics. This Perspective will give a summary of examples of the types of experiments that could be achieved, and it will additionally summarize the laser, pulse shaping, and characterization that would be recommended as standard equipment for time-resolved XFEL beamlines, with an emphasis on ultrafast time-resolved serial femtosecond crystallography.

Optical Control of Crossing the Conical Intersection in β-Carotene

Optical control of dynamic processes has been challenging yet has only been demonstrated in several chemical and biological systems. The control of a reaction passing the widely present conical intersection has not been realized. Here, we modulated the phase of the excitation pulse to control the dynamics of β-carotene through accessing the conical intersection (CI). We observed different dynamics in 110-220 fs into the CI and the consecutive process in 400-600 fs through another CI by various chirped excitation pulses. We successfully controlled those ultrafast wavepacket dynamics passing the CIs on the femtosecond time scales. The method developed here can be used to control a various of ultrafast chemical and biological reactions through the CI(s).

Influence of Mode-Specific Excitation on the Nonadiabatic Dynamics of Methyl Nitrate (CH3ONO2)

The impact of mode-specific vibrational excitations on initial-preparation conditions was studied by examining the excited-state population decay rates in the nonadiabatic dynamics of methyl nitrate (CH₃ONO₂). In particular, exciting a few specific modes by adding a single quantum of energy clearly decelerated the nonadiabatic dynamics population decay rates. The underlying reason for this slower population decay was explained by analyzing the profiles of the excited-state potential energy surfaces in the Franck-Condon regions and the topology of the S₁/S₀ conical intersection. This study not only provides physical insights into the key mechanisms controlling nonadiabatic dynamics but also shows the possibility of controlling nonadiabatic dynamics via mode-specific vibrational excitations.

Time-resolved X-ray and XUV based spectroscopic methods for nonadiabatic processes in photochemistry

The photochemistry of numerous molecular systems is influenced by conical intersections (CIs). These omnipresent nonadiabatic phenomena provide ultra-fast radiationless relaxation channels by creating degeneracies between electronic states and decide over the final photoproducts. In their presence, the Born-Oppenheimer approximation breaks down, and the timescales of the electron and nuclear dynamics become comparable. Due to the ultra-fast dynamics and the complex interplay between nuclear and electronic degrees of freedom, the direct experimental observation of nonadiabatic processes close to CIs remains challenging. In this article, we give a theoretical perspective on novel spectroscopic techniques capable of observing clear signatures of CIs. We discuss methods that are based on ultra-short laser pulses in the extreme ultraviolet and X-ray regime, as their spectral and temporal resolution allow for resolving the ultra-fast dynamics near CIs.

Waveform control of molecular dynamics close to a conical intersection

Conical intersections are ubiquitous in chemical systems but, nevertheless, extraordinary points on the molecular potential energy landscape. They provide ultra-fast radiationless relaxation channels, their topography influences the product branching, and they equalize the timescales of the electron and nuclear dynamics. These properties reveal optical control possibilities in the few femtosecond regime. In this theoretical study, we aim to explore control options that rely on the carrier envelope phase of a few-cycle IR pulse. The laser interaction creates an electronic superposition just before the wave packet reaches the conical intersection. The imprinted phase information is varied by the carrier envelope phase to influence the branching ratio after the conical intersection. We test and analyze this scenario in detail for a model system and show to what extent it is possible to transfer this type of control to a realistic system like uracil.

On Weak-Field (One-Photon) Coherent Control of Photoisomerization

Photochemistry induced by phase-coherent laser light is an intriguing topic. The possibility of weak-field (one-photon) phase-only control of photoisomerization is controversial. Experimental studies on the weak-field coherent control of cis-trans isomerization have led to conflicting results. Here we address this issue by quantum dynamical calculations, focusing on a mechanism where different "phase-shaped" wave packets are quickly stabilized ("dumped") in the trans configuration because of prompt energy dissipation. We systematically investigate different relaxation rates with the system-bath dynamics described within the time-dependent Hartree approximation leading to a friction-type force. We find evidence for phase control of trans-isomer yields (about 5%) in this model with pure energy dissipation given sufficiently strong dampening.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

Uncovering the Functional Role of Coherent Phonons during the Photoinduced Phase Transition in a Molecular Crystal

The atomic motions that make up phonons and molecular vibrations in molecular crystals influence their photophysical and electronic properties, including polaron formation, carrier mobility, and phase transitions. Discriminating between spectator and driving motions is a significant challenge hindering optimization. Unlocking this information and developing fine-tuned controls over actively participating phonon modes would not only lead to a stronger understanding of photochemistry but also provide a significant new tool in controlling solid state chemistry. We present a strategy using rationally designed double pulses to unveil the unique function of specific excited state phonon modes. Using ultrafast spectroscopy, we identified 50 and 90 cm^-1 phonons involved in modulating the photoinduced spin-Peierls melting of potassium tetracyanoquinodimethane crystals. We show that the 50 cm^-1 phonon specifically corresponds to the coherent nuclear wavepacket involved in the charge transfer component of the overall spin-Peierls phase melting process, while the 90 cm^-1 phonon facilitates the phase transition component.

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.

An efficient spectral method for numerical time-dependent perturbation theory

We develop the Fourier-Laplace Inversion of the Perturbation Theory (FLIPT), a novel numerically exact "black box" method to compute perturbative expansions of the density matrix with rigorous convergence conditions. Specifically, the FLIPT method is extremely well-suited to simulate multiphoton pulsed laser experiments with complex pulse shapes. The n-dimensional frequency integrals of the nth order perturbative expansion are evaluated numerically using tensor products. The N-point discretized integrals are computed in O(N²) operations, a significant improvement over the O(Nⁿ) scaling of standard quadrature methods.

Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.

Control of Nuclear Dynamics through Conical Intersections and Electronic Coherences

The effect of nuclear dynamics and conical intersections on electronic coherences is investigated employing a two-state, two-mode linear vibronic coupling model. Exact quantum dynamical calculations are performed using the multiconfiguration time-dependent Hartree method. It is found that the presence of a nonadiabatic coupling close to the Franck-Condon point can preserve electronic coherence to some extent. Additionally, the possibility of steering the nuclear wave packets by imprinting a relative phase between the electronic states during the photoionization process is discussed. It is found that the steering of nuclear wave packets is possible given that a coherent electronic wave packet embodying the phase difference passes through a conical intersection. A conical intersection close to the Franck-Condon point is thus a necessary prerequisite for control, providing a clear path towards attochemistry.

The origin of absorptive features in the two-dimensional electronic spectra of rhodopsin

A three-state three-mode model Hamiltonian reveals the origin of the absorptive features in the two-dimensional electronic spectra of rhodopsin.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

Interfering resonance as an underlying mechanism in the adaptive feedback control of radiationless transitions: Retinal isomerization

Control of molecular processes via adaptive feedback often yields highly structured laser pulses that have eluded physical explanation. By contrast, coherent control approaches propose physically transparent mechanisms but are not readily visible in experimental results. Here, an analysis of a condensed phase adaptive feedback control experiment on retinal isomerization shows that it manifests a quantum interference based coherent control mechanism: control via interfering resonances. The result promises deep insight into the physical basis for the adaptive feedback control of a broad class of bound state processes.

Light penetration-coupled photoisomerization modeling for photodeformation of diarylethene single crystal: upscaling isomerization to macroscopic deformation

Diarylethene is one of the photo-responsive materials that show rapid and reversible changes in their color/electrochemical properties and macroscopic deformations in the crystalline phase by light irradiation. Photoisomerization is the main cause of the photo reactivity of diarylethene, and we established a statistical model based on the density matrix formalism, which predicts quantitative isomerization progress as a population term. The model reflects photo-switching properties of the target molecule, which were characterized by first principle calculations, and external stimulus factors (light irradiation conditions and temperature). By merging light penetration physics with the model, we derived light penetration depth dependent isomerization progress to theoretically investigate photodeformation of single crystal. The model well reproduced in-plane shear deformation under ultraviolet light irradiation which would provide guideline for photoactuator design. In addition, the statistical model addressed crucial findings (primary stimuli and molecular design parameter for increasing the isomerization rate, external stimuli enhancing fluorescence performance) itself.

Probing the Interstate Coupling near a Conical Intersection by Optical Spectroscopy

Conical intersections are points where adiabatic potential energy surfaces cross. The interstate coupling between the potential energy surfaces plays a crucial role in many processes associated with conical intersections. Still no method exists to measure this coupling driving the chemical reactions between the potential energy surfaces involved. In this Letter, using a generic model for photoisomerization, we propose a novel experimental approach to estimate the coupling that mixes the electronic states near a conical intersection. The approach is based on analyzing the vibrational wavepacket of the reactant in the adiabatic ground and excited electronic states. The nuclear wavepacket dynamics are extracted from linear absorption and two-dimensional electronic spectroscopy. Comparing the frequencies of the coupling mode in the adiabatic ground and excited states from models with and without coupling between the potential energy surfaces suggests an experimental tool to determine the interstate coupling.

Communication: Analytical optimal pulse shapes obtained with the aid of genetic algorithms: Controlling the photoisomerization yield of retinal

We recently proposed a Quantum Optimal Control (QOC) method constrained to build pulses from analytical pulse shapes [R. D. Guerrero et al., J. Chem. Phys. 143(12), 124108 (2015)]. This approach was applied to control the dissociation channel yields of the diatomic molecule KH, considering three potential energy curves and one degree of freedom. In this work, we utilized this methodology to study the strong field control of the cis-trans photoisomerization of 11-cis retinal. This more complex system was modeled with a Hamiltonian comprising two potential energy surfaces and two degrees of freedom. The resulting optimal pulse, made of 6 linearly chirped pulses, was capable of controlling the population of the trans isomer on the ground electronic surface for nearly 200 fs. The simplicity of the pulse generated with our QOC approach offers two clear advantages: a direct analysis of the sequence of events occurring during the driven dynamics, and its reproducibility in the laboratory with current laser technologies.

Controlling S2 Population in Cyanine Dyes Using Shaped Femtosecond Pulses

Fast population transfer from higher to lower excited states occurs via internal conversion (IC) and is the basis of Kasha's rule, which states that spontaneous emission takes place from the lowest excited state of the same multiplicity. Photonic control over IC is of interest because it would allow direct influence over intramolecular nonradiative decay processes occurring in condensed phase. Here we tracked the S2 and S1 fluorescence yield for different cyanine dyes in solution as a function of linear chirp. For the cyanine dyes with polar solvation response IR144 and meso-piperidine substituted IR806, increased S2 emission was observed when using transform limited pulses, whereas chirped pulses led to increased S1 emission. The nonpolar solvated cyanine IR806, on the other hand, did not show S2 emission. A theoretical model, based on a nonperturbative solution of the equation of motion for the density matrix, is offered to explain and simulate the anomalous chirp dependence. Our findings, which depend on pulse properties beyond peak intensity, offer a photonic method to control S2 population thereby opening the door for the exploration of photochemical processes initiated from higher excited states.

Quantum Mechanical Wave Packet Dynamics at a Conical Intersection with Strong Vibrational Dissipation

We derive a reduced model for the nonadiabatic quantum dynamics of an electronic wave packet moving through a conical intersection in the presence of strong vibrational damping. Starting from the dissipative two-state two-model model, we transform the tuning and the coupling mode to the bath. The resulting quantum two-state model with two highly structured environments is solved numerically exactly in the regime of strong vibrational damping. We find negative cross peaks in the ultrafast optical 2D spectra as clear signatures of the conical intersection. They arise from secondary excitations of the wave packet after having passed through the photophysical energy funnel. This feature is in agreement with recent transient absorption measurements of rhodopsin.

Accurate nonadiabatic dynamics

This Perspective addresses the use of coupled diabatic potential energy surfaces (PESs) together with rigorous quantum dynamics in full or reduced dimensional coordinate spaces to obtain accurate solutions to problems in nonadiabatic dynamics.

Ab initio study of potential ultrafast internal conversion routes in oxybenzone, caffeic acid, and ferulic acid: implications for sunscreens

Oxybenzone (OB) and ferulic acid (FA) both find use in commercial sunscreens; caffeic acid (CA) differs from FA by virtue of an -OH group in place of a -OCH3 group on the aromatic ring. We report the results of ab initio calculations designed to explore the excited state nonradiative relaxation pathways that provide photostability to these molecules and the photoprotection they offer toward UV-A and UV-B radiation. In the case of OB, internal conversion (IC) is deduced to occur on ultrafast time scales, via a barrierless electron-driven H atom transfer pathway from the S1(1(1)nπ*) state to a conical intersection (CI) with the ground (S0) state potential energy surface (PES). The situation with respect to CA and FA is somewhat less clear-cut, with low energy CIs identified by linking excited states to the S0 state following photoexcitation and subsequent evolution along (i) a ring centered out-of-plane deformation coordinate, (ii) the E/Z isomerism coordinate and, in the case of CA, (iii) an O-H stretch coordinate. Analogy with catechol suggests that the last of these processes (if active) would lead to radical formation (and thus potential phototoxicity), encouraging a suggestion that FA might be superior to CA as a sunscreen ingredient.

Environmental stability of quantum chaotic ratchets

The transitory and stationary behavior of a quantum chaotic ratchet consisting of a biharmonic potential under the effect of different drivings in contact with a thermal environment is studied. For weak forcing and finite ℏ, we identify a strong dependence of the current on the structure of the chaotic region. Moreover, we have determined the robustness of the current against thermal fluctuations in the very weak coupling regime. In the case of strong forcing, the current is determined by the shape of a chaotic attractor. In both cases the temperature quickly stabilizes the ratchet, but in the latter it also destroys the asymmetry responsible for the current generation. Finally, applications to isomerization reactions are discussed.

Using wave-packet interferometry to monitor the external vibrational control of electronic excitation transfer

We investigate the control of electronic energy transfer in molecular dimers through the preparation of specific vibrational coherences prior to electronic excitation, and its observation by nonlinear wave-packet interferometry (nl-WPI). Laser-driven coherent nuclear motion can affect the instantaneous resonance between site-excited electronic states and thereby influence short-time electronic excitation transfer (EET). We first illustrate this control mechanism with calculations on a dimer whose constituent monomers undergo harmonic vibrations. We then consider the use of nl-WPI experiments to monitor the nuclear dynamics accompanying EET in general dimer complexes following impulsive vibrational excitation by a subresonant control pulse (or control pulse sequence). In measurements of this kind, two pairs of polarized phase-related femtosecond pulses following the control pulse generate superpositions of coherent nuclear wave packets in optically accessible electronic states. Interference contributions to the time- and frequency-integrated fluorescence signals due to overlaps among the superposed wave packets provide amplitude-level information on the nuclear and electronic dynamics. We derive the basic expression for a control-pulse-dependent nl-WPI signal. The electronic transition moments of the constituent monomers are assumed to have a fixed relative orientation, while the overall orientation of the complex is distributed isotropically. We include the limiting case of coincident arrival by pulses within each phase-related pair in which control-influenced nl-WPI reduces to a fluorescence-detected pump-probe difference experiment. Numerical calculations of pump-probe signals based on these theoretical expressions are presented in the following paper [J. D. Biggs and J. A. Cina, J. Chem. Phys. 131, 224302 (2009)].

`Making the molecular movie': first frames

Recent advances in high-intensity electron and X-ray pulsed sources now make it possible to directly observe atomic motions as they occur in barrier-crossing processes. These rare events require the structural dynamics to be triggered by femtosecond excitation pulses that prepare the system above the barrier or access new potential energy surfaces that drive the structural changes. In general, the sampling process modifies the system such that the structural probes should ideally have sufficient intensity to fully resolve structures near the single-shot limit for a given time point. New developments in both source intensity and temporal characterization of the pulsed sampling mode have made it possible to make so-called `molecular movies',i.e.measure relative atomic motions faster than collisions can blur information on correlations. Strongly driven phase transitions from thermally propagated melting to optically modified potential energy surfaces leading to ballistic phase transitions and bond stiffening are given as examples of the new insights that can be gained from an atomic level perspective of structural dynamics. The most important impact will likely be made in the fields of chemistry and biology where the central unifying concept of the transition state will come under direct observation and enable a reduction of high-dimensional complex reaction surfaces to the key reactive modes, as long mastered by Mother Nature.

Control of retinal isomerization in bacteriorhodopsin in the high-intensity regime

A learning algorithm was used to manipulate optical pulse shapes and optimize retinal isomerization in bacteriorhodopsin, for excitation levels up to 1.8 x 10(16) photons per square centimeter. Below 1/3 the maximum excitation level, the yield was not sensitive to pulse shape. Above this level the learning algorithm found that a Fourier-transform-limited (TL) pulse maximized the 13-cis population. For this optimal pulse the yield increases linearly with intensity well beyond the saturation of the first excited state. To understand these results we performed systematic searches varying the chirp and energy of the pump pulses while monitoring the isomerization yield. The results are interpreted including the influence of 1-photon and multiphoton transitions. The population dynamics in each intermediate conformation and the final branching ratio between the all-trans and 13-cis isomers are modified by changes in the pulse energy and duration.

A new quantum control scheme for multilevel systems based on effective decomposition by intense laser fields

A new quantum control scheme for general multilevel systems using intense laser fields is proposed. In the present scheme, the target subspace consisting of several quantum levels is effectively isolated by applying intense cw lasers with specific conditions. The formulation is carried out using the Green function with the help of projection operator method. Dynamics of the isolated target subspace is governed by an effective Hamiltonian. The developed scheme is applied to the quantum control of dissipative four- and five-level systems. It is clarified that the present method makes it possible not only to manipulate the coherent population dynamics but also to suppress the dissipative dynamics.

Optimal alignment control of a nonpolar molecule through nonresonant multiphoton transitions

Alignment control of an ensemble of nonpolar molecules is numerically studied by means of optimal control simulation. A nitrogen molecule that is modeled by a quantum rigid rotor is adopted. Controlled rotational wave packets are created through nonresonant optical transitions induced by polarizability coupling. Optimal pulses are designed to achieve the alignment control at a specified time in the absence/presence of external static fields in zero- and finite-temperature cases, as well as to maintain an aligned state. When maintaining an aligned state over a specified time interval is chosen as a target, the control mechanism is primarily attributed to a dynamical one. Multiple optimal solutions that lead to virtually the same control achievement are found, which are consistent with the topology of the quantum control landscape.