Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion.

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. We use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

Dynamic embedding of effective harmonic normal mode vibrations in all-atomistic energy gap fluctuations: Case study of light harvesting 2 complex.

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. We address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics.

Modified Fermi's golden rule rate expressions.

Fermi's golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Simple model calculations illustrate the utility and implications of new rate expressions.

Partially polaron-transformed quantum master equation for exciton and charge transport dynamics.

Polaron-transformed quantum master equation (PQME) offers a unified framework to describe the dynamics of quantum systems in both limits of weak and strong couplings to environmental degrees of freedom. Thus, the PQME serves as an efficient method to describe charge and exciton transfer/transport dynamics for a broad range of parameters in condensed or complex environments. However, in some cases, the polaron transformation (PT) being employed in the formulation invokes an over-relaxation of slow modes and results in premature suppression of important coherence terms. A formal framework to address this issue is developed in the present work by employing a partial PT that has smaller weights for low frequency bath modes. It is shown here that a closed form expression of a second order time-local PQME including all the inhomogeneous terms can be derived for a general form of partial PT, although more complicated than that for the full PT. All the expressions needed for numerical calculation are derived in detail. Applications to a model of a two-level system coupled to a bath of harmonic oscillators, with test calculations focused on those due to homogeneous relaxation terms, demonstrate the feasibility and the utility of the present approach.

A simple generalization of the energy gap law for nonradiative processes.

For more than 50 years, an elegant energy gap (EG) law developed by Englman and Jortner [Mol. Phys. 18, 145 (1970)] has served as a key theory to understand and model the nearly exponential dependence of nonradiative transition rates on the difference of energy between the initial and final states. This work revisits the theory, clarifies the key assumptions involved in the rate expression, and provides a generalization for the cases where the effects of temperature dependence and low-frequency modes cannot be ignored. For a specific example where the low-frequency vibrational and/or solvation responses can be modeled as an Ohmic spectral density, a simple generalization of the EG law is provided. Test calculations demonstrate that this generalized EG law brings significant improvement over the original EG law. Both the original and generalized EG laws are also compared with the stationary phase approximations developed for electron transfer theory, which suggests the possibility of a simple interpolation formula valid for any value of EG.

Theoretical investigation of non-Förster exciton transfer mechanisms in perylene diimide donor, phenylene bridge, and terrylene diimide acceptor systems.

The rates of exciton transfer within dyads of perylene diimide and terrylene diimide connected by oligophenylene bridge units have been shown to deviate significantly from those of Förster's resonance energy transfer theory, according to single molecule spectroscopy experiments. The present work provides a detailed computational and theoretical study investigating the source of such a discrepancy. Electronic spectroscopy data are calculated by time-dependent density functional theory and then compared with experimental results. Electronic couplings between the exciton donor and the acceptor are estimated based on both the transition density cube method and transition dipole approximation. These results confirm that the delocalization of the exciton to the bridge parts contributes to significant enhancement of donor-acceptor electronic coupling. Mechanistic details of exciton transfer are examined by estimating the contributions of the bridge electronic states, vibrational modes of the dyads commonly coupled to both donor and acceptor, inelastic resonance energy transfer mechanism, and dark exciton states. These analyses suggest that the contribution of common vibrational modes serves as the main source of deviation from Förster's spectral overlap expression.

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments.

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge-Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi's Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller-Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

Implications for human odor sensing revealed from the statistics of odorant-receptor interactions.

Binding of odorants to olfactory receptors (ORs) elicits downstream chemical and neural signals, which are further processed to odor perception in the brain. Recently, Mainland and colleagues have measured more than 500 pairs of odorant-OR interaction by a high-throughput screening assay method, opening a new avenue to understanding the principles of human odor coding. Here, using a recently developed minimal model for OR activation kinetics, we characterize the statistics of OR activation by odorants in terms of three empirical parameters: the half-maximum effective concentration EC50, the efficacy, and the basal activity. While the data size of odorants is still limited, the statistics offer meaningful information on the breadth and optimality of the tuning of human ORs to odorants, and allow us to relate the three parameters with the microscopic rate constants and binding affinities that define the OR activation kinetics. Despite the stochastic nature of the response expected at individual OR-odorant level, we assess that the confluence of signals in a neuron released from the multitude of ORs is effectively free of noise and deterministic with respect to changes in odorant concentration. Thus, setting a threshold to the fraction of activated OR copy number for neural spiking binarizes the electrophysiological signal of olfactory sensory neuron, thereby making an information theoretic approach a viable tool in studying the principles of odor perception.

Fourth order expressions for the electronic absorption lineshape of molecular excitons.

The line shape of electronic absorption spectroscopy reflects the information on quantum dynamical processes accompanying the electronic excitation, and its accurate description is an important component for validating theoretical models and assumptions. The present work provides detailed expressions for the absorption line shape of molecular excitons that are valid up to the fourth order of exciton-bath interactions within the quantum master equation approach. These expressions can serve as the basis for developing general and systematic methods to model the line shape for a broad class of molecular exciton systems and environments. For the bath model of linearly coupled harmonic oscillators, more detailed expressions employing the spectral densities of the bath are presented. These expressions are then tested for a linear aggregate of identical chromophores each coupled to the harmonic oscillator bath. Calculation results for a super-Ohmic spectral density with exponential cutoff demonstrate the feasibility of calculations and also offer insights into the utility and difficulty of going beyond the second order approximation.

Implausibility of the vibrational theory of olfaction.

The vibrational theory of olfaction assumes that electron transfer occurs across odorants at the active sites of odorant receptors (ORs), serving as a sensitive measure of odorant vibrational frequencies, ultimately leading to olfactory perception. A previous study reported that human subjects differentiated hydrogen/deuterium isotopomers (isomers with isotopic atoms) of the musk compound cyclopentadecanone as evidence supporting the theory. Here, we find no evidence for such differentiation at the molecular level. In fact, we find that the human musk-recognizing receptor, OR5AN1, identified using a heterologous OR expression system and robustly responding to cyclopentadecanone and muscone, fails to distinguish isotopomers of these compounds in vitro. Furthermore, the mouse (methylthio)methanethiol-recognizing receptor, MOR244-3, as well as other selected human and mouse ORs, responded similarly to normal, deuterated, and (13)C isotopomers of their respective ligands, paralleling our results with the musk receptor OR5AN1. These findings suggest that the proposed vibration theory does not apply to the human musk receptor OR5AN1, mouse thiol receptor MOR244-3, or other ORs examined. Also, contrary to the vibration theory predictions, muscone-d30 lacks the 1,380- to 1,550-cm(-1) IR bands claimed to be essential for musk odor. Furthermore, our theoretical analysis shows that the proposed electron transfer mechanism of the vibrational frequencies of odorants could be easily suppressed by quantum effects of nonodorant molecular vibrational modes. These and other concerns about electron transfer at ORs, together with our extensive experimental data, argue against the plausibility of the vibration theory.

Non-uniqueness of quantum transition state theory and general dividing surfaces in the path integral space.

Despite the fact that quantum mechanical principles do not allow the establishment of an exact quantum analogue of the classical transition state theory (TST), the development of a quantum TST (QTST) with a proper dynamical justification, while recovering the TST in the classical limit, has been a long standing theoretical challenge in chemical physics. One of the most recent efforts of this kind was put forth by Hele and Althorpe (HA) [J. Chem. Phys. 138, 084108 (2013)], which can be specified for any cyclically invariant dividing surface defined in the space of the imaginary time path integral. The present work revisits the issue of the non-uniqueness of QTST and provides a detailed theoretical analysis of HA-QTST for a general class of such path integral dividing surfaces. While we confirm that HA-QTST reproduces the result based on the ring polymer molecular dynamics (RPMD) rate theory for dividing surfaces containing only a quadratic form of low frequency Fourier modes, we find that it produces different results for those containing higher frequency imaginary time paths which accommodate greater quantum fluctuations. This result confirms the assessment made in our previous work [Jang and Voth, J. Chem. Phys. 144, 084110 (2016)] that HA-QTST does not provide a derivation of RPMD-TST in general and points to a new ambiguity of HA-QTST with respect to its justification for general cyclically invariant dividing surfaces defined in the space of imaginary time path integrals. Our analysis also offers new insights into similar path integral based QTST approaches.

Generalized quantum Fokker-Planck equation for photoinduced nonequilibrium processes with positive definiteness condition.

This work provides a detailed derivation of a generalized quantum Fokker-Planck equation (GQFPE) appropriate for photo-induced quantum dynamical processes. The path integral method pioneered by Caldeira and Leggett (CL) [Physica A 121, 587 (1983)] is extended by utilizing a nonequilibrium influence functional applicable to different baths for the ground and the excited electronic states. Both nonequilibrium and non-Markovian effects are accounted for consistently by expanding the paths in the exponents of the influence functional up to the second order with respect to time. This procedure results in approximations involving only single time integrations for the exponents of the influence functional but with additional time dependent boundary terms that have been ignored in previous works. The boundary terms complicate the derivation of a time evolution equation but do not affect position dependent physical observables or the dynamics in the steady state limit. For an effective density operator with the boundary terms factored out, a time evolution equation is derived, through short time expansion of the effective action and Gaussian integration in analytically continued complex domain of space. This leads to a compact form of the GQFPE with time dependent kernels and additional terms, which renders the resulting equation to be in the Dekker form [Phys. Rep. 80, 1 (1981)]. Major terms of the equation are analyzed for the case of Ohmic spectral density with Drude cutoff, which shows that the new GQFPE satisfies the positive definiteness condition in medium to high temperature limit. Steady state limit of the GQFPE is shown to approach the well-known expression derived by CL in the high temperature and Markovian bath limit and also provides additional corrections due to quantum and non-Markovian effects of the bath.

Can quantum transition state theory be defined as an exact t = 0+ limit?

The definition of the classical transition state theory (TST) as a t → 0+ limit of the flux-side time correlation function relies on the assumption that simultaneous measurement of population and flux is a well defined physical process. However, the noncommutativity of the two measurements in quantum mechanics makes the extension of such a concept to the quantum regime impossible. For this reason, quantum TST (QTST) has been generally accepted as any kind of quantum rate theory reproducing the TST in the classical limit, and there has been a broad consensus that no unique QTST retaining all the properties of TST can be defined. Contrary to this widely held view, Hele and Althorpe (HA) [J. Chem. Phys. 138, 084108 (2013)] recently suggested that a true QTST can be defined as the exact t → 0+ limit of a certain kind of quantum flux-side time correlation function and that it is equivalent to the ring polymer molecular dynamics (RPMD) TST. This work seeks to question and clarify certain assumptions underlying these suggestions and their implications. First, the time correlation function used by HA as a starting expression is not related to the kinetic rate constant by virtue of linear response theory, which is the first important step in relating a t = 0+ limit to a physically measurable rate. Second, a theoretical analysis calls into question a key step in HA's proof which appears not to rely on an exact quantum mechanical identity. The correction of this makes the true t = 0+ limit of HA's QTST different from the RPMD-TST rate expression, but rather equal to the well-known path integral quantum transition state theory rate expression for the case of centroid dividing surface. An alternative quantum rate expression is then formulated starting from the linear response theory and by applying a recently developed formalism of real time dynamics of imaginary time path integrals [S. Jang, A. V. Sinitskiy, and G. A. Voth, J. Chem. Phys. 140, 154103 (2014)]. It is shown that the t → 0+ limit of the new rate expression vanishes in the exact quantum limit.

Generalized master equation with non-Markovian multichromophoric Förster resonance energy transfer for modular exciton densities.

A generalized master equation (GME) governing quantum evolution of modular exciton density (MED) is derived for large scale light harvesting systems composed of weakly interacting modules of multiple chromophores. The GME-MED offers a practical framework to incorporate real time coherent quantum dynamics calculations of small length scales into dynamics over large length scales, and also provides a non-Markovian generalization and rigorous derivation of the Pauli master equation employing multichromophoric Förster resonance energy transfer rates. A test of the GME-MED for four sites of the Fenna-Matthews-Olson complex demonstrates how coherent dynamics of excitonic populations over coupled chromophores can be accurately described by transitions between subgroups (modules) of delocalized excitons. Application of the GME-MED to the exciton dynamics between a pair of light harvesting complexes in purple bacteria demonstrates its promise as a computationally efficient tool to investigate large scale exciton dynamics in complex environments.

Kinetic Model for the Desensitization of G Protein-Coupled Receptor.

Desensitization of G-protein-coupled receptors (GPCR) is a general regulatory mechanism adopted by biological organisms against overstimulation of G protein signaling. Although the details of the mechanism are extensively studied, it is not easy to gain an overarching understanding of the process constituted by a multitude of molecular events with vastly differing time scales. To offer a semiquantitative yet predictive understanding of the mechanism, we formulate a kinetic model for the G protein signaling and desensitization by considering essential biochemical steps from ligand binding to receptor internalization. The internalization, followed by receptor depletion from the plasma membrane, attenuates the downstream signal. Together with the kinetic model and its full numerics of the expression derived for the dose-response relation, an approximated form of the expression clarifies the role played by the individual biochemical processes and allows us to identify four distinct regimes for the downregulation that emerge from the balance between phosphorylation, dephosphorylation, and the cellular level of ß-arrestin.

Nonadiabatic Derivative Couplings through Multiple Franck-Condon Modes Dictate the Energy Gap Law for Near and Short-Wave Infrared Dye Molecules.

Near infrared (NIR, 700-1000 nm) and short-wave infrared (SWIR, 1000-2000 nm) dye molecules exhibit significant nonradiative decay rates from the first singlet excited state to the ground state. While these trends can be empirically explained by a simple energy gap law, detailed mechanisms of nearly universal behavior have remained unsettled for many cases. Theoretical and experimental results for two representative NIR/SWIR dye molecules reported here clarify the key mechanism for the observed energy gap law behavior. It is shown that the first derivative nonadiabatic coupling terms serve as major coupling pathways for nonadiabatic decay processes from the first excited singlet state to the ground state for these NIR and SWIR dye molecules and that vibrational modes other than the highest frequency modes also make significant contributions to the rate. This assessment is corroborated by further theoretical comparison with possible alternative mechanisms of intersystem crossing to triplet states and also by comparison with experimental data for deuterated molecules.

Emission lineshapes of the B850 band of light-harvesting 2 (LH2) complex in purple bacteria: a second order time-nonlocal quantum master equation approach.

The emission lineshape of the B850 band in the light harvesting complex 2 of purple bacteria is calculated by extending the approach of 2nd order time-nonlocal quantum master equation [S. Jang and R. J. Silbey, J. Chem. Phys. 118, 9312 (2003)]. The initial condition for the emission process corresponds to the stationary excited state density where exciton states are entangled with the bath modes in equilibrium. This exciton-bath coupling, which is not diagonal in either site excitation or exciton basis, results in a new inhomogeneous term that is absent in the expression for the absorption lineshape. Careful treatment of all the 2nd order terms are made, and explicit expressions are derived for both full 2nd order lineshape expression and the one based on secular approximation that neglects off-diagonal components in the exciton basis. Numerical results are presented for a few representative cases of disorder and temperature. Comparison of emission line shape with the absorption line shape is also made. It is shown that the inhomogeneous term coming from the entanglement of the system and bath degrees of freedom makes significant contributions to the lineshape. It is also found that the perturbative nature of the theory can result in negative portion of lineshape in some situations, which can be removed significantly by inclusion of the inhomogeneous term and completely by using the secular approximation. Comparison of the emission and absorption lineshapes at different temperatures demonstrates the role of thermal population of different exciton states and exciton-phonon couplings.