Work fluctuations in a time-dependent harmonic potential: rigorous results beyond the overdamped limit

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

Heat statistics in the relaxation process of the Edwards-Wilkinson elastic manifold

The stochastic thermodynamics of systems with a few degrees of freedom has been studied extensively so far. We would like to extend the study to systems with more degrees of freedom and even further-continuous fields with infinite degrees of freedom. The simplest case for a continuous stochastic field is the Edwards-Wilkinson elastic manifold. It is an exactly solvable model of which the heat statistics in the relaxation process can be calculated analytically. The cumulants require a cutoff spacing to avoid ultraviolet divergence. The scaling behavior of the heat cumulants with time and the system size as well as the large deviation rate function of the heat statistics in the large size limit is obtained.

Hierarchical structure of fluctuation theorems for a driven system in contact with multiple heat reservoirs

For driven open systems in contact with multiple heat reservoirs, we find the marginal distributions of work or heat do not satisfy any fluctuation theorem, but only the joint distribution of work and heat satisfies a family of fluctuation theorems. A hierarchical structure of these fluctuation theorems is discovered from microreversibility of the dynamics by adopting a step-by-step coarse-graining procedure in both classical and quantum regimes. Thus, we put all fluctuation theorems concerning work and heat into a unified framework. We also propose a general method to calculate the joint statistics of work and heat in the situation of multiple heat reservoirs via the Feynman-Kac equation. For a classical Brownian particle in contact with multiple heat reservoirs, we verify the validity of the fluctuation theorems for the joint distribution of work and heat.

Microscopic theory of the Curzon-Ahlborn heat engine based on a Brownian particle

The Curzon-Ahlborn (CA) efficiency, as the efficiency at the maximum power (EMP) of the endoreversible Carnot engine, has significant impact on finite-time thermodynamics. However, the CA engine is based on many assumptions. In the past few decades, although a lot of efforts have been made, a microscopic theory of the CA engine is still lacking. By adopting the method of the stochastic differential equation of energy, we formulate a microscopic theory of the CA engine realized with a highly underdamped Brownian particle in a class of nonharmonic potentials. This theory gives microscopic interpretation of all assumptions made by Curzon and Ahlborn. In other words, we find a microscopic counterpart of the CA engine in stochastic thermodynamics. Also, based on this theory, we derive the explicit expression of the protocol associated with the maximum power for any given efficiency, and we obtain analytical results of the power and the efficiency statistics for the Brownian CA engine. Our research brings new perspectives to experimental studies of finite-time microscopic heat engines featured with fluctuations.

Perturbative expansion of irreversible work in symmetric and asymmetric processes

The systematic expansion method of the solution of the Fokker-Planck equation is developed by generalizing the formulation proposed in [J. Phys. A: Math. Theor. 50, 325001 (2017)10.1088/1751-8121/aa7af4]. Using this method, we obtain an alternative formula to calculate the mean work perturbatively which is applicable to systems with degeneracy in the eigenvalues of the Fokker-Planck operator. This method enables us to study how the geometrical symmetry affects thermodynamic description of a Brownian particle. To illustrate the application of the derived theory, we consider the Fokker-Planck equation with a two-dimensional harmonic potential. To investigate the effect of symmetry of the potential, we study thermodynamic properties in symmetric and asymmetric deformation processes of the potential: the rotational symmetry of the harmonic potential is held in the former, but it is broken in the latter. Optimized deformations in these processes are defined by minimizing mean work. Comparing these optimized processes, we find that the difference between the symmetric and asymmetric processes is maximized when the deformation time of the potential is given by a critical time which is characterized by the relaxation time of the Fokker-Planck equation. This critical time in the mean work is smaller than that of the change of the mean energy because of the hysteresis effect in the irreversible processes.

Solution to the Fokker-Planck equation for slowly driven Brownian motion: Emergent geometry and a formula for the corresponding thermodynamic metric

Considerable progress has recently been made with geometrical approaches to understanding and controlling small out-of-equilibrium systems, but a mathematically rigorous foundation for these methods has been lacking. Towards this end, we develop a perturbative solution to the Fokker-Planck equation for one-dimensional driven Brownian motion in the overdamped limit enabled by the spectral properties of the corresponding single-particle Schrödinger operator. The perturbation theory is in powers of the inverse characteristic timescale of variation of the fastest varying control parameter, measured in units of the system timescale, which is set by the smallest eigenvalue of the corresponding Schrödinger operator. It applies to any Brownian system for which the Schrödinger operator has a confining potential. We use the theory to rigorously derive an exact formula for a Riemannian "thermodynamic" metric in the space of control parameters of the system. We show that up to second-order terms in the perturbation theory, optimal dissipation-minimizing driving protocols minimize the length defined by this metric. We also show that a previously proposed metric is calculable from our exact formula with corrections that are exponentially suppressed in a characteristic length scale. We illustrate our formula using the two-dimensional example of a harmonic oscillator with time-dependent spring constant in a time-dependent electric field. Lastly, we demonstrate that the Riemannian geometric structure of the optimal control problem is emergent; it derives from the form of the perturbative expansion for the probability density and persists to all orders of the expansion.

Experimental Verification of the Work Fluctuation-Dissipation Relation for Information-to-Work Conversion

We study experimentally work fluctuations in a Szilard engine that extracts work from information encoded as the occupancy of an electron level in a semiconductor quantum dot. We show that as the average work extracted per bit of information increases toward the Landauer limit k_{B}Tln2, the work fluctuations decrease in accordance with the work fluctuation-dissipation relation. We compare the results to a protocol without measurement and feedback and show that when no information is used, the work output and fluctuations vanish simultaneously, contrasting the information-to-energy conversion case where increasing amount of work is produced with decreasing fluctuations. Our study highlights the importance of fluctuations in the design of information-to-work conversion processes.

Stochastic Thermodynamics of a Piezoelectric Energy Harvester Model

We experimentally study a piezoelectric energy harvester driven by broadband random vibrations. We show that a linear model, consisting of an underdamped Langevin equation for the dynamics of the tip mass, electromechanically coupled with a capacitor and a load resistor, can accurately describe the experimental data. In particular, the theoretical model allows us to define fluctuating currents and to study the stochastic thermodynamics of the system, with focus on the distribution of the extracted work over different time intervals. Our analytical and numerical analysis of the linear model is succesfully compared to the experiments.

Work distribution in thermal processes

We find the moment generating function (mgf) of the nonequilibrium work for open systems undergoing a thermal process, i.e., when the stochastic dynamics maps thermal states into time-dependent thermal states. The mgf is given in terms of a temperaturelike scalar satisfying a first-order ordinary differential equation. We apply the result to some paradigmatic situations: a levitated nanoparticle in a breathing optical trap, a Brownian particle in a box with a moving piston, and a two-state system driven by an external field, where the work mgfs are obtained for different timescales and compared with Monte Carlo simulations.

Standard and fractional Ornstein-Uhlenbeck process on a growing domain

We study normal diffusive and subdiffusive processes in a harmonic potential (Ornstein-Uhlenbeck process) on a uniformly growing or contracting domain. Our starting point is a recently derived fractional Fokker-Planck equation, which covers both the case of Brownian diffusion and the case of a subdiffusive continuous-time random walk (CTRW). We find a high sensitivity of the random walk properties to the details of the domain growth rate, which gives rise to a variety of regimes with extremely different behaviors. At the origin of this rich phenomenology is the fact that the walkers still move while they wait to jump, since they are dragged by the deterministic drift arising from the domain growth. Thus, the increasingly long waiting times associated with the aging of the subdiffusive CTRW imply that, in the time interval between two consecutive jumps, the walkers might travel over much longer distances than in the normal diffusive case. This gives rise to seemingly counterintuitive effects. For example, on a static domain, both Brownian diffusion and subdiffusive CTRWs yield a stationary particle distribution with finite width when a harmonic potential is at play, thus indicating a confinement of the diffusing particle. However, for a sufficiently fast growing or contracting domain, this qualitative behavior breaks down, and differences between the Brownian case and the subdiffusive case are found. In the case of Brownian particles, a sufficiently fast exponential domain growth is needed to break the confinement induced by the harmonic force; in contrast, for subdiffusive particles such a breakdown may already take place for a sufficiently fast power-law domain growth. Our analytic and numerical results for both types of diffusion are fully confirmed by random walk simulations.

Brownian motion of a circle swimmer in a harmonic trap

We study the dynamics of a Brownian circle swimmer with a time-dependent self-propulsion velocity in an external temporally varying harmonic potential. For several situations, the noise-free swimming paths, the noise-averaged mean trajectories, and the mean-square displacements are calculated analytically or by computer simulation. Based on our results, we discuss optimal swimming strategies in order to explore a maximum spatial range around the trap center. In particular, we find a resonance situation for the maximum escape distance as a function of the various frequencies in the system. Moreover, the influence of the Brownian noise is analyzed by comparing noise-free trajectories at zero temperature with the corresponding noise-averaged trajectories at finite temperature. The latter reveal various complex self-similar spiral or rosette-like patterns. Our predictions can be tested in experiments on artificial and biological microswimmers under dynamical external confinement.

Stochastic Thermodynamics of a Particle in a Box

The piston system (particles in a box) is the simplest paradigmatic model in traditional thermodynamics. However, the recently established framework of stochastic thermodynamics (ST) fails to apply to this model system due to the embedded singularity in the potential. In this Letter, we study the ST of a particle in a box by adopting a novel coordinate transformation technique. Through comparing with the exact solution of a breathing harmonic oscillator, we obtain analytical results of work distribution for an arbitrary protocol in the linear response regime and verify various predictions of the fluctuation-dissipation relation. When applying to the Brownian Szilard engine model, we obtain the optimal protocol λ_{t}=λ_{0}2^{t/τ} for a given sufficiently long total time τ. Our study not only establishes a paradigm for studying ST of a particle in a box but also bridges the long-standing gap in the development of ST.

Non-Markovian barotropic-type and Hall-type fluctuation relations in crossed electric and magnetic fields

In this paper we derive the non-Markovian barotropic-type and Hall-type fluctuation relations for noninteracting charged Brownian particles embedded in a memory heat bath and under the action of crossed electric and magnetic fields. We first obtain a more general non-Markovian fluctuation relation formulated within the context of a generalized Langevin equation with arbitrary friction memory kernel and under the action of a constant magnetic field and an arbitrary time-dependent electric field. It is shown that this fluctuation relation is related to the total amount of an effective work done on the charged particle as it is driven out of equilibrium by the applied time-dependent electric field. Both non-Markovian barotropic- and Hall-type fluctuation relations are then derived when the electric field is assumed to be also a constant vector pointing along just one axis. In the Markovian limit, we show explicitly that they reduce to the same results reported in the literature.

Non-Markovian work fluctuation theorem in crossed electric and magnetic fields

The validity of the transient work fluctuation theorem for a charged Brownian harmonic oscillator embedded in a non-Markovian heat bath and under the action of crossed electric and magnetic fields is investigated. The aforementioned theorem is verified to be valid within the context of the generalized Langevin equation with an arbitrary memory kernel and arbitrary dragging in the potential minimum. The fluctuation-dissipation relation of the second kind is assumed to be valid and shows that the non-Markovian stochastic dynamics associated with the particle, in the absence of the external time-dependent electric field, reaches an equilibrium state, as is precisely demanded by such a relation. The Jarzynski equality in this problem is also analyzed.

Hidden entropy production by fast variables

We investigate nonequilibrium underdamped Langevin dynamics of Brownian particles that interact through a harmonic potential with coupling constant K and are in thermal contact with two heat baths at different temperatures. The system is characterized by a net heat flow and an entropy production in the steady state. We compare the entropy production of the harmonic system with that of Brownian particles linked with a rigid rod. The harmonic system may be expected to reduce to the rigid rod system in the infinite K limit. However, we find that the harmonic system in the K→∞ limit produces more entropy than the rigid rod system. The harmonic system has the center-of-mass coordinate as a slow variable and the relative coordinate as a fast variable. By identifying the contributions of the degrees of freedom to the total entropy production, we show that the hidden entropy production by the fast variable is responsible for the extra entropy production. We discuss the K dependence of each contribution.

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.

Nonequilibrium fluctuations for a single-particle analog of gas in a soft wall

We investigate the motion of a colloidal particle driven out of equilibrium by a time-varying stiffness of the optical trap that produces persistent nonequilibrium work. Measurements of work production for repeated cycles composed of the compression and expansion processes for the optical potential show huge fluctuations due to thermal motion. Using a precise technique to modulate the stiffness in time, we accurately estimate the probability distributions of work produced for the compression and expansion processes. We confirm the fluctuation theorem from the ratio of the two distributions. We also show that the average values of work for the two processes comply with the Jarzynski equality. This system has an analogy with a gas in a breathing soft wall. We discuss about its applicability to a heat engine and an information engine operated by feedback control.

Work fluctuations for a Brownian particle driven by a correlated external random force

We have considered the underdamped motion of a Brownian particle in the presence of a correlated external random force. The force is modeled by an Ornstein-Uhlenbeck process. We investigate the fluctuations of the work done by the external force on the Brownian particle in a given time interval in the steady state. We calculate the large deviation functions as well as the complete asymptotic form of the probability density function of the performed work. We also discuss the symmetry properties of the large deviation functions for this system. Finally we perform numerical simulations and they are in a very good agreement with the analytic results.

Heat fluctuations and initial ensembles

Time-integrated quantities such as work and heat increase incessantly in time during nonequilibrium processes near steady states. In the long-time limit, the average values of work and heat become asymptotically equivalent to each other, since they only differ by a finite energy change in average. However, the fluctuation theorem (FT) for the heat is found not to hold with the equilibrium initial ensemble, while the FT for the work holds. This reveals an intriguing effect of everlasting initial memory stored in rare events. We revisit the problem of a Brownian particle in a harmonic potential dragged with a constant velocity, which is in contact with a thermal reservoir. The heat and work fluctuations are investigated with initial Boltzmann ensembles at temperatures generally different from the reservoir temperature. We find that, in the infinite-time limit, the FT for the work is fully recovered for arbitrary initial temperatures, while the heat fluctuations significantly deviate from the FT characteristics except for the infinite initial-temperature limit (a uniform initial ensemble). Furthermore, we succeed in calculating finite-time corrections to the heat and work distributions analytically, using the modified saddle point integral method recently developed by us. Interestingly, we find noncommutativity between the infinite-time limit and the infinite-initial-temperature limit for the probability distribution function (PDF) of the heat.