Non-equilibrium spin-boson model: Counting statistics and the heat exchange fluctuation theorem

Nonequilibrium steady state full counting statistics in the noncrossing approximation

Quantum transport is often characterized not just by mean observables like the particle or energy current but by their fluctuations and higher moments, which can act as detailed probes of the physical mechanisms at play. However, relatively few theoretical methods are able to access the full counting statistics (FCS) of transport processes through electronic junctions in strongly correlated regimes. While most experiments are concerned with steady state properties, most accurate theoretical methods rely on computationally expensive propagation from a tractable initial state. Here, we propose a simple approach for computing the FCS through a junction directly at the steady state, utilizing the propagator noncrossing approximation. Compared to time propagation, our method offers reduced computational cost at the same level of approximation, but the idea can also be used within other approximations or as a basis for numerically exact techniques. We demonstrate the method's capabilities by investigating the impact of lead dimensionality on electronic transport in the nonequilibrium Anderson impurity model at the onset of Kondo physics. Our results reveal a distinct signature of one dimensional leads in the noise and Fano factor not present for other dimensionalities, showing the potential of FCS measurements as a probe of the environment surrounding a quantum dot.

Generalized Quantum Fluctuation Theorem for Energy Exchange

The nonequilibrium fluctuation relation is a cornerstone of quantum thermodynamics. It is widely believed that the system-bath heat exchange obeys the famous Jarzynski-Wójcik fluctuation theorem. However, this theorem is established in the Born-Markovian approximation under the weak-coupling condition. Via studying the energy exchange between a harmonic oscillator and its coupled bath in the non-Markovian dynamics, we establish a generalized quantum fluctuation theorem for energy exchange being valid for arbitrary coupling strength. The Jarzynski-Wójcik fluctuation theorem is recovered in the weak-coupling limit. We also find the average energy exchange exhibits rich nonequilibrium characteristics when different numbers of system-bath bound states are formed, which suggests a useful way to control the quantum heat. Deepening our understanding of the fluctuation relation in quantum thermodynamics, our result lays the foundation to design high-efficiency quantum heat engines.

Bath-induced interactions and transient dynamics in open quantum systems at strong coupling: Effective Hamiltonian approach

Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.

Challenges in molecular dynamics simulations of heat exchange statistics

We study heat exchange in temperature-biased metal-molecule-metal molecular junctions by employing the molecular dynamics simulator LAMMPS. Generating the nonequilibrium steady state with Langevin thermostats at the boundaries of the junction, we show that the average heat current across a gold-alkanedithiol-gold nanojunction behaves physically, with the thermal conductance value matching the literature. In contrast, the full probability distribution function for heat exchange, as generated by the simulator, violates the fundamental fluctuation symmetry for entropy production. We trace this failure back to the implementation of the thermostats and the expression used to calculate the heat exchange. To rectify this issue and produce the correct statistics, we introduce single-atom thermostats as an alternative to conventional many-atom thermostats. Once averaging heat exchange over the hot and cold thermostats, this approach successfully generates the correct probability distribution function, which we use to study the behavior of both the average heat current and its noise. We further examine the thermodynamic uncertainty relation in the molecular junction and show that it holds, albeit demonstrating nontrivial trends. Our study points to the need to carefully implement nonequilibrium molecular dynamics solvers in atomistic simulation software tools for future investigations of noise phenomena in thermal transport.

Simulations of heat transport in single-molecule junctions: Investigations of the thermal diode effect

With the objective of understanding microscopic principles governing thermal energy flow in nanojunctions, we study phononic heat transport through metal-molecule-metal junctions using classical molecular dynamics (MD) simulations. Considering a single-molecule gold-alkanedithiol-gold junction, we first focus on aspects of method development and compare two techniques for calculating thermal conductance: (i) The Reverse Nonequilibrium MD (RNEMD) method, where heat is inputted and extracted at a constant rate from opposite metals. In this case, the thermal conductance is calculated from the nonequilibrium temperature profile that is created at the junction. (ii) The Approach-to-Equilibrium MD (AEMD) method, with the thermal conductance of the junction obtained from the equilibration dynamics of the metals. In both methods, simulations of alkane chains of a growing size display an approximate length-independence of the thermal conductance, with calculated values matching computational and experimental studies. The RNEMD and AEMD methods offer different insights, and we discuss their benefits and shortcomings. Assessing the potential application of molecular junctions as thermal diodes, alkane junctions are made spatially asymmetric by modifying their contact regions with the bulk, either by using distinct endgroups or by replacing one of the Au contacts with Ag. Anharmonicity is built into the system within the molecular force-field. We find that, while the temperature profile strongly varies (compared with the gold-alkanedithiol-gold junctions) due to these structural modifications, the thermal diode effect is inconsequential in these systems-unless one goes to very large thermal biases. This finding suggests that one should seek molecules with considerable internal anharmonic effects for developing nonlinear thermal devices.

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.

Strong system-bath coupling effects in quantum absorption refrigerators

We study the performance of three-level quantum absorption refrigerators, paradigmatic autonomous quantum thermal machines, and reveal central impacts of strong couplings between the working system and the thermal baths. Using the reaction coordinate quantum master equation method, which treats system-bath interactions beyond weak coupling, we demonstrate that in a broad range of parameters the cooling window at strong coupling can be captured by a weak-coupling theory, albeit with parameters renormalized by the system-bath coupling energy. As a result, at strong system-bath couplings the window of cooling is significantly reshaped compared to predictions of weak-coupling treatments. We further show that strong coupling admits direct transport pathways between the thermal reservoirs. Such beyond-second-order transport mechanisms are typically detrimental to the performance of quantum thermal machines. Our study reveals that it is inadequate to claim for either a suppression or an enhancement of the cooling performance as one increases system-bath coupling-when analyzed against a single parameter and in a limited domain. Rather, a comprehensive approach should be adopted so as to uncover the reshaping of the operational window.

Full counting statistics of the particle currents through a Kitaev chain and the exchange fluctuation theorem

Exchange fluctuation theorems (XFTs) describe a fundamental symmetry relation for particle and energy exchange between several systems. Here we study the XFTs of a Kitaev chain connected to two reservoirs at the same temperature but different bias. By varying the parameters in the Kitaev chain model, we calculate analytically the full counting statistics of the transport current and formulate the corresponding XFTs for multiple current components. We also demonstrate the XFTs with numerical results. We find that due to the presence of the U(1) symmetry breaking terms in the Hamiltonian of the Kitaev chain, various forms of the XFTs emerge, and they can be interpreted in terms of various well-known transport processes.

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.

Large deviations and fluctuation theorem for the quantum heat current in the spin-boson model

We study the heat current flowing between two baths consisting of harmonic oscillators interacting with a qubit through a spin-boson coupling. An explicit expression for the generating function of the total heat flowing between the right and left baths is derived by evaluating the corresponding Feynman-Vernon path integral by performing the noninteracting blip approximation (NIBA). We recover the known expression, obtained by using the polaron transform. This generating function satisfies the Gallavotti-Cohen fluctuation theorem, both before and after performing the NIBA. We also verify that the heat conductance is proportional to the variance of the heat current, retrieving the well-known fluctuation dissipation relation. Finally, we present numerical results for the heat current.

Cooling condition for multilevel quantum absorption refrigerators

Models for quantum absorption refrigerators serve as test beds for exploring concepts and developing methods in quantum thermodynamics. Here we depart from the minimal, ideal design and consider a generic multilevel model for a quantum absorption refrigerator, which potentially suffers from lossy processes. Based on a full-counting statistics approach, we derive a formal cooling condition for the refrigerator, which can be feasibly evaluated analytically and numerically. We exemplify our approach on a three-level model for a quantum absorption refrigerator that suffers from different forms of nonideality (heat leakage, competition between different cooling pathways) and examine the cooling current with different designs. This study assists in identifying the cooling window of imperfect thermal machines.

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.

Heat transfer statistics in mixed quantum-classical systems

The modelling of quantum heat transfer processes at the nanoscale is crucial for the development of energy harvesting and molecular electronic devices. Herein, we adopt a mixed quantum-classical description of a device, in which the open subsystem of interest is treated quantum mechanically and the surrounding heat baths are treated in a classical-like fashion. By introducing such a mixed quantum-classical description of the composite system, one is able to study the heat transfer between the subsystem and bath from a closed system point of view, thereby avoiding simplifying assumptions related to the bath time scale and subsystem-bath coupling strength. In particular, we adopt the full counting statistics approach to derive a general expression for the moment generating function of heat in systems whose dynamics are described by the quantum-classical Liouville equation (QCLE). From this expression, one can deduce expressions for the dynamics of the average heat and heat current, which may be evaluated using numerical simulations. Due to the approximate nature of the QCLE, we also find that the steady state fluctuation symmetry holds up to order ℏ for systems whose subsystem-bath couplings and baths go beyond bilinear and harmonic, respectively. To demonstrate the approach, we consider the nonequilibrium spin boson model and simulate its time-dependent average heat and heat current under various conditions.

Frequency-dependent current noise in quantum heat transfer: A unified polaron calculation

To investigate frequency-dependent current noise (FDCN) in open quantum systems at steady states, we present a theory which combines Markovian quantum master equations with a finite time full counting statistics. Our formulation of the FDCN generalizes previous zero-frequency expressions and can be viewed as an application of MacDonald's formula for electron transport to heat transfer. As a demonstration, we consider the paradigmatic example of quantum heat transfer in the context of a non-equilibrium spin-boson model. We adopt a recently developed polaron-transformed Redfield equation which allows us to accurately investigate heat transfer with arbitrary system-reservoir coupling strength, arbitrary values of spin bias, and temperature differences. We observe a turn-over of FDCN in the intermediate coupling regimes, similar to the zero-frequency case. We find that the FDCN with varying coupling strengths or bias displays a universal Lorentzian-shape scaling form in the weak coupling regime, and a white noise spectrum emerges with zero bias in the strong coupling regime due to distinctive spin dynamics. We also find that the bias can suppress the FDCN in the strong coupling regime, in contrast to its zero-frequency counterpart which is insensitive to bias changes. Furthermore, we utilize the Saito-Utsumi relation as a benchmark to validate our theory and study the impact of temperature differences at finite frequencies. Together, our results provide detailed dissections of the finite time fluctuation of heat current in open quantum systems.

Current fluctuations in quantum absorption refrigerators

Absorption refrigerators transfer thermal energy from a cold bath to a hot bath without input power by utilizing heat from an additional "work" reservoir. Particularly interesting is a three-level design for a quantum absorption refrigerator, which can be optimized to reach the maximal (Carnot) cooling efficiency. Previous studies of three-level chillers focused on the behavior of the averaged cooling current. Here, we go beyond that and study the full counting statistics of heat exchange in a three-level chiller model. We explain how to obtain the complete cumulant generating function of the refrigerator in a steady state, then derive a partial cumulant generating function, which yields closed-form expressions for both the averaged cooling current and its noise. Our analytical results and simulations are beneficial for the design of nanoscale engines and cooling systems far from equilibrium, with their performance optimized according to different criteria, efficiency, power, fluctuations, and dissipation.

Energy transfer in the nonequilibrium spin-boson model: From weak to strong coupling

To explore energy transfer in the nonequilibrium spin-boson model (NESB) from weak to strong system-bath coupling regimes, we propose a polaron-transformed nonequilibrium Green's function (NEGF) method. By combining the polaron transformation, we are able to treat the system-bath coupling nonperturbatively, thus in direct contrast to conventionally used NEGF methods which take the system-bath coupling as a perturbation. The Majorana-fermion representation is further utilized to evaluate terms in the Dyson series. This method not only allows us to deal with weak as well as strong coupling regimes but also enables an investigation on the role of bias in the energy transfer. As an application of the method, we study an Ohmic NESB. For an unbiased spin system, our energy current result smoothly bridges predictions of two benchmarks, namely, the quantum master equation and the nonequilibrium noninteracting blip approximation, a considerable improvement over existing theories. In case of a biased spin system, we found a bias-induced nonmonotonic behavior of the energy conductance in the intermediate coupling regime, resulting from the resonant character of the energy transfer. This finding may offer a nontrivial quantum control knob over energy transfer at the nanoscale.

Electron transfer across a thermal gradient

Charge transfer is a fundamental process that underlies a multitude of phenomena in chemistry and biology. Recent advances in observing and manipulating charge and heat transport at the nanoscale, and recently developed techniques for monitoring temperature at high temporal and spatial resolution, imply the need for considering electron transfer across thermal gradients. Here, a theory is developed for the rate of electron transfer and the associated heat transport between donor-acceptor pairs located at sites of different temperatures. To this end, through application of a generalized multidimensional transition state theory, the traditional Arrhenius picture of activation energy as a single point on a free energy surface is replaced with a bithermal property that is derived from statistical weighting over all configurations where the reactant and product states are equienergetic. The flow of energy associated with the electron transfer process is also examined, leading to relations between the rate of heat exchange among the donor and acceptor sites as functions of the temperature difference and the electronic driving bias. In particular, we find that an open electron transfer channel contributes to enhanced heat transport between sites even when they are in electronic equilibrium. The presented results provide a unified theory for charge transport and the associated heat conduction between sites at different temperatures.

Strongly Coupled Quantum Heat Machines

Energy conversion of heat into work at the quantum level is modeled by quantum heat machines (QHMs) generally assumed to operate at weak coupling to the baths. This supposition is grounded in the separability principle between systems and allows the derivation of the evolution equation. In the weak coupling regime, the machine's output is limited by the coupling strength, restricting their application. Seeking to overcome this limitation, we analyze QHMs in the virtually unexplored strong coupling regime here, where separability, as well as other standard thermodynamic assumptions, may no longer hold. We show that strongly coupled QHMs may be as efficient as their weakly coupled counterparts. In addition, we find a novel turnover behavior where their output saturates and disappears in the limit of ultrastrong coupling.

Heat transfer in the spin-boson model: a comparative study in the incoherent tunneling regime

We study the transfer of heat in the nonequilibrium spin-boson model with an Ohmic dissipation. In the nonadiabatic limit we derive a formula for the thermal conductance based on a rate equation formalism at the level of the noninteracting blip approximation, valid for temperatures T>T(K), with T(K) as the Kondo temperature. We evaluate this expression analytically assuming either weak or strong couplings, and demonstrate that our results agree with exact relations. Far-from-equilibrium situations are further examined, showing a close correspondence to the linear response limit.

Nonadiabatic effect on the quantum heat flux control

We provide a general formula of quantum transfer that includes the nonadiabatic effect under periodic environmental modulation by using full counting statistics in Hilbert-Schmidt space. Applying the formula to an anharmonic junction model that interacts with two bosonic environments within the Markovian approximation, we find that the quantum transfer is divided into the adiabatic (dynamical and geometrical phases) and nonadiabatic contributions. This extension shows the dependence of quantum transfer on the initial condition of the anharmonic junction just before the modulation, as well as the characteristic environmental parameters such as interaction strength and cut-off frequency of spectral density. We show that the nonadiabatic contribution represents the reminiscent effect of past modulation including the transition from the initial condition of the anharmonic junction to a steady state determined by the very beginning of the modulation. This enables us to tune the frequency range of modulation, whereby we can obtain the quantum flux corresponding to the geometrical phase by setting the initial condition of the anharmonic junction.

Kondo signature in heat transfer via a local two-state system

We study the Kondo effect in heat transport via a local two-state system. This system is described by the spin-boson Hamiltonian with Ohmic dissipation, which can be mapped onto the Kondo model with anisotropic exchange coupling. We calculate thermal conductance by the Monte Carlo method based on the exact formula. Thermal conductance has a scaling form κ=(k(B)(2)T(K)/ħ)f(α,T/T(K)), where T(K) and α indicate the Kondo temperature and dimensionless coupling strength, respectively. Temperature dependence of conductance is classified by the Kondo temperature as κ is proportional to (T/T(K))(3) for T<<T(K) and κ is proportional to (kB)T/ħω(c))(2α-1) for T>>T(K). Similarities to the Kondo signature in electric transport are discussed.

Stochastic thermodynamics, fluctuation theorems and molecular machines

Stochastic thermodynamics as reviewed here systematically provides a framework for extending the notions of classical thermodynamics such as work, heat and entropy production to the level of individual trajectories of well-defined non-equilibrium ensembles. It applies whenever a non-equilibrium process is still coupled to one (or several) heat bath(s) of constant temperature. Paradigmatic systems are single colloidal particles in time-dependent laser traps, polymers in external flow, enzymes and molecular motors in single molecule assays, small biochemical networks and thermoelectric devices involving single electron transport. For such systems, a first-law like energy balance can be identified along fluctuating trajectories. For a basic Markovian dynamics implemented either on the continuum level with Langevin equations or on a discrete set of states as a master equation, thermodynamic consistency imposes a local-detailed balance constraint on noise and rates, respectively. Various integral and detailed fluctuation theorems, which are derived here in a unifying approach from one master theorem, constrain the probability distributions for work, heat and entropy production depending on the nature of the system and the choice of non-equilibrium conditions. For non-equilibrium steady states, particularly strong results hold like a generalized fluctuation-dissipation theorem involving entropy production. Ramifications and applications of these concepts include optimal driving between specified states in finite time, the role of measurement-based feedback processes and the relation between dissipation and irreversibility. Efficiency and, in particular, efficiency at maximum power can be discussed systematically beyond the linear response regime for two classes of molecular machines, isothermal ones such as molecular motors, and heat engines such as thermoelectric devices, using a common framework based on a cycle decomposition of entropy production.

Full-counting statistics of heat transport in harmonic junctions: transient, steady states, and fluctuation theorems

We study the statistics of heat transferred in a given time interval t_{M}, through a finite harmonic chain, called the center, which is connected to two heat baths, the left (L) and the right (R), that are maintained at two temperatures. The center atoms are driven by external time-dependent forces. We calculate the cumulant generating function (CGF) for the heat transferred out of the left lead, Q_{L}, based on the two-time quantum measurement concept and using the nonequilibrium Green's function method. The CGF can be concisely expressed in terms of Green's functions of the center and an argument-shifted self-energy of the lead. The expression of the CGF is valid in both transient and steady-state regimes. We consider three initial conditions for the density operator and show numerically, for a one-atom junction, how their transient behaviors differ from each other but, finally, approach the same steady state, independent of the initial distributions. We also derive the CGF for the joint probability distribution P(Q_{L},Q_{R}), and discuss the correlations between Q_{L} and Q_{R}. We calculate the CGF for total entropy production in the reservoirs. In the steady state we explicitly show that the CGFs obey steady-state fluctuation theorems. We obtain classical results by taking ℏ→0. We also apply our method to the counting of the electron number and electron energy, for which the associated self-energy is obtained from the usual lead self-energy by multiplying a phase and shifting the contour time, respectively.