Heat transfer statistics in mixed quantum-classical systems

Population and Energy Transfer Dynamics in an Open Excitonic Quantum Battery

In a previous study, we proposed an open quantum network model of a quantum battery (QB) that possesses dark states owing to its structural exchange symmetries. While in a dark state, the QB is capable of storing an exciton without any environment-induced population losses. However, when the structural exchange symmetry is broken, the QB begins to discharge the exciton towards its exit site. In this article, we start by demonstrating that this QB is not only loss-free with respect to exciton population during the storage phase, but also with respect to the QB energy. We then explore the exciton population and energy transfer dynamics of the QB during the discharge phase over a wide range of site energies, bath temperatures, and bath reorganization energies. Our results shed light on how to optimize the QB's population and energy transfer dynamics for different purposes.

DECIDE: A Deterministic Mixed Quantum-Classical Dynamics Approach

Mixed quantum-classical dynamics provides an efficient way of simulating the dynamics of quantum subsystems coupled to many-body environments. Many processes, including proton-transfer reactions, electron-transfer reactions, and vibrational energy transport, for example, take place in such open systems. The most accurate algorithms for performing mixed quantum-classical simulations require very large ensembles of trajectories to obtain converged expectation values, which is computationally prohibitive for quantum subsystems containing even a few degrees of freedom. The recently developed “Deterministic evolution of coordinates with initial decoupled equations” (DECIDE) method has demonstrated high accuracy and low computational cost for a host of model systems; however, these applications relied on representing the equations of motion in subsystem and adiabatic energy bases. While these representations are convenient for certain systems, the position representation is convenient for many other systems, including systems undergoing proton- and electron-transfer reactions. Thus, in this review, we provide a step-by-step derivation of the DECIDE approach and demonstrate how to cast the DECIDE equations in a quantum harmonic oscillator position basis for a simple one-dimensional (1D) hydrogen bond model. After integrating the DECIDE equations of motion on this basis, we show that the total energy of the system is conserved for this model and calculate various quantities of interest. Limitations of casting the equations in an incomplete basis are also discussed.

Quantum thermal transport beyond second order with the reaction coordinate mapping

Standard quantum master equation techniques, such as the Redfield or Lindblad equations, are perturbative to second order in the microscopic system–reservoir coupling parameter λ. As a result, the characteristics of dissipative systems, which are beyond second order in λ, are not captured by such tools. Moreover, if the leading order in the studied effect is higher-than-quadratic in λ, a second-order description fundamentally fails even at weak coupling. Here, using the reaction coordinate (RC) quantum master equation framework, we are able to investigate and classify higher-than-second-order transport mechanisms. This technique, which relies on the redefinition of the system–environment boundary, allows for the effects of system–bath coupling to be included to high orders. We study steady-state heat current beyond second-order in two models: The generalized spin-boson model with non-commuting system–bath operators and a three-level ladder system. In the latter model, heat enters in one transition and is extracted from a different one. Crucially, we identify two transport pathways: (i) System’s current, where heat conduction is mediated by transitions in the system, with the heat current scaling as j q ∝ λ2 to the lowest order in λ. (ii) Inter-bath current, with the thermal baths directly exchanging energy between them, facilitated by the bridging quantum system. To the lowest order in λ, this current scales as j q ∝ λ4. These mechanisms are uncovered and examined using numerical and analytical tools. We contend that the RC mapping brings, already at the level of the mapped Hamiltonian, much insight into transport characteristics.

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. We observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system-bath entanglement.

Coherences and the thermodynamic uncertainty relation: Insights from quantum absorption refrigerators

The thermodynamic uncertainty relation, originally derived for classical Markov-jump processes, provides a tradeoff relation between precision and dissipation, deepening our understanding of the performance of quantum thermal machines. Here, we examine the interplay of quantum system coherences and heat current fluctuations on the validity of the thermodynamics uncertainty relation in the quantum regime. To achieve the current statistics, we perform a full counting statistics simulation of the Redfield quantum master equation. We focus on steady-state quantum absorption refrigerators where nonzero coherence between eigenstates can either suppress or enhance the cooling power, compared with the incoherent limit. In either scenario, we find enhanced relative noise of the cooling power (standard deviation of the power over the mean) in the presence of system coherence, thereby corroborating the thermodynamic uncertainty relation. Our results indicate that fluctuations necessitate consideration when assessing the performance of quantum coherent thermal machines.

Quantum bath effects on nonequilibrium heat transport in model molecular junctions

Quantum-classical dynamics simulations enable the study of nonequilibrium heat transport in realistic models of molecules coupled to thermal baths. In these simulations, the initial conditions of the bath degrees of freedom are typically sampled from classical distributions. Herein, we investigate the effects of sampling the initial conditions of the thermal baths from quantum and classical distributions on the steady-state heat current in the nonequilibrium spin-boson model-a prototypical model of a single-molecule junction-in different parameter regimes. For a broad range of parameter regimes considered, we find that the steady-state heat currents are ∼1.3-4.5 times larger with the classical bath sampling than with the quantum bath sampling. Using both types of sampling, the steady-state heat currents exhibit turnovers as a function of the bath reorganization energy, with sharper turnovers in the classical case than in the quantum case and different temperature dependencies of the turnover maxima. As the temperature gap between the hot and cold baths increases, we observe an increasing difference in the steady-state heat currents obtained with the classical and quantum bath sampling. In general, as the bath temperatures are increased, the differences between the results of the classical and quantum bath sampling decrease but remain non-negligible at the high bath temperatures. The differences are attributed to the more pronounced temperature dependence of the classical distribution compared to the quantum one. Moreover, we find that the steady-state fluctuation theorem only holds for this model in the Markovian regime when quantum bath sampling is used. Altogether, our results highlight the importance of quantum bath sampling in quantum-classical dynamics simulations of quantum heat transport.

Heat Transport in a Spin-Boson Model at Low Temperatures: A Multilayer Multiconfiguration Time-Dependent Hartree Study

Extending our previous work, quantum dynamic simulations are performed to study low temperature heat transport in a spin-boson model where a two-level subsystem is coupled to two independent harmonic baths. Multilayer multiconfiguration time-dependent Hartree theory is used to numerically evaluate the thermal flux, for which the bath is represented by hundreds to thousands of modes. The simulation results are compared with the approximate Redfield theory approach, and the physics is analyzed versus different physical parameters.

Proposal of a Computational Approach for Simulating Thermal Bosonic Fields in Phase Space

When a quantum field is in contact with a thermal bath, the vacuum state of the field may be generalized to a thermal vacuum state, which takes into account the thermal noise. In thermo field dynamics, this is realized by doubling the dimensionality of the Fock space of the system. Interestingly, the representation of thermal noise by means of an augmented space is also found in a distinctly different approach based on the Wigner transform of both the field operators and density matrix, which we pursue here. Specifically, the thermal noise is introduced by augmenting the classical-like Wigner phase space by means of Nosé–Hoover chain thermostats, which can be readily simulated on a computer. In this paper, we illustrate how this may be achieved and discuss how non-equilibrium quantum thermal distributions of the field modes can be numerically simulated.

Thermodynamic uncertainty relation in thermal transport

We use the fundamental nonequilibrium steady-state fluctuation symmetry and derive a condition on the validity of the thermodynamic uncertainty relation (TUR) in thermal transport problems, classical and quantum alike. We test this condition and study the breakdown of the TUR in different thermal transport junctions of bosonic and electronic degrees of freedom. We prove that the TUR is valid in harmonic oscillator junctions. In contrast, in the nonequilibrium spin-boson model, which realizes many-body effects, it is satisfied in the Markovian limit, but violations arise as we tune (reduce) the cutoff frequency of the thermal baths, thus observing non-Markovian dynamics. We consider heat transport by noninteracting electrons in a tight-binding chain model. We show that the TUR is feasibly violated by tuning, e.g., the hybridization energy of the chain to the metal leads. These results manifest that the validity of the TUR relies on the statistics of the participating carriers, their interaction, and the nature of their couplings to the macroscopic contacts (metal electrodes and phonon baths).

Nonequilibrium heat transport in a molecular junction: A mixed quantum-classical approach

In a recent study [J. Liu et al., J. Chem. Phys. 149, 224104 (2018)], we developed a general mixed quantum-classical framework for studying heat transport through molecular junctions, in which the junction molecule is treated quantum mechanically and the thermal reservoirs to which the molecule is coupled are treated classically. This framework yields expressions for the transferred heat and steady-state heat current, which could be calculated using a variety of mixed quantum-classical dynamics methods. In this work, we use the recently developed "Deterministic Evolution of Coordinates with Initial Decoupled Equations" (DECIDE) method for calculating the steady-state heat current in the nonequilibrium spin-boson model in a variety of parameter regimes. Our results are compared and contrasted with those obtained using the numerically exact multilayer multiconfiguration time-dependent Hartree approach, and using approximate methods, including mean field theory, Redfield theory, and adiabatic mixed quantum-classical dynamics. Despite some quantitative differences, the DECIDE method performs quite well, is capable of capturing the expected trends in the steady-state heat current, and, overall, outperforms the approximate methods. These results hold promise for DECIDE simulations of nonequilibrium heat transport in realistic models of nanoscale systems.

Path-integral methodology and simulations of quantum thermal transport: Full counting statistics approach

We develop and test a computational framework to study heat exchange in interacting, nonequilibrium open quantum systems. Our iterative full counting statistics path integral (iFCSPI) approach extends a previously well-established influence functional path integral method, by going beyond reduced system dynamics to provide the cumulant generating function of heat exchange. The method is straightforward; we implement it for the nonequilibrium spin boson model to calculate transient and long-time observables, focusing on the steady-state heat current flowing through the system under a temperature difference. Results are compared to perturbative treatments and demonstrate good agreement in the appropriate limits. The challenge of converging nonequilibrium quantities, currents and high order cumulants, is discussed in detail. The iFCSPI, a numerically exact technique, naturally captures strong system-bath coupling and non-Markovian effects of the environment. As such, it is a promising tool for probing fundamental questions in quantum transport and quantum thermodynamics.