On the nature of two-photon transitions for a collection of molecules in a Fabry-Perot cavity.

We investigate the effect of a cavity on nonlinear two-photon transitions of a molecular system and we analyze how such an effect depends on the cavity quality factor, the field enhancement, and the possibility of dephasing. We find that the molecular response to strong light fields in a cavity with a variable quality factor can be understood as arising from a balance between (i) the ability of the cavity to enhance the field of an external probe and promote multiphoton transitions more easily and (ii) the fact that the strict selection rules on multiphoton transitions in a cavity support only one resonant frequency within the excitation range. Although our simulations use a classical level description of the radiation field (i.e., we solve Maxwell-Bloch or Maxwell-Liouville equations within the Ehrenfest approximation for the field-molecule interaction), based on experience with this level of approximation in the past studies of plasmonic and polaritonic systems, we believe that our results are valid over a wide range of external probing.

Methods to Calculate Electronic Excited-State Dynamics for Molecules on Large Metal Clusters with Many States: Ensuring Fast Overlap Calculations and a Robust Choice of Phase.

We present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

Modeling nonadiabatic dynamics with degenerate electronic states, intersystem crossing, and spin separation: A key goal for chemical physics.

We examine the many open questions that arise for nonadiabatic dynamics in the presence of degenerate electronic states, e.g., for singlet-to-triplet intersystem crossing where a minimal Hamiltonian must include four states (two of which are always degenerate). In such circumstances, the standard surface hopping approach is not sufficient as the algorithm does not include Berry force. Yet, we hypothesize that such a Berry force may be crucial as far as creating chiral induced spin separation, which is now a burgeoning field of study. Thus, this Perspective highlights the fact that if one can generate a robust and accurate semiclassical approach for the case of degenerate states, one will take a big step forward toward merging chemical physics with spintronics.

On the proper derivation of the Floquet-based quantum classical Liouville equation and surface hopping describing a molecule or material subject to an external field.

We investigate different approaches to derive the proper Floquet-based quantum-classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis and then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form that is similar to the usual QCLE with two modifications: (1) The electronic degrees of freedom are expanded to infinite dimension and (2) the nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a fictitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., [J. Chem. Phys. 140, 084104 (2014)] who found that transforming to the adiabatic representation must always be the last operation applied, although now we have extended this result to a time-dependent Hamiltonian. This paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.

Nonadiabatic Dynamics in a Laser Field: Using Floquet Fewest Switches Surface Hopping To Calculate Electronic Populations for Slow Nuclear Velocities.

We investigate two well-known approaches for extending the fewest switches surface hopping (FSSH) algorithm to periodic time-dependent couplings. The first formalism acts as if the instantaneous adiabatic electronic states were standard adiabatic states, which just happen to evolve in time. The second formalism replaces the role of the usual adiabatic states by the time-independent adiabatic Floquet states. For a set of modified Tully model problems, the Floquet FSSH (F-FSSH) formalism gives a better estimate for both transmission and reflection probabilities than the instantaneous adiabatic FSSH (IA-FSSH) formalism, especially for slow nuclear velocities. More importantly, only F-FSSH predicts the correct final scattering momentum. Finally, in order to use Floquet theory accurately, we find that it is crucial to account for the interference between wavepackets on different Floquet states. Our results should be of interest to all those interested in laser-induced molecular dynamics.

Predictive Semiclassical Model for Coherent and Incoherent Emission in the Strong Field Regime: The Mollow Triplet Revisited.

We reinvestigate the famous Mollow triplet and show that most of the well-known quantum characteristics of the Mollow triplet-including incoherent emission and a nonstandard dependence of the sidebands on detuning-can be recovered quantitatively using semiclassical dynamics with a classical light field. In fact, by not relying on the rotating wave approximation, a semiclassical model predicts some quantum effects beyond the quantum optical Bloch equation, including higher-order scattering and asymmetric sideband features. This Letter highlights the fact that, with strong intensities, many putatively quantum features of light-matter interactions arise from a simple balance of mean-field electrodynamics and elementary spontaneous emission, which requires minimal computational cost. Our results suggest that the application of semiclassical electrodynamics to problems with strong light-matter coupling in the fields of nanophotonics and superradiance are likely to yield a plethora of new information.

Comparison of Different Classical, Semiclassical, and Quantum Treatments of Light-Matter Interactions: Understanding Energy Conservation.

The optical response of an electronic two-level system (TLS) coupled to an incident continuous wave (cw) electromagnetic (EM) field is simulated explicitly in one dimension by the following five approaches: (i) the coupled Maxwell-Bloch equations, (ii) the optical Bloch equation (OBE), (iii) Ehrenfest dynamics, (iv) the Ehrenfest+R approach, and (v) classical dielectric theory (CDT). Our findings are as follows: (i) standard Ehrenfest dynamics predict the correct optical signals only in the linear response regime where vacuum fluctuations are not important; (ii) both the coupled Maxwell-Bloch equations and CDT predict incorrect features for the optical signals in the linear response regime due to a double-counting of self-interaction; (iii) by exactly balancing the effects of self-interaction versus the effects of quantum fluctuations (and insisting on energy conservation), the Ehrenfest+R approach generates the correct optical signals in the linear regime and slightly beyond, yielding, e.g., the correct ratio between the coherent and incoherent scattering EM fields. As such, Ehrenfest+R dynamics agree with dynamics from the quantum OBE, but whereas the latter is easily applicable only for a single TLS in vacuum, the former should be applicable to large systems in environments with arbitrary dielectrics. Thus, this benchmark study suggests that the Ehrenfest+R approach may be very advantageous for simulating light-matter interactions semiclassically.

Ehrenfest+R dynamics. I. A mixed quantum-classical electrodynamics simulation of spontaneous emission.

The dynamics of an electronic system interacting with an electromagnetic field is investigated within mixed quantum-classical theory. Beyond the classical path approximation (where we ignore all feedback from the electronic system on the photon field), we consider all electron-photon interactions explicitly according to Ehrenfest (i.e., mean-field) dynamics and a set of coupled Maxwell-Liouville equations. Because Ehrenfest dynamics cannot capture certain quantum features of the photon field correctly, we propose a new Ehrenfest+R method that can recover (by construction) spontaneous emission while also distinguishing between electromagnetic fluctuations and coherent emission.

Ehrenfest+R dynamics. II. A semiclassical QED framework for Raman scattering.

In Paper I [Chen et al., J. Chem. Phys. 150, 044102 (2019)], we introduced Ehrenfest+R dynamics for a two-level system and showed how spontaneous emission can be heuristically included such that, after averaging over an ensemble of Ehrenfest+R trajectories, one can recover both coherent and incoherent electromagnetic fields. In the present paper, we now show that Ehrenfest+R dynamics can also correctly describe Raman scattering, whose features are completely absent from standard Ehrenfest dynamics. Ehrenfest+R dynamics appear to be quantitatively accurate both for resonant and off-resonant Raman signals, as compared with Kramers-Heisenberg-Dirac theory.

A Necessary Trade-off for Semiclassical Electrodynamics: Accurate Short-Range Coulomb Interactions versus the Enforcement of Causality?

We investigate two key representative semiclassical approaches for propagating resonant energy transfer (RET) between a pair of electronic two-level systems (donor and acceptor) with coupled Maxwell-Liouville equations. On the one hand, when the electromagnetic (EM) field is treated classically and Coulomb interactions are treated quantum-mechanically, we find that a quantum-classical mismatch leads to a violation of causality, i.e., the acceptor can be excited before the retarded EM field arrives. On the other hand, if we invoke a classical intermolecular Coulomb operator, we find that the energy transfer in the near field loses quantitative accuracy compared with Förster theory, even though causality is strictly obeyed. Thus, our work raises a fundamental paradox when choosing a semiclassical electrodynamics algorithm. Namely, which is more important: Accurate short-range interactions or long-range causality? Apparently, one cannot have one's cake and eat it too.

Inchworm Monte Carlo for exact non-adiabatic dynamics. II. Benchmarks and comparison with established methods.

In this second paper of a two part series, we present extensive benchmark results for two different inchworm Monte Carlo expansions for the spin-boson model. Our results are compared to previously developed numerically exact approaches for this problem. A detailed discussion of convergence and error propagation is presented. Our results and analysis allow for an understanding of the benefits and drawbacks of inchworm Monte Carlo compared to other approaches for exact real-time non-adiabatic quantum dynamics.

Inchworm Monte Carlo for exact non-adiabatic dynamics. I. Theory and algorithms.

In this paper, we provide a detailed description of the inchworm Monte Carlo formalism for the exact study of real-time non-adiabatic dynamics. This method optimally recycles Monte Carlo information from earlier times to greatly suppress the dynamical sign problem. Using the example of the spin-boson model, we formulate the inchworm expansion in two distinct ways: The first with respect to an expansion in the system-bath coupling and the second as an expansion in the diabatic coupling. The latter approach motivates the development of a cumulant version of the inchworm Monte Carlo method, which has the benefit of improved scaling. This paper deals completely with methodology, while Paper II provides a comprehensive comparison of the performance of the inchworm Monte Carlo algorithms to other exact methodologies as well as a discussion of the relative advantages and disadvantages of each.

On the accuracy of the Padé-resummed master equation approach to dissipative quantum dynamics.

Well-defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expected to be applicable. Extensive comparison of Padé-resummed master equations with numerically exact results in the context of the spin-boson model demonstrates that the proposed criteria correctly demarcate the regions of parameter space where the Padé approximation is reliable. The applicability analysis we present is not confined to the specifics of the Hamiltonian under consideration and should provide guidelines for other classes of resummation techniques.

On the accuracy of surface hopping dynamics in condensed phase non-adiabatic problems.

We perform extensive benchmark comparisons of surface hopping dynamics with numerically exact calculations for the spin-boson model over a wide range of energetic and coupling parameters as well as temperature. We find that deviations from golden-rule scaling in the Marcus regime are generally small and depend sensitively on the energetic bias between electronic states. Fewest switches surface hopping (FSSH) is found to be surprisingly accurate over a large swath of parameter space. The inclusion of decoherence corrections via the augmented FSSH algorithm improves the accuracy of dynamical behavior compared to exact simulations, but the effects are generally not dramatic, at least for the case of an environment modeled with the commonly used Debye spectral density.

An unusual function of RON receptor tyrosine kinase as a transcriptional regulator in cooperation with EGFR in human cancer cells.

Homodimerization of RON (MST1R), a receptor tyrosine kinase, usually occurs in cells stimulated by a ligand and leads to the downstream activation of signaling pathways. Here we report that bladder cancer cells, in response to physiological stress, use an alternative mechanism for signaling activation. Time-course studies indicated that RON migrated directly from the membrane to the nucleus of bladder cancer cells in response to serum starvation. Biochemical and genetic studies implied that this nuclear internalization was complexed with epidermal growth factor receptor (EGFR) and required the docking of importins. In vivo analysis confirmed that nuclear RON was present in 38.4% (28/73) of primary bladder tumors. Chromatin immunoprecipitation (ChIP) on microarray analysis further revealed that this internalized complex bound to at least 134 target genes known to participate in three stress-responsive networks: p53, stress-activated protein kinase/c-jun N-terminal kinase and phosphatidylinositol 3-kinase/Akt. These findings suggest that RON, in a complex with EGFR, acts as a transcriptional regulator in response to acute disturbances (e.g. serum starvation) imposed on cancer cells. In an attempt to re-establish homeostasis, these cells bypass regular mechanisms required by ligand stimulation and trigger the RON-directed transcriptional response, which confers a survival advantage.