Statistics of work performed on a forced quantum oscillator

Work distribution for unzipping processes

A simple zipper model is introduced, representing in a simplified way, e.g., the folded DNA double helix or hairpin structures in RNA. The double stranded hairpin is connected to a heat bath at temperature T and subject to an external force f, which couples to the free length L of the unzipped sequence. The leftmost zipped position can be seen as the position of a random walker in a special external field. Increasing the force leads to a zipping-unzipping first-order phase transition at a critical force f_{c}(T) in the thermodynamic limit of a very large chain. We compute analytically, as a function of temperature T and force f, the full distribution P(L) of free lengths in the thermodynamic limit and show that it is qualitatively very different for ff_{c}. Next we consider quasistatic work processes where the force is incremented according to a linear protocol. Having obtained P(L) already allows us to derive an analytical expression for the work distribution P(W) in the zipped phase f Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Peter Werner Institut für Physik, Universität Oldenburg, 26111 Oldenburg, Germany
Alexander K Hartmann Institut für Physik, Universität Oldenburg, 26111 Oldenburg, Germany
Satya N Majumdar Laboratoire de Physique Théorique et Modèles Statistique, Université Paris-Saclay, 91405 Orsay CEDEX, France

Collapse

Quasistatic work processes: When slowness implies certainty

Two approaches are outlined to characterize the fluctuation behavior of work applied to a system by a slow change of a parameter. One approach uses the adiabatic theorems of quantum and classical mechanics, and the other one is based on the behavior of the correlations of the generalized coordinate that is conjugate to the changed parameter. Criteria are obtained under which the work done on small thermally isolated as well as on open systems ceases to fluctuate in a quasistatic process.

Nonequilibrium Green's Function's Approach to the Calculation of Work Statistics

The calculation of work distributions in a quantum many-body system is of significant importance and also of formidable difficulty in the field of nonequilibrium quantum statistical mechanics. To solve this problem, inspired by the Schwinger-Keldysh formalism, we propose the contour-integral formulation for work statistics. Based on this contour integral, we show how to do the perturbation expansion of the characteristic function of work (CFW) and obtain the approximate expression of the CFW to the second order of the work parameter for an arbitrary system under a perturbative protocol. We also demonstrate the validity of fluctuation theorems by utilizing the Kubo-Martin-Schwinger condition. Finally, we use noninteracting identical particles in a forced harmonic potential as an example to demonstrate the powerfulness of our approach.

Work Statistics across a Quantum Phase Transition

We investigate the statistics of the work performed during a quench across a quantum phase transition using the adiabatic perturbation theory when the system is characterized by independent quasiparticles and the "single-excitation" approximation is assumed. It is shown that all the cumulants of work exhibit universal scaling behavior analogous to the Kibble-Zurek scaling for the average density of defects. Two kinds of transformations are considered: quenches between two gapped phases in which a critical point is traversed, and quenches that end near the critical point. In contrast to the scaling behavior of the density of defects, the scaling behavior of the cumulants of work are shown to be qualitatively different for these two kinds of quenches. However, in both cases the corresponding exponents are fully determined by the dimension of the system and the critical exponents of the transition, as in the traditional Kibble-Zurek mechanism (KZM). Thus, our study deepens our understanding about the nonequilibrium dynamics of a quantum phase transition by revealing the imprint of the KZM on the work statistics.

Path-integral approach to the calculation of the characteristic function of work

Work statistics characterizes important features of a nonequilibrium thermodynamic process, but the calculation of the work statistics in an arbitrary nonequilibrium process is usually a cumbersome task. In this work, we study the work statistics in quantum systems by employing Feynman's path-integral approach. We derive the analytical work distributions of two prototype quantum systems. The results are proved to be equivalent to the results obtained based on Schrödinger's formalism. We also calculate the work distributions in their classical counterparts by employing the path-integral approach. Our study demonstrates the effectiveness of the path-integral approach for the calculation of work statistics in both quantum and classical thermodynamics, and brings important insights to the understanding of the trajectory work in quantum systems.

Symmetry and its breaking in a path-integral approach to quantum Brownian motion

We study the Caldeira-Leggett model where a quantum Brownian particle interacts with an environment or a bath consisting of a collection of harmonic oscillators in the path-integral formalism. Compared to the contours that the paths take in the conventional Schwinger-Keldysh formalism, the paths in our study are deformed in the complex time plane as suggested by the recent study by C. Aron, G. Biroli, and L. F. Cugliandolo [SciPost Phys. 4, 008 (2018)10.21468/SciPostPhys.4.1.008]. This is done to investigate the connection between the symmetry properties in the Schwinger-Keldysh action and the equilibrium or nonequilibrium nature of the dynamics in an open quantum system. We derive the influence functional explicitly in this setting, which captures the effect of the coupling to the bath. We show that in equilibrium the action and the influence functional are invariant under a set of transformations of path-integral variables. The fluctuation-dissipation relation is obtained as a consequence of this symmetry. When the system is driven by an external time-dependent protocol, the symmetry is broken. From the terms that break the symmetry, we derive a quantum Jarzynski-like equality for a quantum mechanical worklike quantity given as a function of fluctuating quantum trajectory. In the classical limit, the transformations becomes those used in the functional integral formalism of the classical stochastic thermodynamics to derive the classical fluctuation theorem.

Computing characteristic functions of quantum work in phase space

In phase space, we analytically obtain the characteristic functions (CFs) of a forced harmonic oscillator [Talkner et al., Phys. Rev. E 75, 050102(R) (2007)PLEEE81539-375510.1103/PhysRevE.75.050102], a time-dependent mass and frequency harmonic oscillator [Deffner and Lutz, Phys. Rev. E 77, 021128 (2008)PLEEE81539-375510.1103/PhysRevE.77.021128], and coupled harmonic oscillators under driving forces in a simple and unified way. For general quantum systems, a numerical method that approximates the CFs to ℏ^{2} order is proposed. We exemplify the method with a time-dependent frequency harmonic oscillator and a family of quantum systems with time-dependent even power-law potentials.

Quantum corrections of work statistics in closed quantum systems

We investigate quantum corrections to the classical work characteristic function (CF) as a semiclassical approximation to the full quantum work CF. In addition to explicitly establishing the quantum-classical correspondence of the Feynman-Kac formula, we find that these quantum corrections must be in even powers of ℏ. Exact formulas of the lowest corrections (ℏ^{2}) are proposed, and their physical origins are clarified. We calculate the work CFs for a forced harmonic oscillator and a forced quartic oscillator respectively to illustrate our results.

Path Integral Approach to Quantum Thermodynamics

Work belongs to the most basic notions in thermodynamics but it is not well understood in quantum systems, especially in open quantum systems. By introducing a novel concept of the work functional along an individual Feynman path, we invent a new approach to study thermodynamics in the quantum regime. Using the work functional, we derive a path integral expression for the work statistics. By performing the ℏ expansion, we analytically prove the quantum-classical correspondence of the work statistics. In addition, we obtain the quantum correction to the classical fluctuating work. We can also apply this approach to an open quantum system in the strong coupling regime described by the quantum Brownian motion model. This approach provides an effective way to calculate the work in open quantum systems by utilizing various path integral techniques. As an example, we calculate the work statistics for a dragged harmonic oscillator in both isolated and open quantum systems.

Quantum work for sudden quenches in Gaussian random Hamiltonians

In the context of nonequilibrium quantum thermodynamics, variables like work behave stochastically. A particular definition of the work probability density function (pdf) for coherent quantum processes allows the verification of the quantum version of the celebrated fluctuation theorems, due to Jarzynski and Crooks, that apply when the system is driven away from an initial equilibrium thermal state. Such a particular pdf depends basically on the details of the initial and final Hamiltonians, on the temperature of the initial thermal state, and on how some external parameter is changed during the coherent process. Using random matrix theory we derive a simple analytic expression that describes the general behavior of the work characteristic function G(u), associated with this particular work pdf for sudden quenches, valid for all the traditional Gaussian ensembles of Hamiltonians matrices. This formula well describes the general behavior of G(u) calculated from single draws of the initial and final Hamiltonians in all ranges of temperatures.

Quantum work fluctuations in connection with the Jarzynski equality

A result of great theoretical and experimental interest, the Jarzynski equality predicts a free energy change ΔF of a system at inverse temperature β from an ensemble average of nonequilibrium exponential work, i.e., 〈e^{-βW}〉=e^{-βΔF}. The number of experimental work values needed to reach a given accuracy of ΔF is determined by the variance of e^{-βW}, denoted var(e^{-βW}). We discover in this work that var(e^{-βW}) in both harmonic and anharmonic Hamiltonian systems can systematically diverge in nonadiabatic work protocols, even when the adiabatic protocols do not suffer from such divergence. This divergence may be regarded as a type of dynamically induced phase transition in work fluctuations. For a quantum harmonic oscillator with time-dependent trapping frequency as a working example, any nonadiabatic work protocol is found to yield a diverging var(e^{-βW}) at sufficiently low temperatures, markedly different from the classical behavior. The divergence of var(e^{-βW}) indicates the too-far-from-equilibrium nature of a nonadiabatic work protocol and makes it compulsory to apply designed control fields to suppress the quantum work fluctuations in order to test the Jarzynski equality.

Quantum-to-classical transition in the work distribution for chaotic systems

The work distribution is a fundamental quantity in nonequilibrium thermodynamics mainly due to its connection with fluctuation theorems. Here, we develop a semiclassical approximation to the work distribution for a quench process in chaotic systems that provides a link between the quantum and classical work distributions. The approach is based on the dephasing representation of the quantum Loschmidt echo and on the quantum ergodic conjecture, which states that the Wigner function of a typical eigenstate of a classically chaotic Hamiltonian is equidistributed on the energy shell. Using numerical simulations, we show that our semiclassical approximation accurately describes the quantum distribution as the temperature is increased.

Work distributions for random sudden quantum quenches

The statistics of work performed on a system by a sudden random quench is investigated. Considering systems with finite dimensional Hilbert spaces we model a sudden random quench by randomly choosing elements from a Gaussian unitary ensemble (GUE) consisting of Hermitian matrices with identically, Gaussian distributed matrix elements. A probability density function (pdf) of work in terms of initial and final energy distributions is derived and evaluated for a two-level system. Explicit results are obtained for quenches with a sharply given initial Hamiltonian, while the work pdfs for quenches between Hamiltonians from two independent GUEs can only be determined in explicit form in the limits of zero and infinite temperature. The same work distribution as for a sudden random quench is obtained for an adiabatic, i.e., infinitely slow, protocol connecting the same initial and final Hamiltonians.

Quantum Performance of Thermal Machines over Many Cycles

The performance of quantum heat engines is generally based on the analysis of a single cycle. We challenge this approach by showing that the total work performed by a quantum engine need not be proportional to the number of cycles. Furthermore, optimizing the engine over multiple cycles leads to the identification of scenarios with a quantum enhancement. We demonstrate our findings with a quantum Otto engine based on a two-level system as the working substance that supplies power to an external oscillator.

Quantum work statistics of charged Dirac particles in time-dependent fields

The quantum Jarzynski equality is an important theorem of modern quantum thermodynamics. We show that the Jarzynski equality readily generalizes to relativistic quantum mechanics described by the Dirac equation. After establishing the conceptual framework we solve a pedagogical, yet experimentally relevant, system analytically. As a main result we obtain the exact quantum work distributions for charged particles traveling through a time-dependent vector potential evolving under Schrödinger as well as under Dirac dynamics, and for which the Jarzynski equality is verified. Special emphasis is put on the conceptual and technical subtleties arising from relativistic quantum mechanics.

Jarzynski Equality in PT-Symmetric Quantum Mechanics

We show that the quantum Jarzynski equality generalizes to PT-symmetric quantum mechanics with unbroken PT symmetry. In the regime of broken PT symmetry, the Jarzynski equality does not hold as also the CPT norm is not preserved during the dynamics. These findings are illustrated for an experimentally relevant system-two coupled optical waveguides. It turns out that for these systems the phase transition between the regimes of unbroken and broken PT symmetry is thermodynamically inhibited as the irreversible work diverges at the critical point.

Degenerate optimal paths in thermally isolated systems

We present an analysis of the work performed on a system of interest that is kept thermally isolated during the switching of a control parameter. We show that there exists, for a certain class of systems, a finite-time family of switching protocols for which the work is equal to the quasistatic value. These optimal paths are obtained within linear response for systems initially prepared in a canonical distribution. According to our approach, such protocols are composed of a linear part plus a function which is odd with respect to time reversal. For systems with one degree of freedom, we claim that these optimal paths may also lead to the conservation of the corresponding adiabatic invariant. This points to an interesting connection between work and the conservation of the volume enclosed by the energy shell. To illustrate our findings, we solve analytically the harmonic oscillator and present numerical results for certain anharmonic examples.

Interference of identical particles and the quantum work distribution

Quantum-mechanical particles in a confining potential interfere with each other while undergoing thermodynamic processes far from thermal equilibrium. By evaluating the corresponding transition probabilities between many-particle eigenstates we obtain the quantum work distribution function for identical bosons and fermions, which we compare with the case of distinguishable particles. We find that the quantum work distributions for bosons and fermions significantly differ at low temperatures, while, as expected, at high temperatures the work distributions converge to the classical expression. These findings are illustrated with two analytically solvable examples, namely the time-dependent infinite square well and the parametric harmonic oscillator.

Work and efficiency of quantum Otto cycles in power-law trapping potentials

We study the performance of a quantum Otto cycle operating in trapping potentials of different shapes. We show that, while both the mean work output and the efficiency of two Otto cycles in different trapping potentials can be made equal, the work probability distribution will still be strongly affected by the difference in structure of the energy levels. To exemplify our results, we study the family of potentials of the form V(t)(x) ∼ x(2q). This family of potentials possesses a simple scaling property that allows for analytical insights into the efficiency and work output of the cycle. We perform a comparison of quantum Otto cycles in various physically relevant scenarios and find that in certain instances, the efficiency of the cycle is greater when using potentials with larger values of q, while in other cases, the efficiency is greater with harmonic traps.

Microcanonical work and fluctuation relations for an open system: An exactly solvable model

We calculate the probability distribution of work for an exactly solvable model of a system interacting with its environment. The system of interest is a harmonic oscillator with a time-dependent control parameter, the environment is modeled by N-independent harmonic oscillators with arbitrary frequencies, and the system-environment coupling is bilinear and not necessarily weak. The initial conditions of the combined system and environment are sampled from a microcanonical distribution and the system is driven out of equilibrium by changing the control parameter according to a prescribed protocol. In the limit of infinitely large environment, i.e., N→∞, we recover the nonequilibrium work relation and Crooks's fluctuation relation. Moreover, the microcanonical Crooks relation is verified for finite environments. Finally, we show the equivalence of multitime correlation functions of the system in the infinite environment limit for canonical and microcanonical ensembles.

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.

Quantum-trajectory approach to the stochastic thermodynamics of a forced harmonic oscillator

I formulate a quantum stochastic thermodynamics for the quantum trajectories of a continuously monitored forced harmonic oscillator coupled to a thermal reservoir. Consistent trajectory-dependent definitions are introduced for work, heat, and entropy, through engineering the thermal reservoir from a sequence of two-level systems. Within this formalism the connection between irreversibility and entropy production is analyzed and confirmed by proving a detailed fluctuation theorem for quantum trajectories. Finally, possible experimental verifications are discussed.

Validity of nonequilibrium work relations for the rapidly expanding quantum piston

Recent work by Teifel and Mahler [Eur. Phys. J. B 75, 275 (2010)] raises legitimate concerns regarding the validity of quantum nonequilibrium work relations in processes involving moving hard walls. We study this issue in the context of the rapidly expanding one-dimensional quantum piston. Utilizing exact solutions of the time-dependent Schrödinger equation, we find that the evolution of the wave function can be decomposed into static and dynamic components, which have simple semiclassical interpretations in terms of particle-piston collisions. We show that nonequilibrium work relations remain valid at any finite piston speed, provided both components are included, and we study explicitly the work distribution for this model system.

Nonequilibrium work statistics of an Aharonov-Bohm flux

We investigate the statistics of work performed on a noninteracting electron gas confined in a ring as a threaded magnetic field is turned on. For an electron gas initially prepared in a grand canonical state it is demonstrated that the Jarzynski equality continues to hold in this case, with the free energy replaced by the grand potential. The work distribution displays a marked dependence on the temperature. While in the classical (high-temperature) regime, the work probability density function follows a Gaussian distribution and the free energy difference entering the Jarzynski equality is null, the free energy difference is finite in the quantum regime, and the work probability distribution function becomes multimodal. We point out the dependence of the work statistics on the number of electrons composing the system.

Quantum Bochkov-Kuzovlev work fluctuation theorems

The quantum version of the Bochkov-Kuzovlev identity is derived on the basis of the appropriate definition of work as the difference of the measured internal energies of a quantum system at the beginning and the end of an external action on the system given by a prescribed protocol. According to the spirit of the original Bochkov-Kuzovlev approach, we adopt the 'exclusive' viewpoint, meaning that the coupling to the external work source is not counted as part of the internal energy. The corresponding canonical and microcanonical quantum fluctuation theorems are derived as well, and are compared with the respective theorems obtained within the 'inclusive' approach. The relations between the quantum inclusive work w, the exclusive work w(0) and the dissipated work w(dis), are discussed and clarified. We show by an explicit example that w(0) and w(dis) are distinct stochastic quantities obeying different statistics.

Work distributions for Ising chains in a time-dependent magnetic field

Master equations can be conveniently used to investigate many particle systems driven out of equilibrium by time-dependent external fields. This topic is of vital interest in connection with fluctuation theorems and the associated microscopic work and heat distributions. We present an exact Monte Carlo simulation algorithm, which allows us to study interaction effects on these distributions. The method is applied to an Ising chain with Glauber dynamics. We find that the distributions are characterized by delta and step functions with a smooth part in between. With decreasing sweeping rate of the external field or decreasing interaction strength, the singular part becomes less dominant and the Gaussian fluctuation regime is approached.

Increase of Boltzmann entropy in a quantum forced harmonic oscillator

Recently, a quantum-mechanical proof of the increase of Boltzmann entropy in quantum systems that are coupled to an external classical source of work has been given. Here we illustrate this result by applying it to a forced quantum harmonic oscillator. We show plots of the actual temporal evolution of work and entropy for various forcing protocols. We note that entropy and work can be partially or even fully returned to the source, although both work and entropy balances are non-negative at all times in accordance with the minimal work principle and the Clausius principle, respectively. A necessary condition for the increase of entropy is that the initial distribution is decreasing (e.g., canonical). We show evidence that for a nondecreasing distribution (e.g., microcanonical), the quantum expectation of entropy may decrease slightly. Interestingly, the classical expectation of entropy cannot decrease, irrespective of the initial distribution, in the forced harmonic oscillator.