Suppression of work fluctuations by optimal control: An approach based on Jarzynski's equality

Experimental Test of the Jarzynski Equality in a Single Spin-1 System Using High-Fidelity Single-Shot Readouts

The Jarzynski equality (JE), which connects the equilibrium free energy with nonequilibrium work statistics, plays a crucial role in quantum thermodynamics. Although practical quantum systems are usually multilevel systems, most tests of the JE were executed in two-level systems. A rigorous test of the JE by directly measuring the work distribution of a physical process in a high-dimensional quantum system remains elusive. Here, we report an experimental test of the JE in a single spin-1 system. We realized nondemolition projective measurement of this three-level system via cascading high-fidelity single-shot readouts and directly measured the work distribution utilizing the two-point measurement protocol. The validity of the JE was verified from the nonadiabatic to adiabatic zone and under different effective temperatures. Our work puts the JE on a solid experimental foundation and makes the nitrogen-vacancy (NV) center system a mature toolbox to perform advanced experiments of stochastic quantum thermodynamics.

Calculating the free energy difference by applying the Jarzynski equality to a virtual integrable system

The Jarzynski equality (JE) provides a nonequilibrium method to measure and calculate the free energy difference (FED). Note that if two systems share the same Hamiltonian at two equilibrium states, respectively, they share the same FED between these two equilibrium states as well. Therefore the calculation of the FED of a system may be facilitated by considering instead another virtual system designed to this end. Taking advantage of this flexibility and the JE, we show that by introducing an integrable virtual system, the evolution problem involved in the JE can be solved. As a consequence, FED is expressed in the form of an equilibrium equality, in contrast with the nonequilibrium JE it is based on. Numerically, this result allows FED to be computed by sampling the canonical ensemble directly and the computational cost can be significantly reduced. The effectiveness and efficiency of this scheme are illustrated with numerical studies of several representative model systems.

Phase Transition in Protocols Minimizing Work Fluctuations

For two canonical examples of driven mesoscopic systems-a harmonically trapped Brownian particle and a quantum dot-we numerically determine the finite-time protocols that optimize the compromise between the standard deviation and the mean of the dissipated work. In the case of the oscillator, we observe a collection of protocols that smoothly trade off between average work and its fluctuations. However, for the quantum dot, we find that as we shift the weight of our optimization objective from average work to work standard deviation, there is an analog of a first-order phase transition in protocol space: two distinct protocols exchange global optimality with mixed protocols akin to phase coexistence. As a result, the two types of protocols possess qualitatively different properties and remain distinct even in the infinite duration limit: optimal-work-fluctuation protocols never coalesce with the minimal-work protocols, which therefore never become quasistatic.

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.

Universal Work Fluctuations During Shortcuts to Adiabaticity by Counterdiabatic Driving

Counterdiabatic driving (CD) exploits auxiliary control fields to tailor the nonequilibrium dynamics of a quantum system, making possible the suppression of dissipated work in finite-time thermodynamics and the engineering of optimal thermal machines with no friction. We show that while the mean work done by the auxiliary controls vanishes, CD leads to a broadening of the work distribution. We derive a fundamental inequality that relates nonequilibrium work fluctuations to the operation time and quantifies the thermodynamic cost of CD in both critical and noncritical systems.

Trade-Off Between Speed and Cost in Shortcuts to Adiabaticity

Achieving effectively adiabatic dynamics is a ubiquitous goal in almost all areas of quantum physics. Here, we study the speed with which a quantum system can be driven when employing transitionless quantum driving. As a main result, we establish a rigorous link between this speed, the quantum speed limit, and the (energetic) cost of implementing such a shortcut to adiabaticity. Interestingly, this link elucidates a trade-off between speed and cost, namely, that instantaneous manipulation is impossible as it requires an infinite cost. These findings are illustrated for two experimentally relevant systems-the parametric oscillator and the Landau-Zener model-which reveal that the spectral gap governs the quantum speed limit as well as the cost for realizing the shortcut.

Merits and qualms of work fluctuations in classical fluctuation theorems

Work is one of the most basic notions in statistical mechanics, with work fluctuation theorems being one central topic in nanoscale thermodynamics. With Hamiltonian chaos commonly thought to provide a foundation for classical statistical mechanics, here we present general salient results regarding how (classical) Hamiltonian chaos generically impacts on nonequilibrium work fluctuations. For isolated chaotic systems prepared with a microcanonical distribution, work fluctuations are minimized and vanish altogether in adiabatic work protocols. For isolated chaotic systems prepared at an initial canonical distribution at inverse temperature β, work fluctuations depicted by the variance of e^{-βW} are also minimized by adiabatic work protocols. This general result indicates that, if the variance of e^{-βW} diverges for an adiabatic work protocol, it diverges for all nonadiabatic work protocols sharing the same initial and final Hamiltonians. Such divergence is hence not an isolated event and thus greatly impacts on the efficiency of using Jarzynski's equality to simulate free-energy differences. Theoretical results are illustrated in a Sinai model. Our general insights shall boost studies in nanoscale thermodynamics and are of fundamental importance in designing useful work protocols.

Work fluctuation and total entropy production in nonequilibrium processes

Work fluctuation and total entropy production play crucial roles in small thermodynamic systems subject to large thermal fluctuations. We investigate a trade-off relation between them in a nonequilibrium situation in which a system starts from an arbitrary nonequilibrium state. We apply a variational method to study this problem and find a stationary solution against variations over protocols that describe the time dependence of the Hamiltonian of the system. Using the stationary solution, we find the minimum of the total entropy production for a given amount of work fluctuation. An explicit protocol that achieves this is constructed from an adiabatic process followed by a quasistatic process. The obtained results suggest how one can control the nonequilibrium dynamics of the system while suppressing its work fluctuation and total entropy production.

Elastic Barrier Dynamical Freezing in Free Energy Calculations: A Way To Speed Up Nonequilibrium Molecular Dynamics Simulations by Orders of Magnitude

An important issue concerning computer simulations addressed to free energy estimates via nonequilibrium work theorems, such as the Jarzynski equality [Phys. Rev. Lett. 1997, 78, 2690], is the computational effort required to achieve results with acceptable accuracy. In this respect, the dynamical freezing approach [Phys. Rev. E 2009, 80, 041124] has been shown to improve the efficiency of this kind of simulations, by blocking the dynamics of particles located outside an established mobility region. In this report, we show that dynamical freezing produces a systematic spurious decrease of the particle density inside the mobility region. As a consequence, the requirements to apply nonequilibrium work theorems are only approximately met. Starting from these considerations, we have developed a simulation scheme, called "elastic barrier dynamical freezing", according to which a stiff potential-energy barrier is enforced at the boundaries of the mobility region, preventing the particles from leaving this region of space during the nonequilibrium trajectories. The method, tested on the calculation of the distance-dependent free energy of a dimer immersed into a Lennard-Jones fluid, provides an accuracy comparable to the conventional steered molecular dynamics, with a computational speedup exceeding a few orders of magnitude.

Shortcuts to adiabaticity from linear response theory

A shortcut to adiabaticity is a finite-time process that produces the same final state as would result from infinitely slow driving. We show that such shortcuts can be found for weak perturbations from linear response theory. With the help of phenomenological response functions, a simple expression for the excess work is found-quantifying the nonequilibrium excitations. For two specific examples, i.e., the quantum parametric oscillator and the spin 1/2 in a time-dependent magnetic field, we show that finite-time zeros of the excess work indicate the existence of shortcuts. Finally, we propose a degenerate family of protocols, which facilitates shortcuts to adiabaticity for specific and very short driving times.

Principle of minimal work fluctuations

Understanding and manipulating work fluctuations in microscale and nanoscale systems are of both fundamental and practical interest. For example, in considering the Jarzynski equality 〈e-βW〉=e-βΔF, a change in the fluctuations of e-βW may impact how rapidly the statistical average of e-βW converges towards the theoretical value e-βΔF, where W is the work, β is the inverse temperature, and ΔF is the free energy difference between two equilibrium states. Motivated by our previous study aiming at the suppression of work fluctuations, here we obtain a principle of minimal work fluctuations. In brief, adiabatic processes as treated in quantum and classical adiabatic theorems yield the minimal fluctuations in e-βW. In the quantum domain, if a system initially prepared at thermal equilibrium is subjected to a work protocol but isolated from a bath during the time evolution, then a quantum adiabatic process without energy level crossing (or an assisted adiabatic process reaching the same final states as in a conventional adiabatic process) yields the minimal fluctuations in e-βW, where W is the quantum work defined by two energy measurements at the beginning and at the end of the process. In the classical domain where the classical work protocol is realizable by an adiabatic process, then the classical adiabatic process also yields the minimal fluctuations in e-βW. Numerical experiments based on a Landau-Zener process confirm our theory in the quantum domain, and our theory in the classical domain explains our previous numerical findings regarding the suppression of classical work fluctuations [G. Y. Xiao and J. B. Gong, Phys. Rev. E 90, 052132 (2014)].

Construction and optimization of a quantum analog of the Carnot cycle

The quantum analog of Carnot cycles in few-particle systems consists of two quantum adiabatic steps and two isothermal steps. This construction is formally justified by use of a minimum work principle. It is then shown, using minimal assumptions of work or heat in nanoscale systems, that the heat-to-work efficiency of such quantum heat engine cycles can be further optimized via two conditions regarding the expectation value of some generalized force operators evaluated at equilibrium states. In general the optimized efficiency is system specific, lower than the Carnot efficiency, and dependent upon both temperatures of the cold and hot reservoirs. Simple computational examples are used to illustrate our theory. The results should be an important guide towards the design of favorable working conditions of a realistic quantum heat engine.

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.