Search

In this study, we investigate the crucial role of measurement time ($t_m$), information gain and energy consumption in information engines (IEs) utilizing a von-Neumann measurement model. These important measurement parameters allow us to analyze the efficiency and power output of these devices. As the measurement time increases, the information gain and subsequently the extracted work also increase. However, there is a corresponding increase in the energetic cost. The efficiency of converting information into free energy diminishes as $t_m$ approaches both 0 and infinity, peaking at intermediate values of $t_m$. The power output (work extracted per times) also reaches a maximum at specific operational time regimes. By considering the product of efficiency and power as a performance metric, we can identify the optimal operating conditions for the IE.

Molecular chirality has traditionally been viewed as a binary property where a molecule is classified as either chiral or achiral, yet in the recent decades mathematical methods for quantifying chirality have been explored. Here we use toy molecular systems to systematically compare the performance of two state of the art chirality measures (1) the Continuous Chirality Measure (CCM) and (2) the Chirality Characteristic ($\chi$). We find that both methods exhibit qualitatively similar behavior when applied to simple molecular systems such as a four-site molecule or the polymer double-helix, but we show that the CCM may be more suitable for evaluating the chirality of arbitrary molecules or abstract structures such as normal vibrational modes. We discuss a range of considerations for applying these methods to molecular systems in general, and we provide links to user-friendly codes for both methods. We aim this paper to serve as a concise resource for scientists attempting to familiarize themselves with these chirality measures or attempting to implement chirality measures in their own work.

The lack of quantitative methods for studying chirality has stunted our theoretical understanding of chirality-induced physical phenomena. The chiral phonon has been proposed to be involved in various physical phenomena, but it is not yet well defined mathematically. Here we examine two approaches for assigning and quantifying the chirality of molecular normal modes in double-helical molecular wires with various levels of twist. First, associating with each normal mode a structure obtained by imposing the corresponding motion on a common atomic origin, we apply the Continuous Chirality Measure (CCM) to quantitatively assess the relationship between the chirality-weighted normal mode spectrum and the chirality of the underlying molecular structure. We find that increasing the amount of twist in the double helix shifts the mean normal mode CCM to drastically higher values, implying that the chirality of molecular normal modes is strongly correlated with that of the underlying molecular structure. Second, we assign to each normal mode a pseudoscaler defined as the product of atomic linear and angular momentum summed over all atoms, and we analyze the handedness of the normal mode spectrum with respect to this quantity. We find that twisting the double-chain structure introduces asymmetry between right and left-handed normal modes so that in twisted structures different frequency bands are characterized by distinct handedness. This may give rise to global phenomena such as thermal chirality.

We apply the Kubo-Anderson stochastic theory of molecular spectral lineshape to the case of polaritons formed in the collective strong coupling regime. We investigate both the fast and slow limits of the random frequency modulation of the emitter as well as the intermediate regime and show how the interplay between the characteristic timescales of the cavity and the molecular disorder is expressed in the observed polaritons lineshapes. The analytical solution obtained for the slow limit is valid for any ratio between the inhomogeneous broadening of the molecules and the Rabi splitting, especially relevant for molecular polaritons where these two quantities can be of the same order of magnitude.

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (non-zero) in-plane wavevectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. We find that the transport of the molecular excitons as induced by the external light field is synchronized with the evolution of the enhanced and localized electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps to better understand how experimental observations of polaritons should be interpreted., Comment: 22 Pages, 6 figures

Despite the desirability of polymers for use in many products due to their flexibility, light weight, and durability, their status as thermal insulators has precluded their use in applications where thermal conductors are required. However, recent results suggest that the thermal conductance of polymers can be enhanced and that their heat transport behaviors may be highly sensitive to nanoscale control. Here we use non-equilibrium molecular dynamics (MD) simulations to study the effect of mechanical twist on the steady-state thermal conductance across multi-stranded polyethylene wires. We find that a highly twisted double-helical polyethylene wire can display a thermal conductance up to three times that of its untwisted form, an effect which can be attributed to a structural transition in the strands of the double helix. We also find that in thicker wires composed of many parallel strands, adding just one twist can increase its thermal conductance by over 30%. However, we find that unlike stretching a polymer wire, which causes a monotonic increase in thermal conductance, the effect of twist is highly non-monotonic, and certain amounts of twist can actually decrease the thermal conductance. Finally, we apply the Continuous Chirality Measure (CCM) in an attempt to explore the correlation between heat conductance and chirality. The CCM is found to correlate with twist as expected, but we attribute the observed heat transport behaviors to structural factors other than chirality.

We investigate the effect of a cavity on nonlinear two-photon transitions of a molecular system and 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 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 past studies of plasmonic and polaritonic systems, we believe that our results are valid over a wide range of external probing., Comment: 19 pages, 6 figures

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking with quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical mechanics reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose-Einstein statistics into an effective temperature in the form of modified Langevin equation. We show that the results from such a quasi-classical effective temperature (QCET) MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions, to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

We consider discrete and continuous representations of a thermodynamic process in which a random walker (e.g. a molecular motor on a molecular track) uses a periodically pumped energy (work) to pass $N$ sites and move energetically downhill while dissipating heat. Interestingly, we find that, starting from a discrete model, the limit in which the motion becomes continuous in space and time ($N\to \infty$) is not unique and depends on what physical observables are assumed to be unchanged in the process. In particular, one may (as usually done) choose to keep the speed and diffusion coefficient fixed during this limiting process, in which case the entropy production is affected. In addition, we study also processes in which the entropy production is kept constant as $N\to \infty$ at the cost of modified speed or diffusion coefficient. Furthermore, we also combine this dynamics with work against an opposing force, which makes it possible to study the effect of discretization of the process on the thermodynamic efficiency of transferring power input to power output. Interestingly, we find that the efficiency is increased in the limit of $N\to\infty$. Finally, we investigate the same process when transitions between sites can only happen at finite time intervals and study the impact of this time discretization on the thermodynamic variables as the continuous limit is approached.

A first-principles approach to describe electron dynamics in open quantum systems driven far from equilibrium via external time-dependent stimuli is introduced. Within this approach, the driven Liouville von Neumann methodology is used to impose open boundary conditions on finite model systems, whose dynamics is described using time-dependent density functional theory. As proof of concept, the developed methodology is applied to simple spin-compensated model systems including a hydrogen chain and a graphitic molecular junction. Good agreement between steady-state total currents obtained via direct propagation and from the self-consistent solution of the corresponding Sylvester equation indicates the validity of the implementation. The capability of the new computational approach to analyze, from first principles, non-equilibrium dynamics of open quantum systems in terms of temporally and spatially resolved current densities is demonstrated. Future extensions of the approach towards the description of dynamical magnetization and decoherence effects are briefly discussed.

The realization of single-molecule thermal conductance measurements has driven the need for theoretical tools to describe conduction processes that occur over atomistic length scales. In macroscale systems, the principle that is typically used to understand thermal conductivity is Fourier's law. At molecular length scales, however, deviations from Fourier's law are common in part because microscale thermal transport properties typically depend on the complex interplay between multiple heat conduction mechanisms. Here, the thermal transport properties that arise from electron transfer across a thermal gradient in a molecular conduction junction are examined theoretically. We illustrate how transport in a model junction is affected by varying the electronic structure and length of the molecular bridge in the junction as well as the strength of the coupling between the bridge and its surrounding environment. Three findings are of note: First, the transport properties can vary significantly depending on the characteristics of the molecular bridge and its environment; second, the system's thermal conductance commonly deviates from Fourier's law; and third, in properly engineered systems, the magnitude of electron hopping thermal conductance is similar to what has been measured in single-molecule devices.

The rigorous quantum mechanical description of the collective interaction of many molecules with the radiation field is usually considered numerically intractable, and approximation schemes must be employed. Standard spectroscopy usually contains some levels of perturbation theory, but under strong coupling conditions, other approximations are used. A common approximation is the 1-exciton model in which processes involving weak excitations are described using a basis comprising of the ground state and singly excited states of the molecule cavity-mode system. In another frequently used approximation in numerical investigations, the electromagnetic field is described classically, and the quantum molecular subsystem is treated in the mean-field Hartree approximation with its wavefunction assumed to be a product of single molecules' wavefunctions. The former disregards states that take long time to populate, and is therefore essentially a short time approximation. The latter is not limited in this way, but by its nature, disregards some intermolecular and molecule-field correlations. In this work, we directly compare results obtained from these approximations when applied to several prototype problems involving the optical response of molecule-in-optical cavities systems. In particular, we show that our recent model investigation (J. Chem. Phys. 157, 114108 (2022)) of the interplay between the electronic strong coupling and molecular nuclear dynamics using the truncated 1-exciton approximation agrees very well with the semiclassical mean-field calculation.

Like a free particle, the initial growth of a broad (relative to lattice spacing) wavepacket placed on an ordered lattice is slow (zero initial slope) and becomes linear in $t$ at long time. On a disordered lattice, the growth is inhibited at long time (Anderson localization). We consider site disorder with nearest-neighbor hopping on 1- and 2-dimensional systems, and show via numerical simulations supported by the analytical study that the short time growth of the particle distribution is faster on the disordered lattice than on the ordered one. Such faster spread takes place on time and length scale that may be relevant to the exciton motion in disordered systems.

The interplay between collective optical response and molecular static and dynamic disorder is studied using simple effective Hamiltonians for an ensemble of two-level emitters inside and outside a single-mode cavity. We model environmental disorder by randomly modulating the molecular transition frequencies and the coupling between the emitters and the electromagnetic field. We also consider effects of intermolecular interactions and orientational disorder. We investigate how these effects lead to new features in the steady state absorption (outside the cavity), transmission spectra (inside the cavity) and the yield of local molecular processes such as a unimolecular reaction. Outside the cavity, the collective behavior is manifested in the linewidth of the steady state absorption, the emission spectrum, and the local chemical yield. Inside the cavity, however, the collective behavior primarily determines the Rabi splitting. Effects of intermolecular interactions under orientational disorder are also studied. For the most part, for all types of disorder, if we increase disorder, we find a reduction in the collective nature of the molecular response (smaller effective N), and therefore, the Rabi splitting contraction occurs with orientational disorder. Moreover, we find that static disorder is more destructive to collective behavior than dynamic disorder., Comment: 32 pages, 15 figures

We determine the zero-frequency charge current noise in a metal-molecule-metal junction embedded in a thermal environment, e.g., a solvent, dominated by sequential charge transmission described by a classical master equation, and study its dependence of specific model parameters, i.e., the environmental reorganization energy and relaxation behavior. Interestingly, the classical current noise term has the same structure as its quantum analog which reflects a charge correlation due to the bridging molecule. We further determine the thermodynamic uncertainty relation (TUR) which defines a bound on the relationship between the average charge current, its fluctuation and the entropy production in an electrochemical junction in the Marcus regime. In a second part, we use the same methodology to calculate the current noise and the TUR for a protoype photovoltaic cell in order to predict its upper bound for the efficiency of energy conversion into useful work., Comment: (in press)

We consider an ensemble of diatomic molecules resonantly coupled to an optical cavity under strong coupling conditions at normal incidence. Photodissociation dynamics is examined via direct numerical integration of the coupled Maxwell-Schrodinger equations with molecular ro-vibrational degrees of freedom explicitly taken into account. It is shown that the dissociation is significantly affected (slowed down) when the system is driven at its polaritonic frequencies. The observed effect is demonstrated to be of transient nature and has no classical analog. An intuitive explanation of the dissociation slowdown at polaritonic frequencies is proposed.

A model designed to mimic the implications of the collective optical response of molecular ensembles in optical cavities on molecular vibronic dynamics is investigated. Strong molecule-radiation field coupling is often reached when a large number N of molecules respond collectively to the radiation field. In electronic strong coupling, molecular nuclear dynamics following polariton excitation reflects (a) the timescale separation between the fast electronic and photonic dynamics and the slow nuclear motion on one hand, and (b) the interplay between the collective nature of the molecule-field coupling and the local nature of the molecules nuclear response on the other. The first implies that the electronic excitation takes place, in the spirit of the Born approximation, at an approximately fixed nuclear configuration. The second can be rephrased as the intriguing question, can the collective nature of the optical excitation lead to collective nuclear motion following polariton formation, resulting in so-called polaron decoupled dynamics. We address this issue by studying the dynamical properties of a simplified Holstein-Tavis-Cummings type model, in which boson modes representing molecular vibrations are replaced by two-level systems while the boson frequency and the vibronic coupling are represented by the coupling between these levels (that induces Rabi oscillations between them) and electronic state dependence of this coupling. We investigate the short-time behavior of this model following polariton excitation as well as its response to CW driving and its density of states spectrum. We find that, while some aspects of the dynamical behavior appear to adhere to the polaron decoupling picture, the observed dynamics mostly reflect the local nature of the nuclear configuration of the electronic polariton rather than this picture.

The interplay between static and dynamic disorder and collective optical response in molecular ensembles is an important characteristic of nanoplasmonic and nanophotonic molecular systems. Here we investigate the cooperative superradiant response of a molecular ensemble of quantum emitters under the influence of environmental disorder, including inhomogeneous broadening (as induced by static random distribution of the molecular transition frequencies) and motional narrowing (as induced by stochastic modulation of these excitation energies). The effect of inhomogeneous broadening is to destroy the coherence of the collective molecular excitation and suppress superradiant emission. However, fast stochastic modulation of the molecular excitation energy can effectively restore the coherence of the quantum emitters and lead to a recovery of superradiant emission, which is an unexpected manifestation of motional narrowing. For a light scattering process as induced by an off-resonant incident pulse, stochastic modulation leads to inelastic fluorescence emission at the average excitation energy at long times and suggests that dynamic disorder effects can actually lead to collective excitation of the molecular ensemble.

In standard molecular junctions, a molecular structure is placed between and connected to metal leads. Understanding how mechanical tuning in such molecular junctions can change heat conductance has interesting applications in nanoscale energy transport. In this work, we use nonequilibrium molecular dynamics simulations to address the effect of stretching on the phononic contribution to the heat conduction of molecular junctions consisting of single long-chain alkanes and various metal leads such as Ag, Au, Cu, Ni, and Pt. The thermal conductance of such junctions is found to be much smaller than the intrinsic thermal conductance of the polymer and significantly depends on the nature of metal leads as expressed by the metal-molecule coupling and metal vibrational density of states. This behavior is expected and reflects the mismatch of phonon spectra at the metal molecule interfaces. As a function of stretching, we find a behavior similar to what was observed earlier [J. Chem. Phys. 153, 164903 (2020)] for pure polymeric structures. At relatively short electrode distances, where the polyethylene chains are compressed, it is found that the thermal conductances of the molecular junctions remain almost constant as one stretches the polymer chains. At critical electrode distances, the thermal conductances start to increase, reaching the values of the fully-extended molecular junctions. Similar behaviors are observed for junctions in which several long-chain alkanes are sandwiched between various metal leads. These findings indicate that this behavior under stretching is an intrinsic property of the polymer chain and not significantly associated with the interfacial structures., Comment: 23 pages, 7 figures, 3 tables

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)

We consider a molecular junction immersed in a solvent where the electron transfer is dominated by Marcus-type steps. However, the successive nature of the charge transfer through the junction does not imply that the solvent reach thermal equilibrium throughout the transport. In our previous work \cite{kir2020} we have determined the nonequilibrium distribution of the solvent where its dynamics, expressed by a friction, is considered in two limiting regimes of fast and slow solvent relaxation. In dependence of the nonequilibrium solvent dynamics, we investigate now the electrical, thermal and thermoelectric properties of the molecular junction. We show that by suitable tuning the friction, we can reduce the heat dissipation into the solvent and enhance the heat transfer between the electrodes. Interestingly, we find that the Seebeck coefficient grows significantly by adapting the solvent friction in both regimes.

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

The interaction between excited states of a molecule and excited states of metal nanostructure (e.g. plasmons) leads to hybrid states with modified optical properties. When plasmon resonance is swept through molecular transition frequency an avoided crossing may be observed, which is often regarded as a signature of strong coupling between plasmons and molecules. Such strong coupling is expected to be realized when $2|U|/{\hbar\Gamma}>1$, where $U$ and ${\Gamma}$ are the molecule-plasmon coupling and the spectral width of the optical transition respectively. Because both $U$ and ${\Gamma}$ strongly increase with decreasing distance between a molecule and a plasmonic structure it is not obvious that this condition can be satisfied for any molecule-metal surface distance. In this work we investigate the behavior of $U$ and ${\Gamma}$ for several geometries. Surprisingly, we find that if the only contributions to ${\Gamma}$ are lifetime broadenings associated with the radiative and nonradiative relaxation of a single molecular vibronic transition, including effects on molecular radiative and nonradiative lifetimes induced by the metal, the criterion $2|U|/{\hbar\Gamma}>1$ is easily satisfied by many configurations irrespective of the metal-molecule distance. This implies that the Rabi splitting can be observed in such structures if other sources of broadening are suppressed. Additionally, when the molecule-metal surface distance is varied keeping all other molecular and metal parameters constant, this behavior is mitigated due to the spectral shift associated with the same molecule-plasmon interaction, making the observation of Rabi splitting more challenging.

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.

Molecular conduction operating in dielectric solvent environments are often described using kinetic rates based on Marcus theory of electron transfer at a molecule-metal electrode interface. However, the successive nature of charge transfer in such system implies that the solvent does not necessarily reach equilibrium in such process. Here we generalize the theory to account for solvent nonequilibrium and consider a molecular junction consisting of an electronic donor-acceptor system coupled to two metallic electrodes and placed in a polarizable solvent. We determine the nonequilbrium distribution of the solvent by solving diffusion equations in the strong- and weak-friction limits and calculate the charge current and its fluctuating behavior. In extreme limits: the absence of the solvent or fast solvent relaxation, the charge transfer statistics is Poissonian, while it becomes correlated by the dynamic solvent in between these limits. A Kramers-like turnover of the nonequilibrium current as a function of the solvent damping is found. Finally, we propose a way to tune the solvent-induced damping using geometrical control of the solvent dielectric response in nanostructured solvent channels.

We present a consistent theory of energy balance and conversion in a single-molecule junction with strong interactions between electrons on the molecular linker (dot) and phonons in the nuclear environment where the Marcus-type electron hopping processes predominate in the electron transport. It is shown that the environmental reorganization and relaxation that accompany electron hopping energy exchange between the electrodes and the nuclear (molecular and solvent) environment may bring a moderate local cooling of the latter in biased systems. The effect of a periodically driven dot level on the heat transport and power generated in the system is analyzed and energy conservation is demonstrated both within and beyond the quasistatic regime. Finally, a simple model of atomic scale engine based on a Marcus single-molecule junction with a driven electron level is suggested and discussed.

We present a scattering approach for the study of the transport and thermodynamics of quantum systems strongly coupled to their thermal environment(s). This formalism recovers the standard non-equilibrium Green's function expressions for quantum transport and reproduces recently obtained results for the quantum thermodynamic of slowly driven systems. Using this approach, new results have been obtained. First, we derived of a general explicit expression for non-equilibrium steady state density matrix of a system compromised of multiple infinite baths coupled through a general interaction. Then, we obtained a general expression for the dissipated power for the driven non-interacting resonant level to first order in the driving speeds, where both the dot energy level and its couplings are changing, without invoking the wide band approximation. In addition, we also showed that the symmetric splitting of system bath interaction, employed for the case of a system coupled to one bath to determine the effective system Hamiltonian [Phys. Rev. B 93, 115318 (2016)] is valid for the multiple baths case as well. Finally, we demonstrated an equivalence of our method to the Landauer-Buttiker formalism and its extension to slowly driven systems developed by von Oppen and co-workers [Phys. Rev. Lett. 120, 107701 (2018)]. To demonstrate the use of this formalism we analyze the operation a device in which the dot is driven cyclically between two leads under strong coupling conditions. We also generalize the previously obtained expression for entropy production in such driven processes to the many-bath case., Comment: 86 pages, 2 figures

Non-equilibrium thermodynamics of the driven resonant level model is studied using numerical simulations based on the driven Liouville von-Neumann formalism. The approach is first validated against recently obtained analytical results for quasi-static level shifts and the corresponding first order corrections. The numerical approach is then used to study far-from-equilibrium thermodynamic properties of the system under finite level shift rates. The proposed methodology allows the study of unexplored non-equilibrium thermodynamic regimes in open quantum systems.

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

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. 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 Floquet states. Our results should be of interest to all those interested in laser induced molecular dynamics., Comment: 34 pages, 6 figures, typo in Eq. (46) corrected

The Wiedemann-Franz (WF) law is a fundamental result in solid-state physics that relates the thermal and electrical conductivity of a metal. It is derived from the predominant origin of energy conversion in metals: the motion of quasi-free charge-carrying particles. Here, an equivalent WF relationship is developed for molecular systems in which charge carriers are moving not as free particles but instead hop between redox sites. We derive a concise analytical relationship between the electrical and thermal conductivity generated by electron hopping in molecular systems and find that the linear temperature dependence of their ratio as expressed in the standard WF law is replaced by a linear dependence on the nuclear reorganization energy associated with the electron hopping process. The robustness of the molecular WF relation is confirmed by examining the conductance properties of a paradigmatic molecular junction. This result opens a new way to analyze conductivity in molecular systems, with possible applications advancing the design of molecular technologies that derive their function from electrical and/or thermal conductance.

Understanding phononic heat transport processes in molecular junctions is a central issue in the developing field of nanoscale heat conduction and manipulation. Here we present a Stochastic Nonequlibrium Molecular Dynamics simulation framework to investigate heat transport processes in molecular junctions in and beyond the linear response regime. We use extended molecular models which filter Markovian heat reservoirs through an intermediate substrate region, to provide a realistic and controllable effective bath spectral density. The results obtained for alkanedithol molecules connecting gold substrates agree with previous nonequilibrium Green's function calculations in frequency domain, and match recent experimental measurements (e.g. thermal conductance around 20 pW/K for alkanedithiols in single molecular junctions) Classical MD simulations using the full molecular forcefield and quantum Landauer-type calculations based on the harmonic part of the same forcefield are compared, and the similarity of the results indicate that heat transport is dominated by modes in the lower frequency range. Heat conductance simulations on polyynes of different lengths illuminates the effects of molecular conjugation on thermal transport.

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

A theory is developed to describe the coupled transport of energy and charge in networks of electron donor-acceptor sites which are seated in a thermally heterogeneous environment, where the transfer kinetics are dominated by Marcus-type hopping rates. It is found that the coupling of heat and charge transfer in such systems gives rise to exotic transport phenomena which are absent in thermally homogeneous systems and cannot be described by standard thermoelectric relations. Specifically, the directionality and extent of thermal transistor amplification and cyclical electronic currents in a given network can be controlled by tuning the underlying temperature gradient in the system. The application of these findings toward optimal control of multithermal currents is illustrated on a paradigmatic nanostructure.

Controlling the direction and magnitude of both heat and electronic currents using rectifiers has significant implications for the advancement of molecular circuit design. In order to facilitate the implementation of new transport phenomena in such molecular structures, we examine thermal and thermoelectric rectification effects that are induced by an electron transfer process that occurs across a temperature gradient between molecules. Historically, the only known heat conduction mechanism able to generate thermal rectification in purely molecular environments is phononic heat transport. Here, we show that electron transfer between molecular sites with different local temperatures can also generate a thermal rectification effect and that electron hopping through molecular bridges connecting metal leads at different temperatures gives rise to asymmetric Seebeck effects, that is, thermoelectric rectification, in molecular junctions.

The rate of electron transfer between a molecular species and a metal, each at a different local temperature, is examined theoretically through implementation of a bithermal (characterized by two temperatures) Marcus formalism. Expressions for the rate constant and the electronic contribution to a heat transfer mechanism which is induced by the temperature gradient between molecule and metal are constructed. The system of coupled dynamical equations describing the electronic and thermal currents are derived and examined over diverse ranges of reaction geometries and temperature gradients. It is shown that electron transfer across the molecule-metal interface is associated with heat transfer and that the electron exchange between metal and molecule makes a distinct contribution to the interfacial heat conduction even when the net electronic current vanishes.

Statistical properties of Brownian motion that arise by analyzing, separately, trajectories over which the system energy increases (upside) or decreases (downside) with respect to a threshold energy level, are derived. This selective analysis is applied to examine transport properties of a nonequilibrium Brownian process that is coupled to multiple thermal sources characterized by different temperatures. Distributions, moments, and correlation functions of a free particle that occur during upside and downside events are investigated for energy activation and energy relaxation processes, and also for positive and negative energy fluctuations from the average energy. The presented results are sufficiently general and can be applied without modification to standard Brownian motion. This article focuses on the mathematical basis of this selective analysis. In subsequent articles in this series we apply this general formalism to processes in which heat transfer between thermal reservoirs is mediated by activated rate processes that take place in a system bridging them.

The energy partitioning during activation and relaxation events under steady-state conditions for a Brownian particle driven by multiple thermal reservoirs of different local temperatures is investigated. Specifically, we apply the formalism derived in a previous article [G. T. Craven and A. Nitzan, J. Chem. Phys. 148, 044101 (2018)] to examine the thermal transport properties of two sub-ensembles of Brownian processes, distinguished at any given time by the specification that all the trajectories in each group have, at that time, energy either above (upside) or below (downside) a preselected energy threshold. Dynamical properties describing energy accumulation and release during activation/relaxation events and relations for upside/downside energy partitioning between thermal reservoirs are derived. The implications for heat transport induced by upside and downside events are discussed.

The physics of collective optical response of molecular assemblies, pioneered by Dicke in 1954, has long been at the center of theoretical and experimental scrutiny. The influence of the environment on such phenomena is also of great interest due to various important applications in e.g. energy conversion devices. In this manuscript we demonstrate both experimentally and theoretically the spatial modulations of the collective decay rates of molecules placed in proximity to a metal interface. We show in a very simple framework how the cooperative optical response can be analyzed in terms of intermolecular correlations causing interference between the response of different molecules and the polarization induced on a nearby metallic boundary and predict similar collective interference phenomena in excitation energy transfer between molecular aggregates.

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.

The interaction between molecular (atomic) electron(s) and the vacuum field of a reflective cavity generates a significant interest thanks to the rapid developments in nanophotonics. Such interaction which lies within the realm of cavity quantum electrodynamic can substantially affect transport properties of molecular systems. In this work we consider non-adiabatic electron transfer process in the presence of a cavity mode. We present a generalized framework for the interaction between a charged molecular system and a quantized electromagnetic field of a cavity and apply it to the problem of electron transfer between a donor and an acceptor placed in a confined vacuum electromagnetic field. The effective system Hamiltonian corresponds to a unified Rabi and spin-boson model which includes a self-dipole energy term. Two limiting cases are considered: one where the electron is assumed much faster than the cavity mode and another in which the electron tunneling time is significantly larger than the mode period. In both cases a significant rate enhancement can be produced by coupling to the cavity mode in the Marcus inverted region. The results of this work offer new possibilities for controlling electron transfer processes using visible and infrared plasmonics, Comment: 33 pages, 4 figures

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

Methodological aspects of using the driven Liouville-von Neumann (DLvN) approach for simulating dynamical properties of molecular junctions are discussed. As a model system we consider a non-interacting resonant level uniformly coupled to a single Fermionic bath. We demonstrate how a finite system can mimic the depopulation dynamics of the dot into an infinite band bath of continuous and uniform density of states. We further show how the effects of spurious energy resolved currents, appearing due to the approximate nature of the equilibrium state obtained in DLvN calculations, can be avoided. Several ways to approach the wide band limit that is often adopted in analytical treatments, using a finite numerical model system are discussed including brute-force increase of the lead model bandwidth as well as efficient cancellation or direct subtraction of finite-bandwidth effect. These methodological considerations may be relevant also for other numerical schemes that aim to study non-equilibrium thermodynamics via simulations of open quantum systems., Comment: 30 pages, 6 figures