## Comprehensive Study of Dynamical Phase Transitions in One-Dimensional XY Model

Bao-Ming Xu ∗

Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China

## Xiu-Ying Zhang

School of Physics and Optoelectronic Engineering, Ludong University, YanTai 264025, China (Dated: Submitted September 11, 2024)

Short-time transient dynamical phase transition (SDQPT), serving as the theoretical framework on far-from-equilibrium physics of quantum many-body systems, has been experimentally observed recently. However, the existence of different types of SDQPTs and the associated conditions are yet to be elusive. In this paper, we comprehensively classify SDQPTs by considering the onedimensional transverse field quantum XY model initialized in the ground state. We find that there are only six types of SDQPTs, and three of which are completely new that had not yet been identified previously. The associated conditions for different SDQPTs are discussed in detail, and based on this, the dynamical phase diagrams are given. Our analytical method used to classify SDQPTs in one-dimensional transverse field quantum XY model can be applied to any other two-band models in one- and two-spatial dimensions.

## I. INTRODUCTION

The past decades have witnessed the flourishing of dynamical quantum phase transitions (DQPTs), both in experimental and theoretical research, among condensed matter physicists, because quantum simulators have nowadays achieved experimental access to the realtime dynamics of closed quantum many-body systems; see Refs. [1-3] for the recent reviews. Such quantum simulators have been realized on various experimental platforms, such as ultra-cold atoms [4-7] or trapped ions [8-11]. DQPTs are the sudden or non-analytic changes in the behaviour of a large quantum system during its time evolution, and generally refer to two largely independent concepts: long-time relaxation DQPT and shorttime transient DQPT. It should be noted that the abbreviation of DQPT refers to the dynamical quantum phase transition in general, and does not distinguish whether it belongs to a long-time relaxation dynamical phase transition or a short-time transient dynamical phase transition. Some papers represent long-time relaxation DQPT as DQPT-I, and short-time transient dynamical phase transition as DQPT-II. In this paper, for the convenience of discussion, we use LDQPT to represent long-time relaxation DQPT and SDQPT to represent short-time transient DQPT. LDQPT describes the dynamical relaxation behaviors of two-point correlators [12-19], which is related to the thermalization of closed quantum systems [20]. The difference of the correlators at a time t from their respective steady state values, reached asymptotically, scales as a power law in time, t -µ [14]. The transition is reflected as the change of the exponent of this power-law relaxation behavior. The long-time limit of physically relevant observables, namely the order parameters, are finite for one phase, whereas they vanish for another, resembling the equilibrium Landau theory [3]. SDQPT generalizes the nonanalytic behavior of the free energy at a phase transition in the thermodynamic limit to the out-of-equilibrium case [2, 21]. To this end, the Loschmidt echo, measuring the return probability of the time-evolving system onto its initial state, has been introduced. SDQPT is defined by the zeros in the Loschmidt echo or the nonanalyticities in its rate function obtained by large deviation principle. It is analogous to the zeros in the partition function and the nonanalyticities in the free energy in equilibrium phase transitions, and therefore the rate function of Loschmidt echo is also called the dynamical free energy per particle. Some SDQPTs obey the dynamical scaling defined by a corresponding outof-equilibrium analog of the universality class [22, 23]. Both kinds of DQPTs have aroused great interests, and this paper focuses exclusively on SDQPT.

∗ Electronic address: xbmv@bit.edu.cn

SDQPTs were originally introduced in the nearestneighbor transverse-field Ising model [21] and was generalized to mixed state [24-26], finite temperature [27], Floquet systems [28-39], slow quench process [40-43] and Kitaev honeycomb model [44]. Even SDQPT can also be generated by changing the boundary conditions [45]. Efforts have also been made to generalize this concept to non-Hermitian systems [46-51], nonintegrable models [52-58] and disorder systems [23, 59-61]. Experimentally, SDQPTs have been observed in various systems, ranging from trapped ions [62], ultracold atoms [63] to superconducting qubits [64], and some other platforms [65-67]. The creation of different phases is not immediately apparent from the nature of the effect, but some dynamical order parameters, such as longitudinal magnetization [21, 68], string order parameter [58, 69], the winding number [40, 69, 70], vortex loops [25, 71] and some others [7274] have been identified so far. Entanglement, as a key concept for the characterization of equilibrium quantum phases and criticality [75, 76], is naturally expected to be generalized to DQPTs. In pursuit of this goal, the connections between entanglement dynamics and SDQPTs

were explored [32, 42, 62, 77-83], although the principles underlying these connections are unclear on a general level. In addition, quantum coherence has been shown to be a valuable resource for SDQPTs [84, 85].

Although the links between the SDQPTs and many physical observables have been established, a thorough understanding of this transition still calls for more studies. To the best of our knowledge, the general relation between SDQPTs and the underlying equilibrium quantum critical point is unclear. The seminal paper [21] suggests that SDQPT occurs when the quenching process is ramped through an equilibrium critical point, seemingly implying that SDQPT is closely related to equilibrium phase transitions [68]. However, it was later demonstrated that this connection lacks rigor, meaning that SDQPT can occur for the quench within the same phases, and it may not occur for the quench crossing the equilibrium critical point [86]. Although SDQPT has been disentangled from equilibrium phase transitions, the conditions for it remain hazy. Beyond that, some SDQPTs lead to topological numbers with integer jumps [29, 39, 40, 46-49, 69, 70, 85, 87-89], while others cause topological numbers with half-integer jumps [47, 48, 89, 90]. The existence of other types of SDQPTs and the associated conditions require further investigation. This paper conducts a comprehensive study of SDQPTs by taking into account the one-dimensional twoband noninteracting translationally invariant Hamiltonian, identifying six distinct types of SDQPTs and delineating their respective conditions. Two-band model is simple and has analytical results, but captures some of the essential physics, such as symmetry protected topological phases [91]. Two-band models are the effective model of spin chains, including the transverse field Ising and XY chains. In this study, we examine the onedimensional transverse field XY model, chosen as a prototypical representation of the broader class of two-band models. Such two-band systems are widely studied in the literature, and several of them have been realized experimentally, in condensed matter or with analog quantum simulators using ultracold atoms in optical lattices [9295].

This paper is organized as follows: In the next section, we introduce the quenched one-dimensional XY model, a prototypical two-band model, and briefly review the crucial method for the analytical solution, which will be used in the following. Sec. III is devoted to study the central objects of SDQPT, including Loschmidt echo, the rate function, Fisher zeros and the topological order parameter of the winding number. In Sec. IV, a relatively comprehensive classification of SDQPTs is provided, along with the associated criteria. According to the conditions under which different SDQPTs occur, the dynamical phase diagrams are given in Sec. V. Finally, Sec. VI closes the paper with some concluding remarks.

## II. QUENCHED ONE-DIMENSIONAL XY MODEL

The Hamiltonian of the one-dimensional XY model is [96]

$$\underset { \text{rough} } { \text{many} } \quad \hat { H } ( \lambda, \gamma ) = - \frac { J } { 2 } \sum _ { j = 1 } ^ { N } \left [ \frac { 1 + \gamma } { 2 } \hat { \sigma } _ { j } ^ { x } \hat { \sigma } _ { j + 1 } ^ { x } + \frac { 1 - \gamma } { 2 } \hat { \sigma } _ { j } ^ { y } \hat { \sigma } _ { j + 1 } ^ { y } + \lambda \hat { \sigma } _ { j } ^ { z } \right ], \\ \text{lation} \quad \text{where $1$ is non-until-min-min} \text{ $\pi$} \text{ $\sigma$} \text{ over} \text{s}.$$

where J is longitudinal spin-spin coupling, γ governs the anisotropic coupling between spins along the x and y directions, λ is a dimensionless parameter measuring the strength of the transverse field with respect to the longitudinal spin-spin coupling. In this work, we set J = 1 as the overall energy scale without loss of generality. ˆ σ α j ( α = x, y, z ) is the spin-1/2 Pauli operator at lattice site j and the periodic boundary conditions are imposed as ˆ σ α N +1 = ˆ σ α 1 . Here we only consider that N is even. The model exhibits competitions between anisotropic and magnetic couplings, which results in the existence of multiple phases. The quantum phase transition from the ferromagnetic phase (FM) to the paramagnetic phase driven by the transverse field λ is called the Ising transition with the quantum critical point λ c = 1. On the other hand, the quantum phase transition between two FMs, with magnetic ordering in the x -direction and the y -direction, respectively, driven by the anisotropy parameter γ , is called the anisotropic transition with the critical point γ c = 0. In fact, in the absence of the transverse field, the ground state of the system is in the Luttinger liquid phase at γ c = 0.

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 and the Fourier transformation, the Hamiltonain converts from spin operators into spinless fermionic operators as [97, 98]

$$\inf _ { \substack { \text{a pro-} \\ \vdash \text{band} \\ \text{in the} \\ \text{exper-} \\ \text{antum} } } \hat { H } ( \lambda, \gamma ) = \sum _ { k } \left ( \hat { c } _ { - k } ^ { \dagger } \hat { c } _ { k } \right ) \begin{pmatrix} \lambda + \cos k & - i \gamma \sin k \\ i \gamma \sin k & - \lambda - \cos k \end{pmatrix} \begin{pmatrix} \hat { c } _ { - k } \\ \hat { c } _ { k } ^ { \dagger } \end{pmatrix},$$

where ˆ c k and ˆ c † k are respectively fermion annihilation and creation operators for mode k = (2 n -1) π/N with n = 1 · · · N/ 2, corresponding to antiperiodic boundary conditions for N is even. Each ˆ H k ( λ, γ ) acts on a twodimensional Hilbert space generated by { ˆ c † -k ˆ c † k | 0 〉 , | 0 〉} , where | 0 〉 is the vacuum of the Jordan-Wigner fermions ˆ c k and ˆ c -k , and can be represented in that basis by a 2 × 2 matrix

$$\hat { H } _ { k } ( \lambda, \gamma ) = \mathbf d _ { k } \cdot \hat { \sigma }$$

with

$$\mathbf d _ { k } = ( 0, \gamma \sin k, \lambda + \cos k ). \text{ \quad \ \ } ( 4 )$$

The instantaneous eigenvalues are d ± k = ± d k with

$$d _ { k } = \sqrt { ( \lambda + \cos k ) ^ { 2 } + \gamma ^ { 2 } \sin ^ { 2 } k }. \quad \quad ( 5 )$$

The corresponding eigenvectors are

$$| d _ { k } ^ { + } \rangle & = \left [ \frac { 1 + d _ { k } ^ { z } / d _ { k } } { \sqrt$$

and

$$| d _ { k } ^ { - } \rangle & = \left [ i \frac { d _ { k } ^ { y } / d _ { k } } { \sqrt {$$

respectively.

At time t ≤ 0, the system is prepared in the ground state

$$| G \rangle = \bigotimes _ { k } | d _ { k } ^ { - } \rangle. \quad \quad ( 8 ) \quad \quad \quad$$

Then at time t = 0, the transverse field is suddenly changed from λ to λ ′ or the anisotropic parameter is suddenly changed from γ to γ ′ . We assume that this process is so sudden that the system state has no time to change. After quenching, the dynamics of the system is governed by the evolution operator ( /planckover2pi1 = 1)

$$\hat { U } ( t ) = e ^ { - i \hat { H } ( \lambda ^ { \prime }, \gamma ^ { \prime } ) t } = \bigot$$

so that the state | ψ t ( ) 〉 at a time t after the quench is given by

$$| \psi ( t ) \rangle = \bigotimes _ { k } e ^ { - i \hat { H } _ { k } ( \lambda ^ { \prime }, \$$

## III. THE FUNDAMENTAL THEORY OF SDQPT

The central object within the theory of SDQPTs is the Loschmidt overlap amplitude quantifying the deviation of the time-evolved state from the initial condition. Here, the Loschmidt overlap amplitude is generalized to general density matrices defined as

$$\mathcal { G } ( t ) = \langle G | \hat { U } ( t ) | G \rangle. \text{ \quad \ \ } ( 1 1 ) \text{ \ \ } \text{ be}$$

Because different modes are independent with each other, the Loschmidt overlap amplitude can be written as

$$\mathcal { G } ( t ) = \prod _ { k } \mathcal { G } _ { k } ( t ) = \prod _ { k } \langle d _ { k } ^ { - } | e ^ { - i \hat { H } _ { k } ( \lambda ^ { \prime }, \gamma ^ { \prime } ) t } | d _ { k } ^ { - } \rangle. \quad ( 1 2 ) \quad \text{num} \quad \text{is def} \\ \text{--} \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \end{cases}$$

The behavior of the momentum-wise Loschmidt overlap amplitude can be visualized by the vector introduced by Ding [90]:

$$\vec { R } _ { k } = ( x _ { k }, y _ { k } ) = | \mathcal { G } _ { k } ( t ) | e ^ { i \phi _ { k } ^ { G } ( t ) }, \quad \ \ ( 1 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where

$$\phi _ { k } ^ { G } ( t ) = \phi _ { k } ( t ) - \phi _ { k } ^ { D } ( t ) \quad \quad ( 1 4 ) \quad \text{the}$$

FIG. 1: (Color online) The schematic drawing of the Bloch sphere.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-b4df7310-2783-4f84-b2e1-d7c497b7e5d1.jpg)

TABLE I: The conditions for SDQPTs.

/negationslash

/negationslash

| SDQPTs   | conditions                                                  | conditions                                                                          | conditions   | conditions           |
|----------|-------------------------------------------------------------|-------------------------------------------------------------------------------------|--------------|----------------------|
| SDQPTs   | γ → γ ′                                                     | λ → λ ′                                                                             | ∆            | γγ ′ + λλ ′ 1 - γγ ′ |
| SDQPT-1  | γγ ′ = 0 γ = 0, γ ′ = 0 γ = 0, γ ′ = 0 γγ ′ ( γγ ′ - 1) = 0 | ( | λ | - 1)( | λ ′ | - 1) < 0 | λ ′ | < 1 | λ | < 1, λ ′ = λ λ = λ ′ , | λ ′ | = 1 | -            | - - - < 1            |
| SDQPT-2  | γγ ′ ( γγ ′ - 1) = 0                                        | ( | λ | - 1)( | λ ′ | - 1) > 0                                                      | > 0          | < 1                  |
| SDQPT-3  | γγ ′ ( γγ ′ - 1) = 0                                        | -                                                                                   | =0           | < 1                  |
| SDQPT-4  | γγ ′ = 1 γ = 0, γ ′ = 0 γγ ′ ( γγ ′ - 1) = 0                | | λ | = 1, λ ′ = λ                                                                  | - - > 0      | - - > 1              |
| SDQPT-5  | γγ ′ ( γγ ′ - 1) = 0                                        | | λ | = 1, λ ′ = λ                                                                  | > 0          | < 1                  |
| SDQPT-6  | γγ ′ = 1                                                    | | λ | = 1, λ ′ = - λ                                                                | -            | -                    |

/negationslash

/negationslash

/negationslash is the geometric phase of the Loschmidt overlap amplitude with

$$\phi _ { k } ( t ) = \arg \mathcal { G } _ { k } ( t ) \text{ \quad \ \ } ( 1 5 )$$

being the total phase of G k ( ) and t

$$\phi _ { k } ^ { D } ( t ) = - \langle d _ { k } ^ { - } | \hat { H } _ { k } ( \gamma ^ { \prime }, \lambda ^ { \prime } ) | d _ { k } ^ { - } \rangle t. \quad \ \ ( 1 6 )$$

being the dynamical phase. At a given time, the winding number of the geometric phase in the first Brillouin zone is defined as [69]

$$\nu ( t ) = \frac { 1 } { 2 \pi } \int _ { 0 } ^ { \pi } \frac { \partial \phi _ { k } ^ { G } ( t ) } { \partial k } d k. \quad \quad ( 1 7 ) \\.. \. \. \. \. \. \. \.$$

In general, it takes a quantized jump whenever the evolution of the system passes through a critical time of SDQPT; therefore, it is generally regarded as the dynamical topological order parameter of SDQPT. It is obvious that ν t ( ) can be understood as the winding number of the trajectory /vector R k around the origin.

/negationslash

/negationslash

/negationslash

/negationslash

/negationslash

/negationslash

/negationslash

/negationslash

FIG. 2: (Color online) Lines of Fisher zeros (a1 and a2), the time evolution of rate function r t ( ) (blue curves in b1 and b2), the winding number ν t ( ) (red curves in b1 and b2) and the trajectory of vector /vector R k = ( x , y k k ) (c1 and c2) for the quenches from λ = 0 5 to . λ ′ = 1 5 (a1-c1) and . λ ′ = -1 5 (a2-c2). . Trajectories are plotted before (red bubbles), at (blue bubbles) and after (green bubbles) the first critical time t ∗ 1 . The other parameters are γ = 0 5, . γ ′ = 2 and n = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-4c7f40ca-a790-473d-bd7e-fa86e273b629.jpg)

FIG. 3: (Color online) Lines of Fisher zeros (a), the time evolution of rate function r t ( ) (blue curves in b), the winding number ν t ( ) (red curves in b) and the trajectory of vector /vector R k = ( x , y k k ) (c) for the quenches from λ = 1 1 to . λ ′ = 1 3. . The first and the second critical times are t ∗ 1 ≈ 2 07 and . t ∗ 2 ≈ 3 95. . The other parameters are γ = 1, γ ′ = -0 5 and . n = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-4a62e2d6-ab06-4878-8c67-67dc9976de91.jpg)

The rate function of the Loschmidt echo, measuring the return probability |G ( ) t | 2 , is defined as

$$r ( t ) = - \lim _ { N \to \infty } \frac { 1 } { N } \ln | \mathcal { G } ( t ) | ^ { 2 }. \quad \quad ( 1 8 ) \quad \stackrel { \text{aler} } { \langle d _ { k } ^ { - } } \\ \text{the}$$

It is the dynamical analogues of the equilibrium free energy, namely the dynamical free energy. In the thermodynamic limit one can derive an exact result for r t ( ):

the values of z that make Z z ( ) = 0. Due to the multiplicative properties of Eq. (20), a zero in Z z ( ) is equivalent to finding at least one mode k and one z making 〈 d -k | e -zH ˆ k ( λ ,γ ′ ′ ) | d -k 〉 = 0. In the thermodynamic limit the zeros of the boundary partition function in the complex plane coalesce to a family of lines labeled by a number n ∈ Z

$$r ( t ) = - \frac { 1 } { \pi } \text{Re} \left [ \int _ { 0 } ^ { \pi } d k \ln \mathcal { G } _ { k } ( t ) \right ]. \quad \text{(19)} \\ \text{O} \, \text{O} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$} \, \text{$\dots$}$$

A SDQPT will occur when r t ( ) exhibits a nonanalyticity at some times t ∗ , namely critical times. A powerful method to analyse the nonanalyticity of r t ( ) is Fisher zeros [21, 99]. In order to apply this concept, we should expend time t into the complex plane and focus on the boundary partition function Z z ( ) = 〈 G e | -zH λ ,γ ˆ ( ′ ′ ) | G 〉 . The boundary partition function can be expressed as

$$Z ( z ) = \prod _ { k } \langle d _ { k } ^ { - } | e ^ { - z \hat { H } _ { k } ( \lambda ^ { \prime }, \gamma ^ { \prime } ) } | d _ { k } ^ { - } \rangle \quad \ \ ( 2 0 ) \quad \ \ \\ \circ \, \ast \, \leq \, \cap \, \text{For information} \, \ast \, - \, i t \text{ this inet} \, \text{aecoribase} \quad \text{In } c$$

where z ∈ C . For imaginary z = it this just describes the overlap amplitude of Eq. (11). The Fisher zeros are

$$& \quad z _ { n } ( k ) = \frac { 1 } { 2 d _ { k } ^ { \prime } } \left [ \ln \frac { 1 - \frac { \frac { \frac { \Delta _ { k } } { d _ { k } } } \cdot \frac { \frac { \Delta _ { k } ^ { \prime } } { d _ { k } ^ { \prime } } } { 1 + \frac { \frac { \Delta _ { k } } { d _ { k } } } \cdot \frac { \frac { \Delta _ { k } ^ { \prime } } { d _ { k } ^ { \prime } } } { d _ { k } ^ { \prime } } } + i ( 2 n + 1 ) \pi \right ]. \quad ( 2 1 ) \\ &.$$

SDQPTs, manifested as the non-analytic behaviors of the dynamical free energy, are therefore contingent upon the existence of purely imaginary zeros of the boundary partition function, corresponding to real-time zeros of the Loschmidt overlap amplitude. From Eq. (21), we can obtain the well known condition of SDQPT d k · d ′ k = 0, i.e.,

$$( \lambda + \cos k ) ( \lambda ^ { \prime } + \cos k ) + \gamma \gamma ^ { \prime } \sin ^ { 2 } k = 0. \quad ( 2 2 )$$

In other words, given a set of parameters ( λ, γ, λ , γ ′ ′ ), a SDQPT can only occur if at least one mode satisfy Eq.

(22) (or d k · d ′ k = 0), this mode is called critical mode k ∗ . Conversely, solving Eq. (22) yields the critical mode k ∗ , and discussing the existence of the solution yields the concrete parameters ( λ, γ, λ , γ ′ ′ ) for SDQPT. Then, there is a series of critical times t ∗ n

$$t _ { n } ^ { * } = \frac { 1 } { 2 d _ { k ^ { * } } ^ { \prime } } ( 2 n + 1 ) \pi \, \quad \ \ ( 2 3 ) \, \quad \ \text{will}$$

/negationslash depending on the corresponding critical mode k ∗ . Fisher zeros are the key basis for determining whether SDQPT occurs, which have been used to discuss the relation between SDQPTs and equilibrium phase transitions [86]. Beyond that, Fisher zeros are also the key basis for determining the type of SDQPT, which will be used to achieve the main goal of this paper: classifying SDQPTs and discussing the classification criteria. In this paper, the modes k = 0 and k = π is named as the edge modes, and the others ( k = 0 or k = π ) are named as the edge modes. The critical mode can be either a body or an edge mode. Based on the number and attribution of the critical mode, as well as the properties of Fisher zeros, SDQPTs can be classified into six types, which will be discussed in the next section in detail. SDQPTs in the two-band mode have been classified by the ground state topology of both the initial and final Hamiltonians [100]. In contrast, our classification is based on the attribution and number of critical modes, and each resulting type of SDQPTs has its own topological properties.

/negationslash

Physically, the condition of SDQPT can be understood with the help of the Bloch sphere. According to the rate function, SDQPT occurring is equivalent to making 〈 d -k ∗ | e -iH ˆ k ∗ ( λ ,γ ′ ′ ) t | d -k ∗ 〉 = 0, which means a flip of the vector of | d -k ∗ 〉 in the Bloch sphere (see Fig. 1). It is conceivable that the axis of rotation must be perpendicular to a vector in order to flip it (see the green axis in Fig. 1). Therefore, SDQPT is only possible if the axis corresponding to the evolution operator ˆ U k ∗ ( ) t = e -iH ˆ k ∗ ( λ ,γ ′ ′ ) t is perpendicular to the initial state | d -k ∗ 〉 , otherwise it is impossible (see the blue axis in Fig. 1). The axis corresponding to the evolution operator can be understood as d ′ k ∗ and the vector of initial state can be understood as d k ∗ , therefore the condition of SDQPT can be demonstrated as d k ∗ · d ′ k ∗ = 0.

/negationslash

/negationslash

However, it must be noted that d k ∗ · d ′ k ∗ = 0 is not always effective in causing SDQPT. Specifically, the discussions above have acquiesced d k ∗ = 0 and d ′ k ∗ = 0, i.e., the energy levels of the pre- and post-quench systems are both non-degenerated. Breaking this acquiescence, i.e., d k ∗ = 0 or d ′ k ∗ = 0, d k ∗ · d ′ k ∗ = 0 does not always imply the perpendicularity, because the orientation of zero vector is uncertain and can no longer be used to represent the orientation of the initial state or the axis of evolution. In this subtle case, we use whether the initial state can be flipped as a criterion to judge whether SDQPT occurs, which is a more fundamental and thoroughgoing approach to determine SDQPT. If d ′ k ∗ = 0, the initial state | d -k ∗ 〉 would not evolve at all, let along

/negationslash undergo SDQPT, because ˆ U k ∗ ( ) = t e -iH ˆ k ∗ ( λ ,γ ′ ′ ) t = I ( I is an identity matrix). If d k ∗ = 0 but d ′ k ∗ = 0, the orientation of the evolution operator would still be represented by d ′ k ∗ , and as long as it is perpendicular to the orientation of the initial state (not d k ∗ ), the initial state would be flipped and SDQPT would occur. The analysis above will be concretized in the following sections.

## IV. THE CLASSIFICATION OF DYNAMICAL PHASE TRANSITIONS

In this section, we mainly classify SDQPTs and discuss the parameter spaces of different SDQPTs based on Eq. (22). First, we give the solution to Eq. (22). If γγ ′ = 1, Eq. (22) tells us that

$$\cos k = - \frac { 1 + \lambda \lambda ^ { \prime } } { \lambda + \lambda ^ { \prime } } \, \text{ \quad \ \ } ( 2 4 )$$

/negationslash with the condition λ + λ ′ = 0. If γγ ′ = 1, the solutions of Eq. (22) are

/negationslash

$$\cos k _ { \pm } = \frac { 1 } { 2 ( 1 - \gamma \gamma ^ { \prime } ) } \left [ - ( \lambda + \lambda ^ { \prime } ) \pm \sqrt { \Delta } \right ] \quad ( 2 5 ) \\..................$$

with the condition

$$\ell _ { 0 } \Delta = ( \lambda + \lambda ^ { \prime } ) ^ { 2 } - 4 ( 1 - \gamma \gamma ^ { \prime } ) ( \gamma \gamma ^ { \prime } + \lambda \lambda ^ { \prime } ) \geq 0. \quad ( 2 6 )$$

Specially, if γγ ′ = 0, these two solutions would be simplified as

$$\cos k _ { 1 } = - \lambda, \ \cos k _ { 2 } = - \lambda ^ { \prime }. \quad \quad ( 2 7 )$$

And if λ = ± 1 or λ ′ = ± 1, the solutions (25) would be reduced to

$$\cos k _ { 1 } & = \mp 1, \ \cos k _ { 2 } = \mp \frac { \lambda \lambda ^ { \prime } + \gamma \gamma ^ { \prime } } { 1 - \gamma \gamma ^ { \prime } }. \quad \ ( 2 8 ) \\ \cdot \,. \,. \quad. \quad. \quad. \quad. \quad. \quad.$$

These solutions allow for single, two, or even infinite critical modes to cause SDQPT.

## A. SDQPT-1

SDQPT-1 is caused by a single critical body mode, wherein the line of Fisher zeros cuts Im( z ) axis. First, we discuss the parameter space of SDQPT-1.

For γγ ′ = 1, a critical body mode k ∗ = arccos[ -1+ λλ ′ λ + λ ′ ] can be found if | 1+ λλ ′ λ + λ ′ | &lt; 1. And we find that ( d k&lt;k ∗ · d ′ k&lt;k ∗ )( d k&gt;k ∗ · d ′ k&gt;k ∗ ) &lt; 0, which means that the line of Fisher zeros cuts Im( z ) axis at this critical body mode, thereby implying the occurrence of SDQPT-1. The condition | 1+ λλ ′ λ + λ ′ | &lt; 1 means that

$$( | \lambda | - 1 ) ( | \lambda ^ { \prime } | - 1 ) < 0. \text{ \quad \ \ } ( 2 9 )$$

FIG. 4: (Color online) Lines of Fisher zeros (a), the time evolution of rate function r t ( ) (blue curves in b), the winding number ν t ( ) (red curves in b) and the trajectory of vector /vector R k = ( x , y k k ) (c) for the quenches from λ = 1 1 to . λ ′ = (1 -2 γγ ′ ) λ -√ 4 γγ ′ (1 -γγ ′ )(1 -λ 2 ). The first critical times is t ∗ 1 ≈ 2 47. . The other parameters are γ = 1, γ ′ = -0 5 and . n = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-d891f61e-bfc0-4a46-ac9e-eb2945522fc9.jpg)

FIG. 5: (Color online) Lines of Fisher zeros (a1 and a2), the time evolution of rate function r t ( ) (blue curves in b1 and b2), the winding number ν t ( ) (red curves in b1 and b2) and the trajectory of vector /vector R k = ( x , y k k ) (c1 and c2) for the quenches from λ = 1 to λ ′ = 1 5 (a1-c1) and . λ ′ = -1 5 (a2-c2). . Trajectories are plotted before (red bubbles), at (blue bubbles) and after (green bubbles) the first critical time t ∗ 1 . The other parameters are γ = 0 5, . γ ′ = 2 and n = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-b201efdb-86b6-43a8-9d18-153f12dd3d4d.jpg)

In other words, if the transverse field is quenched across the quantum critical point | λ c | = 1, i.e., | λ &lt; | 1 →| λ ′ | &gt; 1 or | λ &gt; | 1 →| λ ′ | &lt; 1, SDQPT-1 would occur.

/negationslash

/negationslash

If γγ ′ = 0, two critical body modes k ∗ 1 = arccos( -λ ) and k ∗ 2 = arccos( -λ ′ ) could be found for | λ | &lt; 1 and | λ ′ | &lt; 1. However, whether SDQPT-1 occurs depends on the quench protocol of anisotropic parameter. If γ = 0 → γ ′ = 0, the initial state would be | d -k 〉 = 0 , which could be flipped by the evolution op-| 〉 erator ˆ U k ∗ 2 ( ) t = e -id ′ k ∗ 2 σ t y but could not by the evolution operator ˆ U k ∗ 1 ( ) t = e -iH ˆ k ∗ 1 ( λ ,γ ′ ′ ) t , therefore SDQPT only occurs at the critical body mode k ∗ 2 = arccos( -λ ′ ). This SDQPT belongs to SDQPT-1 because the line of Fisher zeros cuts Im( ) axis z at this critical mode. If γ = 0 → γ ′ = 0, SDQPT-1 would occur at the critical mode k ∗ 1 = arccos( -λ ), but would not at the critical mode k ∗ 2 = arccos( -λ ′ ) because d ′ k ∗ 2 = 0, which prevents the initial state from evolving at all.

/negationslash

Given γγ ′ ( γγ ′ -1) = 0, one critical mode k ∗ = arccos(cos k ± ) ( | cos k ± | &lt; | cos k ∓ | ) could be found if ( | cos k + | -1)( | cos k -| -1) &lt; 0, implying that SDQPT-1 occurs. The condition ( cos | k + | -1)( | cos k -| -1) &lt; 0

also means ( | λ | -1)( | λ ′ | -1) &lt; 0 that the transverse field is quenched crossing the quantum critical point | λ c | = 1. It can be proven that ( | λ | -1)( | λ ′ | -1) &lt; 0 implies ∆ &gt; 0, therefore only ( | λ | -1)( | λ ′ | -1) &lt; 0 deserves consideration in the parameter space of SDQPT1. In addition, for | λ ′ | = 1, two critical modes k ∗ 1 = π and k ∗ 2 = arccos -λλ ′ + γγ ′ 1 -γγ ′ ( λ ′ = 1) or k ∗ 1 = 0 and k ∗ 2 = arccos λλ ′ + γγ ′ 1 -γγ ′ ( λ ′ = -1) could be found if

$$\left | \frac { \lambda \lambda ^ { \prime } + \gamma \gamma ^ { \prime } } { 1 - \gamma \gamma ^ { \prime } } \right | < 1. \quad \text{(30)} \\ \text{dge mode $k^{*}_{1}=\pi$ or $0$ does not cause any}$$

All the elements consisting the parameter space of SDQPT-1 are investigated above. To sum up, they include (see Tab. I)

∣ ∣ The critical edge mode k ∗ 1 = π or 0 does not cause any SDQPT because d ′ k ∗ 1 = 0, but the critical body mode k ∗ 2 = arccos -λλ ′ + γγ ′ 1 -γγ ′ ( λ ′ = 1) or k ∗ 2 = arccos λλ ′ + γγ ′ 1 -γγ ′ ( λ ′ = -1) can cause SDQPT-1.

(i) γγ ′ = 0, ( | λ | -1)( | λ ′ | -1) &lt; 0;

/negationslash

(iii) γ = 0, γ ′ = 0, | λ &lt; | 1, λ ′ = λ ;

/negationslash

/negationslash

(ii) γ = 0, γ ′ = 0, | λ ′ | &lt; 1;

/negationslash

/negationslash

(iv) γγ ′ ( γγ ′ -1) = 0, λ = λ ′ , | λ ′ | = 1, and | λλ ′ + γγ ′ 1 -γγ ′ | &lt; 1. Now, we discuss the properties of SDQPT-1. Consider γ = 0 5 . → γ ′ = 2 and λ = 0 5 . → λ ′ = ± 1 5 as an exam-. ple, we plot the lines of Fisher zeros, the time evolution of rate function r t ( ), the winding number ν t ( ) and the trajectory of vector in Fig. 2. It can be seen that the line of Fisher zeros cuts Im( z ) axis at the critical mode [see Figs. 2(a1) and (a2)], giving rise to nonanalytic behavior (cusp singularity) of the rate function of the Loschmidt echo [see the blue curves in Figs. 2(b1) and (b2)]. The trajectory of vector /vector R k forms a closed loop anticounterclockwise [see Fig. 2(c1)] or counterclockwise for λ ′ &lt; 0 [see Fig. 2(c2)], as k varies from 0 to π . The trajectory does not encircle the origin until it firstly crosses it at the first critical time [see Figs. 2(c1) and (c2)], and then the winding number jumps up into 1 [see the red curves in Fig. 2(b1)] or jumps down into -1 [see the red curves in Fig. 2(b2)]. After the first critical time, the trajectory will encircle the origin until the second critical time, when the trajectory secondly crosses the origin and the winding number jumps up into 2 or jumps down into -2, and so on. It should be noted that SDQPT-1 is just the result discovered by the seminal paper [21].

/negationslash

## B. SDQPT-2

SDQPT-2 is caused by two critical body modes. These two critical body modes can only be k ∗ + = arccos(cos k + ) and k ∗ -= arccos(cos k -) derived from Eq. (25). The presence of them requires ( | cos k + | -1)( | cos k -| -1) &gt; 0 and | cos k + cos k -| &lt; 1. ( | cos k + | -1)( | cos k -| -1) &gt; 0 implies that

$$( | \lambda | - 1 ) ( | \lambda ^ { \prime } | - 1 ) > 0, \quad \ \ ( 3 1 ) \quad \stackrel { \overline { \nu } _ { \text{i.e} } } { \text{$\overline{tuc}$} }$$

/negationslash i.e., the transverse field is quenched within the same phases. | cos k + cos k -| &lt; 1 is equivalent to | γγ ′ + λλ ′ 1 -γγ ′ | &lt; 1 [Eq. (30)]. In summary, the parameter space of SDQPT2 is given by the conditions γγ ′ ( γγ ′ -1) = 0, ( | λ | -1)( | λ ′ | -1) &gt; 0, ∆ &gt; 0 and | γγ ′ + λλ ′ 1 -γγ ′ | &lt; 1 (see Tab. I). Considering γ = 1 → γ ′ = -0 5 and . λ = 1 1 . → λ ′ = 1 3 . as an example, we can find that the line of Fisher zeros cuts Im( z ) axis twice at two critical modes [see Figs. 3(a)], giving rise to nonanalytic behavior (cusp singularity) of the rate function which implies SDQPT occurring [see the blue curves in Figs. 3(b)]. The trajectory of vector /vector R k forms a closed loop anticounterclockwise as k varies from 0 to π . The trajectory does not encircle the origin until it firstly crosses it at the first critical time t ∗ 1 ≈ 2 07 [see Figs. . 3(c)], and then the winding number jumps up into 1 [see the red curves in Fig. 3(b)]. After the first critical time, the trajectory will encircle the origin until the second critical time t ∗ 2 ≈ 3 95, when it puts . the origin out the loop, resulting in the winding number jumping down to 0. SDQPT-2 has been discussed in Ref. [69, 86], and has been shown to occur with the aid of quantum coherence [85].

## C. SDQPT-3

/negationslash

Notably, two Fisher zeros cutting Im( z ) axis [see Figs. 3(a)] can approach each other until they merge, and then leave Im( ) z axis, by adjusting the external field and anisotropic parameter. The merged Fisher zero means that there is only one solution to Eq. (22), i.e., cos k + = cos k -= -λ + λ ′ 2(1 -γγ ′ ) , which requires ∆ = 0. One critical body mode k ∗ = arccos[ -λ + λ ′ 2(1 -γγ ′ ) ] can be obtained if | λ + λ ′ 2(1 -γγ ′ ) | = | γγ ′ + λλ ′ 1 -γγ ′ | &lt; 1. And it should be noted that ( d k&lt;k ∗ · d ′ k&lt;k ∗ )( d k&gt;k ∗ · d ′ k&gt;k ∗ ) &gt; 0, which implies that the line of Fisher zeros only touches but not cross Im( ) axis z at this critical body mode. The SDQPT caused by this kind of critical body mode is named SDQPT-3, an entirely new SDQPT, and the parameter space of SDQPT-3 is determined by the conditions γγ ′ ( γγ ′ -1) = 0, ∆ = 0 and | γγ ′ + λλ ′ 1 -γγ ′ | &lt; 1 (see Tab. I). Considering γ = 1 → γ ′ = -0 5 . and λ = 1 1 . → λ ′ = (1 -2 γγ ′ ) λ -√ 4 γγ ′ (1 -γγ ′ )(1 -λ 2 ) as an example, we plot Fisher zeros, the rate function, winding number and the trajectory in Fig. 4. The line of Fisher zeros touches Im( z ) axis at the critical mode [see Figs. 4(a)], but the rate function does not show the nonanalytic behavior (cusp singularity) [see the blue curves in Figs. 4(b)]. How to describe SDQPT in this case? The trajectory of vector /vector R k forms a closed loop anticounterclockwise for λ &gt; ′ 0 [see Fig. 2(c)], as k varies from 0 to π . This loop becomes larger and larger as time goes on from t = 0, and its trajectory will cross the origin when the first critical time is reached. After the critical time, the loop leaves the origin rather than encircle it [see Figs. 2(c)]. Mathematically, the argument of the trajectory /vector R k (i.e., the geometric phase of the Loschmidt overlap amplitude) will suddenly changes with π when the trajectory crosses the origin, resulting in singularities of the winding number. On the other hand, since the trajectories before and after the critical time do not encircle the origin, the winding numbers before and after the critical time are both 0. In other words, the appearance of the singularities of the winding number indicates that the trajectory crosses the origin, which means SDQPT occurring.

## D. SDQPT-4

/negationslash

SDQPT-4 is caused by one edge mode. This edge mode requires that one of the solutions of Eq. (22) should be cos k = ± 1, implying | λ | = 1 or | λ ′ | = 1. If the transverse field is quenched to the quantum critical point λ ′ = ± 1, the critical edge mode would not cause SDQPT because d ′ k ∗ = 0 and the initial state does not evolve at all. The transverse field quenched from the quantum critical point λ = ± 1 to others λ ′ = λ is a prerequisite of SDQPT-4. In this case, the energy levels at critical edge mode k ∗ will be degenerate with d k ∗ = 0. However, whether the SDQPT occurs also depends on the anisotropy param-

FIG. 6: (Color online) Lines of Fisher zeros (a), the time evolution of rate function r t ( ) (blue curves in b), the winding number ν t ( ) (red curves in b) and the trajectory of vector /vector R k = ( x , y k k ) (c) for the quenches from λ = 1 to λ ′ = 0 1. . The first and the second critical times are t ∗ 1 ≈ 1 73 and . t ∗ 2 ≈ 2 94. . The other parameters are γ = 1, γ ′ = 0 3 and . n = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-8812175b-817a-4bbe-9f80-0b4fb2e7c94f.jpg)

FIG. 7: (Color online) Lines of Fisher zeros (a), the time evolution of rate function r t ( ) (blue curves in b), the winding number ν t ( ) (red curves in b) and the trajectory of vector /vector R k = ( x , y k k ) (c) for the quenches from λ = 1 to λ ′ = -1. The other parameters are γ = 0 5, . γ ′ = 2 and n = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-ef3840ed-ce00-4559-a31d-ee97d6b0d1cf.jpg)

/negationslash

/negationslash

/negationslash

(c2)]. This half-loop becomes larger and larger as time goes on from t = 0 by moving its end, and its end reaches the origin at the first critical time. After the critical time, the end of half-loop moves to the other side of the origin. In this case, the argument of the trajectory /vector R k (i.e., the geometric phase of the Loschmidt overlap amplitude) will suddenly changes with π , resulting in the jumping of the winding number with half-integer. This result is consistent with the observation in Ref. [89, 90]. The trajectory of vector /vector R k swipes half a circle anticounterclockwise for λ ′ &gt; 0 [see Fig. 5(c1)] and counterclockwise for λ ′ &lt; 0 [see Fig. 5(c2)], then the winding number jumps up [see the red curves in Fig. 5(b1)] or jumps down [see the red curves in Fig. 5(b2)] with half integer at the critical times.

/negationslash

/negationslash

$$| \frac { \gamma \gamma ^ { \prime } + \lambda \lambda ^ { \prime } } { 1 - \gamma \gamma ^ { \prime } } | > \\ \tilde { C o n e i d o }$$

eter γ . If γ = 0, the initial state would be | d -k 〉 = 0 | 〉 and the evolution operator at critical edge mode would be ˆ U k ∗ ( ) = t e -id ′ k ∗ σ t z , so that the initial state can not be flipped not by the evolution operator. Therefore, γ = 0 is also a prerequisite of SDQPT-4. If γ = 0, the ground state (initial state) would continuously become the superposition of | d -k ∗ 〉 = √ 1 / 2( i | 1 〉 + 0 ), but the evolution | 〉 operator would continuously become e -i ( λ ′ ∓ 1)ˆ σ t z , as k approaches into critical edge mode k ∗ . The initial state can be flipped by the evolution operator, so SDQPT-4 occurs. It should be noted that, in addition to this critical edge mode there could be another critical body mode, if γγ ′ ( γγ ′ -1) = 0, ∆ &gt; 0 and | γγ ′ + λλ ′ 1 -γγ ′ | &lt; 1, which would cause SDQPT-5, discussed in the next section. To sum up, the parameter spaces of SDQPT-4 are (see Tab. I) (i) γγ ′ = 1, | λ | = 1, λ ′ = λ ; (ii) γ = 0, γ ′ = 0, | λ | = 1, λ ′ = λ ; (iii) γγ ′ ( γγ ′ -1) = 0, | λ | = 1, λ = λ ′ , ∆ &gt; 0, γγ ′ + λλ ′ ′ &gt; 1.

/negationslash

/negationslash

/negationslash

Consider γ = 0 5 . → γ ′ = 2 and λ = 1 → λ ′ = ± 1 5 . as an example, we plot the lines of Fisher zeros, the time evolution of rate function r t ( ), the winding number ν t ( ) and the trajectory of vector in Fig. 5. It can be seen that the Fisher zeros coalesce into a continuous curve that touches Im( z ) axis at a critical edge mode k ∗ [see Fig. 5(a1) and (a2)], giving rise to nonanalytic behavior (cusp singularity) of the rate function of the Loschmidt echo which implies SDQPT occurring [see the blue curves in Fig. 5(b1) and (b2)]. It should be noted that the trajectory of /vector R k only swipes half a circle [see Fig. 5(c1) and

## E. SDQPT-5

/negationslash

/negationslash

SDQPT-5 is caused by one critical edge mode and one critical body mode. These two critical modes are derived from Eq. (28) for | γγ ′ + λλ ′ 1 -γγ ′ | &lt; 1. In summary, the parameter space satisfies the following conditions γγ ′ ( γγ ′ -1) = 0, | λ | = 1, λ = λ ′ , ∆ &gt; 0 and | γγ ′ + λλ ′ 1 -γγ ′ | &lt; 1 (see Tab. I). Consider γ = 1 → γ ′ = 0 3 . and λ = 1 → λ ′ = 0 1 as example, . we plot SDQPT5 in Fig. 6. From Fig. 6 we can see that the line of Fisher zeros cuts Im( z ) axis at the critical mode k ∗ 2 and touches it at the critical mode k ∗ 1 [see Figs. 6(a)], giving rise to nonanalytic behavior of the rate function of the

FIG. 8: (Color online) The dynamical phase diagrams. Note that the purple and green lines represent the same DQPTs as the purple and green areas. We have been darkened the colors to emphasize the boundaries.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00597v2-figure-dba63b9a-2d66-455a-8d04-6dd736c87c17.jpg)

Loschmidt echo [see the blue curves in Figs. 6(b)]. As k varies from 0 to π , the coordinate of /vector R k moves first clockwise and then counterclockwise, forms an open trajectory [see Figs. 6(b)]. The end of trajectory crosses the origin at the first critical time, and then the winding number jumps down into -1 / 2 [see the red curves in Fig. 6(b)]. After the first critical time, the trajectory will half-encircle the origin until the second critical time, when the trajectory secondly cross the origin counterclockwise. Notably, it is not the end of trajectory that crosses the origin secondly. As a result, the winding number will jumps up from -1 / 2 to 1 / 2 [see the red curves in Fig. 6(b)]. SDQPT-5 is also an entirely new SDQPT that is not discovered before.

/negationslash

/negationslash

## F. SDQPT-6

It is worth noting that, if the transverse field is quenched from one critical point λ = ± 1 to another one λ ′ = ∓ 1 and the anisotropic parameters satisfy γγ ′ = 1, Eq. (22) is an equality independent of k . In other words, all the Fisher zeros are located on Im( z ) axis [see Fig. 7(a)], implying a very anomalous DQPT that the critical time is no longer periodic. The trajectory of /vector R k is always along the imaginary axis [see Fig. 7(c)], in this case, the winding number is always 0 [see Fig. 7(b)]. This anomalous SDQPT is also an entirely new SDQPT, and is referred to in this paper as SDQPT-6.

In summary, we demonstrate that there are only six possible SDQPTs in the one-dimensional XY model, three of which were previously discovered, and the others are entirely new SDQPTs that have not been discovered before. Furthermore, the associated conditions for different SDQPTs are shown in Table I. It is important to emphasize that this classification of SDQPT and the associate criteria are universal, holding for any other twoband models in one and two dimensions, including SSH model, Kitaev-chain, Haldane model, p + ip superconductor, etc. To be specific, there are two types of the critical mode: the critical body mode ( k ∗ = 0 or k ∗ = π ) and the critical edge mode ( k ∗ = 0 or k ∗ = π ). If only one critical body mode is generated by controlling the parameters, and the line of Fisher zeros cuts Im(z) axis at this mode, SDQPT-1 will take place in which the winding number jumps up or down with an integer step. On the contrary, if two critical body modes are generated SDQPT-2 will occur where the winding number jumps up and down alternately with the integer step. If these two critical body modes merge into one by adjusting the parameters, SDQPT-3 occurs with the singularity of the winding number. A single critical edge mode generates SDQPT-4 where the winding number jumps up or down with the half-integer step. A critical body mode and a critical edge mode commonly cause SDQPT-5 whose topological property is the alternative jumps of the winding number with integer and half-integer steps. If all the modes are critical, SDQPT-6 will be caused. The exact values of the initial and final parameters λ , γ , λ ′ , γ ′ in given classification criteria do not modify qualitatively the behavior of the corresponding types of SDQPT as a sign of some kind of universality.

## V. DYNAMICAL PHASE DIAGRAM

In one-dimensional XY model, the ground state phase diagram has already been well known [14, 96-98], but the dynamical phase diagram is still unclear. Now we give the dynamical phase diagram according to the classification of SDQPT and the associated conditions (see Table I). Although the conditions under which SDQPTs occur are listed in Table 1, the dynamical phase diagram allows for a more clear and intuitive understanding.

/negationslash

/negationslash

/negationslash

We should consider six cases: γγ ′ &lt; 0, γγ ′ = 0 including γ = 0 → γ ′ = 0 and γ = 0 → γ ′ = 0, 0 &lt; γγ ′ &lt; 1, γγ ′ = 1 and γγ ′ &gt; 1 to investigate the dynamical phase diagrams, and plot them in Fig. 8. It can be seen that for γγ ′ = 0, i.e., the anisotropic parameters ( γ and γ ′ ) before and after quench are both not at the critical point γ c = 0, SDQPT-1 always happen if the transverse field is quenched across quantum critical point | λ c | = 1 [see the purple areas in Fig. 8(a)-(d)]. However, this does not mean SDQPT must depend on the quantum phase transition because SDQPT-2 (the yellow areas) and SDQPT-3 (the red curves) can happen when the transverse field is quenched within the same phases. To be specific, SDQPT-2 and SDQPT-3 can happen at the ferromagnetic phase if 0 &lt; γγ ′ &lt; 1, but can happen at the paramagnetic phase if γγ ′ &gt; 1 [see Fig. 8(b) and (c)]. Beyond that, SDQPT-2 and SDQPT-3 can happen at both ferromagnetic and paramagnetic phases, if the anisotropic parameter is quenched across the critical point γ c = 0 [see Fig. 8(d)]. It is worth noting that for a given suitable field λ = λ ′ , only the quench of the anisotropic parameter across the critical point γ c = 0 can still cause SDQPT-2 and SDQPT-3, and even SDQPT-1.

/negationslash

At the critical boundaries, the situations are more complicated. If the transverse field is quenched from the quantum critical point | λ c | = 1, SDQPT can always occur, either SDQPT-4 or SDQPT-5 (see the blue and cyan lines in Fig. 8), depending on λ ′ and γγ ′ , as shown in Table I. On the other hand, if the transverse field is quenched to the quantum critical point | λ ′ | = 1, only SDQPT-1 can occur, but not always, depending on λ ′ and γγ ′ (see the purple lines in Fig. 8). If the anisotropic parameter is quenched from the critical point γ c = 0, SDQPT-1 can always occur when the transverse field is quench to the ferromagnetic phase | λ ′ | &lt; 1, independent of the pre-quenched field λ [see the purple areas in Fig. 8(e)]; on the contrary, if the anisotropic parameter is quenched to the critical point γ c = 0, SDQPT-1 can always occur when the transverse field is quench from the ferromagnetically ordered phase | λ ′ | &lt; 1, independent of the post-quenched field λ ′ = λ [see the purple areas in Fig. 8(f)]. These results furthermore demonstrate that SDQPT is independent of quantum phase transition, which was first pointed out in Ref. [86].

In addition to the Ising transition and anisotropic transition, there is another type of transition, namely the commensurate-incommensurate transition in onedimensional XY model [14, 96-98]. This transition causes

LDQPT (i.e., long-time relaxation dynamical phase transition). To be specific, two distinct power-law relaxation behaviors following a sudden quench can be observed, depending upon whether the quenched Hamiltonian lies in the commensurate or the incommensurate phase [14]. A natural question arises, does the commensurateincommensurate transition has any effect on SDQPT? The answer is no, because SDQPT can occur regardless of whether the commensurate-incommensurate transition occurs or not. This result demonstrate that SDQPT is of fundamentally different origin from LDQPT.

## VI. CONCLUSIONS

In this paper, we have comprehensively investigated SDQPTs by focusing the Fisher zeros, the rate function of Loschmidt echo and the winding number after quenching the strength of transverse field and the anisotropic parameter in the one-dimensional XY model. Accordingly, SDQPTs have been classified into six types based on the attribution and number of critical modes, and each type has its own topological properties. Specifically, SDQPT1 is caused by a single critical body mode at which the line of Fisher zeros cuts the imagery axis, and its topological parameter order jumps up or down with an integer unit at the critical times. SDQPT-2 happens when two critical body modes exist. In SDQPT-2, the winding number jumps up at one set of critical times determined by one critical mode, but jumps down at another set of critical times determined by another critical mode, the jump heights are all integer units. By adjusting the external field and anisotropic parameter, two critical body mode above can approach each other until they merge. In this case, SDQPT-3 occurs. SDQPT-3 can not be described by the rate function of the Loschmidt echo, but can be described by the singularities of the winding number. SDQPT-4 is caused by a single critical edge mode, wherein the winding number jumps up or down with a half-integer unit. SDQPT-5 is commonly caused by two critical modes: one critical body mode and one critical edge mode. The critical body mode makes the winding number jump with the integer unit, whereas the critical edge mode makes it jump with the half-integer unit. SDQPT-6 is very anomalous that all the modes locate at the imagery axis. Although SDQPT-1, SDQPT-2 and SDQPT-4 have been previously discovered in other papers, SDQPT-3, SDQPT-5 and SDQPT-6 are entirely new and had not been discussed before. Furthermore, the parameter spaces of these six SDQPTs have been determined and the phase diagrams have been given. It is important to emphasize that this classification of SDQPT and the associate criteria are universal, holding for any other two-band models in one- and two-dimensions (e.g. SSH model, Kitaev-chain, Haldane model, p+ip superconductor, etc.).

In equilibrium, continuous phase transitions exhibit scaling and universality, connecting to renormalization

group theory and the associated fixed points. To what extent the concepts of scaling and universality can be applied to SDQPTs is of fundamental importance for the theoretical understanding of the out-of-equilibrium dynamics of quantum many-body systems. In Ref. [22], Heyl proved that these concepts can be extended to SDQPTs both in one- and two-dimensions Ising models, showing that SDQPTs are critical points associated with unstable fixed points of equilibrium Ising models. Therefore, these SDQPTs obey scaling and universality. It is important to note that the SDQPTs in Ising mode belongs to SDQPT-1. In other words, SDQPT-1 obeys scaling and universality. Investigating to what extent universality and scaling generalize to the other types of SDQPT

- [1] M. Heyl, Dynamical quantum phase transitions: A brief survey, Europhys. Lett. 125 , 26001 (2019).
- [2] M. Heyl, Dynamical quantum phase transitions: a review, Rep. Prog. Phys. 81 , 054001 (2018).
- [3] J. Marino, M. Eckstein, M. S. Foster and A. M. Rey, Dynamical phase transitions in the collisionless prethermal states of isolated quantum systems: theory and experiments, Rep. Prog. Phys. 85 , 116001 (2022).
- [4] K. Levin, A. Fetter, and D. Stamper-Kurn, Ultracold Bosonic and Fermionic Gases, Contemporary Concepts of Condensed Matter Science (Elsevier Science, New York, 2012).
- [5] V. I. Yukalov, Equilibration and thermalization in finite quantum systems, Laser Phys. Lett. 8 , 485 (2011).
- [6] I. Bloch, J. Dalibard, and W. Zwerger, Many-body physics with ultracold gases, Rev. Mod. Phys. 80 , 885 (2008).
- [7] M. Greiner, O. Mandel, T. W. H¨ ansch, and I. Bloch, Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms, Nature (London) 419 , 51 (2002).
- [8] D. Porras and J. I. Cirac, Effective quantum spin systems with trapped ions, Phys. Rev. Lett. 92 , 207901 (2004).
- [9] K. Kim, M.-S. Chang, R. Islam, S. Korenblit, L.-M. Duan, and C. Monroe, Entanglement and tunable spinspin couplings between trapped ions using multiple transverse modes, Phys. Rev. Lett. 103 , 120502 (2009).
- [10] P. Jurcevic, B. P. Lanyon, P. Hauke, C. Hempel, P. Zoller, R. Blatt, and C. F. Roos, Quasiparticle engineering and entanglement propagation in a quantum manybody system, Nature (London) 511 , 202 (2014).
- [11] J. Eisert, M. Friesdorf and C. Gogolin, Quantum manybody systems out of equilibrium, Nature Phys. 11 , 124 (2015).
- [12] A. Sen, S. Nandy, and K. Sengupta, Entanglement generation in periodically driven integrable systems: Dynamical phase transitions and steady state, Phys. Rev. B 94 , 214301 (2016).
- [13] S. Nandy, K. Sengupta, and A. Sen, Periodically driven integrable systems with long-range pair potentials, J. Phys. A: Math. Theor. 51 , 334002 (2018).
- [14] A. A. Makki, S. Bandyopadhyay, S. Maity, and A. Dutta, Dynamical crossover behavior in the relaxation

(especially those independent of quantum phase transition), and whether genuine nonequilibrium fixed points can also appear, which are not related to the equilibrium criticality, is a pertinent task for future work on the theory of SDQPTs.

## Acknowledgement

This work was supported by the National Natural Science Foundation of China (Grant No. 11705099) and the Talent Introduction Project of Dezhou University of China (Grant Nos. 2020xjpy03 and 2019xgrc38).

- of quenched quantum many-body systems, Phys. Rev. B 105 , 054301 (2022).
- [15] Y.-T. Zou and C. Ding, Dynamical relaxation behaviors of a critical quench, Phys. Rev. B 108 , 014303 (2023).
- [16] S. Aditya, S. Samanta, A. Sen, K. Sengupta, and D. Sen, Dynamical relaxation of correlators in periodically driven integrable quantum systems, Phys. Rev. B 105 , 104303 (2022).
- [17] S. Aditya, D. Sen, Dynamical localization and slow thermalization in a class of disorder-free periodically driven one-dimensional interacting systems, SciPost Phys. Core 6 , 083 (2023).
- [18] K. Cao, Y. Hu, P. Tong, and G. Yang, Dynamical relaxation behavior of an extended XY chain with a gapless phase following a quantum quench, Phys. Rev. B 109 , 024303 (2024).
- [19] Y.-H. Huang, Y.-T. Zou, and C. Ding, Dynamical relaxation of a long-range Kitaev chain, Phys. Rev. B 109 , 094309 (2024).
- [20] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, Colloquium: Nonequilibrium dynamics of closed interacting quantum systems, Rev. Mod. Phys. 83 , 863 (2011).
- [21] M. Heyl, A. Polkovnikov, S. Kehrein, Dynamical Quantum Phase Transitions in the Transverse Field Ising Model, Phys. Rev. Lett. 110 , 135704 (2013).
- [22] M. Heyl, Scaling and universality at dynamical quantum phase transitions, Phys. Rev. Lett. 115 , 140602 (2015).
- [23] D. Trapin, J. C. Halimeh, and M. Heyl, Unconventional critical exponents at dynamical quantum phase transitions in a random ising chain, Phys. Rev. B 104 , 115159 (2021).
- [24] U. Bhattacharya, S. Bandyopadhyay, and A. Dutta, Mixed state dynamical quantum phase transitions, Phys. Rev. B 96 , 180303(R) (2017).
- [25] M. Heyl and J. C. Budich, Dynamical topological quantum phase transitions for mixed states, Phys. Rev. B 96 , 180304(R) (2017).
- [26] H. Lang, Y. Chen, Q. Hong, and H. Fan, Dynamical quantum phase transition for mixed states in open systems, Phys. Rev. B 98 , 134310 (2018).
- [27] B. Mera, C. Vlachou, N. Paunkovi´, V. R. Vieira, and O. c Viyuela, Dynamical phase transitions at finite temperature from fidelity and interferometric Loschmidt echo

induced metrics, Phys. Rev. B 97 , 094110 (2018).

- [28] S. Sharma, A. Russomanno, G. E. Santoro, and A. Dutta, Loschmidt echo and dynamical fidelity in periodically driven quantum systems, Europhys. Lett. 106 , 67003 (2014).
- [29] K. Yang, L. Zhou, W. Ma, X. Kong, P. Wang, X. Qin, X. Rong, Y. Wang, F. Shi, J. Gong, and J. Du, Floquet dynamical quantum phase transitions, Phys. Rev. B 100 , 085308 (2019).
- [30] S. Zamani, R. Jafari, and A. Langari, Floquet dynamical quantum phase transition in the extended XY model: Nonadiabatic to adiabatic topological transition, Phys. Rev. B 102 , 144306 (2020).
- [31] L. Zhou and Q. Du, Floquet dynamical quantum phase transitions in periodically quenched systems, J. Phys.: Condens. Matter 33 , 345403 (2021).
- [32] R. Jafari and A. Akbari, Floquet dynamical phase transition and entanglement spectrum, Phys. Rev. A 103 , 012204 (2021).
- [33] R. Hamazaki, Exceptional dynamical quantum phase transitions in periodically driven systems, Nature Commun. 12 , 5108 (2021).
- [34] S. Zamani, R. Jafari, and A. Langari, Out-of-time-order correlations and Floquet dynamical quantum phase transition, Phys. Rev. B 105 , 094304 (2022).
- [35] L.-N. Luan, M.-Y. Zhang, and L. Wang, Floquet dynamical quantum phase transitions of the XY spinchain under periodic quenching, Physica A 604 , 127866 (2022).
- [36] R. Jafari, A. Akbari, U. Mishra, and H. Johannesson, Floquet dynamical quantum phase transitions under synchronized periodic driving, Phys. Rev. B 105 , 094311 (2022).
- [37] A. Kosior and K. Sacha, Dynamical quantum phase transitions in discrete time crystals, Phys. Rev. A 97 , 053621 (2018).
- [38] A. Kosior, A. Syrwid, and K. Sacha, Dynamical quantum phase transitions in systems with broken continuous time and space translation symmetries, Phys. Rev. A 98 , 023612 (2018).
- [39] J. Naji, R. Jafari, L. Zhou, and A. Langari, Engineering Floquet dynamical quantum phase transitions, Phys. Rev. B 106 , 094314 (2022).
- [40] S. Sharma, U. Divakaran, A. Polkovnikov, and A. Dutta, Slow quenches in a quantum Ising chain: Dynamical phase transitions and topology, Phys. Rev. B 93 , 144306 (2016).
- [41] N. Sedlmayr, P. Jaeger, M. Maiti, and J. Sirker, Bulkboundary correspondence for dynamical phase transitions in one-dimensional topological insulators and superconductors, Phys. Rev. B 97 , 064304 (2018).
- [42] T. Mas/suppressowski l and N. Sedlmayr, Quasiperiodic dynamical quantum phase transitions in multiband topological insulators and connections with entanglement entropy and fidelity susceptibility, Phys. Rev. B 101 , 014301 (2020).
- [43] S. Zamani, J. Naji, R. Jafari and A. Langari, Scaling and universality at ramped quench dynamical quantum phase transitions, J. Phys.: Condens. Matter 36 355401 (2024).
- [44] M. Schmitt and S. Kehrein, Dynamical quantum phase transitions in the Kitaev honeycomb model, Phys. Rev. B 92 , 075114 (2015).
- [45] L. Rossi and F. Dolcini, Nonlinear current and dy-
- namical quantum phase transitions in the flux-quenched Su-Schrieffer-Heeger model, Phys. Rev. B 106 , 045410 (2022).
- [46] L. Zhou, Q.-h. Wang, H. Wang, and J. Gong, Dynamical quantum phase transitions in non-Hermitian lattices, Phys. Rev. A 98 , 022129 (2018).
- [47] L. Zhou and Q. Du, Non-Hermitian topological phases and dynamical quantum phase transitions: A generic connection, New J. Phys. 23 , 063041 (2021).
- [48] D. Mondal and T. Nag, Anomaly in the dynamical quantum phase transition in a non-Hermitian system with extended gapless phases, Phys. Rev. B 106 , 054308 (2022)
- [49] D. Mondal and T. Nag, Finite temperature dynamical quantum phase transition in a non-Hermitian system, Phys. Rev. B 107 , 184311 (2023).
- [50] D. Mondal and T. Nag, Persistent anomaly in dynamical quantum phase transition in long-range non-Hermitian p-wave Kitaev chain. Eur. Phys. J. B 97, 59 (2024).
- [51] Y. Jing, J.-J. Dong, Y.-Y. Zhang, and Z.-X. Hu, Biorthogonal Dynamical Quantum Phase Transitions in Non-Hermitian Systems, Phys. Rev. Lett. 132 , 220402 (2024).
- [52] C. Karrasch and D. Schuricht, Dynamical phase transitions after quenches in nonintegrable models, Phys. Rev. B 87 , 195104 (2013).
- [53] S. Sharma, S. Suzuki, and A. Dutta, Quenches and dynamical phase transitions in a nonintegrable quantum Ising model, Phys. Rev. B 92 , 104306 (2015).
- [54] H. Cheraghi, S. Mahdavifar, Dynamical Quantum Phase Transitions in the 1D Nonintegrable Spin-1/2 Transverse Field XZZ Model, Ann. Phys. (Berlin) 533 , 2000542 (2021).
- [55] D. M. Kennes, D. Schuricht, and C. Karrasch, Controlling dynamical quantum phase transitions, Phys. Rev. B 97 , 184302 (2018).
- [56] F. Andraschko and J. Sirker, Dynamical quantum phase transitions and the Loschmidt echo: A transfer matrix approach, Phys. Rev. B 89 , 125120 (2014).
- [57] J. C. Halimeh, M. V. Damme, L. Guo, J. Lang, and P. Hauke, Dynamical phase transitions in quantum spin models with antiferromagnetic long-range interactions, Phys. Rev. B 104 , 115133 (2021).
- [58] I. Hagym´ asi, C. Hubig, O. Legeza, and U. Schollw¨ck, ¨ o Dynamical Topological Quantum Phase Transitions in Nonintegrable Models, Phys. Rev. Lett. 122 , 250601 (2019).
- [59] K. Cao, W. Li, M. Zhong, and P. Tong, Influence of weak disorder on the dynamical quantum phase transitions in the anisotropic xy chain, Phys. Rev. B 102 , 014207 (2020).
- [60] O. N. Kuliashov, A. A. Markov, and A. N. Rubtsov, Dynamical quantum phase transition without an order parameter, Phys. Rev. B 107 , 094304 (2023).
- [61] R. Jafari, A. Langari, S. Eggert, and H. Johannesson, Dynamical quantum phase transitions following a noisy quench, Phys. Rev. B 109 , L180303 (2024).
- [62] P. Jurcevic, H. Shen, P. Hauke, C. Maier, T. Brydges, C. Hempel, B. P. Lanyon, M. Heyl, R. Blatt, and C. F. Roos, Direct Observation of Dynamical Quantum Phase Transitions in an Interacting Many-Body System, Phys. Rev. Lett. 119 , 080501 (2017).
- [63] N. Fl¨schner, D. Vogel, M. Tarnowski, B. S. Rem, D.-S. a

- L¨ uhmann, M. Heyl, J. C. Budich, L. Mathey, K. Sengstock, and C. Weitenberg, Observation of dynamical vortices after quenches in a system with topology, Nat. Phys. 14 , 265 (2018).
- [64] X.-Y. Guo, C. Yang, Y. Zeng, Y. Peng, H.-K. Li, H. Deng, Y.-R. Jin, S. Chen, D. Zheng, and H. Fan, Observation of a Dynamical Quantum Phase Transition by a Superconducting Qubit Simulation, Phys. Rev. Applied 11 , 044080 (2019).
- [65] T. Tian, Y. Ke, L. Zhang, S. Lin, Z. Shi, P. Huang, C. Lee, and J. Du, Observation of dynamical phase transitions in a topological nanomechanical system, Phys. Rev. B 100 , 024310 (2019).
- [66] K. Wang, X. Qiu, L. Xiao, X. Zhan, Z. Bian, W. Yi, and P. Xue, Simulating Dynamic Quantum Phase Transitions in Photonic Quantum Walks, Phys. Rev. Lett. 122 , 020501 (2019).
- [67] X.-Y. Xu, Q.-Q. Wang, M. Heyl, J. C. Budich, W.-W. Pan, Z. Chen, M. Jan, K. Sun, J.-S. Xu, Y.-J. Han, C.-F. Li, and G.-C. Guo, Measuring a dynamical topological order parameter in quantum walks, Light Sci. Appl. 9 , 7 (2020).
- [68] M. Heyl, Dynamical Quantum Phase Transitions in Systems with Broken-Symmetry Phases, Phys. Rev. Lett. 113 , 205701 (2014).
- [69] J. C. Budich and M. Heyl, Dynamical topological order parameters far from equilibrium, Phys. Rev. B 93 , 085416 (2016).
- [70] U. Bhattacharya and A. Dutta, Emergent topology and dynamical quantum phase transitions in twodimensional closed quantum systems, Phys. Rev. B 96 , 014302 (2017).
- [71] A. Kosior and M. Heyl, Vortex loop dynamics and dynamical quantum phase transitions in three-dimensional fermion matter, Phys. Rev. B 109 , L140303 (2024).
- [72] W. C. Yu, P. D. Sacramento, Y. C. Li, and H.-Q. Lin, Correlations and dynamical quantum phase transitions in an interacting topological insulator, Phys. Rev. B 104 , 085104 (2021).
- [73] J. C. Halimeh, D. Trapin, M. Van Damme, and M. Heyl, Local measures of dynamical quantum phase transitions, Phys. Rev. B 104 , 075130 (2021).
- [74] C.-X. Li, S. Yang, J.-B. Xu, and H.-Q. Lin, Exploring dynamical quantum phase transitions in a spin model with deconfined critical point via the quantum steering ellipsoid, Phys. Rev. B 107 , 085130 (2023).
- [75] L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Entanglement in many-body systems, Rev. Mod. Phys. 80 , 517 (2008).
- [76] N. Laflorencie, Quantum entanglement in condensed matter systems, Phys. Rep. 646 , 1 (2016).
- [77] G. Torlai, L. Tagliacozzo and G. De Chiara, Dynamics of the entanglement spectrum in spin chains, J. Stat. Mech. (2014) P06001.
- [78] E. Canovi, E. Ercolessi, P. Naldesi, L. Taddia, and D. Vodola, Dynamics of entanglement entropy and entanglement spectrum crossing a quantum phase transition, Phys. Rev. B 89 , 104303 (2014).
- [79] M. Schmitt, M. Heyl, Quantum dynamics in transversefield Ising models from classical networks, SciPost Phys. 4 , 013 (2018).
- [80] J. Surace, L. Tagliacozzo, and E. Tonni, Operator content of entanglement spectra in the transverse field Ising chain after global quenches, Phys. Rev. B 101 ,
- 241107(R) (2020).
- [81] S. De Nicola, A. A. Michailidis, and M. Serbyn, Entanglement View of Dynamical Quantum Phase Transitions, Phys. Rev. Lett. 126 , 040602 (2021).
- [82] R. Modak and D. Rakshit, Many-body dynamical phase transition in a quasiperiodic potential, Phys. Rev. B 103 , 224310 (2021).
- [83] S. De Nicola, A. A. Michailidis, and M. Serbyn, Entanglement and precession in two-dimensional dynamical quantum phase transitions, Phys. Rev. B 105 , 165149 (2022).
- [84] B.-M. Xu and C.-Q. Wang, The singularities of the rate function of quantum coherent work in one-dimensional transverse field Ising model, New J. Phys. 25 073020 (2023).
- [85] B.-M. Xu, Quantum coherence assisted dynamical phase transition, arXiv:2305.08400 (2023).
- [86] S. Vajna and B. D´ra, Disentangling dynamical phase o transitions from equilibrium phase transitions, Phys. Rev. B 89 , 161105(R) (2014).
- [87] S. Bhattacharjee and A. Dutta, Dynamical quantum phase transitions in extended transverse Ising models, Phys. Rev. B 97 , 134306 (2018).
- [88] X. Qiu, T.-S. Deng, G.-C. Guo, and W. Yi, Dynamical topological invariants and reduced rate functions for dynamical quantum phase transitions in two dimensions, Phys. Rev. A 98 , 021601(R) (2018).
- [89] C. Y. Wong, H. Cheraghi, and W. C. Yu, Quantum spin fluctuations in dynamical quantum phase transitions, Phys. Rev. B 108 , 064305 (2023).
- [90] C. Ding, Dynamical quantum phase transition from a critical quantum quench, Phys. Rev. B 102 , 060409(R) (2020).
- [91] C.-K. Chiu, J. C.-Y. Teo, A. P. Schnyder, and S. Ryu, Classification of topological quantum matter with symmetries, Rev. Mod. Phys. 88 , 035005 (2016).
- [92] M. Atala, M. Aidelsburger, J. T. Barreiro, D. Abanin, T. Kitagawa, E. Demler, and I. Bloch, Direct measurement of the Zak phase in topological Bloch bands, Nature Phys. 9 , 795 (2013).
- [93] G. Jotzu, M. Messer, R. Desbuquois, M. Lebrat, T. Uehlinger, D. Greif, and T. Esslinger, Experimental realization of the topological Haldane model with ultracold fermions, Nature (London) 515 , 237 (2014).
- [94] V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, and L. P. Kouwenhoven, Signatures of Majorana Fermions in Hybrid SuperconductorSemiconductor Nanowire Devices, Science 336 , 1003 (2012).
- [95] S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. H. MacDonald, B. A. Bernevig, and A. Yazdani, Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor, Science 346 , 602 (2014).
- [96] A. Dutta, G. Aeppli, B. K. Chakrabarti, U. Divakaran, T. F. Rosenbaum, and D. Sen, Quantum Phase Transitions in Transverse Field Spin Models: From Statistical Physics to Quantum Information (Cambridge University Press 2015).
- [97] V. Mukherjee, U. Divakaran, A. Dutta, and D. Sen, Quenching dynamics of a quantum XY spin1 2 chain in a transverse field, Phys. Rev. B 76, 174303 (2007).
- [98] J. E. Bunder and R. H. McKenzie, Effect of disorder on quantum phase transitions in anisotropic XY spin

chains in a transverse field, Phys. Rev. B 60 , 344 (1999).

- [99] M. E. Fisher, in Boulder Lectures in Theoretical Physics (University of Colorado, Boulder, 1965), Vol. 7.
- [100] S. Vajna and B. D´ra, Topological classification of dyo

namical phase transitions, Phys. Rev. B 91 , 155127 (2015).
Bao-Ming Xu Xiu-Ying Zhang
Comprehensive Study of Dynamical Phase Transitions in One-Dimensional XY Model
Short-time transient dynamical phase transition (SDQPT), serving as the
theoretical framework on far-from-equilibrium physics of quantum many-body
systems, has been experimentally observed recently. However, the existence of
different types of SDQPTs and the associated conditions are yet to be elusive.
In this paper, we comprehensively classify SDQPTs by considering the
one-dimensional transverse field quantum XY model initialized in the ground
state. We find that there are only six types of SDQPTs, and three of which are
completely new that had not yet been identified previously. The associated
conditions for different SDQPTs are discussed in detail, and based on this, the
dynamical phase diagrams are given. Our analytical method used to classify
SDQPTs in one-dimensional transverse field quantum XY model can be applied to
any other two-band models in one- and two-spatial dimensions.