Exceptional Dynamical Quantum Phase Transitions in Periodically Driven Systems

Phase transition landau2013statistical ; sachdev2007quantum is one of the most fundamental collective phenomena in macroscopic systems. Recent experiments on artificial quantum many-body systems motivate researchers to understand phases and their transitions in systems out of equilibrium. Various nonequilibrium phases are proposed including e.g., many-body localized phases Schreiber15 ; Smith16 , Floquet topological phases lindner2011floquet ; rechtsman2013photonic , and discrete time crystals Choi17 ; Zhang17 ; Moessner17 .

Recently, dynamical quantum phase transitions (DQPTs) particularly gather great attention as a nonequilibrium counterpart of equilibrium phase transition, which occurs for transient times of quantum relaxation Heyl13 ; Heyl18 . Defined as the singularity of the so-called dynamical free energy (especially at critical times), which is calculated from the overlap between the time-evolved and reference states, the DQPT has been actively studied theoretically Karrasch13 ; Heyl14 ; Andraschko14 ; Canovi14 ; Heyl15 ; Sharma15 ; Budich16 ; Sharma16 ; Halimeh17 ; Zauner-Stauber17 ; Zunkovic18 ; Bhattacharyya20 ; Bandyopadhyay20 ; Halimeh20 and experimentally Jurcevic17 ; Flaschner18 .

Despite extensive studies, the nature of DQPTs is yet to be elusive. One of the important problems is what mechanism leads to DQPTs. Several studies find that some DQPTs are associated with equilibrium/steady-state phase transition Heyl14 ; Zunkovic18 . On the other hand, DQPTs without such relations may also exist Andraschko14 ; Halimeh17 , which indicates that DQPTs can be caused by an unconventional mechanism unique to finite-time (high-frequency) regime of quantum relaxation. Another open problem is universality and criticality of DQPTs. Although typical DQPTs are accompanied by cusps of dynamical free energies Heyl13 ; Heyl18 , several works report DQPTs with different types of singularities Heyl15 ; Bhattacharyya20 . However, a clear understanding of universality and criticality of DQPTs is far from complete.

In this work, we find universal DQPTs in periodically driven unitary dynamics caused by the spontaneous “antiunitary symmetry (AUS)” breaking. While spontaneous symmetry breaking is a fundamental mechanism for conventional phase transitions, several distinct features appear in our results. First, the AUS breaking in our model occurs uniquely at finite times and cannot be captured by conventional equilibrium or steady-state phases. Second, the AUS appears as a symmetry of a hidden nonunitary transfer operator, which is obtained by switching the role of space and time. Consequently, the universality and criticality found in the unitary dynamics are characterized by those of the exceptional point, which recently gathers great attention in non-Hermitian physics ElGanainy18 ; Miri19 ; thus we call the transition the exceptional DQPT. To demonstrate our discovery, we particularly use a stroboscopic chaotic Ising chain and show that the derivative of dynamical free energy defined at finite times can diverge through changes of a parameter (Fig. 1a,b). Using the recently developed technique called the spacetime duality Akila16 ; Bertini18 ; Bertini19 ; Bertini19EC and determining the hidden nonunitary operator, we discuss several properties of the exceptional DQPT besides the divergence of the dynamical free energy (Fig. 1c). For example, instead of the long-range order associated with conventional symmetry breaking, we show that the generalized correlation function has the divergent correlation length at transition and exhibits oscillatory long-range order after antiunitary symmetry breaking. Finally, we demonstrate that the signatures of the exceptional DQPTs are observed through quasi-local observables that are accessible by state-of-the-art experiments Zhang17 ; Tan21 . Notably, we argue that the signature of the exceptional DQPTs is easier to observe than that of the normal DQPTs because of their strong singularity. Our results make an important step toward understanding the nature of phase transitions occurring in a short-time regime, which goes beyond conventional phase transitions since time serves as another crucial parameter here, with their hidden connection to nonunitary physics.

Results

Stroboscopic Ising chains and dynamical free energy. To demonstrate our finding, we introduce a one-dimensional quantum stroboscopic spin model Prosen02 ; Akila16 ; Bertini18 composed of Ising interaction and subsequent global rotation. This model is a prototypical model for quantum chaotic dynamics and can be realized in experiments of e.g., trapped ions Zhang17 . Its unitary time evolution for a single step can be written as

where we impose a periodic boundary condition.

Let us consider a time-evolved state $U^{T}\ket{\psi_{i}}$ after $T$ steps from an initial state $\ket{\psi_{i}}$. To characterize this nonequilibrium state, we focus on the overlap with another state $\ket{\psi_{f}}$, i.e., $\braket{\psi_{f}}{U^{T}}{\psi_{i}}$. The logarithm of the absolute value of this overlap per system size, $F_{L,T}$, is dubbed as the dynamical free energy density Heyl18 . We here consider three types of dynamical free energy density. The first one is to take $\ket{\psi_{i}}=\ket{\psi_{f}}=\ket{\psi}$ and average the overlap over $\ket{\psi}$ randomly taken from the unitary Haar measure before taking the absolute value and the logarithm. Then, the (modified) dynamical free energy density reads

We note that $F_{L,T}^{\mathrm{Tr}}$ is the logarithm of the two-point spectral measure through $|\mathrm{Tr}[U^{T}]|^{2}=\sum_{a,b}e^{iT(z_{a}-z_{b})}$, where $e^{iz_{a}}$ are the eigenvalues for $U$. Since the appearance of trace simplifies the discussion, we mainly use this quantity to show our results.

The second one is to take $\ket{\psi_{i}}=\bigotimes_{j=1}^{L}\ket{\uparrow_{j}}$ and $\ket{\psi_{f}}=\bigotimes_{j=1}^{L}\ket{\downarrow_{j}}$, where $\ket{\uparrow_{j}}$/$\ket{\downarrow_{j}}$ is the eigenstate of $\sigma_{j}^{z}$ with an eigenvalue $+1$/$-1$. In this case, we have

The third one is to take $\ket{\psi_{i}}=\ket{\psi_{f}}=\bigotimes_{j=1}^{L}\ket{\uparrow_{j}}$, leading to $F^{\uparrow\uparrow}_{L,T}(b,J,h)=-\frac{1}{L}\log|\braket{\uparrow\cdots\uparrow}{U^{T}}{\uparrow\cdots\uparrow}|$.

The derivative of $F_{L,T}$ gives the (imaginary part of) so-called generalized expectation values. For example, we have

where $\tilde{\rho}=U^{T}/\mathrm{Tr}[U^{T}]$ and we have used translation invariance. Importantly, the dynamical free energy density and the generalized expectation values can be in principle measured with an interferometric experiment Canovi14 ; Heyl18 .

We seek for singularities of $F_{\infty,T}$ when some continuous parameter is varied. In Ref. Heyl13 , $F_{\infty,T}$ exhibits singularity at critical times for continuous-time models. Since $T$ is discrete in our model, instead of changing $T$, we consider continuously changing other parameters (such as $b$) for fixed $T$.

Dynamical phases and their transitions. As a prime example that highlights our discovery, we show in Fig. 2 the (real-part of) dynamical free energy density $F^{\mathrm{Tr}}_{\infty,T(=6)}$ and $\mathrm{Im}[\braket{\sigma_{1}^{x}}_{\mathrm{gexp}}]$ as a function of the rotation angle $b$ for $J=-\pi/4$ and $h=3.0$ (see Supplementary Note 1 for the data with other parameters and initial/final states). This is calculated from the eigenvalue with the largest modulus of the space-time dual operator, as detailed later. We find different singular behaviors for $F_{\infty,T}^{\mathrm{Tr}}$, signaling distinct DQPTs at critical parameters. Many cusps of $F^{\mathrm{Tr}}_{\infty,T}$ with varying $b$ are are analogous to (continuous time) DQPTs studied previously, where $\mathrm{Im}[\braket{\sigma_{1}^{x}}_{\mathrm{gexp}}]$ exhibits a finite jump.

Notably, we find a distinct singularity at $b=b_{c}\simeq 0.0257$ for $T=6$, where the derivative diverges as $\frac{\mathrm{d}F^{\mathrm{Tr}}_{\infty,T}}{\mathrm{d}b}\propto\mathrm{Im}[\braket{\sigma_{1}^{x}}_{\mathrm{gexp}}]\sim|b_{c}-b|^{-1/2}$ for $b\lesssim b_{c}$. Such a strong singularity is prohibited for equilibrium free energy density since the thermal expectation value of a local observable cannot diverge. We call this transition an exceptional DQPT, as it turns out to originate from the occurrence of an exceptional point of a nonunitary operator that is dual to $U$. As shown below, an exceptional DQPT can occur for $F^{\{\mathrm{Tr}/\downarrow\uparrow\}}_{\infty,T}$ with $J=\frac{\pi}{4}+\frac{n\pi}{2}\>(n\in\mathbb{Z})$ and even/odd $T$ and is robust under certain weak perturbation (such as $h$), which is deeply related to the hidden symmetry of our setup. We note that the value of $b_{c}$ itself depends on the parameters, such as $T$. We also note that, while the divergence of the derivative of dynamical free energy was recently found in Ref. Bhattacharyya20 for an integrable system, the connection to the underlying symmetry was not discussed.

The exceptional DQPT occurs at a different point from the self-dual points, which are $J=\frac{\pi}{4}+\frac{n\pi}{2}$ and $b=\frac{\pi}{4}+\frac{m\pi}{2}\>(n,m\in\mathbb{Z})$ and known in the context of quantum many-body chaos Bertini18 ; Bertini19 . As discussed in Supplementary Note 5, we find that crossing self-dual points entails DQPT universally for $F^{\mathrm{Tr}/\uparrow\uparrow/\downarrow\uparrow}_{\infty,T}$ with any $T$ and $h$, whose criticality is analogous to that for the conventional DQPT (see Fig. 2).

We stress that DQPTs in our model do not appear as infinite-time averages of expectation values of local observables (see Supplementary Note 4), in contrast with the observation in Ref. Zunkovic18 . Indeed, our DQPTs occur at nonintegrable points, where the infinite-time averages of expectation values trivially thermalize because of the Floquet eigenstate thermalization hypothesis Kim14 . This means that our DQPTs are unique to finite-time regimes, in which time serves as an important parameter in stark contrast with conventional phase transitions.

Spacetime duality and hidden symmetries. To understand the above behaviors, we employ the space-time duality Akila16 of our Floquet operator. This is an exact method to switch the role of time and space and rewrite $U^{T}$ with $L$ product of a space-time-dual transfer matrix $\tilde{U}$, which involves $T$ spins. Using this method, we can rewrite the dynamical free energy as

where the nonunitary operator $\tilde{U}$ depends on the type of $F_{L,T}$. For example, we have

with the periodic boundary condition for $F^{\mathrm{Tr}}_{L,T}$ Akila16 ; Bertini18 ; Bertini19 (see Supplementary Note 2 for the proof and the similar construction for ${\tilde{U}}_{\uparrow\uparrow/\downarrow\uparrow}$, which corresponds to $F^{\uparrow\uparrow/\downarrow\uparrow}_{L,T}$). Here, $\tilde{b}=-\pi/4-i\log(\tan J)/2$, $\tilde{J}=-\pi/4-i\log(\tan b)/2$ and $C=(\sin 2b/\sin 2\tilde{b})^{T/2}/2$.

Let $\lambda_{\mathrm{M,\alpha}}=|\lambda_{\mathrm{M}}|e^{i\theta_{\alpha}}$ be eigenvalues of $\tilde{U}$ whose modulus gives the largest one among all eigenvalues. Here, $\alpha\>(=1,\cdots,n_{\mathrm{deg}})$ is the label of the degeneracy, where $n_{\mathrm{deg}}$ is the number of eigenvalues giving the maximum modulus. For large $L$, $F_{L,T}$ is dominated by these largest eigenvalues, i.e.,

In the thermodynamic limit, the second term vanishes.

Similar to the discussion noted in Ref. Andraschko14 , DQPTs occur when the eigenstate that gives the largest eigenvalue switches. For typical cases, conventional DQPTs occur when maximum of two eigenvalues with different $\theta_{\alpha}$ switches accidentally, where $n_{\mathrm{deg}}=1$ for each phase and $n_{\mathrm{deg}}=2$ at transition (Fig. 3a).

In contrast, hitherto unknown dynamical phases and transitions can appear when $\tilde{U}$ possesses AUS Gong18T ; Kawabata18T ; Kawabata19S . In nonunitary physics, the operator $\tilde{U}$ is said to have the AUS when some unitary operator $V$ and $\phi\in\mathbb{R}$ exist and $V\tilde{U}^{*}V^{\dagger}=e^{i\phi}\tilde{U}$ is satisfied (see Table I). As detailed in the Methods section, nonunitary operator $\tilde{U}$ is called Class A if $\tilde{U}$ does not have the AUS, Class AI when the AUS exists and the corresponding $V$ satisfies $VV^{*}=\mathbb{I}$, and Class AII when the AUS exists and the corresponding $V$ satisfies $VV^{*}=\mathbb{I}$. A particularly important class is Class AI, where the spectral transition unique to nonunitarity, i.e., spontaneous AUS breaking, occurs with the change of parameters. In this case, the eigenstates do and do not respect the AUS for each phase separated at the critical point, which is called the exceptional point. Through the transition, two eigenvalues are attracted, degenerated (at the exceptional point), and repelled in a singular manner (see Fig. 3b).

We find that some of our Floquet operators $U$ can have hidden AUS of $\tilde{U}$ for $J=\frac{\pi}{4}+\frac{n\pi}{2}\>(n\in\mathbb{Z})$ (see Table II and Methods section). Particularly, $\tilde{U}_{\mathrm{Tr}}$ belongs to Class AI for even $T$ (and AII for odd $T$), and $\tilde{U}_{\downarrow\uparrow}$ belongs to Class AI for odd $T$ (and AII for even $T$) as long as $J=\frac{\pi}{4}+\frac{n\pi}{2}\>(n\in\mathbb{Z})$. In contrast, $\tilde{U}_{\uparrow\uparrow}$ does not have AUS and belongs to Class A in general.

The above symmetries clearly explain the origin of the exceptional DQPT: as shown in Fig. 3b, this transition occurs when eigenvalues with the largest and the second-largest modulus collide under Class AI AUS, i.e., at the many-body exceptional point Ashida17 ; Hamazaki19N ; Luitz19 for $\lambda_{\mathrm{M}}$. It is known that this (second-order) exceptional point entails a universal spectral singularity, where the gap between two eigenvalues behave like $|b-b_{c}|^{1/2}$. This leads to the previously-mentioned notable divergence of the generalized expectation value $\sim(b_{c}-b)^{-1/2}$ for $b<b_{c}$, where $-1/2$ is also known to be a universal critical exponent.

For $F^{\mathrm{Tr}}_{\mathrm{\infty},{T=6}}$, the phases with $b<b_{c}\simeq 0.0257\pi$ and $b>b_{c}$ correspond to hidden AUS-unbroken and AUS-broken phases, respectively. This is highlighted by the generalized correlation function, $C(r)=|\braket{\sigma_{1}^{z}\sigma_{r+1}^{z}}_{\mathrm{gexp}}-\braket{\sigma_{1}^{z}}_{\mathrm{gexp}}\braket{\sigma_{r+1}^{z}}_{\mathrm{gexp}}|$ (see Fig. 4 and the Methods section). While $C(r)$ decays exponentially as $\sim e^{-r/\xi_{\mathrm{cor}}}$ in the AUS-unbroken phase, the correlation length diverges as $\xi_{\mathrm{cor}}\sim(b_{c}-b)^{-1/2}$ as it approaches the exceptional DQPT point. At AUS-broken phases, $\xi_{\mathrm{cor}}$ diverges and long-range order appears. Notably, we find that $C(r)$ oscillates with the oscillation length $\xi_{\mathrm{osc}}$, which also diverges near the exceptional DQPT $\xi_{\mathrm{osc}}\sim(b-b_{c})^{-1/2}$. We remark that the qualitative signature of the transition can be captured by the existence of the long-range order even for relatively small systems, which are relevant for experiments (see Supplementary Note 6).

Here, we comment on the relation with the seminal work by Lee, Yang Lee52 ; Yang52 and Fisher Fisher78 , who investigated thermodynamic phase transitions by non-Hermitian operators. While our motivation is to investigate DQPTs occurring at finite times, which is different from their motivation, there exists some mathematical analogy. In fact, the exceptional DQPT can be regarded as the realization of the edge singularity of the partition-function zeros at physical (i.e., real) parameters, as discussed in the Methods section.

Hidden Class AI AUS also enables us to discuss conditions for having exceptional DQPTs. In our prototypical stroboscopic Ising model, we can observe the exceptional DQPT by considering $F^{\mathrm{Tr}}_{\infty,T}$ with even $T$ and $F^{\downarrow\uparrow}_{\infty,T}$ with odd $T$ under the condition $J=\pi/4+n\pi/2\>(n\in\mathbb{Z})$ (see Supplementary Note 3 for the example of $F^{\downarrow\uparrow}_{L,T}$). Note that this transition is robust even if the value of $h$ is slightly perturbed since the transition is protected by AUS. We also stress that $J$ cannot be generic in our anlysis: $J=\pi/4+n\pi/2\>(n\in\mathbb{Z})$ is important for the exceptional DQPT because it ensures the antiunitary symmetry for $\tilde{U}$. Investigation of the exceptional DQPT for other values of $J$ is a future problem.

Signature through quasi-local observables. Next, we show that the signature of our DQPTs is accessible through the expectation values of quasi-local observables, which are more experimentally friendly than the overlap itself (in other words, the DQPT affects the behavior of the expectation values of the quasi-local observables). We also demonstrate that the exceptional DQPT is easier to measure with finite-size scaling analysis than the conventional DQPT, thanks to its strong singularity. We here explain this fact by focusing on $F^{\downarrow\uparrow}_{L,T}$ in Eq. (3), instead of $F^{\mathrm{Tr}}_{L,T}$, since its operational meaning in experimental situations is more direct. We note that $F^{\downarrow\uparrow}_{\infty,T}$ shows the exceptional DQPT for $b=b_{c}\simeq 0.446\pi$ with $h=1.3$, $T=5$ and $J=-\pi/4$, where the AUS is broken for $b<b_{c}$ and unbroken for $b>b_{c}$ (this is opposite to the case for $F^{\mathrm{Tr}}_{\infty,T}$).

To see our argument, we introduce the following quantity

where $P_{f}^{(l)}=\bigotimes_{i=1}^{l}\ket{\downarrow}_{i}\bra{\downarrow}_{i}$ and $P_{f}^{(l)}(T)=U^{-T}P_{f}^{(l)}U^{T}$ is the Heisenberg representation. While $P_{f}^{(l=L)}=\bigotimes_{i=1}^{L}\ket{\downarrow}_{i}\bra{\downarrow}_{i}=\ket{\psi_{f}}\bra{\psi_{f}}$ and Eq. (8) reduces to $F^{\downarrow\uparrow}_{L,T}$ for $l=L$, $P_{f}^{(l)}$ becomes quasi-local when $l=\mathrm{O}(L^{0})\ll L$ Halimeh20 ; Bandyopadhyay20 . For the latter case, Eq. (8) is represented by the standard expectation value of the quasi-local observable, which describes the presence of consecutive spin-down domain at size $l$, at time $T$. Note that such spin domains have been measured in ion experiments using single-site imaging Zhang17 ; Tan21 .

We argue that the signature of the exceptional DQPT can be captured by $F^{\downarrow\uparrow(l)}_{L,T}$ and its derivative even for relatively small $l$, which is more experimentally friendly than the dynamical free energy density itself. Figure 5 shows the $b$-dependence of $F^{\downarrow\uparrow(l)}_{L,T}$ and $\partial F^{\downarrow\uparrow(l)}_{L,T}/\partial b$ for different $l\>(=2,3,4,5,6,\infty)$. We find that the peak develops even for small $l$ around the exceptional DQPT ($b\simeq 0.44\pi$). Particularly, the peaks for the derivative become rapidly sharper as increasing $l$, reflecting the divergence for $l=L\rightarrow\infty$. Our results physically mean that, in this setting, large spin-down domains are rapidly suppressed toward the exceptional DQPT critical point.

We also note that the sharp peaks indicate the experimental advantage of considering the exceptional DQPT compared with the conventional DQPT. Indeed, as shown in Fig. 5, we cannot find sharp peaks for $l\leq 6$ for the conventional DQPT ($b\simeq 0.33\pi$). This indicates that the exceptional DQPT is easier to detect even with small $l$ than the conventional DQPT because of its unique singularity, which is another advantage for our analysis.

Discussions

Although we have demonstrated the singularity of the dynamical free energy and the oscillatory long-range order for the spontaneous antiunitary symmetry breaking, one may wonder whether we can define an order parameter that is nonzero only for the symmetry-breaking phase. As detailed in Supplementary Note 7, we show that an order parameter can be explicitly constructed using different-time generalized observables. This indicates that antiunitary symmetry breaking physically means the symmetry breaking of the dynamical interference patterns, which cannot be diagnosed by the usual single-time expectation values.

The exceptional DQPT appears in other situations as well as the above situation. When we change $h$ instead of $b$, AUS of $\tilde{U}_{\mathrm{Tr}/\downarrow\uparrow}$ is preserved and the exceptional DQPT appears for even/odd $T$, meaning that $\braket{\sigma^{z}_{i}}_{\mathrm{gexp}}$ diverges. We also stress that the exceptional DQPT is not restricted to the stroboscopic Ising model but occurs for a broader class of Floquet systems, as shown in Supplementary Note 8.

To conclude, we have shown that the spontaneous antiunitary symmetry breaking leads to the unconventional universal DQPT, i.e., the exceptional DQPT, uniquely at finite times in Floquet quantum many-body systems. The appearance of finite-time phase transitions related to nonunitary physics can be understood from the spacetime duality. We have also demonstrated that the signatures of the exceptional DQPTs are observed through quasi-local observables that are accessible by state-of-the-art experiments Zhang17 ; Tan21 . Notably, the signature of the exceptional DQPTs is easier to observe than that of the normal DQPTs because of their strong singularity.

Our result paves the way to study completely unknown phases in short-time regimes, where time is regarded as a crucial parameter. As demonstrated in this work, our method via spacetime duality is useful for investigating unconventional finite-time phase transitions for quantum many-body unitary dynamics through the scope of nonunitary many-body physics. One of the promising directions is to classify such dynamical phases by the symmetries of the spacetime-dual operator in light of non-Hermitian symmetries, which are completely classified only recently Gong18T ; Kawabata19S .

Methods

Antiunitary symmetry of $\tilde{U}$. Let us assume that a nonunitary operator $\tilde{U}$ satisfies $V\tilde{U}^{*}V^{\dagger}=e^{i\phi}\tilde{U}$ for some unitary operator $V$ and $\phi\in\mathbb{R}$. According to the recent classification of non-Hermitian systems Kawabata19S , $\tilde{U}$ is called Class A without AUS, Class AI when $V$ with $VV^{*}=\mathbb{I}$ exists, and Class AII when $V$ with $VV^{*}=-\mathbb{I}$ exists. If we consider $\phi=0$ without loss of generality, the eigenvalues of $\tilde{U}$ in Class AI are either real or form complex conjugate pairs. Furthermore, at certain parameters, two real eigenvalues collide and form a complex conjugate pair, which can be called spontaneous AUS-breaking transition. In fact, while eigenstates $\ket{\phi}$ are symmetric under AUS in the AUS-unbroken phase, i.e., $V\ket{\phi}^{*}=\ket{\phi}$, $V\ket{\phi}^{*}$ and $\ket{\phi}$ are different int the AUS-broken phase. At the transition point, known as the exceptional point, two eigenstates become equivalent, which offers a unique feature for nonnormal matrices. In Class AII, on the other hand, eigenvalues generically form complex conjugate pairs and are not real in the presence of the level repulsion Hamazaki20U .

Our Floquet operators $U$ can have such hidden antiunitary symmetries of $\tilde{U}$ for $J=\frac{\pi}{4}+\frac{n\pi}{2}\>(n\in\mathbb{Z})$: indeed, we find $V=\prod_{\tau=1}^{T}e^{i\frac{\pi}{2}\sigma_{\tau}^{y}}$ for $\tilde{U}_{\mathrm{Tr}}$ and $V=\mathcal{P}\prod_{\tau=1}^{T-1}e^{i\frac{\pi}{2}\sigma_{\tau}^{y}}$ for $\tilde{U}_{\downarrow\uparrow}$, where $\mathcal{P}$ is the parity operator exchanging $\tau$ and $T-\tau$ (see Supplementary Note 3 for the detailed calculation). Since $VV^{*}$ takes either $+\mathbb{I}$ or $-\mathbb{I}$ depending on $T$, we find that $\tilde{U}_{\mathrm{Tr}}$ belongs to Class AI for even $T$ and AII for odd $T$, and that $\tilde{U}_{\downarrow\uparrow}$ belongs to Class AI for odd $T$ and AII for even $T$ as long as $J=\frac{\pi}{4}+\frac{n\pi}{2}\>(n\in\mathbb{Z})$. On the other hand, $\tilde{U}_{\uparrow\uparrow}$ does not have AUS and belongs to Class A in general.

Generalized correlation function. To calculate the generalized correlation function, we first note the dual representation

Here, we choose the time point $\tau$ for the dual spin $\sigma^{z}_{\tau}$ as $\tau=1$. Inserting the eigenstate decomposition of $\tilde{U}=\sum_{\alpha}\lambda_{\alpha}\ket{\phi_{\alpha}}\bra{\chi_{\alpha}}$, we have

for $b\lesssim b_{c}$ and large $L$. Here, $0$ and $1$ respectively indicate the labels of eigenvalues with the largest and the second-largest modulus. From this the generalized correlation length is obtained as $\xi_{\mathrm{cor}}=-(\ln\lambda_{1}/\lambda_{0})^{-1}\simeq\lambda_{0}/(\lambda_{0}-\lambda_{1})$ and behaves as $\sim(b_{c}-b)^{-1/2}$ near the exceptional DQPT.

For $b>b_{c}$, $C(r)$ contains a term $e^{-ir\Delta}$ even in the thermodynamic limit, where $\Delta\>(<\pi)$ is the difference between angles of two complex-conjugate eigenvalues. Thus the oscillation length becomes $\xi_{\mathrm{osc}}=\frac{2\pi}{\Delta}$ and behaves as $\sim(b-b_{c})^{-1/2}$ near the exceptional DQPT.

Partition-function zeros. Phase transitions occur when the zeros of the partition function $e^{-LF_{L,T}}$, whose parameter (especially $b$ in our context) regime is extended to a complex one, accumulate at real values in the thermodynamic limit Fisher78 ; Heyl18 . Accumulation points of the partition-function zeros are thus read out from the points where maximum eigenvalues switch when we add proper perturbation $\delta b\>(\in\mathbb{C})$ whose magnitude is infinitesimal Andraschko14 . Notably, the partition-function zeros accumulate along the real axis when the complex-conjugate pair contributes to maximum eigenvalues with $n_{\mathrm{deg}}=2$ owing to AUS of $\tilde{U}$. This is because one of the eigenvalues that form the complex conjugate at $b\in\mathbb{R}$ becomes larger and smaller than the other for $b+\delta b$ and $b-\delta b\>(\delta b\in i\mathbb{R})$, respectively. Moreover, we find that these zeros on the real axis (say $b\geq b_{c}$) terminate at the exceptional DQPT ($b=b_{c}$). This means that the exceptional DQPT corresponds to the realization of the edge singularity of the partition-function zeros at physical parameters on the real axis.

Acknowledgments We are grateful to Kohei Kawabata, Nobuyuki Yoshioka, Keiji Saito, and Mamiko Tatsuta for fruitful comments. We thank Jad C. Halimeh, Amit Dutta, and Subinay Dasgupta for notifying us of related papers. The numerical calculations were carried out with the help of QUSPIN Weinberg19 .