Linear response approximations are central to our understanding and simulations of nonequilibrium statistical mechanics. Despite the success of these approaches in predicting nonequilibrium dynamics, open questions remain. Laird and Thompson [J. Chem. Phys. 126, 211104 (2007)] previously formalized, in the context of solvation dynamics, the connection between the static linear-response approximation and the assumption of Gaussian statistics. The Gaussian statistics perspective is useful in understanding why linear response approximations are still accurate for perturbations much larger than thermal energies. In this paper, we use this approach to address three outstanding issues in the context of the "dipole-flip" model, which is known to exhibit nonlinear response. First, we demonstrate how non-Gaussian statistics can be predicted from purely equilibrium molecular dynamics (MD) simulations (i.e., without resort to a full nonequilibrium MD as is the current practice). Second, we show that the Gaussian statistics approximation may also be used to identify the physical origins of nonlinear response residing in a small number of coordinates. Third, we explore an approach for correcting the Gaussian statistics approximation for nonlinear response effects using the same equilibrium simulation. The results are discussed in the context of several other examples of nonlinear responses throughout the literature.

We present a framework for simulating relaxation dynamics through a conical intersection of an open quantum system that combines methods to approximate the motion of degrees of freedom with disparate time and energy scales. In the vicinity of a conical intersection, a few degrees of freedom render the nuclear dynamics nonadiabatic with respect to the electronic degrees of freedom. We treat these strongly coupled modes by evolving their wavepacket dynamics in the absence of additional coupling exactly. The remaining weakly coupled nuclear degrees of freedom are partitioned into modes that are fast relative to the nonadiabatic coupling and those that are slow. The fast degrees of freedom can be traced out and treated with second-order perturbation theory in the form of the time-convolutionless master equation. The slow degrees of freedom are assumed to be frozen over the ultrafast relaxation, and treated as sources of static disorder. In this way, we adopt the recently developed frozenmode extension to second-order quantum master equations. We benchmark this approach to numerically exact results in models of pyrazine internal conversion and rhodopsin photoisomerization. We use this framework to study the dependence of the quantum yield on the reorganization energy and the characteristic timescale of the bath, in a two-mode model of photoisomerization. We find that the yield is monotonically increasing with reorganization energy for a Markovian bath, but monotonically decreasing with reorganization energy for a non-Markovian bath. This reflects the subtle interplay between dissipation and decoherence in conical intersection dynamics in the condensed phase. arXiv:1905.00029v2 [cond-mat.stat-mech]

We present a method to study rare nonadiabatic dynamics in open quantum systems using transition path sampling and quantum jump trajectories. As with applications of transition path sampling to classical dynamics, the method does not rely on prior knowledge of transition states or reactive pathways, and thus can provide mechanistic insight into ultrafast relaxation processes in addition to their associated rates. In particular, we formulate a quantum path ensemble using the stochastic realizations of an unravelled quantum master equation, which results in trajectories that can be conditioned on starting and ending in particular quantum states. Because the dynamics rigorously obeys detailed balance, rate constants can be evaluated from reversible work calculations in this conditioned ensemble, allowing for branching ratios and yields to be computed in an unbiased manner. We illustrate the utility of this method with three examples: energy transfer in a donor-bridge-acceptor model, and models of photo-induced proton-coupled electron transfer and thermally activated electron transfer. These examples demonstrate the efficacy of path ensemble methods and pave the way for their use in studying of complex reactive quantum dynamics.