Intrinsic Jump Character of the First-Order Quantum Phase Transitions Qiang Luo,1, ∗ Jize Zhao,2, † and Xiaoqun Wang3, 4, 5, ‡ 1Department of Physics, Renmin University of China, Beijing 100872, China 2School of Physical Science and Technology & Key Laboratory for Magnetism and Magnetic Materials of the MoE, Lanzhou University, Lanzhou 730000, China 3Key Laboratory of Artificial Structures and Quantum Control (Ministry of Education), School of Physics and Astronomy, Tsung-Dao Lee Institute, Shanghai Jiao Tong University, Shanghai 200240, China 4Collaborative Innovation Center for Advanced Microstructures, Nanjing 210093, China 5Beijing Computational Science Research Center, Beijing 100084, China (Dated: August 20, 2019) We find that the first-order quantum phase transitions (QPTs) are characterized by intrinsic jumps of relevant operators while the continuous ones are not. Based on such an observation, we propose a bond reversal method where a quantity D, the difference of bond strength (DBS), is introduced to judge whether a QPT is of first order or not. This method is firstly applied to an exactly solvable spin-1/2 XXZ Heisenberg chain and a quantum Ising chain with longitudinal field where distinct jumps of D appear at the first-order transition points for both cases. We then use it to study the topological QPT of a cross-coupled (J×) spin ladder where the Haldane–rung-singlet transition switches from being continuous to exhibiting a first-order character at J×,I ≃ 0.30(2). Finally, we study a recently proposed one-dimensional analogy of deconfined quantum critical point connecting two ordered phases in a spin-1/2 chain. We rule out the possibility of weakly first-order QPT because the DBS is smooth when crossing the transition point. Moreover, we affirm that such transition belongs to the Gaussian universality class with the central charge c = 1. Introduction.– Understanding how strongly correlated a first-order QPT called a photon-blockade breakdown systems order into different phases as well as the tran- was observed experimentally in a driven circuit quantum sitions among them remains one of the most fundamen- electrodynamics system20. The finite-size scaling of gap tal and significant problems in modern condensed mat- and various probes borrowed from quantum information ter physics1–3. In particular, the quantum phase tran- sciences near the transition points of the first-order QPTs sitions (QPTs), which occurr at zero temperature, are have been discussed until recently21–25. Therefore, it is of omnipresent phenomena and could in general be classi- vital importance to devise appropriate tools for a proper fied into two types. One is continuous when the ground characterization of their dominating features. state of the system changes continuously at the transi- Let us consider a Hamiltonian of the form1 tion point, accompanied by a diverging correlation length and vanishing energy gap. The other is, instead, of (λ)= 0 + λ I (1) first order when the order parameter and other rele- H H H vant observables display discontinuity across the tran- where λ is a driving parameter. If 0 and I com- H H sition point. While traditional continuous QPTs are well mute, then both of them could be simultaneously diag- described by the Landau-Ginzburg-Wilson (LGW) the- onalized and eigenfunctions are independent of λ. This ory, recent years have witnessed some exceptions such means that the spectra of 0 and I are irrelevant of 4–6 H H as topological QPT and deconfined quantum criti- λ, while the total ground-state energy Eg(λ) = (λ) hH i cal point (DQCP)7–10 that beyond the scope of LGW could vary linearly with λ. Consequently, there can be paradigm. The topological QPT may take place between a level-crossing point λt where the ground-state energy two disordered phases, thus it can not be detected by per site eg(λ) = (1/L)Eg(λ) = e0 + λeI as a function of local order parameters11,12. What’s more, it could be ei- λ exhibits nonanalyticity. Here, L is the number of lat- arXiv:1906.06553v2 [cond-mat.str-el] 19 Aug 2019 ther continuous or of first order upon a fine-tuned interac- tice sites. It should be aware that continuous QPT could tion strength13–16. The DQCP was proposed by Senthil also occur and we move the discussion to the supplemen- 7,8 26 et. al. whereby a continuous QPT occurs between two tal material (SM) . Taking the derivative of eg(λ) with spontaneously symmetry breaking (SSB) phases and the respect to λ, a jump of eI at λt will appear, indicating critical point implies an emergent symmetry. The J- of a first-order QPT. The jump also reflects the struc- Q model9 is such an example where extensive numeri- tural change of ground-state wave function. Because of cal studies provide evidences for a continuous (or weakly the continuity of energy eg(λ), a similar jump for e0 is first order) transition between a N´eel phase and a valence also expected to eliminate the singularity. To detect the 17–19 bond solid (VBS) phase . transition point λt, a quantity dubbed the difference of bond strength (DBS) is introducedD to magnify the jump In contrast to the continuous QPTs, the first-order behaviors. It is defined as QPTs are less studied so far despite they appear fre- quently in quantum many-body systems. Remarkably, = e sgn(λ )e , (2) D 0 − t I 2 where the minus sign reflects the spirit of the bond re- years later37. According to Eq. (2), the DBS is defined versal method. Nevertheless, in most systems the two as = Sz Sz + 2 Sx Sx . DL h L/2 L/2+1i h L/2 L/2+1i terms 0 and I do not commute, resulting in a cum- At ∆ = 1 the Hamiltonian Eq. (3) posses a hidden H H − bersome expression of eg(λ) versus λ. However, the main FM SU(2) symmetry, resulting in a ground-state mani- spirit remains unchanged in that the jump character of fold of degenerate SU(2) multiplet corresponding to the D faithfully inherits the discontinuity of first-order QPTs. largest total spin. The model is not conformal invari- In what follows we will firstly illustrate the bond re- ant and dramatic changes of its entanglement behaviors versal method in two different but thoroughly studied occur38–41. In Fig. 1(a) we show the ground-state energy one-dimensional (1D) spin models that possess first-order eg around the transition point ∆ = 1. A vivid cusp of QPTs: (i) a celebrated spin-1/2 Heisenberg XXZ chain the energy curve could be spotted at− ∆ = 1, while it for which all the energy as well as the DBS can be cal- turns to be a jump of in Fig. 1(b). This− gives clearly D L culated analytically and (ii) a quantum Ising chain with a first glance of the jumpD character of in the first-order both longitudinal and transverse fields which does not QPT26. D host exact solution but the transition line is well-known. Having established the cornerstone of our method, we -0.2 0.3 then apply it to (iii) a topological QPT of a cross-coupled spin ladder and (iv) a recently proposed spin-1/2 chain 0.1 with DQCP, both of which are beyond the scope of the -0.3 conventional LGW paradigm. All the models are studied by the density-matrix renormalization group (DMRG) -0.1 method27–30, which is a powerful tool for dealing with quantum-mechanical problems in 1D systems. We uti- -0.4 lize the periodic boundary condition (PBC) for the first -0.3 two cases to have a better comparison with analytical results. For the latter cases, however, we turn to the -0.5 -0.5 open boundary condition (OBC) which is beneficial to -2 -1 0 -2 -1 0 large-scale numerical calculations. XXZ chain.– The 1D spin-1/2 Heisenberg XXZ chain has long served as the workhorse for the study of quan- FIG. 1: (a) The ground-state energy eg for L = 64 (red rhom- tum magnetism31. Its Hamiltonian is given by bus), 128 (blue circle), and TDL (black line). (b) The same setup as (a) for DBS DL. L 1 = S+S− + S−S+ + ∆SzSz (3) H 2 i i+1 i i+1 i i+1 Quantum Ising chain with longitudinal field.– The 1D i=1 X quantum Ising chain is integrable and its exact solution 42 ± x y was presented by Pfeuty in 1970 . It owns a continuous where Si = Si + iSi is the raising/lowering operator at QPT of the Ising universality class separating a FM and site i and ∆ is the anisotropic parameter. In particular, a paramagnetic phase at the critical value of transverse in the region 1 < ∆ 1 the ground state is a Lut- 43 − ≤ field hz( 0) . What’s more, an emergent E8 symmetry tinger liquid (LL) with a gapless excitation spectrum. was experimentally≥ verified around the critical point44. The ground-state energy per site eg in the thermody- By introducing a longitudinal field h the model is no namic limit (TDL) L can be calculated as32–34 x → ∞ longer integrable except for a specially fine-tuned weak 45 ∞ longitudinal field . The total Hamiltonian is thus given LL ∆ sin πυ tanh υx eg (∆) = 1 dx (4) by 4 − π 0 − tanh x Z L with ∆ = cos πυ. Beyond the critical region it presents = σxσx + h σz + h σx (5) H − i i+1 z i x i a long-range ferromagnetic (FM) or antiferromagnetic i=1 (AFM) order, exhibiting in correspondence to the FM X whereσ ˆ = (σx, σy, σz) are the Pauli matrices.