Effects Induced by the Initial Condition in the Quantum Kibble–zurek Scaling for Changing the Symmetry-breaking Field

The Kibble–Zurek scaling describes the driven critical dynamics starting with an equilibrium state far away from the critical point. Recently, it has been shown that scaling behaviors also exist when the fluctuation term changes starting near the critical point. In this case, the relevant initial conditions should be included in the scaling theory as additional scaling variables. Here, we study the driven quantum critical dynamics in which a symmetry-breaking field is linearly changed starting from the vicinity of the critical point. We find that, similar to the case of changing the fluctuation term, scaling behaviors in the driven dynamics can be described by the Kibble–Zurek scaling with the initial symmetry-breaking field being included as its additional scaling variable. Both the cases of zero and finite temperatures are considered, and the scaling forms of the order parameter and the entanglement entropy are obtained. We numerically verify the scaling theory by taking the quantum Ising model as an example.


Introduction
Recent experimental progress made in cold-atom systems has led to many studies on the nonequilibrium behaviors in quantum many-body systems [1][2][3].Among them, the Kibble-Zurek mechanism, which was first proposed by Kibble in cosmology [4], and then by Zurek in condensed matter physics [5], has been generalized to describe the driven quantum critical dynamics starting with a state far away from the critical point [6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23].According to the Kibble-Zurek mechanism [1][2][3], the whole driven process can be divided into three stages: one impulse stage sandwiched by two adiabatic stages.At the first stage, the energy gap is very large, the excitation induced by the external driving can be ignored, and thus the system evolves adiabatically along the equilibrium state; then, when the system is near its critical point, the gap becomes smaller and cannot suppress the excitations induced by the external driving, and the system is in an impulse stage in which the external driving dominates.After this stage, the system enters another adiabatic stage.At the boundary between this latter adiabatic stage and the impulse stage, the Kibble-Zurek mechanism shows that the number of topological defects can be scaled by the driving rate [1][2][3].Recently, the Kibble-Zurek mechanism has been verified experimentally in trapped-ion systems and Bose-Einstein condensates [24][25][26][27].Furthermore, the full scaling forms have been proposed to explain the scaling behaviors in the whole process [18][19][20][21][22][23]28,29].Additionally, these scaling forms have been employed to numerically estimate the critical point and critical exponents in both classical and quantum phase transitions [28][29][30][31][32].
However, when the driven dynamics are started up in the vicinity of the critical point, the initial effects will be remembered for a macroscopic time scale because the relaxation time near the critical point is very large.As a result, the initial condition will affect the driven dynamics.This inspires many studies on how to generalize the Kibble-Zurek mechanism to take the effects induced by the initial conditions into account in both classical and quantum phase transitions [33][34][35][36][37][38].
Since thermal effects cannot be perfectly excluded in real experiments, the study on the quantum critical dynamics at finite temperatures [35][36][37][38][39] is another important issue.For the driven dynamics in a closed system, the thermal effects can only be involved in the initial condition [35][36][37][38].It has been demonstrated that for a closed system, the thermal effects involved in the driven critical dynamics are closely related to the initial distance to the critical point [38].
In contrast to previous studies, in which the fluctuation term is changed from the vicinity of the critical point [35][36][37][38], here we study the quantum driven critical dynamics by changing a symmetry-breaking field starting with an equilibrium state near the critical point.We will show that the scaling behaviors also exist.Similar to the case of changing the fluctuation term, the initial symmetry-breaking field should be included in the full scaling form as its additional scaling variable.The general scaling forms will be obtained for both the zero temperature and finite temperature cases.In particular, we will consider the scaling forms of the order parameter and the entanglement entropy.The scaling theory will be verified in the one-dimensional (1D) quantum Ising model.
The rest of the paper is organized as follows.In Section 2, we briefly review the Kibble-Zurek scaling (KZS) and its generalization.Then, in Section 3, we show the scaling theory of driven critical dynamics for changing the symmetry-breaking field starting near the critical point.Section 4 will introduce the 1D quantum Ising model and the numerical method.Then, the numerical results to verify the scaling theory are shown in Section 5. Finally, we will give a summary in Section 6.

Kibble-Zurek Scaling and Its Generalization
Here we give a brief review on the KZS.The original KZS often considers the case of changing the fluctuation field starting from a state far from the critical point and only focuses on the scaling of the number of topological defects at the boundary between the impulse and adiabatic region, while the dynamics in the impulse region are assumed to be frozen [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17].However, this is apparently not the real situation.Finite-time scaling (FTS) [28,29] provides a theory to understand the dynamics in the impulse region and build a full scaling theory of the KZS.The FTS shows that an effective time scale, which is induced by the external driving, dominates the dynamics in the impulse region, since in this region this effective time scale is the shortest time scale compared with other macroscopic time scales, such as the reaction time scale of the system.This is similar to the finite-size scaling theory, which states that the lattice size dominates the critical behavior when it is smaller than the correlation length.
We start from the usual finite-size scaling theory.Near the critical point, the scale transformation of the order parameter M reads as [39] in which g is the distance of the fluctuation term to the critical point and L is the lattice size, β is defined as M ∼ (−g) β for g < 0 [39], and ν is defined as ξ ∼|g| −ν with ξ being the correlation length.The finite-size scaling theory shows that when L ξ, the scaling behavior is dominated by L. Accordingly, M should satisfy the finite-size scaling form [39].
in which f 1 is an analytic scaling function and the effect induced by g is only a perturbation.In contrast, when L In this case, g dominates the scaling behavior and L is a perturbation.Similar to the finite-size scaling theory, the FTS theory explains the driven dynamic scaling behavior by noting that the external driving will introduce an effective time scale ζ e , which competes against the relaxation time ζ s ∼ ξ z [28,29], in which z is the dynamic.To be explicit, we consider the case for changing the fluctuation term g as g = R g t + g 0 with g 0 being far from the critical point.In this case, the driving-induced effective time scale is ζ e ∼ R g −z/r g , in which r g = z + 1/ν.Similar to Equation ( 1), the evolution of the order parameter M under the driving should satisfy In the FTS region, in which ζ e ζ s , ζ e controls the dynamics and the external driving dominates.By setting R g b r g = 1 in Equation ( 1), one can obtain the scaling form of M, Compared with the original KZS, one finds that the FTS region is just the impulse region in the original KZS.When ζ e ∼ ζ s , gR g −1/νr g is a constant and Equation ( 5) gives M ∝ R g β/νr g , which is consistent with the original KZS.Moreover, the FTS shows that the system in the impulse region does not cease to evolve.Instead, its evolution time scale is given by ζ e .
When g 0 is near the critical point, it has been shown that the driven dynamics of M still exhibit scaling behaviors, which can be described by [38] where g 0 has been included in the scaling form.Since g 0 brings another time scale into the dynamics, i.e., the initial time scale ζ i ∼ g 0 −νz , which is larger than ζ e ∼ R g −z/r g for a small g 0 , there is an initial crossover stage, bridging the region dominated by ζ i and the FTS region.However, the critical slowing down makes this crossover stage macroscopically long.Thus, g 0 affects the driven dynamics during the whole process.Another generalization of the KZS is that the thermal effects can be taken into account by including the temperature T in the scaling theory [21,[35][36][37][38].However, it has been shown that for a closed quantum system, the status of temperature T is closely related to the initial distance to the critical point [38].To see the reason, one should at first note that for a closed system, the temperature involved in the driven dynamics can only be from the initial condition.Only when g 0 is near the critical point is the initial gap small compared with the temperature, and the number of thermal excitations can be macroscopically large enough to affect the following dynamics.Accordingly, the evolution of M starting with a thermal state near the critical point satisfies [38]

Effects Induced by the Initial Condition in the Kibble-Zurek Scaling for Changing the Symmetry-Breaking Field
Besides the case of changing g, it has been shown that the dynamics of changing the symmetry-breaking field h as h = R h t + h 0 for a large h 0 also shows a scaling behavior [22,23,28,29,34].In this case, h 0 is irrelevant.According to the FTS theory [28,29,34], the effective time scale induced by the external driving is ζ e ∼ R h −z/r h , in which r h = z + βδ/ν, with δ being defined as M ∼ h 1/δ at g = 0.
Similar to the case of a small g 0 , when h 0 is small, h 0 should also be included in the scaling theory.Accordingly, the order parameter satisfies in which we have set g = 0 for simplicity.
In Equation ( 8), when h 0 is large, the initial stage is an adiabatic stage.In this case, h 0 is irrelevant and Equation ( 8) restores the usual scaling form of FTS [22,23,28].However, when h 0 is small, in the initial stage, the characterized time scale is ζ i ∼ h 0 −νz/βδ , and M ∼ h 0 1/δ .Since ζ i ζ e , the system will cross over to the FTS region.The cooperation between these crossover dynamics and the driven dynamics will make the dynamic behavior different from the usual KZS.For example, at the impulse-adiabatic boundary, at which hR h −1/νr h is a constant, M depends on the driving rate R h and h 0 .However, for the usual FTS, M only depends on R h .When the initial state is an equilibrium state at finite T, the scaling theory should also include the temperature T, Similar to the case of changing g [38], here T only has nontrivial effects when h 0 is small.When h 0 is large, the number of initial thermal excitations is tiny and thus has no effect on the following dynamics.At T = 0, the entanglement entropy S also demonstrates scaling behaviors affected by the initial condition [40][41][42][43][44].The entanglement entropy S is defined as S = −Tr(ρlogρ), in which the base of the logarithm is 2 throughout and ρ is the reduced density matrix of half of the system.For the 1D system near the critical point, the entanglement entropy scales as S = (c/6)logξ, where c is the central charge [40][41][42][43][44].
For h = R h t + h 0 , the scaling form of ξ reads Therefore, one can obtain the scaling form of the entanglement entropy S, The scaling behavior of the entanglement entropy has been studied in the driven dynamics of the changing g [45].It has been shown that for a large g 0 and fixed gR g −1/νr g , S ∼ −(c/6r g )logR g .Here, we demonstrate that similar scaling behavior also exists for changing the symmetry-breaking field.

Model and Numerical Method
The 1D quantum Ising model in a longitudinal field will be employed to illustrate the scaling theory.The Hamiltonian is [39] in which σ z n and σ x n are Pauli matrices in the z and x directions, respectively, at site n; w x is the transverse field which is the fluctuation term; and h is the longitudinal field which breaks the Z 2 symmetry.The order parameter M is defined as M ≡< σ z n >, where the angle brackets denote the expectation value of the operator at site n.The critical point for Equation ( 12) is w xc = 1 (and thus g ≡ w x − 1) and h c = 0.The exact critical exponents are β = 1, ν = 1, δ = 15, z = 1 [39], r h = 23/8 [28,29], and the central charge c = 1/2 [40,41].In the following, we will consider the case for g = 0.
We will employ the infinite time-evolving block decimation (iTEBD) algorithm [45] to simulate the driven dynamics.This algorithm is a variant of the time-dependent density-matrix renormalization group method [46][47][48].For the zero temperature case, two steps must be implemented: first, calculating the ground state at h 0 ; second, calculating the evolution of M. To realize the first step, one should at first prepare a random state, and represent it in a matrix product state form via Vidal's decomposition [46], and then the ground state can be determined by imposing the imaginary time evolution operator on the matrix product state.The second step can be realized by a similar procedure: from the ground state, time dynamics can be calculated by imposing the evolution operator on the matrix product state, and at every time step, h should be changed by a rate R h .
For the driven dynamics starting from a thermal equilibrium state, we also need two steps.The first step is to calculate the thermal equilibrium state at h 0 , i.e., ρ 0 = exp[−H(h 0 )/T]/Z, in which Z is the partition function.To do this, one should at first purify the identity matrix I into a pure state |ϕ I by introducing an auxiliary system [47,48], which is maximally entangled to the physical system.After representing |ϕ I in the matrix product state, one can obtain a "thermal" pure state |ϕ T ≡ exp[−H/(2T)]|ϕ I by imposing the imaginary time evolution operator on |ϕ I [47,48].The second step is to calculate the real time evolution state |ϕ(t) T , which can be directly calculated by imposing the real time evolution operator on |ϕ T .It has been shown that the dynamics of the order parameter M can be directly calculated from |ϕ(t) T [47,48].
In the following calculations, the time interval is chosen to be 0.005 and 100 states are kept for both the zero temperature and finite temperature cases.

Numerical Results
We consider the driven dynamics of Equation ( 12) for changing the symmetry-breaking field h as h = R h t + h 0 .For simplicity, we fix g = 0, i.e., w x = 1.In Figure 1, we compare the driven dynamics starting in the vicinity of the critical point and the driven dynamics starting far from the critical point.One finds that the curves for these two initial conditions are quite different.In Figure 1a, one finds that the hysteresis effect of M for the evolution starting from the vicinity of the critical point is weaker than that for the evolution starting far from the critical point.In contrast, in Figure 1b, we find the hysteresis for the dynamics starting from the vicinity of the critical point is stronger than that for the dynamics starting far from the critical point.
To verify Equation ( 8), we plot M versus h for a small fixed h 0 R h −1/νr h in Figure 2.After rescaling M versus h by R h in Figure 2b, we find that the rescaled curves collapse onto each other, verifying Equation (8).Therefore, we confirm that although the driven dynamics are affected by the initial condition, scaling behaviors also exist and the initial symmetry-breaking field h 0 should be included in the scaling form.
Entropy 2016, 18, 446 6 of 9    The scaling theory including 0 h can also be verified according to the scaling behavior of the entanglement entropy S .According to Equation ( 11), we calculate ( / 6 )log 3 shows that the curves for different driving rates h R collapse onto each other after rescaling.These results confirm Equation (11).Then, we verify the scaling theory at the finite temperature.For a small fixed TR − , we find the evolution curves for different driving rates match each other after rescaling in Figure 4.This result confirms Equation (9).Since all the experiments are implemented The scaling theory including h 0 can also be verified according to the scaling behavior of the entanglement entropy S. According to Equation ( 11), we calculate S + (c/6r h )logR h for a small fixed h 0 R h −1/νr h .Figure 3 shows that the curves for different driving rates R h collapse onto each other after rescaling.These results confirm Equation (11).Then, we verify the scaling theory at the finite temperature.For a small fixed  Then, we verify the scaling theory at the finite temperature.For a small fixed h 0 R h −1/νr h and a fixed TR h −z/r h , we find the evolution curves for different driving rates match each other after rescaling in Figure 4.This result confirms Equation ( 9).Since all the experiments are implemented at finite temperatures, this result also provides the experimental feasibility to examine the present scaling theory.
at finite temperatures, this result also provides the experimental feasibility to examine the present scaling theory.

Summary and Discussion
In this paper, we have investigated the effects induced by the initial condition in the driven critical dynamics by changing the symmetry-breaking field.We have shown that when the initial symmetry-breaking field is small, the driven dynamics are quite different from the dynamics starting with a large symmetry-breaking field, which is required in the original finite-time scaling.In spite of this, we have demonstrated that the dynamics also exhibit scaling behaviors.By using the finite-time scaling theory, we obtain the scaling forms for the driven dynamics.In these scaling forms, the initial symmetry-breaking field should be included as an additional scaling variable.Furthermore, we also consider the finite temperature case.So this scaling theory can be examined in experiments.We have numerically verified the scaling theories by taking the 1D quantum Ising model as an example.We emphasize that for some experiments in which the usual Kibble-Zurek scaling cannot work, the present scaling theory may be a possible candidate to modify the scaling analysis.We also note that the Kibble-Zurek scaling has been generalized to the inhomogeneous cases, in which the relevant parameter depends on both the space and the time [49].It has been shown that scaling behaviors in these systems are quite different from the original Kibble-Zurek mechanism [49].So, it is appealing to explore the finite-time scaling and the initial effects therein.

Summary and Discussion
In this paper, we have investigated the effects induced by the initial condition in the driven critical dynamics by changing the symmetry-breaking field.We have shown that when the initial symmetry-breaking field is small, the driven dynamics are quite different from the dynamics starting with a large symmetry-breaking field, which is required in the original finite-time scaling.In spite of this, we have demonstrated that the dynamics also exhibit scaling behaviors.By using the finite-time scaling theory, we obtain the scaling forms for the driven dynamics.In these scaling forms, the initial symmetry-breaking field should be included as an additional scaling variable.Furthermore, we also consider the finite temperature case.So this scaling theory can be examined in experiments.We have numerically verified the scaling theories by taking the 1D quantum Ising model as an example.We emphasize that for some experiments in which the usual Kibble-Zurek scaling cannot work, the present scaling theory may be a possible candidate to modify the scaling analysis.We also note that the Kibble-Zurek scaling has been generalized to the inhomogeneous cases, in which the relevant parameter depends on both the space and the time [49].It has been shown that scaling behaviors in these systems are quite different from the original Kibble-Zurek mechanism [49].So, it is appealing to explore the finite-time scaling and the initial effects therein.

Figure 1 .
Figure 1.Comparison of the evolution of (a) the order parameter M and (b) the entanglement entropy S under decreasing h with different initial conditions but the identical driving rate 0.009 h R =

Figure 1 .
Figure 1.Comparison of the evolution of (a) the order parameter M and (b) the entanglement entropy S under decreasing h with different initial conditions but the identical driving rate R h = 0.009.Arrows show the driven direction.

Figure 1 ..
Figure 1.Comparison of the evolution of (a) the order parameter M and (b) the entanglement entropy S under decreasing h with different initial conditions but the identical driving rate 0.009 h R =

Figure 3 .
Figure 3.At 0 T = , the evolution of S under decreasing h with fixed 1/ 0

Figure 2 .
Figure 2. At T = 0 the evolution of M under decreasing h with fixed h 0 R h −1/νr h = 0.002159 for different three driving rates indicated.The curves before and after being rescaled are shown in (a) and (b), respectively.

Figure 1 ..
Figure 1.Comparison of the evolution of (a) the order parameter M and (b) the entanglement entropy S under decreasing h with different initial conditions but the identical driving rate 0.009 h R =

Figure 3 .
Figure 3.At 0 T = , the evolution of S under decreasing h with fixed 1/ 0 − , we find the evolution curves for different driving rates match each other after rescaling in Figure4.This result confirms Equation(9).Since all the experiments are implemented

Figure 3 .
Figure 3.At T = 0, the evolution of S under decreasing h with fixed h 0 R h −1/νr h = 0.002159 for different three driving rates indicated is shown in (a).The curves for S + (c/6r h )logR h versus the rescaled h are shown in (b).

Figure 4 .
Figure 4.The evolution of M under decreasing h with fixed 1/ 0

Figure 4 .
Figure 4.The evolution of M under decreasing h with fixed h 0 R −1/νr h = 0.002159 and TR h −z/r h = 1.583 for four T indicated.The curves before and after rescaling are shown in (a) and (b), respectively.