• Abstract and Figures
• Public Full-text

Toward conjecture: Heating is faster than cooling Tan Van Vu1, 2, ∗ and Yoshihiko Hasegawa1, y 1Department of Information and Communication Engineering, Graduate School of Information Science and Technology, The University of Tokyo, Tokyo 113-8656, Japan 2Current Address: Department of Physics, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan (Dated: July 27, 2021) An asymmetry in thermal relaxation toward equilibrium has been uncovered for Langevin systems near stable minima [Phys. Rev. Lett. 125, 110602 (2020)]. It is conjectured that, given the same degree of nonequilibrium of the initial distributions, relaxation from a lower temperature state (heating) is faster than that from a higher temperature state (cooling). In this study, we elucidate the conjecture for discrete-state Markovian systems described by the master equation. We rigorously prove that heating is faster than cooling for arbitrary two-state systems, whereas for systems with more than two distinct energy levels, the conjecture is, in general, invalid. Furthermore, for systems whose energy levels degenerate into two energy states, we find that there exist critical thresholds of the energy gap. Depending on the magnitude of this energy gap, heating can be faster or slower than cooling, irrespective of the transition rates between states. Our results clarify the conjecture for discrete-state systems and reveal several hidden features inherent in thermal relaxation. I. INTRODUCTION via the convergence rate of the system state toward the equilibrium state [18]. Given two identical systems initi- Systems attached to thermal reservoirs will relax to- ated in thermal states, one at a lower and the other at a ward a stationary state. Such thermal relaxation pro- higher temperature than the given temperature, a natu- cesses are ubiquitous in nature and possess rich prop- ral question arises: which one relaxes faster? Recently, erties from both dynamic and thermodynamic perspec- this question has been addressed by Lapolla and Godec tives. One of the counterintuitive behaviors is the [19] for continuous-state Langevin systems. By consider- Mpemba effect [1], where cooling a hot system is faster ing a pair of thermodynamically equidistant temperature than cooling a cold system. Such nonmonotonic relax- quenches (which have the same nonequilibrium free en- ation phenomena have been observed in various systems ergy difference), they unveiled an unforeseen asymmetry [2{6] and theoretically analyzed for microscopic dynamics in thermal relaxation; i.e., relaxation from a lower tem- [7{9]. In addition, it was found that cooling a system be- perature is faster than that from a higher temperature. fore heating it could lead to exponentially fast relaxation Roughly speaking, it implies that heating up cold ob- [10]. From the perspective of thermodynamics, thermal jects is faster than cooling down hot objects. This phe- relaxation processes exhibit universal relations regard- nomenon has been proven for quenches of dynamics near ing irreversibility, which is quantified by irreversible en- stable minima, whereas the case of general overdamped tropy production [11]. Notably, it has been shown that dynamics remains a conjecture. Considering the average irreversible entropy production during thermal relaxation energy, it can be observed that the energy of a thermal is lower-bounded by information-theoretical [12{14] and state at a higher temperature is more than that at a lower geometrical [15] distances between the initial and final temperature. This allows us to say that, from the ener- states in both classical and quantum regimes. These re- getic perspective, uphill relaxation is faster than downhill lations imply stronger inequalities than the conventional relaxation, which is counterintuitive to an extent. More- second law of thermodynamics and impose geometrical over, relaxation speed cannot be characterized solely by constraints on the possible relaxation path. Since ther- thermodynamic quantities such as dissipation or frenesy mal relaxation is important in condensed matter [16] and [20]. Therefore, it is highly nontrivial that free energy heat engines [17], deepening our understanding of ther- plays an essential role as a quantifier of nonequilibrium mal relaxation would benefit research in these areas. degree in equidistant temperature quenches. Consider preparing a thermal state corresponding to a In this study, we elucidate the conjecture for discrete- arXiv:2102.07429v2 [cond-mat.stat-mech] 24 Jul 2021 given temperature via thermal relaxation; i.e., attaching state systems modeled by Markov jump processes, thus a system to a single reservoir and allowing it to relax to- improving our understanding of thermal relaxation. ward equilibrium. In this setting, the relaxation time is a First, we prove that for arbitrary two-state systems, heat- quantity of interest and can be approximately estimated ing is faster than cooling, which partially validates the conjecture. However, we find that it is not the case for general systems with at least three distinct energy lev- els. By analytically constructing counterexamples, we ∗ [email protected] demonstrate that heating can be faster or slower than y [email protected] cooling, depending on the transition rates, which implies 2 that the conjecture does not universally hold for discrete- time-evolution distribution corresponding to the initial state systems. Nevertheless, restricted to a particular state jπii, then heating is said to be faster (slower) than c f h f class of systems, some universal results on relaxation cooling if D(pt kπ ) < (>)D(pt kπ ) in the long time. In asymmetry are obtained. We show that when the en- what follows, we explain in detail how to determine the ergy levels of the system are two-state degenerate to two relaxation speed in the long-time regime. energy states, there exist two critical energy gap thresh- Let 0 = λ1 > λ2 > ··· > λN be the eigenvalues of the olds. Depending on whether an energy gap is larger or transition matrix and fjvnign be the set of corresponding smaller than these thresholds, it can be concluded with eigenvectors, certainty that heating is faster or slower than cooling. Rjvni = λnjvni: (2) Notably, all eigenvalues are real numbers since the ma- II. SETUP trix R satisfies the detailed balance condition [22]. The N eigenvectors fjvnig form a basis for the space R with f We consider the thermal relaxation process of an open jv1i = jπ i, and jvni is a traceless vector for all n ≥ 2 system with N states. The system is coupled to a ther- (i.e., h1jvni = 0 since h1jR = h0j). Here, j0i (j1i) denotes −1 mal reservoir at the inverse temperature βf = (kBTf ) , the N-dimensional vector with all zero (one) elements. i where kB is the Boltzmann constant. Owing to inter- Therefore, the initial distribution jπ i can be expressed action with the thermal reservoir, stochastic transitions as a linear combination of fjvnig as follows: between states are induced. The dynamics of the system N are governed by the master equation, i f X i jπ i = jπ i + γnjvni; (3) n=2 @tjpti = Rjpti; (1) where γi 's are real numbers. Consequently, the proba- > n where jpti := [p1(t); : : : ; pN (t)] denotes the probabil- bility distribution at time t can be analytically written ity distribution of the system at time t; the matrix R = in the following form: N×N [Rmn] 2 R is time-independent with Rmn ≥ 0 de- noting the transition rate from state n to state m (6= n), N i f X i λnt P jpti = jπ i + γne jvni: (4) and m Rmn = 0. Without loss of generality, we assume n=2 that E1 ≥ · · · ≥ EN , where En denotes the energy of state n. The transition rates satisfy the detailed balance i In the long-time limit, the probability distribution jpti −β En −β Em condition, Rmne f = Rnme f , which is the suf- can be approximated up to the second-order term as ficient condition such that the system always relaxes to f i f i λ2t the thermal Gibbs state jπ i after a sufficiently long time, jp (t)i ' jπ i + γ2e jv2i: (5) irrespective of the initial state. Here, πf = e−βf En =Z n βf Thus, the relaxation speed can be quantified via the value PN −βf En and Zβ = e . The transition rates can be i f n=1 of jγ2j [7]. Accordingly, heating is faster (slower) than −βf (Bmn−En) c h expressed as Rmn = Γe for m 6= n, where cooling if jγ2j < (>)jγ2 j (see AppendixA for proof). i Bmn = Bnm are the barrier coefficients, and Γ is a posi- A closed form of γ2 can be obtained analytically [8]. tive constant. The transition rate matrix R can be transformed to a Now, let us formulate the problem. We consider re- symmetric matrix R as follows: laxation that initiates from a thermal state jπii associ- 1=2 −1=2 −1 R = F RF ; (6) ated with the inverse temperature βi = (kBTi) . This can be regarded as a temperature quench Ti ! Tf at N×N βf En − where F = [Fmn] 2 R with Fmn = e δmn. time t = 0 . Given a pair of cold and hot tempera- Notably, the matrix R has the same eigenvalues as R, tures, Tc and Th, satisfying Tc < Tf < Th, we investigate and its eigenvectors fjfnig are related to those of R as the relaxation speed depending on the quench direction, 1=2 jfni = F jvni. Moreover, these eigenvectors are mu- Ti = Tc " Tf (heating) and Ti = Th # Tf (cooling). tually orthogonal, hfmjfni = hvnjFjvmiδmn. Multiplying The degree of nonequilibrium (or free energy) of each hf jF1=2 on both sides of Eq. (3), we can show that γi initial state is the same, i.e., D(πckπf ) = D(πhkπf ), 2 2 P is proportional to the inner product between the initial where D(pkq) := n pn ln(pn=qn) is the relative entropy distribution jπii and the vector jf i, given by between two distributions jpi and jqi.

Load more

  6

Copy Link