Quantum dynamics of multi-dimensional rhodopsin photoisomerization models: Approximate versus accurate treatment of the secondary modes

Tensor-Train Split-Operator KSL (TT-SOKSL) Method for Quantum Dynamics Simulations

Numerically exact simulations of quantum reaction dynamics, including nonadiabatic effects in excited electronic states, are essential to gain fundamental insights into ultrafast chemical reactivity and rigorous interpretations of molecular spectroscopy. Here, we introduce the tensor-train split-operator KSL (TT-SOKSL) method for quantum simulations in tensor-train (TT)/matrix product state (MPS) representations. TT-SOKSL propagates the quantum state as a tensor train using the Trotter expansion of the time-evolution operator, as in the tensor-train split-operator Fourier transform (TT-SOFT) method. However, the exponential operators of the Trotter expansion are applied using a rank-adaptive TT-KSL scheme instead of using the scaling and squaring approach as in TT-SOFT. We demonstrate the accuracy and efficiency of TT-SOKSL as applied to simulations of the photoisomerization of the retinal chromophore in rhodopsin, including nonadiabatic dynamics at a conical intersection of potential energy surfaces. The quantum evolution is described in full dimensionality by a time-dependent wavepacket evolving according to a two-state 25-dimensional model Hamiltonian. We find that TT-SOKSL converges faster than TT-SOFT with respect to the maximally allowed memory requirement of the tensor-train representation and better preserves the norm of the time-evolving state. When compared to the corresponding simulations based on the TT-KSL method, TT-SOKSL has the advantage of avoiding the need to construct the matrix product state Laplacian by exploiting the linear scaling of multidimensional tensor-train Fourier transforms.

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.

Ultrafast pulse shaping modulates perceived visual brightness in living animals

Vision is usually assumed to be sensitive to the light intensity and spectrum but not to its spectral phase. However, experiments performed on retinal proteins in solution showed that the first step of vision consists in an ultrafast photoisomerization that can be coherently controlled by shaping the phase of femtosecond laser pulses, especially in the multiphoton interaction regime. The link between these experiments in solution and the biological process allowing vision was not demonstrated. Here, we measure the electric signals fired from the retina of living mice upon femtosecond multipulse and single-pulse light stimulation. Our results show that the electrophysiological signaling is sensitive to the manipulation of the light excitation on a femtosecond time scale. The mechanism relies on multiple interactions with the light pulses close to the conical intersection, like pump-dump (photoisomerization interruption) and pump-repump (reverse isomerization) processes. This interpretation is supported both experimentally and by dynamics simulations.

Computational approaches to efficient generation of the stationary state for incoherent light excitation

Light harvesting processes are often computationally studied from a time-dependent viewpoint, in line with ultrafast coherent spectroscopy experiments. Yet, natural processes take place in the presence of incoherent light, which induces a stationary state. Such stationary states can be described using the eigenbasis of the molecular Hamiltonian, but for realistic systems, a full diagonalization is prohibitively expensive. We propose three efficient computational approaches to obtain the stationary state that circumvents system Hamiltonian diagonalization. The connection between the incoherent perturbations, decoherence, and Kraus operators is established.

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.