Quantum Carnot thermal machines reexamined: Definition of efficiency and the effects of strong coupling.

Whether the strong coupling to thermal baths can improve the performance of quantum thermal machines remains an open issue under active debate. Here we revisit quantum thermal machines operating with the quasistatic Carnot cycle and aim to unveil the role of strong coupling in maximum efficiency. Our analysis builds upon definitions of excess work and heat derived from an exact formulation of the first law of thermodynamics for the working substance, which captures the non-Gibbsian thermal equilibrium state that emerges at strong couplings during quasistatic isothermal processes. These excess definitions differ from conventional ones by an energetic cost for maintaining the non-Gibbsian characteristics. With this distinction, we point out that one can introduce two different yet thermodynamically allowed definitions for efficiency of both the heat engine and refrigerator modes. We dub them excess and hybrid definitions, which differ in the way of defining the gain for the thermal machines at strong couplings by either just analyzing the energetics of the working substance or instead evaluating the performance from an external system upon which the thermal machine acts, respectively. We analytically demonstrate that the excess definition predicts that the Carnot limit remains the upper bound for both operation modes at strong couplings, whereas the hybrid one reveals that strong coupling can suppress the maximum efficiency rendering the Carnot limit unattainable. These seemingly incompatible predictions thus indicate that it is imperative to first gauge the definition for efficiency before elucidating the exact role of strong coupling, thereby shedding light on the ongoing investigation on strong-coupling quantum thermal machines.

Electron transfer at electrode interfaces via a straightforward quasiclassical fermionic mapping approach.

Electron transfer at electrode interfaces to molecules in solution or at the electrode surface plays a vital role in numerous technological processes. However, treating these processes requires a unified and accurate treatment of the fermionic states of the electrode and their coupling to the molecule being oxidized or reduced in the electrochemical processes and, in turn, the way the molecular energy levels are modulated by the bosonic nuclear modes of the molecule and solvent. Here we present a physically transparent quasiclassical scheme to treat these electrochemical electron transfer processes in the presence of molecular vibrations by using an appropriately chosen mapping of the fermionic variables. We demonstrate that this approach, which is exact in the limit of non-interacting fermions in the absence of coupling to vibrations, is able to accurately capture the electron transfer dynamics from the electrode even when the process is coupled to vibrational motions in the regimes of weak coupling. This approach, thus, provides a scalable strategy to explicitly treat electron transfer from electrode interfaces in condensed-phase molecular systems.

Quantum Kinetic Rates within the Nonequilibrium Steady State.

The nonequilibrium steady state (NESS) of a quantum network is central to a host of physical and biological scenarios. Examples include natural processes such as vision and photosynthesis as well as technical devices such as photocells, both activated by incoherent light (e.g., sunlight) and leading to quantum transport. Assessing time scales of the relevant chemical processes in the steady state is thus of utmost interest and is our goal in this paper. Here, a completely general approach to defining components of a quantum network in the NESS and obtaining rates of processes between these components is provided. Quantum effects are explicitly included throughout, both in (a) defining network components via projection operators and (b) determining the role of coherences in rate processes. As examples, the methodology is applied to model cases, two versions of the V-level system, and to the spin-boson model, wherein the roles of the environment and of internal system properties in determining the rates are examined. In addition, the role of Markovian vs non-Markovian contributions is quantified, exposing conditions under which NESS rates can be obtained by perturbing the nonequilibrium steady state.

2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions.

The third-order response lies at the heart of simulating and interpreting nonlinear spectroscopies ranging from two-dimensional infrared (2D-IR) to 2D electronic (2D-ES), and 2D sum frequency generation (2D-SFG). The extra time and frequency dimensions in these spectroscopic techniques provide access to rich information on the electronic and vibrational states present, the coupling between them, and the resulting rates at which they exchange energy that are obscured in linear spectroscopy, particularly for condensed phase systems that usually contain many overlapping features. While the exact quantum expression for the third-order response is well established, it is incompatible with the methods that are practical for calculating the atomistic dynamics of large condensed phase systems. These methods, which include both classical mechanics and quantum dynamics methods that retain quantum statistical properties while obeying the symmetries of classical dynamics, such as LSC-IVR, centroid molecular dynamics, and Ring Polymer Molecular Dynamics (RPMD), naturally provide short-time approximations to the multi-time symmetrized Kubo transformed correlation function. Here, we show how the third-order response can be formulated in terms of equilibrium symmetrized Kubo transformed correlation functions. We demonstrate the utility and accuracy of our approach by showing how it can be used to obtain the third-order response of a series of model systems using both classical dynamics and RPMD. In particular, we show that this approach captures features such as anharmonically induced vertical splittings and peak shifts while providing a physically transparent framework for understanding multidimensional spectroscopies.

Optimal linear cyclic quantum heat engines cannot benefit from strong coupling.

Uncovering whether strong system-bath coupling can be an advantageous operation resource for energy conversion can facilitate the development of efficient quantum heat engines (QHEs). Yet, a consensus on this ongoing debate is still lacking owing to challenges arising from treating strong couplings. Here, we conclude the debate for optimal linear cyclic QHEs operated under a small temperature difference by revealing the detrimental role of strong system-bath coupling in their optimal operations. We analytically demonstrate that both the efficiency at maximum power and maximum efficiency of strong-coupling linear cyclic QHEs are upper bounded by their weak-coupling counterparts with the same degree of time-reversal symmetry breaking. Under strong time-reversal symmetry breaking, we further reveal a quadratic suppression of the optimal efficiencies relative to the Carnot limit when away from the weak-coupling regime, along with a quadratic enhancement of the mean entropy production rate.

Optically Induced Anisotropy in Time-Resolved Scattering: Imaging Molecular-Scale Structure and Dynamics in Disordered Media with Experiment and Theory.

Time-resolved scattering experiments enable imaging of materials at the molecular scale with femtosecond time resolution. However, in disordered media they provide access to just one radial dimension thus limiting the study of orientational structure and dynamics. Here we introduce a rigorous and practical theoretical framework for predicting and interpreting experiments combining optically induced anisotropy and time-resolved scattering. Using impulsive nuclear Raman and ultrafast x-ray scattering experiments of chloroform and simulations, we demonstrate that this framework can accurately predict and elucidate both the spatial and temporal features of these experiments.

Periodically Driven Quantum Thermal Machines from Warming up to Limit Cycle.

Theoretical treatments of periodically driven quantum thermal machines (PD-QTMs) are largely focused on the limit-cycle stage of operation characterized by a periodic state of the system. Yet, this regime is not immediately accessible for experimental verification. Here, we present a general thermodynamic framework that handles the performance of PD-QTMs both before and during the limit-cycle stage of operation. It is achieved by observing that periodicity may break down at the ensemble average level, even in the limit-cycle phase. With this observation, and using conventional thermodynamic expressions for work and heat, we find that a complete description of the first law of thermodynamics for PD-QTMs requires a new contribution, which vanishes only in the limit-cycle phase under rather weak system-bath couplings. Significantly, this contribution is substantial at strong couplings even at limit cycle, thus largely affecting the behavior of the thermodynamic efficiency. We demonstrate our framework by simulating a quantum Otto engine building upon a driven resonant level model. Our results provide new insights towards a complete description of PD-QTMs, from turn-on to the limit-cycle stage and, particularly, shed light on the development of quantum thermodynamics at strong coupling.

Ring-polymer, centroid, and mean-field approximations to multi-time Matsubara dynamics.

Based on a recently developed generalization of Matsubara dynamics to the multi-time realm, we present a formal derivation of multi-time generalizations of ring-polymer molecular dynamics, thermostatted ring-polymer molecular dynamics (TRPMD), centroid molecular dynamics (CMD), and mean-field Matsubara dynamics. Additionally, we analyze the short-time accuracy of each methodology. We find that for multi-time correlation functions of linear operators, (T)RPMD is accurate up to order t3, while CMD is only correct up to t, indicating a degradation in the accuracy of these methodologies with respect to the single-time counterparts. The present work provides a firm justification for the use of path-integral-based approximations for the calculation of multi-time correlation functions.

Energy transfer under natural incoherent light: Effects of asymmetry on efficiency.

The non-equilibrium stationary coherences that form in donor-acceptor systems are investigated to determine their relationship to the efficiency of energy transfer to a neighboring reaction center. It is found that the effects of asymmetry in the dimer are generally detrimental to the transfer of energy. Four types of systems are examined, arising from combinations of localized trapping, delocalized (Forster) trapping, eigenstate dephasing, and site basis dephasing. In the cases of site basis dephasing, the interplay between the energy gap of the excited dimer states and the environment is shown to give rise to a turnover effect in the efficiency under weak dimer coupling conditions. Furthermore, the nature of the coherences and associated flux is interpreted in terms of pathway interference effects. In addition, regardless of the cases considered, the ratio of the real part and the imaginary part of the coherences in the energy-eigenbasis tends to a constant value in the steady state limit.

Two-dimensional Raman spectroscopy of Lennard-Jones liquids via ring-polymer molecular dynamics.

The simulation of multidimensional vibrational spectroscopy of condensed-phase systems including nuclear quantum effects is challenging since full quantum-mechanical calculations are still intractable for large systems comprising many degrees of freedom. Here, we apply the recently developed double Kubo transform (DKT) methodology in combination with ring-polymer molecular dynamics (RPMD) for evaluating multi-time correlation functions [K. A. Jung et al., J. Chem. Phys. 148, 244105 (2018)], providing a practical method for incorporating nuclear quantum effects in nonlinear spectroscopy of condensed-phase systems. We showcase the DKT approach in the simulation of the fifth-order two-dimensional (2D) Raman spectroscopy of Lennard-Jones liquids as a prototypical example, which involves nontrivial nonlinear spectroscopic observables of systems described by anharmonic potentials. Our results show that the DKT can faithfully reproduce the 2D Raman response of liquid xenon at high temperatures, where the system behaves classically. In contrast, liquid neon at low temperatures exhibits moderate but discernible nuclear quantum effects in the 2D Raman response compared to the responses obtained with classical molecular dynamics approaches. Thus, the DKT formalism in combination with RPMD simulations enables simulations of multidimensional optical spectroscopy of condensed-phase systems that partially account for nuclear quantum effects.

Multi-time formulation of Matsubara dynamics.

Matsubara dynamics has recently emerged as the most general form of a quantum-Boltzmann-conserving classical dynamics theory for the calculation of single-time correlation functions. Here, we present a generalization of Matsubara dynamics for the evaluation of multitime correlation functions. We first show that the Matsubara approximation can also be used to approximate the two-time symmetrized double Kubo transformed correlation function. By a straightforward extension of these ideas to the multitime realm, a multitime Matsubara dynamics approximation can be obtained for the multitime fully symmetrized Kubo transformed correlation function. Although not a practical method, due to the presence of a phase-term, this multitime formulation of Matsubara dynamics represents a benchmark theory for future development of Boltzmann preserving semiclassical approximations to general higher order multitime correlation functions. It also reveals a connection between imaginary time-ordering in the path integral and the classical dynamics of multitime correlation functions.

Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy.

The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.

Electron Transfer Assisted by Vibronic Coupling from Multiple Modes.

Understanding the effect of vibronic coupling on electron transfer (ET) rates is a challenge common to a wide range of applications, from electrochemical synthesis and catalysis to biochemical reactions and solar energy conversion. The Marcus-Jortner-Levich (MJL) theory offers a model of ET rates based on a simple analytic expression with a few adjustable parameters. However, the MJL equation in conjunction with density functional theory (DFT) has yet to be established as a predictive first-principles methodology. A framework is presented for calculating transfer rates modulated by molecular vibrations, that circumvents the steep computational cost which has previously necessitated approximations such as condensing the vibrational manifold into a single empirical frequency. Our DFT-MJL approach provides robust and accurate predictions of ET rates spanning over 4 orders of magnitude in the 106-1010 s-1 range. We evaluate the full MJL equation with a Monte Carlo sampling of the entire active space of thermally accessible vibrational modes, while using no empirical parameters. The contribution to the rate of individual modes is illustrated, providing insight into the interplay between vibrational degrees of freedom and changes in electronic state. The reported findings are valuable for understanding ET rates modulated by multiple vibrational modes, relevant to a broad range of systems within the chemical sciences.