Search

We present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear-electronic orbital (NEO) theory and real-time time-dependent density functional theory (RT-TDDFT) for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase., Comment: 30 pages 2 figures

Atomic-scale simulations have progressed tremendously over the past decade, largely due to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques, thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities.

Vibrational polaritons form in a planar Fabry-Perot microcavity when a vibrational mode of a layer of molecules is near resonant with an infrared cavity mode. Herein, dispersion relations of vibrational polaritons are studied when the molecular density distribution breaks the macroscopic translational symmetry along the cavity mirror plane. Both perturbative theory and numerical calculations show that, if a homogeneous in-plane molecular distribution is modulated by sinusoidal fluctuations, in addition to a pair of upper and lower polariton branches, a discrete number of side polariton branches may emerge in the polariton dispersion relation. Moreover, for a periodic Gaussian molecular in-plane density distribution, only two, yet significantly broadened polariton branches exist in the spectra. This polariton linewidth broadening is caused by the scattering between cavity modes at neighboring in-plane frequencies due to the symmetry breaking, which is distinguished from known origins of polariton broadening such as the homogeneous broadening of molecules, the cavity loss, or the large energetic disorder of molecules. Associated with the broadened polariton branches, under the periodic Gaussian in-plane inhomogeneity, a significant number of the VSC eigenstates contain a non-zero contribution from the cavity photon mode at zero in-plane frequency, blurring the distinction between the bright and the dark modes. Looking forward, our theoretical investigation should facilitate the experimental exploration of vibrational polaritons with patterned in-plane molecular density distributions.

Developing theoretical frameworks for vibrational strong coupling (VSC) beyond the single-mode approximation is crucial for a comprehensive understanding of experiments with planar Fabry-P\'erot cavities. Herein, a generalized cavity molecular dynamics (CavMD) scheme is developed to simulate VSC of a large ensemble of realistic molecules coupled to an arbitrary 1D or 2D photonic environment. This approach is built upon the Power-Zienau-Woolley Hamiltonian in the normal mode basis and uses a grid representation of the molecular ensembles to reduce the computational cost. When simulating the polariton dispersion relation for a homogeneous distribution of molecules in planar Fabry-P\'erot cavities, our data highlight the importance of preserving the in-plane translational symmetry of the molecular distribution. In this homogeneous limit, CavMD yields the consistent polariton dispersion relation as analytic theory, i.e., incorporating many cavity modes with varying in-plane wave vectors ($k_{\parallel}$) produces the same spectrum as the system with a single cavity mode. Furthermore, CavMD reveals that the validity of the single mode approximation is challenged when nonequilibrium polariton dynamics are considered, as polariton-polariton scattering occurs between modes with nearest neighbor $k_{\parallel}$. The procedure for numerically approaching the macroscopic limit is also demonstrated with CavMD by increasing the system size. Looking forward, our generalized CavMD approach may facilitate understanding vibrational polariton transport and condensation.

A remote energy transfer pathway from electronic to vibrational degrees of freedom is identified inside an infrared optical cavity under vibrational strong coupling conditions. This mechanism relies on the dynamical Casimir effect, whereby real infrared photons are generated due to a sudden electronic transition of anisotropic molecules. Moreover, the formation of vibrational polaritons enables the excited photon energy to be transferred to the vibrational degrees of freedom before any dissipation occurs. Both analytic solutions and numerical simulations reveal that the magnitude of this electronic to vibrational energy transfer depends quadratically on the number of molecules and resonantly on the vibration-cavity detuning. During this "supervibronic" transition process, because the vibrational energy gain per molecule can be meaningful in the macroscopic limit, this process may potentially be observed using conventional vibrational strong coupling devices.

The coupled quantum dynamics of electrons and protons is ubiquitous in many dynamical processes involving light-matter interaction, such as solar energy conversion in chemical systems and photosynthesis. A first-principles description of such nuclear-electronic quantum dynamics requires not only the time-dependent treatment of nonequilibrium electron dynamics but also that of quantum protons. Quantum mechanical correlation between electrons and protons adds further complexity to such coupled dynamics. Here we extend real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) to periodic systems and perform first-principles simulations of coupled quantum dynamics of electrons and protons in complex heterogeneous systems. The process studied is electronically excited state intramolecular proton transfer of o-hydroxybenzaldehyde in water and at a silicon (111) semiconductor-molecule interface. These simulations illustrate how environments such as hydrogen-bonding water molecules and an extended material surface impact the dynamical process on the atomistic level. Depending on how the molecule is chemisorbed on the surface, excited state electron transfer from the molecule to the semiconductor surface can inhibit ultrafast proton transfer within the molecule. This work elucidates how heterogeneous environments influence the balance between the quantum mechanical proton transfer and excited electron dynamics. The periodic RT-NEO-TDDFT approach is applicable to a wide range of other photoinduced heterogeneous processes.

Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born--Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.

Vibrational strong coupling (VSC) provides a novel means to modify chemical reactions and energy transfer pathways. To efficiently model chemical dynamics under VSC in the collective regime, herein a hybrid quantum mechanical/molecular mechanical (QM/MM) cavity molecular dynamics (CavMD) scheme is developed and applied to an experimentally studied chemical system. This approach can achieve linear scaling with respect to the number of molecules for a dilute solution under VSC by assuming that each QM solute molecule is surrounded by an independent MM solvent bath. Application of this approach to a dilute solution of Fe(CO)$_5$ in n-dodecane under VSC demonstrates polariton dephasing to the dark modes and polariton-enhanced molecular nonlinear absorption. These simulations predict that strongly exciting the lower polariton may provide an energy transfer pathway that selectively excites the equatorial CO vibrations rather than the axial CO vibrations. Moreover, these simulations also directly probe the cavity effect on the dynamics of the Fe(CO)$_5$ Berry pseudorotation reaction for comparison to recent two-dimensional infrared spectroscopy experiments. This theoretical approach is applicable to a wide range of other polaritonic systems and provides a tool for exploring the use of VSC for selective infrared photochemistry., Comment: 19 pages, 5 figures in the main text

Molecular polaritons have become an emerging platform for remotely controlling molecular properties through strong light-matter interactions. Herein, a semiclassical approach is developed for describing molecular polaritons by self-consistently propagating the real-time dynamics of classical cavity modes and a quantum molecular subsystem described by the nuclear-electronic orbital (NEO) method, where electrons and specified nuclei are treated quantum mechanically on the same level. This semiclassical real-time NEO approach provides a unified description of electronic and vibrational strong couplings and describes the impact of the cavity on coupled nuclear-electronic dynamics while including nuclear quantum effects. For a single o-hydroxybenzaldehyde molecule under electronic strong coupling, this approach shows that the cavity suppression of excited state intramolecular proton transfer is influenced not only by the polaritonic potential energy surface but also by the timescale of the chemical reaction. This work provides the foundation for exploring collective strong coupling in nuclear-electronic quantum dynamical systems within optical cavities., Comment: manuscript (33 pages) + SI (13 pages) + TOC figure (1 page)

A quantum simulation of vibrational strong coupling (VSC) in the collective regime via thermostatted ring-polymer molecular dynamics (TRPMD) is reported. For a collection of liquid-phase water molecules resonantly coupled to a single lossless cavity mode, the simulation shows that, as compared with a fully classical calculation, the inclusion of nuclear and photonic quantum effects does not lead to a change in the Rabi splitting but does broaden polaritonic linewidths roughly by a factor of two. Moreover, under thermal equilibrium, both quantum and classical simulations predict that the static dielectric constant of liquid water is largely unchanged inside versus outside the cavity. This result disagrees with a recent experiment demonstrating that the static dielectric constant of liquid water can be resonantly enhanced under VSC, suggesting either limitations of our approach or perhaps other experimental factors that have not yet been explored., Comment: manuscript (18 pages) + supporting information (8 pages)

Chemical manifestations of strong light-matter coupling have been recently subject of intense experimental and theoretical studies. Here we review the present status of this field. Section 1 is an introduction to molecular polaritonics and to collective response aspects of light-matter interactions. Section 2 provides an overview of the key experimental observations of these effects while Section 3 describe our current theoretical understanding of the effect of strong light-matter coupling on chemical dynamics. A brief outline of applications to energy conversion processes in given in Section 4. Pending technical issues in the constructions of theoretical approach are briefly described in Section 5. The summary in Section 6 also outlines the paths ahead in this exciting endeavor.

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [J. Chem. Phys., 154, 094124, (2021)], we systematically study the mechanism of polariton relaxation for liquid CO2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations show that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi's golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. We also summarize recent progress on probing nonreactive VSC dynamics with CavMD., Comment: 16 pages, 6 figures

Selectively exciting target molecules to high vibrational states is inefficient in the liquid phase, which restricts the use of IR pumping to catalyze ground-state chemical reactions. Here, we demonstrate that this inefficiency can sometimes be solved by confining the liquid to an optical cavity under vibrational strong coupling conditions. For a liquid solution of 13CO2 solute in a 12CO2 solvent, cavity molecular dynamics simulations show that exciting a polariton (hybrid light-matter state) of the solvent with an intense laser pulse, under suitable resonant conditions, may lead to a very strong (> 3 quanta) and ultrafast (< 1 ps) excitation of the solute, even though the solvent ends up being barely excited. By contrast, outside a cavity the same input pulse fluence can excite the solute by only half a vibrational quantum and the selectivity of excitation is low. Our finding is robust under different cavity volumes, which may lead to observable cavity enhancement on IR photochemical reactions in Fabry-P\'erot cavities., Comment: 10 pages of manuscript + 23 pages of SI

For a small fraction of hot CO2 molecules immersed in a liquid-phase CO2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO2 molecules and a cavity mode accelerates hot-molecule relaxation. The physical mechanism underlying this acceleration is the fact that polaritons, especially the lower polariton, can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates (i) resonantly depend on the cavity mode detuning, (ii) cooperatively depend on molecular concentration or Rabi splitting, and (iii) collectively scale with the number of hot molecules, which is similar to Dicke's superradiance. For larger cavity volumes, due to a balance between this superradiant-like behavior and a smaller light-matter coupling, the total VSC effect on relaxation rates can scale slower than $1/N$, and the average VSC effect per molecule can remain meaningful for up to $N \sim10^4$ molecules forming VSC. Moreover, we find that the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than equally distributing it to all thermal molecules. Finally, we highlight the similarities of parameter dependence between the current finding with VSC catalysis observed in Fabry-Perot microcavities., Comment: 10 pages, 7 figures

Recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry-Perot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics (CavMD) simulations to numerically investigate the linear and nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light-matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption which is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes -- whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states -- under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.

We simulate vibrational strong (VSC) and ultrastrong coupling (V-USC) for liquid water with classical molecular dynamics simulations. When the cavity modes are resonantly coupled to the O-H stretch mode of liquid water, the infrared spectrum shows asymmetric Rabi splitting. The lower polariton (LP) may be suppressed or enhanced relative to the upper polariton (UP) depending on the frequency of the cavity mode. Moreover, although the static properties and the translational diffusion of water are not changed under VSC or V-USC, we do find the modification of the orientational autocorrelation function of H$_2$O molecules especially under V-USC, which could play a role in ground-state chemistry.

Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this "vacuum-field catalysis" remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition configurations. Here, we examine the VSC effect in a cavity environment on chemical rates as calculated by transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual cavities used to examine VSC situations.

We model a collection of $N$ two-level systems (TLSs) coupled to a multimode cavity via Meyer-Miller-Stock-Thoss (MMST) dynamics, sampling both electronic and photonic zero-point energies (ZPEs) and propagating independent trajectories in Wigner phase space. By investigating the ground state stability of a single TLS, we use MMST dynamics to separately study both electronic ZPE effects (which would naively lead to the breakdown of the electronic ground state) as well as photonic ZPE effects (which would naively lead to spontaneous absorption). By contrast, including both effects (i.e., sampling both electronic and photonic ZPEs) leads to the dynamical stability of the electronic ground state. Therefore, MMST dynamics provide a practical way to identify the contributions of self-interaction and vacuum fluctuations. More importantly, we find that MMST dynamics can predict accurate quantum dynamics for both electronic populations and electromagnetic field intensity in the high saturation limit. For a single TLS in a cavity, MMST dynamics correctly predict the initial exponential decay of spontaneous emission, Poincar\'e recurrences, and the positional dependence of a spontaneous emission rate. For an array of $N$ equally spaced TLSs with only one TLS excited initially, MMST dynamics correctly predict the modification of spontaneous emission rate as a function of the spacing between TLSs. Finally, MMST dynamics also correctly model Dicke's superradiance and subradiance (i.e., the dynamics when all TLSs are excited initially) including the correct quantum statistics for the delay time (as found by counting trajectories, for which a full quantum simulation is hard to achieve). Therefore, this work raises the possibility of simulating large-scale collective light-matter interactions with methods beyond mean-field theory., Comment: 18 pages, 13 figures

Semiclassical electrodynamics is an appealing approach for studying light-matter interactions, especially for realistic molecular systems. However, there is no unique semiclassical scheme. On the one hand, intermolecular interactions can be described instantaneously by static two-body interactions connecting different molecules plus a classical transverse E-field; we will call this Hamiltonian #I. On the other hand, intermolecular interactions can also be described as effects that are mediated exclusively through a classical one-body E-field without any quantum effects at all (assuming we ignore electronic exchange); we will call this Hamiltonian #II. Moreover, one can also mix these two Hamiltonians into a third, hybrid Hamiltonian, which preserves quantum electron-electron correlations for lower excitations but describes higher excitations in a mean-field way. To investigate which semiclassical scheme is most reliable for practical use, here we study the real-time dynamics of a pair of identical two-level systems (TLSs) undergoing either resonance energy transfer (RET) or collectively driven dynamics. While all approaches perform reasonably well when there is no strong external excitation, we find that no single approach is perfect for all conditions. Each method has its own distinct problems: Hamiltonian #I performs best for RET but behaves in a complicated manner for driven dynamics. Hamiltonian #II is always stable, but obviously fails for RET at short distances. One key finding is that, under externally driving, a full configuration interaction description of Hamiltonian #I strongly overestimates the long-time electronic energy, highlighting the not obvious fact that, if one plans to merge quantum molecules with classical light, a full, exact treatment of electron-electron correlations can actually lead to worse results than a simple mean-field treatment., Comment: 17 pages, 7 figures

We investigate detailed balance for a quantum system interacting with thermal radiation within mixed quantum-classical theory. For a two-level system coupled to classical radiation fields, three semiclassical methods are benchmarked: (1) Ehrenfest dynamics over-estimate the excited state population at equilibrium due to the failure of capturing vacuum fluctuations. (2) The coupled Maxwell-Bloch equations, which supplement Ehrenfest dynamics by damping at the full golden rule rate, under-estimate the excited state population due to double-counting of the self-interaction effect. (3) Ehrenfest+R dynamics recover detailed balance and the correct thermal equilibrium by enforcing the correct balance between the optical excitation and spontaneous emission of the quantum system. These results highlight the fact that, when properly designed, mixed quantum-classical electrodynamics can simulate thermal equilibrium in the field of nanoplasmonics.

We re-investigate the famous Mollow triplet and show that most of the well-known quantum characteristics of the Mollow triplet--including incoherent emission and a non-standard 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., Comment: 3 figures

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., Comment: 3 figures

In a previous paper, 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 (KHD) theory.

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.

We investigate two key representative semiclassical approaches for propagating resonant energy transfer 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\"orster 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., Comment: 2 figures, 21 pages

The dynamics of an electronic two-level system coupled to an electromagnetic field are simulated explicitly for one and three dimensional systems through semiclassical propagation of the Maxwell-Liouville equations. We consider three flavors of mixed quantum-classical dynamics: the classical path approximation (CPA), Ehrenfest dynamics, and symmetrical quantum-classical (SQC) dynamics. The CPA fails to recover a consistent description of spontaneous emission. A consistent "spontaneous" emission can be obtained from Ehrenfest dynamics--provided that one starts in an electronic superposition state. Spontaneous emission is always obtained using SQC dynamics. Using the SQC and Ehrenfest frameworks, we further calculate the dynamics following an incoming pulse, but here we find very different responses: SQC and Ehrenfest dynamics deviate sometimes strongly in the calculated rate of decay of the transient excited state. Nevertheless, our work confirms the earlier observations by W. Miller [J. Chem. Phys. 69, 2188-2195, 1978] that Ehrenfest dynamics can effectively describe some aspects of spontaneous emission and highlights new possibilities for studying light-matter interactions with semiclassical mechanics., Comment: 15 figures