BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Joseph E. Subotnik

Joseph E. Subotnik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phcond-mat.mes-hallphysics.opticsquant-phphysics.atom-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

A Phase-Space Semiclassical Approach for Modeling Nonadiabatic Nuclear Dynamics with Electronic Spin

Chemical relaxation phenomena, including photochemistry and electron transfer processes, form a vigorous area of research in which nonadiabatic dynamics plays a fundamental role. Here, we show that for nonadiabatic dynamics with two electronic states and a complex-valued Hamiltonian that does not obey time-reversal symmetry, the optimal semiclassical approach is to run surface hopping dynamics on a set of phase-space adiabatic surfaces. In order to generate such phase-adiabats, one must isolate a proper set of diabats and apply a phase gauge transformation, before eventually diagonalizing the total Hamiltonian (which is now parameterized by both R and P). The resulting algorithm is valid in both the adiabatic and nonadiabatic limits, incorporates all Berry curvature effects, and allows for the study of semiclassical nonadiabatic dynamics in the presence of spin-orbit coupling and/or external magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Electron Transfer and Spin-Orbit Coupling: How Strong are Berry Force Effects In and Out of Equilibrium in the Presence of Nuclear Friction?

We investigate a spin-boson inspired model of electron transfer, where the diabatic coupling is given by a position-dependent phase, exp(iWx). We consider both equilibrium and nonequilibrium initial conditions. We show that, for this model, all equilibrium results are completely invariant to the sign of W(to infinite order). However, the nonequilibrium results do depend on the sign of W, suggesting that photo-induced electron transfer dynamics are meaningfully affected by Berry forces even in the presence of nuclear friction; furthermore, whenever there is spin-orbit coupling, electronic spin polarization can emerge (at least for some time).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Incorporating Berry Force Effects into The Fewest Switches Surface Hopping Algorithm: Intersystem Crossing and The Case of Electronic Degeneracy

We present a preliminary surface-hopping approach for modeling intersystem crossing (ISC) dynamics between four electronic states: one singlet and one (triply degenerate) triplet. In order to incorporate all Berry force effects, the algorithm requires that, when moving along an adiabatic surface associated with the triplet manifold, \mycomment{one must also keep track of a quasi-diabatic index (akin to a "$m_s$" quantum number) for each trajectory. For a simple model problem, we find that a great deal of new physics can be captured by our algorithm, setting the stage for larger, more realistic (or perhaps even {\em ab initio}) simulations in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modeling Spin-Dependent Nonadiabatic Dynamics with Electronic Degeneracy: A Phase-Space Surface-Hopping Method

Nuclear Berry curvature effects emerge from electronic spin degeneracy and canlead to non-trivial spin-dependent (nonadiabatic) nuclear dynamics. However, such effects are completely neglected in all current mixed quantum-classical methods such as fewest switches surface-hopping. In this work, we present a phase-space surface-hopping (PSSH) approach to simulate singlet-triplet intersystem crossing dynamics. We show that with a simple pseudo-diabatic ansatz, a PSSH algorithm can capture the relevant Berry curvature effects and make predictions in agreement with exact quantum dynamics for a simple singlet-triplet model Hamiltonian. Thus, this approach represents an important step towards simulating photochemical and spin processes concomitantly, as relevant to intersystem crossing and spin-lattice relaxation dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the Meaning of Berry Force For Unrestricted Systems Treated With Mean-Field Electronic Structure

We show that the Berry force as computed by an approximate, mean-field electronic structure can be meaningful if properly interpreted. In particular, for a model Hamiltonian representing a molecular system with an even number of electrons interacting via a two-body (Hubbard) interaction and a spin-orbit coupling, we show that a meaningful nonzero Berry force emerges whenever there is spin unrestriction--even though the Hamiltonian is real-valued and formally the on-diagonal single-surface Berry force must be zero. Moreover, if properly applied, this mean-field Berry force yields roughly the correct asymptotic motion for scattering through an avoided crossing. That being said, within the context of a ground-state calculation, several nuances do arise as far interpreting the Berry force correctly, and as a practical matter, the Berry force diverges near the Coulson-Fisher point (which can lead to numerical instabilities). We do not address magnetic fields here.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi's golden rule rate

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Simulations of Vibrational Strong Coupling via Path Integrals

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Effect of Duschinskii Rotations on Spin-Dependent Electron Transfer Dynamics

We investigate spin-dependent electron transfer in the presence of a Duschinskii rotation. In particular, we propagate dynamics for a two-level model system for which spin-orbit coupling introduces an interstate coupling of the form $e^{iWx}$, which is both position(x)-dependent and complex-valued. We demonstrate that two-level systems coupled to Brownian oscillators with Duschinskii rotations (and thus entangled normal modes) can produce marked increases in transient spin polarization relative to two-level systems coupled to simple shifted harmonic oscillators. These conclusions should have significant relevance for modeling the effect of nuclear motion on chiral induced spin selectivity.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Cavity molecular dynamics simulations of vibrational polariton enhanced molecular nonlinear absorption

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Protocol for Spectroscopists to Isolate The Effect of Berry Geometric Magnetic Forces on Molecular Dynamics

We propose a novel means to isolate and quantify the effects of Berry force on molecular dynamics using two reasonably strong continuous wave (CW) laser fields with frequencies $ω$ and $2ω$. For molecules or materials with three frequency-matching bright transitions ($|{0}\rangle\rightarrow|{1}\rangle$, $|{1}\rangle\rightarrow|{2}\rangle$, $|{0}\rangle\rightarrow|{2}\rangle$) at frequencies ($ω$, $ω$, $2ω$) respectively, the effects of Berry curvature can be isolated by varying the phase between the two laser fields ($Δϕ$) and monitoring the dynamics. Moreover, we find that the resulting chemical dynamics can depend critically on the sign of $Δϕ$; in other words, the effects of Berry curvature can be enormous. Thus, this manuscript represents an unusual step forward towards using light-matter interactions to affect chemical dynamics, suggesting that topological concepts usually invoked in adiabatic quantum optics and condensed matter can be directly applied to non-adiabatic chemical excited state dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analytic Gradients and Derivative Couplings for Configuration Interaction with All Single Excitations and One Double Excitation -- En Route to Nonadiabatic Dynamics

We present analytic gradients and derivative couplings for the simplest possible multireference configuration interaction method, CIS-1D, an electronic structure ansatz that includes all single excitations and one lone double excitation on top of a Hartree-Fock reference state. We show that the resulting equations are numerically stable and require the evaluation of a similar number of integrals as compared to standard CIS theory; one can easily differentiate the required frontier orbitals ($\mathtt{h}$ and $\mathtt{l}$) with minimal cost. The resulting algorithm has been implemented within the Q-Chem electronic structure package and should be immediately useful for understanding photochemistry with $S_{0}$-$S_{1}$ crossings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cavity molecular dynamics simulations of liquid water under vibrational ultrastrong coupling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Near Total Electronic Spin Separation as Caused by Nuclear Dynamics: Perturbing a Real-Valued Conical Intersection with Complex-Valued Spin-Orbit Coupling

We investigate the nuclear dynamics near a real-valued conical intersection that is perturbed by a complex-valued spin-orbit coupling. For a model Hamiltonian with two outgoing channels, we find that even a small spin-orbit coupling can dramatically affect the pathway selection on account of Berry force, leading to extremely large spin selectivity (sometime as large as 100%). Thus, this Letter opens the door for organic chemists to start designing spintronic devices that use nuclear motion and conical intersections (combined with standard spin-orbit coupling) in order to achieve spin selection. Vice versa, for physical chemists, this Letter also emphasizes that future semiclassical simulations of intersystem crossing (which have heretofore ignored Berry force) should be corrected to account for the spin polarization that inevitably arises when dynamics pass near conical intersections.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Origin of Ground-State Vacuum-Field Catalysis: Equilibrium Consideration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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 two 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, then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form which is similar to the usual QCLE with two modifications: 1. The electronic degrees of freedom are expanded to infinite dimension. 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 factitious 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., who found that transforming to the adiabatic representation must always be the last operation applied, though 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Configuration interaction approaches for solving quantum impurity models

We develop several configuration interaction approaches for characterizing the electronic structure of an adsorbate on a metal surface (at least in model form). When one can separate adsorbate from substrate, these methods can achieve a reasonable description of adsorbate on-site electron-electron correlation in the presence of a continuum of states. While the present paper is restricted to the Anderson impurity model, there is hope that these methods can be extended to ab initio Hamiltonians, and provide insight into the structure and dynamics of molecule-metal surface interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quasiclassical Modeling of Cavity Quantum Electrodynamics

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é 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.

0 citations0 reviewsNo code linkedOpen access