Probing the Possibilities of Ergodicity in the 1D Spin-1/2 XY Chain with Quench Dynamics

Ergodicity sits at the heart of the connection between statistical mechanics and dynamics of a physical system. By fixing the initial state of the system into the ground state of the Hamiltonian at zero temperature and tuning a control parameter, we consider the occurrence of the ergodicity with quench dynamics in the one-dimensional (1D) spin-1/2 XY model in a transverse magnetic field. The ground-state phase diagram consists of two ferromagnetic and paramagnetic phases. It is known the magnetization in this spin system is non-ergodic. We set up two different experiments as we call them single and double quenches and test the dynamics of the magnetization along the Z-axis and the spin-spin correlation function along the X-axis which are the order parameters of the zero-temperature phases . Our exact results reveal that for single quenches at zero-temperature, the ergodicity depends on the initial state and the order parameter. In single quenches for a given order parameter, ergodicity will be observed with an ergodic-region for quenches from another phase, non-correspond to the phase of the order parameter, into itself. In addition, a quench from a ground-state phase point corresponding to the order parameter into or very close to the quantum critical point, hc = 1.0, discloses an ergodic behavior. Otherwise, for all other single quenches, the system behaves non-ergodic. Interestingly on the other setup, a double quench on a cyclic path, ergodicity is completely broken for starting from the phase corresponding to the order parameter. Otherwise, it depends on the first quenched point, and the quench time T when the model spent before a second quench in the way back which gives an ability to controlling the ergodicity in the system. Therefore, and contrary to expectations, in the mentioned model the ergodicity can be observed with probing quench dynamics at zero-temperature. Our results provide further insight into the zero-temperature dynamical behavior of quantum systems and their connections to the ergodicity phenomenon.

What connects a physical system to the real world is the unitary time evolution of it. A typical scenario in this context is the quantum quench i.e., driven out-of-equilibrium of a system by abruptly changing a control parameter 21,22 where the behavior of the system basically is not susceptible to general principles of equilibrium system 23,24 . The issue is exciting when a quench is on 25,26 or crossed from critical points 27,28 where the system undergoes a non-analytic change of its properties. There are some approaches to understanding quench dynamics in many-body systems such as the Kibble-Zurek mechanism 29 and measurement quench 30,31 .
Here, we are going to study quantum ergodicity at zero-temperature with the use of quench dynamics of the two quantities, the magnetization along the Z-axis, and the spin-spin correlation function along the X-axis which reveals the magnetization along the X-axis. We apply two conditions throughout our study: (i) the initial ground state of the Hamiltonian will be chosen as the initial state of the system (ii) the system will be examined at zero temperature. It should note that recently using these conditions, a new approach of quench dynamics was introduced known as the dynamical quantum phase transition 32 that has been investigated by using different analytical and numerical techniques [33][34][35][36][37][38][39][40][41] , and experimental point of view [42][43][44] . Anyway, using these conditions, we consider 1D spin-1/2 XY model in a transverse field which its ground state phase diagram includes two phases, ferromagnetic (FM) (with the non-zero value of the magnetization along the X-axis) and paramagnetic (PM) (with the almost saturated value of the magnetization along the Z-axis). Single and double quenches are considered. The long-time run of the dynamical quantities are compared with their zero-temperature equilibrium values. Our exact results demonstrate the possibility of the occurrence of ergodicity in the mentioned model. We find that, in a single quench, the ergodicity depends on the starting point and its ordering at zero temperature. For a given order parameter the ergodicity can be observed in two situations, when a quench is done: (i) from another phase, non-corresponding to the phase of the order parameter, into itself which leads to emerge of an ergodic-region, and (ii) from the phase corresponding to the order parameter into or near to the quantum critical point, h h f c 1 = . On the other hand, for all other single quenches, non-ergodic behaviors arise in the system. We also discuss the essential effect of the excited states on the creation of the ergodic behavior in the system. Different behaviors are found for cyclic quenches. In a cyclic quench, the system starts from point i and goes to point f 1 by spending a time T. After timing spend, the system is quenched into the starting point f 2 = i. Based on our results, ergodicity at zero-temperature will be broken for starting from the phase which is corresponding to the order parameter. In contrast, for starting from an initial state in a phase irrelevant to the order parameter, both ergodicity and non-ergodicity can appear depending on the first queched point f 1 and the timing spend in this phase point. This outcome gives surprising flexibility in controlling ergodicity by performing double quenches. Consequently, this paper highlights the physical conditions under which the dynamical quantum system may disclose ergodic behaviors at zero-temperature.

the Model
The Hamiltonian of the 1D spin-1/2 XY model in the presence of a transverse field is given by where σ n is the Pauli spin operator on the n-th site, δ and h are the anisotropy parameter and the homogeneous external transverse magnetic field, respectively. N is the system size (or number of spins). The system is considered in the thermodynamic limit, N → ∞. The quantum phase transition from the ferromagnetically ordered phase to the paramagnetic phase driven by the transverse field h is called the Ising transition. On the other hand, the quantum phase transition between two ferromagnetically ordered phases, with magnetic ordering in the X-direction and the Y-direction, respectively, driven by the anisotropy parameter δ, is called the anisotropic transition. In fact, in the absence of the transverse field, the ground state of the system is in the Luttinger liquid phase at δ = 0. Ferromagnetic ordered phase is found in the presence of anisotropy, 0 < δ ≤ 1. The quantum phase transition into the paramagnetic phase is happened at the quantum critical point h c (δ) = 1 45,46 . This Hamiltonian conserves the parity of the particle number and acts differently on the even (Neveu-Schwarz) and odd (Ramond) subspaces. In the fermionic Fock space, the Hamiltonian in the two subspaces is formally the same if one imposes antiperiodic boundary condition for the even and periodic boundary condition for the odd subspace that in wave-number space these boundary conditions translate to different quantization as momentum quantization in half-integer and in integer multiples of N 2π respectively. In the thermodynamic limit the ground states of the odd and even subspaces become degenerate and one recovers the two fully polarized ferromagnetic ground states.
The Hamiltonian is integrable and can be mapped to a system of free fermions and therefore be solved exactly. By applying the Jordan-Wigner transformation 47 , the Hamiltonain converts from spin operators into spinless fermionic operators as n N n n n n n n n n n n 1 The zero-temperature expectation values of the magnetization along the Z-axis and the spin-spin correlation function along the X-axis at t = 0 as a function of the transverse field are shown in Fig. 1(a,b). It should be noted that these two quantities are defined as where the brackets symbolize the expectation value on the ground state of the system in the thermodynamic limit. As is seen in Fig. 1(a,b), the quantum phase transition occurs in the critical transverse field h c = 1. Only at δ = 0, the magnetization will be saturated ( Fig. 1(a)), and the spin-spin correlation function will be vanished ( Fig. 1(b)) at the critical field.

Setup and ergodicity
It is known that, ergodicity in the systems depends on the path where systems evolve with time. Here, we setup two strategies for studying ergodicity in the system. First, in the view of single quench as the tendency of dynamics of the system to match itself to the zero-temperature value of the final quenched point. Accordingly, the system needs to be investigated in a long-time run where it goes or fluctuates around a stable situation. Second, we apply the idea of double quenches to consider ergodicity in the system. For this case, we do: (1) A quench from an initial state in the phase A to a final phase point in the phase B as i to f 1 , (2) A quench from an initial state in the phase B to a final phase point in the phase A as f 1 to f 2 = i.
It means, the starting and end points are the same way (a cyclic quench). The setup is schematically represented in Fig. 1(c). In this way, our instrumentations are the magnetization along the Z -axis and the spin-spin correlation function along the X -axis. The time-dependent Hamiltonian H(t) that models a double quantum quench is that are the initial state at t = 0 and quenched state at t = T, respectively. We investigate our problem with two conditions: (i) Fixing the initial state of the system into the ground state of the Hamiltonian.
(ii) Considering the system at zero temperature.
It should note it was demonstrated that by the use of a double quench one can control the dynamical quantum phase transitions in the 1D spin-1/2 ITF model 48 . This controlling can be done simply by tuning the time between the first and the second quench. As a result, the system can exhibit all four combinations of absence or presence of non-analyticities before and after the second quench in the rate function of the return probability, respectively. In the following, as will be seen, the same situation for controlling ergodicity in the system will be achievable by controlling the quench time T in the dynamical behavior of the order parameters.

ergodicity in a single quench
Based on the quantum ergodic theorem, for an arbitrary initial state, the expectation value of the magnetization and the spin-spin correlation function will ultimately evolve in time to its value predicted by the theory of ensemble, and thereafter will exhibit only small fluctuations around equilibrium value. To consider the time evolution of a closed quantum system, one can use of quench dynamics, putting the system in an equilibrium state described with the Hamiltonian i k k i   = ∑ and the initial state | ⟩ Ψ 0 , afterward, suddenly changing the control parameters from their initial values to their final values. Final Hamiltonian and its time-evolved state will be as = ∑ , respectively. We choose the ground state of the initial Hamiltonian as the initial state, and study the system at zero temperature.
In the following, we focus on the dynamics of the magnetization along the transverse field and the spin-spin correlation function along the X-axis. The magnetization along the Z-axis is the order parameter of the PM phase. On the other side, the spin-spin correlation function along the X-axis reveals the ordering of the FM phase which is analogous to the magnetization along the X-axis. For a single quench where t ≤ T, we have obtained the time-dependent of these quantities in the thermodynamic limit as where k − is the difference between the Bogoliubov angles diagonalizing the pre-quench and post-quench Hamiltonians, respectively.
In the first step, we have focused on the dynamical behavior of the magnetization. We have considered a quench from an initial state in the FM phase. Results are illustrated in Fig. 2 for different values of δ = 1.0, 0.5 and h i = 0.5. We know that the zero-temperature state of the system is in the FM phase at h i = 0.5 and at the value of field h > h c = 1.0 goes into the PM phase where all spins are aligned along the transverse magnetic field. As is apparent from Fig. 2(a,b), for a single quench from the mentioned initial state in the FM phase into the PM region, the value of the magnetization for most times in the long run is very far from the zero-temperature saturated phase. Thus the spin-1/2 XY chain for the mentioned initial state behaves as a non-ergodic system in quenching from the FM into the PM regions. For a quench from the FM into the same FM region, we also have calculated the magnetization and results are presented in Fig. 2(c,d) for h i = 0.5. As is obvious, the system shows ergodic behavior for selected values of the transverse magnetic field. In fact, the magnetization ultimately evolves in time to its value predicted by ground state of quenched Hamiltonian, and thereafter will exhibit small fluctuations around zero-temperature value. We did the same quenches starting from an initial state in the FM region into all phase points. Results for the long-time run dynamics of the magnetization are displayed in Fig. 2(e,f) for δ = 1.0, 0.5 and h i = 0.2. The phenomenon of broken ergodicity is manifestly seen in this figure. In fact, by starting from a selected initial state in the FM phase and quenching into the same FM region, the system does not break the ergodicity in a wide region. Approaching the quenched point to the quantum critical region, deviation from the ergodic behavior starts to reveal. This phenomenon is obviously shown for quenching into a phase point in the PM region.  Fig. 3(a,b). As is seen, the width of the ergodic-region displays a non-monotonic behavior and will be maximized at h . This may be attributed to the fact that at this initial point, the system has a large coherence and behaves as a quantum memory which is capable of storing the quantum information [49][50][51] . Furthermore, approaching the quantum critical region, the width of the ergodic region decreases and will be minimized at the quantum critical point, h i = h c = 1. In other words, its width at the quantum critical point is approximately zero that indicates at quenching from h c into both the FM and the PM phase, system behaves non-ergodic.
In the second step, a quench from an initial state in the PM phase into different regions is considered. Results on the long-time run of magnetization are indicated in Fig. 4 for different values of δ = 1.0, 0.5 and h i = 4.0. It is explicitly shown that, only in a quench into a phase point very close to the quantum critical point, the system unveils an ergodic behavior. We have to mention that there is not any overlap between curves in the high magnetic region. One conclusion that can be drawn from the discussion up to now is that, the magnetization along the transverse magnetic field can be exhibited an ergodic behavior if a quench is done from a region where the magnetization: (1) is not the order parameter into the same region, (2) is the order parameter into the quantum critical point Let us more explain the role of the excited states emerging from the dynamics of the system to create ergodicity. The magnetization along Z-axis as a function of t can be rewritten as   www.nature.com/scientificreports www.nature.com/scientificreports/ In principle, at t → ∞, the Eq. ((8)) transvers to that indicates if a system behaves as a non-ergodic system and if = | | C 1 0 , thereby, the excited states are ineffective in the ergodicity of the system (| | C n must be zero for all n ≠ 0). In Fig. 5 the coefficient | | C 0 is plotted for δ = 1.0, 0.5 that clarity shows |C 0 | = 1 only happens when no quench is done, otherwise, | | < C 1 0 . In the other side, because it should be ∑ = = C 1 n N n 0 2 , consequently, the ergodic behavior arises if at least there is a ≠ | | C 0 n for n ≠ 0 in such a way that Eq. (10) can satisfy Eq. (9). As a result, it reveals for quenching into the PM phase, it is ever zero the existence of the excited states which lead to the behavior of ergodicity in the system. It also expresses for quenching into the FM phase, only the excited states in the ergodic-region have the main role in the ergodicity of the system.
Another result that lies in Fig. 2(e,d) is quenching into the PM phase where for a very large value of the magnetic field it leads to vanishing the long-time run of the magnetization. In this case, since in the PM phase the spins of the system behave as a non-interacting many-particle system, the quenched Hamiltonian can be approximately considered as the only Zeeman interaction,  Fig. 1(a)), which show that by quenching into the PM region, a small density of spin wave excitations will be created. On the other hand, the initial state with the saturated FM ordering along the X-axis is a superposition of all spin wave excitations with the same probability. Therefore, the system behaves non-ergodic when a quench into the PM phase will be done. Now, let us calculate the particular modes in the ergodic-region. Another expression of the Eq. (9) is Thus, from Eq. (7) for the ergodic-region one can obtain It should be stressed, comparing these modes with the particular modes which lead to the dynamical quantum phase transitions 32 discloses that they are different.
In addition to the magnetization along the transverse field, we have also studied the dynamical behavior of the spin-spin correlation function along the X-axis. Results are exhibited in Fig. 6 for the long-time run dynamics versus the transverse magnetic field for single quenches started from (a & b) the FM region with h i = 0.2, and (c & d) the PM region with h i = 1.5, into all points of the quantum states for δ = 1.0 and δ = 0.5 respectively. In general, this quantity behaves inversely with respect to the magnetization along the transverse field. Starting from an initial state in the FM region, no ergodic behavior is seen except in a quench into a final state very close to the quantum critical point, h c = 1.0. But, when we put the system initially in the PM phase, an ergodic region is observed for quenching into a final state in the same PM phase.

ergodicity in a cyclic quench
As we have mentioned, a cyclic quench is defined as a linked sequence of quenches where eventually returns the system to its initial starting point. The magnetization along the transverse field and the spin-spin correlation function along the X-axis at zero temperature as a function of time in a double quench (t ≥ T) are obtained as  www.nature.com/scientificreports www.nature.com/scientificreports/ depends on the value of h f 1 . For example in Fig. 7(a,b), for the first ergodicity time T, we have < = .
T T . On the other note, when the system is initially located in the PM phase, different behavior is observed. Fig. 7(c,d), shows that system behaves non-ergodic in all cyclic quenches independent of the anisotropy parameter, the middle quenched point, and spending time T.
We have also investigated the dynamics of the spin-spin correlation function along the X-axis for cyclic quenches. for δ = 1.0 and δ = 0.5 respectively. As we mentioned, in general, this quantity behaves conversely rather than the magnetization along the transverse field. Furthermore, in one state ergodicity comes out. This is when a quench is started from the PM phase and the first quenched point h f 1 will be as an appropriate phase point where results in the emerge of ergodicity with an ability to control it with quench time T. For other quenches, including quench to the critical point (for not very short time T), the system behaves non-ergodic.

conclusion
The phenomenon of the ergodicity is known as vital importance to the study of quantum many-body systems. A system is ergodic, if a macroscopic variable ultimately evolves in time to its value predicted by the theory of ensemble, and thereafter exhibit small fluctuations around an equilibrium value. Here, we tested the ergodicity analytically on the 1D spin-1/2 anisotropic XY model in a transverse magnetic field by focusing on quench dynamics throughout the zero-temperature phase space of the system.
It is known the model shows a quantum phase transition from an ordered FM phase to the PM phase at a critical point, h c . We let the model experiences quenches in two setups, single quench, and cyclic quench, with two conditions, choosing the ground state of the system as the initial state, and fixing the system at zero temperature. The Hamiltonian is integrable and the results are exact. We considered dynamics of the magnetization along the Z-axis and the spin-spin correlation function along the X-axis for both setups. As a consequence we clearly showed these two order parameters gave the reversely results. Our results displayed that for single quenches the ergodicity obviously depends on both the initial state and the order parameter. For a certain order parameter, ergodicity can be observed with an ergodic-region for quenches from a phase, non-correspond to the phase of the order parameter, into itself. Notably, for quenching from a phase point corresponding to the order parameter into or very close to the quantum critical point, h c = 1.0, the system behaves ergodic. Otherwise, for all other single quenches, the non-ergodic behavior of the system appears. Moreover, we also discussed the inevitable role of the excited states emerging from the dynamics of the system to create ergodicity. Remarkable results obtained for a cyclic quench. In this case, the ergodicity is broken for quenches started from a phase corresponding to the order parameter, independent of the anisotropy parameter, the first quenched point h f 1 , and the quench time T. With a double quench started from a ground-state phase, non-corresponding to the order parameter, by tuning the elapsing time at a given intermediate phase, one can achieve control over ergodicity in the system. Consequently, our results show for a quantum system at zero-temperature, both ergodicity, and non-ergodicity phenomena can appear dependence of the initial state, the middle quenched point, and also the quench time T. It should note we believe these results have a relationship with the conception of heating 53 . Therefore, we give a suggestion one can consider the problem from this view. We hope the presented results provide further evidence for unveiling a relationship between ergodicity, order parameters, and the non-equilibrium dynamics of quantum systems.