Next Article in Journal
Quantum Stream Cipher Based on Holevo–Yuen Theory
Next Article in Special Issue
Critical Quantum Metrology in the Non-Linear Quantum Rabi Model
Previous Article in Journal
Evolution of Brains and Computers: The Roads Not Taken
Previous Article in Special Issue
The Classical–Quantum Passage: A van der Waals Description
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kibble–Zurek Scaling from Linear Response Theory

by
Pierre Nazé
1,
Marcus V. S. Bonança
1 and
Sebastian Deffner
1,2,*
1
Instituto de Física ‘Gleb Wataghin’, Universidade Estadual de Campinas, Campinas 13083-859, Brazil
2
Department of Physics, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(5), 666; https://doi.org/10.3390/e24050666
Submission received: 23 March 2022 / Revised: 27 April 2022 / Accepted: 7 May 2022 / Published: 10 May 2022
(This article belongs to the Special Issue Current Trends in Quantum Phase Transitions)

Abstract

:
While quantum phase transitions share many characteristics with thermodynamic phase transitions, they are also markedly different as they occur at zero temperature. Hence, it is not immediately clear whether tools and frameworks that capture the properties of thermodynamic phase transitions also apply in the quantum case. Concerning the crossing of thermodynamic critical points and describing its non-equilibrium dynamics, the Kibble–Zurek mechanism and linear response theory have been demonstrated to be among the very successful approaches. In the present work, we show that these two approaches are also consistent in the description of quantum phase transitions, and that linear response theory can even inform arguments of the Kibble–Zurek mechanism. In particular, we show that the relaxation time provided by linear response theory gives a rigorous argument for why to identify the “gap” as a relaxation rate, and we verify that the excess work computed from linear response theory exhibits Kibble–Zurek scaling.

1. Introduction

In thermodynamics a phase transition describes the dramatic change of the macroscopically observable physical properties of matter [1]. At the microscopic scale, such a transition requires the fundamental re-ordering and structuring (or lack thereof) of the system’s constituents. Realizing the complexity of the microscopic properties of a system approaching and passing through a phase transition, it is almost obvious to recognize that around the transition the response to external perturbations is strongly inhibited. In renormalization group theory, this insight is formalized as the universal behavior of response functions [2].
All real processes occur at a finite time and are accompanied by the inevitable production of nonequilibrium excitations. If the rate of driving is much slower than the inverse of the relaxation time, effectively quasistatic, equilibrium processes can be facilitated. However, close to critical points the relaxation time diverges (as does the response), and hence any real driving through a phase transition will always exhibit nonequilibrium characteristics. This observation is at the core of the Kibble–Zurek mechanism [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17], which predicts the size of finite domains to be fully determined by the critical exponents and the rate of driving.
Whereas the arguments of the Kibble–Zurek mechanism can be phrased rather intuitively for thermodynamic phase transitions, the situation is more involved for quantum phase transitions [18]. A quantum system undergoes a quantum phase transition, if its macroscopically observable physical properties of the ground state change according to an external field [18]. It has then been argued that the energy difference between the ground and excited state, the so-called “gap”, plays the role of a relaxation rate, and thus, the Kibble–Zurek mechanism can be generalized to the fully quantum domain [19,20,21,22].
Both, the classical and the quantum Kibble–Zurek mechanism describe nonequilibrium excitations in terms of the critical exponents of the underlying equilibrium phase transition. Hence, it appears somewhat natural to assume that the mechanism itself is valid “close enough to equilibrium”. However, like all phenomenological approaches the range of validity cannot be fully determined from within the approach. On the other hand, “close to equilibrium” is the domain of linear response theory [23,24]. Therefore, the natural question arises whether the Kibble–Zurek mechanics can be phrased as a consequence of linear response, or whether the mechanism goes beyond the theory. In previous work, we have found some clear evidence that the Kibble–Zurek mechanism does in fact describe the physics outside the range of validity of linear response [25,26], but also that for slow enough driving, the two approaches are consistent [27].
In the present work, we further investigate to what extent insight from and about the Kibble–Zurek mechanism can be extracted from linear response theory. To this end, we focus on the less intuitive case and analyze the quantum phase transition of the Ising model in the transverse field [28]. Since this model can be solved analytically [29,30], it has become the paradigmatic case study for phase transitions in quantum systems [19,31,32,33,34,35,36]. As a first result, we elucidate the interpretation of the “gap” as a relaxation rate. To this end, we compute the relaxation time directly from the response function, and we find that the quantum phase transition indeed exhibits “critical slowing down”. This insight can then be used to compute the excess work, which quantifies the “amount” of diabatic excitations and which can be computed relatively easily by means of linear response theory [37,38,39,40,41,42,43,44,45,46]. We find that this excess work exhibits exactly the polynomial behavior as a function of the driving time predicted by the Kibble–Zurek mechanism. Finally, benchmarking our results from linear response theory against exact numerics, we obtain a good characterization of the range of validity of linear response theory around quantum phase transitions.

2. Preliminaries

We begin by establishing notions and notations. To this end, we briefly review some elements of the Kibble–Zurek mechanism, as well as how to compute the excess work from linear response theory. For specificity, we phrase our analysis entirely in terms of the quantum Ising chain in the transverse field,
H = J i = 1 N σ i x σ i + 1 x Γ i = 1 N σ i z .
where σ I z and σ i x are the Pauli matrices of the ith spin, J is the coupling energy, and Γ is the transverse magnetic field. For our purposes we choose N to be even, and we work with periodic boundary condition.

2.1. Kibble–Zurek Mechanism

The Kibble–Zurek mechanism is a phenomenological theory that can be used to describe the non-equilibrium dynamics of the Ising chain (1) when crossing its critical point, Γ = J . Renormalization group theory predicts [2,18] that the “relaxation time” diverges polynomially governed by the corresponding critical exponents. In quantum phase transitions the energy gap, Δ , plays the role of the relaxation rate [19], and we can write
τ R ( t ) ħ Δ ( t ) ,
where τ R is the relaxation time. For large systems, N 1 it is a simple exercise to show that
Δ ( t ) 2 | J Γ ( t ) | .
For simplicity and without loss of generality [47] we now assume that the magnetic field Γ changes linearly as a function of time,
Γ ( t ) = J 1 t τ ,
where τ is the duration of the process. The resulting τ R ( t ) is illustrated in Figure 1. Zurek recognized that this “critical freezing out” of the response has crucial ramifications for the nonequilibrium behavior [4]. Far from the critical point, the relaxation dynamics are fast and all nonequlibrium excitations can be mitigated or “healed”. Close to the critical point, this is no longer possible and the nonequilibrium shattering of the order parameter is imprinted onto the system. Thus, the regions far from the critical point are called adiabatic and close to the critical point the system undergoes the impulse regime.
The transition from adiabatic to impulse behavior occurs when the relaxation time becomes equal to the driving time t ^ = τ R ( t ^ ) , which can be solved for ± t ^ . We have
t ^ = ± ħ τ 2 J ,
which is governed by the driving rate 1 / τ , with which the system crosses the critical point.

2.2. Excess work in Linear Response Theory

In the following, we will investigate how much of the Kibble–Zurek mechanism is encoded in linear response theory. To this end, it will be instructive to write the Hamiltonian (1) as
H ( t ) = H 0 + A λ ( t ) ,
where A is some “observable” and | λ ( t ) | 1 . We will be particularly interested in the excess work W ex , i.e., the amount of energy above the ground state that is injected due to the finite time driving. In linear response theory W ex can be written as [25,26,27,37,38,39,40,41,42,45,46,48]
W ex = 1 2 0 τ 0 τ d t d t Ψ ( t t ) λ ˙ ( t ) λ ˙ ( t ) .
where Ψ ( t t ) is the relaxation function. See Ref. [44] for a brief review on W ex and linear response theory. This can be determined from the response function,
ϕ ( t ) = 1 i ħ [ A ( 0 ) , A ( t ) ] 0 ,
and ϕ ( t ) = Ψ ˙ ( t ) . The average is taken over an initial, equilibrium state, here over the ground state wave function, and A ( t ) is evolved according to the Heisenberg equation of motion for the unperturbed Hamiltonian H 0 .

2.3. Excess Work from Kibble–Zurek Arguments

In Ref. [27] it was argued that the behavior of W ex can be predicted with arguments from the Kibble–Zurek mechanism. To this end, it is instructive to recognize that only driving in the impulse regime will appreciably contribute to W ex , and hence the integrals in Equation (7) are evaluated up to t ^ and not τ . Note that strictly speaking, Ref. [27] verified the claim only for thermodynamic phase transitions, and more specifically noise-induced phase transitions. That similar arguments hold for quantum phase transitions is at best a sophisticated guess.
However, if one simply works with the expression of the relaxation function from renormalization group theory for the quantum Ising model, it is easy to show that [27]
W ex τ γ KZ , γ KZ = Λ 2 z ν + 1
where Λ is the critical exponent corresponding to the variation of an external parameter, z the dynamical critical exponent and ν the spatial critical exponent. In the present case, the driven Ising chain, we have Λ = 0 for the magnetic fields, z = 1 , and ν = 1 , and hence γ KZ = 1 , which is consistent with numerical findings [25,32]. However, the question remains whether this is a coincidence or a deep conceptual fact. Quantum phase transitions occur in the ground state and in unitary dynamics. Hence, notions such as “relaxation” are borrowed at best, and must not be taken too literally. Hence, a more thorough analysis of the relaxation function for the quantum Ising chain appears instrumental to elucidate how the Kibble–Zurek mechanism arises from the equilibrium properties of isolated quantum systems.

3. The Relaxation Function

We now need to analyze the relaxation function, Ψ ( t ) , more thoroughly and determine the corresponding relaxation time (within the framework of linear response theory). In Appendix A, we show that Ψ ( t ) for the quantum Ising chain (1) can be written as
Ψ ( t ) = 16 J n = 1 N / 2 J 3 ϵ n 3 sin 2 ( 2 n 1 ) π N cos 2 ϵ n t ħ ,
where we have introduced the eigenenergies
ϵ n = 2 J 2 + Γ 0 2 2 J Γ 0 cos ( 2 n 1 ) π N ,
and Γ 0 Γ ( t = 0 ) . Observe that Ψ ( t ) is a highly oscillatory function, which is expected for an isolated quantum system evolving under unitary dynamics. Moreover, the expression describing the relaxation behavior is governed by the initial value of the transverse magnetic field, which is a consequence of linear response theory. Thus, already at this point we recognize that the Kibble–Zurek mechanism goes beyond linear response theory, as its arguments address the simultaneous response of the system to the external driving. We will see shortly in Section 4.1 that this does not lead to a major complication within the range of validity of linear response theory.

3.1. Large N Limit

Phase transitions and their corresponding singularities are observed strictly only for infinitely large systems N 1 . In this limit, the discrete eigenvalue spectrum (11) becomes continuous and the quantum numbers can be expressed in terms of the wavenumber k = 2 π n / N . Thus, we write,
ψ ( t ) 8 J 2 π 0 π d k sin 2 ( k ) ϵ 3 ( k ) cos 2 ϵ ( k ) t ħ
and the eigenenergies (11) become
ϵ ( k ) = 2 J 2 + Γ 0 2 2 J Γ 0 cos ( k ) .
Note that the ground state n = 0 now corresponds to the zero mode, k = 0 .

Ferromagnetic and Paramagnetic Phases

It is instructive to first inspect the relaxation function far from the critical point. For Γ 0 / J 1 the quantum Ising model (1) assumes ferromagnetic ordering. In this case, the relaxation function (12) can be expanded and the leading order is,
ψ F ( t ) = 1 2 J cos 4 J t ħ .
Such a relaxation function is characteristic for single spins, which is a good description of macroscopic spin ordering. Moreover, observe that this ferromagnetic relaxation function is independent of the external field Γ 0
In the opposite limit, Γ 0 / J 1 , the Ising chain becomes paramagnetic. The corresponding expansion of Ψ ( t ) gives in leading order
ψ P ( t ) = J 2 2 Γ 0 3 cos 4 Γ 0 t ħ ,
which expresses the fact that paramagnetic systems are highly susceptible to external fields. The stark contrast in the response of the ferromagnetic and paramagnetic phases to external driving is indicative of the “dramatic” change that occurs at the phase transition.

Divergence at the Critical Point

It is then easy to see that Equation (12) exhibits a critical divergence if the Ising chain (1) is driven through its phase transition at Γ = J . To this end, we introduce the amplitude density A ( k ) as well as the characteristic frequency Ω ( k ) , with which we can write
ψ ( t ) = 0 π A ( k ) cos Ω ( k ) t .
Now assuming that the chain starts close to the critical point, Γ 0 J , we obtain
A ( k ) = sin 2 k 2 2 J π ( 1 cos k ) 3 / 2
and
Ω ( k ) = 4 2 ħ ( 1 cos k ) ,
for which Ψ ( t ) clearly diverges in the limit k 0 . Moreover, note that in this limit cos Ω ( k ) t becomes a constant as a function of time, which is the characteristic “freezing” of the response around the critical point.

Variance of the Magnetic Moment Per Spin

For time-independent problems, and for quasistatic driving the relaxation function becomes identical to the magnetic susceptibility, χ , [24]. Thus, we now evaluate χ = Ψ ( t = 0 ) for systems prepared in the zero mode, k = 0 by directly integrating Equation (16). We obtain,
χ = Γ 0 2 + J K 4 J Γ 0 ( Γ 0 + J ) 2 ( Γ 0 + J ) 2 E 4 J Γ 0 ( Γ 0 + J ) 2 π Γ 0 2 | Γ 0 + J | ,
where K and E are the complete elliptic integral of first and second kind [49]. Equation (19) is depicted in Figure 2. We observe that, as expected, at the critical point χ diverges, and that decays polynomially into the ferro- and paramagnetic phases.
This establishes that the relaxation function (12) in the limit N 1 exhibits important properties of a thermodynamic system undergoing a phase transition. Next, we will show how a corresponding relaxation time can be derived from Ψ ( t ) .

3.2. Relaxation Time

In linear response theory, the relaxation time, τ R , can be determined directly from the relaxation function [24]. We have
τ R = 0 d t Ψ ( t ) Ψ ( 0 ) ,
which we can now evaluate for the quantum Ising chain with Equation (12). Note, however, that for isolated quantum systems the relaxation function (12) is oscillatory, and hence Equation (20) is an indeterminate integral. Therefore, in Appendix B, we compute the upper envelop of the integral in Equation (20), for which we obtain
τ R UB = ħ ( J + Γ 0 ) 2 J 2 + Γ 0 2 8 J 2 Γ 0 2 1 | J Γ 0 | .
Equation (21) is plotted in Figure 3, which closely resembles Figure 1.
Remarkably, the relaxation time determined by means of linear response in Equation (21) is governed by the gap and we can write
τ R UB | J Γ 0 | 1 .
Consequently, the critical exponent ν = 1 , and more importantly τ R UB , gives a more transparent justification for the identification of the energy gap with a relaxation rate (2). Equation (21) constitutes our first main result. Rather than having to rely on plausibility arguments, the relaxation time in isolated quantum systems can be determined directly from the relaxation function of linear response theory.
It appears plausible that this finding holds generally for any many-body system exhibiting a quantum phase transition. However, a more sophisticated analysis or at least a numerical verification may be required to explore whether the inverse gap is generally related to the relaxation time from linear response theory.

4. Kibble–Zurek Scaling of the Excess Work

Now that we have established that both relaxation function as well as the corresponding relaxation time behave properly, it is tempting to directly compute the excess work (7). However, to guarantee that our comparison with predictions from the Kibble–Zurek mechanism are sound, we first need to more carefully analyze the range of validity of linear response theory around the critical point.

4.1. Range of Validity

To this end, we computed the exact excess work by solving the corresponding time-dependent Schrödinger equation using a standard Runge–Kutta method. The excess work can be written as
W ex = ψ ( τ ) | H ( τ ) | ψ ( τ ) ψ ( 0 ) | H ( 0 ) | ψ ( 0 ) Δ E
where Δ E is the exergy [27,50], which reduces to the energy difference of initial and final groundstates. Expressions for Δ E can be found in the literature [51]. The numerically exact results can then be compared with the expression from linear response theory (7) for the relaxation function in Equation (10). In Figure 4, we plot our findings for a range of system sizes and “perturbation strengths”, and for processes starting in the ferromagnetic, Γ 0 > 1 , as well as the paramagnetic, Γ 0 < 1 , phases.
Intuitively, we would expect linear response theory to be accurate as long as the quantum Ising chain remains close to its ground state. Thus, a natural parameter to quantify the range of validity can be chosen to be δ Γ / ϵ 1 1 , where δ Γ denotes the “strength” of the driving and ϵ 1 is the energy difference between ground and excited state, i.e., the gap. Note that ϵ 1 0 for N . Thus, one would expect a failure of linear response theory for large systems, which means in the limit of “proper” phase transitions. In fact, in Figure 4, we observe very good agreement between the prediction of linear response theory and the exact numerics for small enough δ Γ / ϵ 1 . However, we also observe that for large δ Γ / ϵ 1 linear response theory still qualitatively captures the behavior of the excess work as a function of the external driving.
Note that the critical point is only crossed in Figure 4g–i. However, also for such processes we find regimes in which linear response theory accurately predicts the excess work, and in all other cases we have at least qualitatively accurate results. Thus, we can now continue to analyze the scaling behavior of W ex (7).

4.2. Kibble–Zurek Scaling from Linear Response Theory

Based on our understanding for when linear response theory is accurate, we can now verify the expected Kibble–Zurek scaling. To this end, we consider a case of N = 10 5 and a process that drives through the critical point at a constant rate (4). In complete analogy to Ref. [27] we consider only the excess work accumulated in the impulse regime,
W ex Im = J 2 τ 2 t ^ t ^ t ^ t d t d t Ψ ( t t ) .
Note that for each τ we have a corresponding value of t ^ (5), and that we choose Γ 0 = Γ ( t ^ ) . This is a fair analysis as the Kibble–Zurek arguments only depend on the rate of driving, and not on the initial values of the external field. The resulting values of W ex Im are plotted on a log-log scale as a function of the driving time τ in Figure 5. We observe polynomial behavior over three orders of magnitude, and the numerical Kibble–Zurek exponent γ KZ 1 . This is in full agreement with the aforementioned expectation, and we are now comfortable to conclude that the framework developed in Ref. [27] indeed also applies to quantum phase transitions.

5. Concluding Remarks

In the present analysis, we analyzed the consistency and interplay of two phenomenological frameworks to describe quantum phase transitions, namely the Kibble–Zurek mechanism and linear response theory. We found that while the Kibble–Zurek mechanism does go beyond the range of validity of linear response theory, additional insight can be obtained by studying both frameworks. A key finding of our analysis is that the relaxation time determined from linear response theory gives solid and rigorous justification for the plausibility argument that identifies the “gap” as a relaxation rate. Moreover, we found that the excess work computed from linear response theory exhibits the scaling properties that are predicted by the Kibble–Zurek arguments.

Author Contributions

Conceptualization, M.V.S.B. and S.D.; formal analysis, P.N.; writing—original draft preparation, P.N. and S.D.; writing—review and editing, M.V.S.B. and S.D.; supervision, S.D.; funding acquisition, P.N., M.V.S.B. and S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) (Brazil) (Grant No. 2018/06365-4, No. 2018/21285-7, and No. 2020/02170-4) and by CNPq (Conselho Nacional de Desenvolvimento Científico e Pesquisa) (Brazil) (Grant No. 141018/2017-8). S.D. acknowledges support from the U.S. National Science Foundation under Grant No. DMR-2010127.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Relaxation Function for the Quantum Ising Model

In this appendix, we briefly summarize the derivation of the relaxation function (10). Generally, the response function ϕ ( t ) is defined by
ϕ ( t ) = 1 i ħ [ Γ H ( 0 ) , Γ H ( t ) ] 0 ,
where ( ) 0 is the average under the canonical ensemble and the time evolution is given by Heisenberg equations for the Hamiltonian of the initial ground state. We now express the operators in a basis where the Hamiltonian is diagonal. Following the procedure outlined in Ref. [18], we first use the Jordan-Wigner transformation, which maps the spin chain onto an equivalent system of spinless fermions
σ j x = ( c j + c j ) i < j ( 1 2 c i c i ) , σ j z = 1 2 c j c j ,
where c j and c j are the creation and annihilator fermionic operators. The Hamiltonian becomes
H = J j = 1 N ( c j c j + 1 + c j c j + 1 + H . c . ) Γ j = 1 N ( 1 2 c j c j ) ,
where c N + 1 = c 1 , given the periodic boundary conditions. The next step is applying a Fourier transform to the fermionic operators
c j = e i π / 4 N k K c k e i k j ,
where
K = { ± ( 2 n 1 ) π / N , n = 1 , , N / 2 } .
Therefore, the Hamiltonian can be written as
H = k K [ ( Γ J cos k ) ( c k c k c k c k ) + J ( c k c k + c k c k ) sin k ] .
It is convenient to express the pair ( k , k ) by means of one number k only. Thus, we can write
H = 2 k K + [ ( Γ J cos k ) ( c k c k c k c k ) + J ( c k c k + c k c k ) sin k ]
where
K + = { ( 2 n 1 ) π / N , n = 1 , , N / 2 } .
The diagonalization of the Hamiltonian is performed using the Bogoliubov transformation in each one of the k K + modes. These Bogoliubov transformations U k are unitary transformations, given by
U k = u k * v k * v k u k ,
with
u k = cos θ k 2 , v k = sin θ k 2 ,
where
sin θ k = sin k J 2 + Γ 2 2 J Γ cos k ,
and
cos θ k = Γ cos k J 2 + Γ 2 2 J Γ cos k .
Hence, the Hamiltonian can now be written as
H = k K + ϵ k ( γ k γ k + γ k γ k 1 ) ,
where
ϵ k = 2 J 2 + Γ 2 2 J Γ cos k ,
and
γ k γ k = U k c k c k ,
which are fermionic operators as well. The solutions of the Heisenberg equations of the operators γ k and γ k are then given by
γ k ( t ) = γ k e i ϵ k t ħ , γ k ( t ) = γ k e i ϵ k t ħ .
To calculate the response function, we calculate the derivative of the Hamiltonian using Equation (A7), since its Bogolibouv transformation depends implicitly on the magnetic fields on the fermionic operators γ k and γ k . We have
Γ H = 2 k K + ( c k c k c k c k )
The crucial step now is to observe that the response function is invariant if the operators involved are transformed by Bogoliubov transformations. In particular, the derivative of the Hamiltonian becomes
Γ H = k K + [ 4 | u k | 2 γ k γ k + 4 | v k | 2 γ k γ k 4 u k v k ( γ k γ k + γ k γ k ) 2 ]
Finally, we can express the response function in terms of a sum in each mode k K +
ϕ ( t ) = 1 i ħ k K + [ Γ H k ( 0 ) , Γ H k ( t ) ] k ,
where
H = k K + H k ,
and ( ) k denotes a thermal average with ρ k = exp β H k / tr exp β H k , where β is the inverse temperature. Collecting expressions we finally arrive at
ϕ ( t ) = 32 ħ n = 1 N / 2 J 2 ϵ n 2 sin 2 ( 2 n 1 ) π N × sin 2 ϵ n t ħ tan β 2 ϵ k ,
which is an odd function and ϕ ( 0 ) = 0 . In the zero temperature limit β , the response function becomes
ϕ ( t ) = 32 ħ n = 1 N / 2 J 2 ϵ n 2 sin 2 ( 2 n 1 ) π N sin 2 ϵ n t ħ .
whose derivative is the desired expression (10).

Appendix B. The Upper Envelop for the Relaxation Time

Finally, we show how to compute the relaxation function from linear response theory. Generally, the relaxation time is given by
τ R = 1 χ 0 0 π d k d t A ( k ) cos ( Ω ( k ) t )
To calculate the integral, we consider a finite time T, first,
L = 8 J 2 π 0 T 0 π d k d t sin 2 ( k ) ϵ 3 ( k ) cos 2 ϵ ( k ) t ħ .
This can be evaluated and we obtain
L = 4 ħ J 2 π 0 π d k sin 2 ( k ) ϵ 4 ( k ) sin 2 ϵ ( k ) T ħ .
An envelop is then readily given by the trigonometric inequality
L L = 4 ħ J 2 π 0 π d k sin 2 ( k ) ϵ 4 ( k ) .
and we have
L ħ | δ Γ | δ Γ 2 + 2 J 2 2 δ Γ J δ Γ 2 | 2 J δ Γ | 16 δ Γ 2 ( J δ Γ ) 2 | 2 J δ Γ |
which holds for any T. Finally using, ϵ ( k ) 2 | J Γ 0 | , we can write
1 χ 8 | J + Γ 0 | 3 J 2 ,
which leads to Equation (21).

References

  1. Callen, H.B. Thermodynamics and an Introduction to Thermostatistics, 2nd ed.; John Wiley & Sons: New York, NY, USA, 1985. [Google Scholar]
  2. Fisher, M.E. The renormalization group theory of critical behavior. Rev. Mod. Phys. 1974, 46, 597. [Google Scholar] [CrossRef]
  3. Kibble, T.W. Topology of cosmic domains and strings. J. Phys. A Math. Gen. 1976, 9, 1387. [Google Scholar] [CrossRef]
  4. Zurek, W.H. Cosmological experiments in superfluid helium? Nature 1985, 317, 505–508. [Google Scholar] [CrossRef] [Green Version]
  5. Zurek, W.H. Cosmological experiments in condensed matter systems. Phys. Rep. 1996, 276, 177. [Google Scholar] [CrossRef] [Green Version]
  6. Laguna, P.; Zurek, W.H. Density of Kinks after a Quench: When Symmetry Breaks, How Big are the Pieces? Phys. Rev. Lett. 1997, 78, 2519. [Google Scholar] [CrossRef] [Green Version]
  7. Biroli, G.; Cugliandolo, L.F.; Sicilia, A. Kibble-Zurek mechanism and infinitely slow annealing through critical points. Phys. Rev. E 2010, 81, 050101. [Google Scholar] [CrossRef] [Green Version]
  8. Chandran, A.; Erez, A.; Gubser, S.S.; Sondhi, S.L. Kibble-Zurek problem: Universality and the scaling limit. Phys. Rev. B 2012, 86, 064304. [Google Scholar] [CrossRef] [Green Version]
  9. Chandran, A.; Burnell, F.J.; Khemani, V.; Sondhi, S.L. Kibble-Zurek scaling and string-net coarsening in topologically ordered systems. J. Phys. Condens. Matt. 2013, 25, 404214. [Google Scholar] [CrossRef] [Green Version]
  10. Ulm, S.; Roßnagel, J.; Jacob, G.; Degünther, C.; Dawkins, S.T.; Poschinger, U.G.; Nigmatullin, R.; Retzker, A.; Plenio, M.B.; Schmidt-Kaler, F.; et al. Observation of the Kibble-Zurek scaling law for defect formation in ion crystals. Nat. Commun. 2013, 4, 2290. [Google Scholar] [CrossRef]
  11. Del Campo, A.; Kibble, T.W.B.; Zurek, W.H. Causality and non-equilibrium second-order phase transitions in inhomogeneous systems. J. Phys. Condens. Matter 2013, 25, 404210. [Google Scholar] [CrossRef] [Green Version]
  12. Partner, H.L.; Nigmatullin, R.; Burgermeister, T.; Pyka, K.; Keller, J.; Retzker, A.; Plenio, M.B.; Mehlstäubler, T.E. Dynamics of topological defects in ion Coulomb crystals. New J. Phys. 2013, 15, 102013. [Google Scholar] [CrossRef]
  13. del Campo, A.; Zurek, W.H. Universality of phase transitions: Toplological defects from symmetry breaking. Int. J. Mod. Phys. A 2014, 29, 1430018. [Google Scholar] [CrossRef] [Green Version]
  14. Deutschländer, S.; Dillmann, P.; Maret, G.; Keim, P. Kibble-Zurek mechanism in colloidal monolayers. Proc. Natl. Acad. Sci. USA 2015, 112, 6925–6930. [Google Scholar] [CrossRef] [Green Version]
  15. Keesling, A.; Omran, A.; Levine, H.; Bernien, H.; Pichler, H.; Choi, S.; Samajdar, R.; Schwartz, S.; Silvi, P.; Sachdev, S.; et al. Quantum Kibble–Zurek mechanism and critical dynamics on a programmable Rydberg simulator. Nature 2019, 568, 207–211. [Google Scholar] [CrossRef] [PubMed]
  16. Zamora, A.; Dagvadorj, G.; Comaron, P.; Carusotto, I.; Proukakis, N.P.; Szymańska, M.H. Kibble-Zurek Mechanism in Driven Dissipative Systems Crossing a Nonequilibrium Phase Transition. Phys. Rev. Lett. 2020, 125, 095301. [Google Scholar] [CrossRef] [PubMed]
  17. Cui, J.M.; Gómez-Ruiz, F.J.; Huang, Y.F.; Li, C.F.; Guo, G.C.; del Campo, A. Experimentally testing quantum critical dynamics beyond the Kibble–Zurek mechanism. Commun. Phys. 2020, 3, 44. [Google Scholar] [CrossRef]
  18. Sachdev, S. Quantum phase transitions. In Handbook of Magnetism and Advanced Magnetic Materials; John Wiley & Sons: New York, NY, USA, 2007. [Google Scholar]
  19. Zurek, W.H.; Dorner, U.; Zoller, P. Dynamics of a Quantum Phase Transition. Phys. Rev. Lett. 2005, 95, 105701. [Google Scholar] [CrossRef] [Green Version]
  20. Scherer, D.R.; Weiler, C.N.; Neely, T.W.; Anderson, B.P. Vortex formation by merging of multiple trapped Bose-Einstein condensates. Phys. Rev. Lett. 2007, 98, 110402. [Google Scholar] [CrossRef] [Green Version]
  21. Weiler, C.N.; Neely, T.W.; Scherer, D.R.; Bradley, A.S.; Davis, M.J.; Anderson, B.P. Spontaneous vortices in the formation of Bose-Einstein condensates. Nature 2008, 455, 14. [Google Scholar] [CrossRef] [Green Version]
  22. Gardas, B.; Dziarmaga, J.; Zurek, W.H. Dynamics of the quantum phase transition in the one-dimensional Bose-Hubbard model: Excitations and correlations induced by a quench. Phys. Rev. B 2017, 95, 104306. [Google Scholar] [CrossRef] [Green Version]
  23. Kubo, R. Statistical-mechanical theory of irreversible processes. I. General theory and simple applications to magnetic and conduction problems. J. Phys. Soc. Jpn. 1957, 12, 570–586. [Google Scholar] [CrossRef]
  24. Kubo, R.; Toda, M.; Hashitsume, N. Statistical Physics II: Nonequilibrium Statistical Mechanics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  25. Soriani, A.; Nazé, P.; Bonança, M.V.S.; Gardas, B.; Deffner, S. Three phases of quantum annealing: Fast, slow, and very slow. Phys. Rev. A 2022, 105, 042423. [Google Scholar] [CrossRef]
  26. Soriani, A.; Nazé, P.; Bonança, M.V.; Gardas, B.; Deffner, S. Assessing performance of quantum annealing with non-linear driving. arXiv 2022, arXiv:2203.17009. [Google Scholar]
  27. Deffner, S. Kibble-Zurek scaling of the irreversible entropy production. Phys. Rev. E 2017, 96, 052125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Pfeuty, P. The one-dimensional Ising model with a transverse field. Ann. Phys. 1970, 57, 79–90. [Google Scholar] [CrossRef]
  29. Dziarmaga, J. Dynamics of a Quantum Phase Transition: Exact Solution of the Quantum Ising Model. Phys. Rev. Lett. 2005, 95, 245701. [Google Scholar] [CrossRef] [Green Version]
  30. Mbeng, G.B.; Russomanno, A.; Santoro, G.E. The quantum Ising chain for beginners. arXiv 2020, arXiv:2009.09208. [Google Scholar]
  31. Fusco, L.; Pigeon, S.; Apollaro, T.J.G.; Xuereb, A.; Mazzola, L.; Campisi, M.; Ferraro, A.; Paternostro, M.; De Chiara, G. Assessing the Nonequilibrium Thermodynamics in a Quenched Quantum Many-Body System via Single Projective Measurements. Phys. Rev. X 2014, 4, 031029. [Google Scholar] [CrossRef] [Green Version]
  32. Francuz, A.; Dziarmaga, J.; Gardas, B.; Zurek, W.H. Space and time renormalization in phase transition dynamics. Phys. Rev. B 2016, 93, 075134. [Google Scholar] [CrossRef] [Green Version]
  33. Gardas, B.; Dziarmaga, J.; Zurek, W.H.; Zwolak, M. Defects in Quantum Computers. Sci. Rep. 2018, 8, 4539. [Google Scholar] [CrossRef] [Green Version]
  34. Piccitto, G.; Silva, A. Dynamical phase transition in the transverse field Ising chain characterized by the transverse magnetization spectral function. Phys. Rev. B 2019, 100, 134311. [Google Scholar] [CrossRef] [Green Version]
  35. Puebla, R.; Deffner, S.; Campbell, S. Kibble-Zurek scaling in quantum speed limits for shortcuts to adiabaticity. Phys. Rev. Res. 2020, 2, 032020. [Google Scholar] [CrossRef]
  36. Carolan, E.; Kiely, A.; Campbell, S. Counterdiabatic control in the impulse regime. Phys. Rev. A 2022, 105, 012605. [Google Scholar] [CrossRef]
  37. Sivak, D.A.; Crooks, G.E. Thermodynamic Metrics and Optimal Paths. Phys. Rev. Lett. 2012, 108, 190602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Zulkowski, P.; Sivak, D.A.; Crooks, G.E.; DeWeese, M.E. Geometry of thermodynamic control. Phys. Rev. E 2012, 86, 041108. [Google Scholar] [CrossRef] [Green Version]
  39. Bonança, M.V.S.; Deffner, S. Optimal driving of isothermal processes close to equilibrium. J. Chem. Phys. 2014, 140, 244119. [Google Scholar] [CrossRef]
  40. Acconcia, T.V.; Bonança, M.V.S. Degenerate optimal paths in thermally isolated systems. Phys. Rev. E 2015, 91, 042141. [Google Scholar] [CrossRef] [Green Version]
  41. Acconcia, T.V.; Bonança, M.V.; Deffner, S. Shortcuts to adiabaticity from linear response theory. Phys. Rev. E 2015, 92, 042148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Bonança, M.V.S.; Deffner, S. Minimal dissipation in processes far from equilibrium. Phys. Rev. E 2018, 98, 042103. [Google Scholar] [CrossRef] [Green Version]
  43. Nazé, P.; Bonança, M.V. Compatibility of linear-response theory with the second law of thermodynamics and the emergence of negative entropy production rates. J. Stat. Mech. Theo. Exp. 2020, 2020, 013206. [Google Scholar] [CrossRef] [Green Version]
  44. Deffner, S.; Bonança, M.V.S. Thermodynamic control—An old paradigm with new applications. EPL (Europhys. Lett.) 2020, 131, 20001. [Google Scholar] [CrossRef]
  45. Bonança, M.V.S.; Nazé, P.; Deffner, S. Negative entropy production rates in Drude-Sommerfeld metals. Phys. Rev. E 2021, 103, 012109. [Google Scholar] [CrossRef] [PubMed]
  46. Bonança, M.V.S.; Deffner, S. Fluctuation theorem for irreversible entropy production in electrical conduction. Phys. Rev. E 2022, 105, L012105. [Google Scholar] [CrossRef] [PubMed]
  47. Polkovnikov, A. Universal adiabatic dynamics in the vicinity of a quantum critical point. Phys. Rev. B 2005, 72, 161201. [Google Scholar] [CrossRef] [Green Version]
  48. Mandal, D.; Jarzynski, C. Analysis of slow transitions between nonequilibrium steady states. J. Stat. Mech. 2016, 2016, 063204. [Google Scholar] [CrossRef] [Green Version]
  49. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; US Government Printing Office: Washington, DC, USA, 1964.
  50. Deffner, S.; Campbell, S. Quantum Thermodynamics: An Introduction to the Thermodynamics of Quantum Information; Morgan & Claypool Publishers: San Rafael, CA, USA, 2019. [Google Scholar]
  51. Katsura, S. Statistical Mechanics of the Anisotropic Linear Heisenberg Model. Phys. Rev. 1962, 127, 1508–1518. [Google Scholar] [CrossRef]
Figure 1. Illustration of the Kibble–Zurek mechanism. Far from the critical point, the dynamics of the system is essentially adiabatic, meaning that the system recovers from the defects of the driving faster than the inverse of the driving rate. Close to the critical point the situation changes dramatically. The healing capacity is lost and finite-size domains are “frozen” into the system.
Figure 1. Illustration of the Kibble–Zurek mechanism. Far from the critical point, the dynamics of the system is essentially adiabatic, meaning that the system recovers from the defects of the driving faster than the inverse of the driving rate. Close to the critical point the situation changes dramatically. The healing capacity is lost and finite-size domains are “frozen” into the system.
Entropy 24 00666 g001
Figure 2. Magnetic susceptibility (19) as a function of the external field Γ 0 for J = 1 .
Figure 2. Magnetic susceptibility (19) as a function of the external field Γ 0 for J = 1 .
Entropy 24 00666 g002
Figure 3. Effective relaxation time (21) for J = 1 .
Figure 3. Effective relaxation time (21) for J = 1 .
Entropy 24 00666 g003
Figure 4. Excess work (7) computed from linear response theory and exact numerics for protocols driving in the ferromagnetic (ac) and paramagnetic (df) phase, and crossing the critical point (gi). Figures (a,d,g) depict situations in which linear response theory and the exact result perfectly match. Figures (b,e,h) depict situations with large N. Figures (c,f,i) depict situations with strong driving.
Figure 4. Excess work (7) computed from linear response theory and exact numerics for protocols driving in the ferromagnetic (ac) and paramagnetic (df) phase, and crossing the critical point (gi). Figures (a,d,g) depict situations in which linear response theory and the exact result perfectly match. Figures (b,e,h) depict situations with large N. Figures (c,f,i) depict situations with strong driving.
Entropy 24 00666 g004
Figure 5. Comparison between Kibble–Zurek scaling of the excess work (7) from exact dynamics and linear response theory.
Figure 5. Comparison between Kibble–Zurek scaling of the excess work (7) from exact dynamics and linear response theory.
Entropy 24 00666 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nazé, P.; Bonança, M.V.S.; Deffner, S. Kibble–Zurek Scaling from Linear Response Theory. Entropy 2022, 24, 666. https://doi.org/10.3390/e24050666

AMA Style

Nazé P, Bonança MVS, Deffner S. Kibble–Zurek Scaling from Linear Response Theory. Entropy. 2022; 24(5):666. https://doi.org/10.3390/e24050666

Chicago/Turabian Style

Nazé, Pierre, Marcus V. S. Bonança, and Sebastian Deffner. 2022. "Kibble–Zurek Scaling from Linear Response Theory" Entropy 24, no. 5: 666. https://doi.org/10.3390/e24050666

APA Style

Nazé, P., Bonança, M. V. S., & Deffner, S. (2022). Kibble–Zurek Scaling from Linear Response Theory. Entropy, 24(5), 666. https://doi.org/10.3390/e24050666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop