Ground-state information geometry and quantum criticality in an inhomogeneous spin model*
Ma Yu-Quan†
School of Applied Science, Beijing Information Science and Technology University, Beijing 100192, China

Corresponding author. E-mail: yqma@bistu.edu.cn

*Project supported by the National Natural Science Foundation of China (Grant Nos. 11404023 and 11347131).

Abstract

We investigate the ground-state Riemannian metric and the cyclic quantum distance of an inhomogeneous quantum spin-1/2 chain in a transverse field. This model can be diagonalized by using a general canonical transformation to the fermionic Hamiltonian mapped from the spin system. The ground-state Riemannian metric is derived exactly on a parameter manifold ring S1, which is introduced by performing a gauge transformation to the spin Hamiltonian through a twist operator. The cyclic ground-state quantum distance and the second derivative of the ground-state energy are studied in different exchange coupling parameter regions. Particularly, we show that, in the case of exchange coupling parameter Ja = Jb, the quantum ferromagnetic phase can be characterized by an invariant quantum distance and this distance will decay to zero rapidly in the paramagnetic phase.

PACS: 03.65.Vf; 75.10.Pq; 73.43.Nq; 05.70.Jk
Keyword: quantum geometry tensor; topological order; quantum phase transition
1. Introduction

Quantum phase transitions (QPTs) are driven purely by the quantum fluctuations when a parameter of the Hamiltonian describing the system varies.[13] Traditionally, QPTs can be well understood in the framework of the Landau– Ginzburg– Wilson paradigm by resorting to the notions of local-order parameter, long-range correlations and the symmetry breaking. In the past few years, a lot of efforts have been devoted into understanding the QPTs from the information-geometry perspectives, [4, 5] such as quantum entanglement, [69] entanglement entropy, [10, 11] quantum discord, [1214] quantum fidelity and fidelity susceptibility, [1521] Berry phase[2232] and the quantum geometric tensor (QGT).[3341]

Generally, QPTs can be witnessed by some qualitative changes of the ground-state properties when some parameters of the Hamiltonian move across the quantum critical point (QCP). The underlying physical mechanism lies in the fact that the different phases are unconnected by the adiabatic evolution of the ground state. In the vicinity of the QCP, the ground state driven by the parameters of the Hamiltonian will lead to an avoided energy-level crossing between the ground state and the first excited state, where the adiabatic evolution can be destroyed by a vanishing energy gap as the system size tends to infinity. From the perspective of the information geometry, a monopole as a gapless point in the Hamiltonian parameters space will generate interesting effects on the ground-state local or topological properties. The local properties can be captured by some local quantities, i.e., the fidelity susceptibility and the Berry curvature. The topological properties will be captured by some topological quantum numbers, i.e., the Chern number, [42, 43]Z2 number, [4447] and recently the Euler number of the Bloch states manifold.[4851]

Recently, the concept of quantum geometric tensor based on the information-geometry[52, 53] of the ground state has been introduced to understand the QPTs in quantum many-body systems.[3741] What is surprising is that the two approaches of the ground-state Berry curvature and the fidelity susceptibility as some witness to QPTs are unified. Specifically, the real part of the QGT is a Riemannian metric defined in the parameter manifold which has been recognized as the essential part of the fidelity susceptibility. Meanwhile its imaginary part is the Berry curvature whose flux gives rise to the Berry phase. Generally speaking, the Riemannian metric and the Berry curvature will exhibit some singularity or scaling behavior in the quantum critical region under the thermodynamic limit. Particularly, a scaling analysis of the ground-state quantum geometric tensor in the vicinity of the critical points has been performed. So far, these approaches have been applied to detect the phase boundaries in various systems.

In this work, we propose a cyclic quantum distance of the ground state to detect the QPTs in a transverse field inhomogeneous Ising spin-1/2 chain, [54, 55] in which the nearest-neighbor exchange interactions will take alternating parameters between the neighbor sites. This model can be solved exactly by introducing a general canonical transformation to diagonalize the fermionic Hamiltonian mapped from the spin Hamiltonian by the Jordan– Wigner transformation. In our scheme, an extra local gauge transformation is performed to the spin system by a twist operator, which endows the Hamiltonian of the system with a topology of a ring S1 without changing its energy spectrum. We obtain the exact expression of the ground-state Riemannian metric and study the cyclic quantum distance of the ground state on the parameter manifold. We study extensively the ground-state Riemannian metric in different parameter regions of the inhomogeneous Ising chain. Particularly, we show that the quantum ferromagnetic phase in the uniform Ising chain can be marked by an invariant cyclic quantum distance of the ground state, and meanwhile, this distance will decay to zero rapidly in the paramagnetic phase.

2. Model

Let us consider an inhomogeneous two-sublattice Ising spin-1/2 chain in an external magnetic field, which consists of N cells with two sites in each cell. The Hamiltonian reads

where are the local Pauli operators, Ja (Jb) is the exchange coupling, h is the external field and the periodic boundary condition (PBC) has been assumed. In order to investigate the model in a two-dimensional quasi-momentum space, we put this model on an external parameter φ -ring S1. This can be achieved by subjecting the system to a local gauge transformation by a twist operator , which in fact makes the system rotate on the spin along the z direction. It can be verified that is π periodic in the parameter φ . Considering the unitarity of the twist operator , the energy spectrum and critical behavior of the system are obviously independent with the parameter φ . The spin Hamiltonian can be mapped exactly on a spinless fermion model through the Jordan– Wigner transformation, i.e., and , where are the spin ladder operators, and Cl, a(b) are the fermion operators. Hamiltonian is transformed into

Note that the PBC on the spin degrees of freedom implies that CN+ 1, m = eiπ NFC1, m, where denotes the total fermion number. Thus the boundary conditions on the fermionic system will obey PBC or anti-PBC depending on whether NF is even or odd. However, the differences between the two boundary conditions are negligible in the thermodynamic limit. Without loss of generality, we take the PBC on the fermionic system, which means CN+ 1, m = C1, m. Apply the following Fourier transformation and to the Hamiltonian H(φ ), where k = (2π /2aN)n, (n = 0, 1, 2, … , N − 1 and Rla (Rlb = Rla + a) is the coordinate of site a (b) in the l-th cell with the lattice parameter 2a. We can rewrite Hamiltonian in the quasi-momentum space

This Hamiltonian H(φ ) can be exactly diagonalized as

by using the following canonical transformation:

where

and the quasi-particle energy spectrums are

3. Ground-state Riemannian metric and cyclic quantum distance

Now, we focus on the geometric properties of the ground state. The Hamiltonian H(φ ) in Eq.  (4) has been diagonalized in the set of quasi-particle number operators, which allows us to determine all the eigenvalues and eigenvectors. We note that the energy spectrums Λ η k ≥ 0, Λ ν k ≥ 0 and Λ γ k ≤ 0, Λ μ k ≤ 0. The ground state, denoted as | g(φ )〉 , corresponds to the state with the lowest energy, which consists of a state with no η and ν fermions occupied but with γ and μ fermions occupied. Explicitly, the ground state can be constructed as follows:

where is the normalized factor, and | 0〉 is the vacuum states of fermionic operators ak and bk, respectively. It is easy to check that η k| g(φ )〉 = 0, ν k| g(φ )〉 = 0 and , for all k. The corresponding ground-state energy Eg is

Now, we introduce the notion of the quantum geometric tensor of the ground state on the Hamiltonian parameter φ manifold. It can be verified that the quantum geometric tensor can be derived from a gauge invariant distance between two ground states on the U(1) line bundle induced by the quantum adiabatic evolution of the ground state in the parameter φ space. The quantum distance dS between two ground states | g(φ )〉 and | g(φ + δ φ )〉 is given by dS2 = 〈 φ g(φ )dφ | φ g(φ )dφ 〉 . Note that the term | φ g(φ )〉 can be decomposed as , where is the projection operator and is the covariant derivative of | g(φ )〉 on the U(1) line bundle. Under the condition of the quantum adiabatic evolution, the evolution from | g(φ )〉 to | g(φ + δ φ )〉 will undergo a parallel transport in the sense of Levi– Cività from φ to φ + δ φ on the parameter manifold, and hence we can get | Dφ g(φ )〉 = 0. Finally, we can obtain the quantum distance as dS2 = Qφ φ dφ 2, where Qφ φ is defined as the quantum geometric tensor

The concrete expression of Qφ φ can be derived by substituting Eq.  (8) into Eq.  (10). Obviously, the straightforward calculation is tedious, however, it can be derived concisely from the following consideration. Notice that the term 〈 g(φ )| φ g(φ )〉 is nothing but the Berry– Simon connection of the ground state, which can be calculated by considering that each term as and will yield the same result as − 2i cos2(θ k/2), and will yield the result as − 2i cos2(β k/2), and meanwhile, all the other terms are zero. Finally, we can get

In order to calculate the term of 〈 φ g(φ )| φ g(φ )〉 , we note that because φ g(φ )| φ g(φ ) 〉 = 0 (see Eq.  (11)), and so we have

Substituting Eqs.  (11) and (12) into Eq.  (10), we can finally obtain the result as

Now, let us focus on the characterization of the geometric properties for the ground state. In our approach, the ground state | g(φ )〉 is defined in a U(1) line bundle located over a one-dimensional parameter manifold S1, and hence the Riemannian metric as the real part of the QGT is just Qφ φ itself. As discussed above, the ground-state Riemannian metric provide a gauge invariant distance measurement of the ground state on the parameter φ manifold. The quantum distance l between two ground states | g(φ A)〉 and | g(φ B)〉 is given by

To have an explicit view of the dependence of the Riemannian metric on the system size, we can perform a scaling transform to the Qφ φ and denote the Riemannian metric as 𝒢 φ φ = Qφ φ /Ld = [∑ k> 0 (sin2θ k + sin2β k)]/N, where Ld = N as the number of the site and here d = 1 is the dimension of the system. To study the quantum criticality, we are interested in the properties under the thermodynamic limit as N → ∞ , and we have the Riemannian metric

where the summation (∑ k> 0)/N has been replaced by the integral . Obviously, the quantum distance l(0, π ) for a cyclical evolution from | g(0)〉 to | g(π )〉 is given by

In Fig.  1(a), we plot the cyclic quantum distance l(0, π ) as a function of h and α = Jb/Ja with the system size N → ∞ . As a comparison, we also provide numerical results of the second derivative of the ground-state energy with respect to h, as a function of Jb/Ja and h with the system size N = 1001 (see Fig.  1(b)). The cyclic quantum distance l(0, π ) as a function of the external field h with the fixed parameters Ja = 1 and Jb = 1.5 is shown in Fig.  2. In the region of the inhomogeneous coupling JaJb, the cyclic quantum distance and the Riemannian metric have similar trends with the second derivative of the ground-state energy.

Fig.  1. (a) The cyclical quantum distance l(0, π ) as a function of Jb/Ja and h with the fixed parameter Ja = 1, and the system size N → ∞ ; (b) the second derivative of the ground-state energy Eg/N with respect to h, as a function of Jb/Ja and h with the fixed parameter Ja = 1, and the system size N = 1001.

Fig.  2. The cyclical quantum distance l(0, π ) as a function h with the fixed parameters Ja = 1 and Jb = 1.5, and with different system sizes.

It is worth reminding that, in the uniform exchange coupling case Ja = Jb = 1 and the system size N → ∞ , the ground-state Riemannian metric 𝒢 φ φ can be exactly solved as

which leads to an invariant ground-state cyclical distance in the ferromagnetic phase, and in the paramagnetic phase. In Fig.  3, we show the properties of the Riemannian metric in the vicinity of the critical points with fixed parameters Ja = Jb = 1 and different system sizes N. As shown in Fig.  3, the Riemannian metric of the ground state in the ferromagnetic phase is in close to the constant 0.5 with the increase of the system size N. In the thermodynamic limit N → ∞ , the first derivative of the cyclic distance and the metric are discontinuous in the critical point.

Fig.  3. The Riemannian metric g as a function of h with the fixed parameters Ja = Ja = 1, and with different system sizes.

For a measurable consequence of the Riemannian metric, it has been recently shown that the quantum geometric tensor of the Bloch bands can be measured by the current noise spectrum of the band insulators. The current noise spectrum can be expressed as (Eq.  (13) in Ref.  [40])

where Ω BZ denotes the volume of the Brillouin zone and d denotes the dimension of the quasi-momentum space. The current noise spectrum is just defined as the Fourier transform of the current– current correlation function

On the other hand, Marzari and Vanderbilt[49] have pointed out that the integral of the trace of the Riemannian metric over the Brillouin zone can be regarded as a gauge invariant measure

of the delocalization or spread of the Wannier functions.

It is worth noting that the metric component Qφ φ of the Bloch band is closely related to the ground-state quantum distance in the parameter φ space. In fact, the ground state | g(φ )〉 is π periodic in the parameter φ . In the condition of the large sites limit N → ∞ , a cyclic ground-state distance lGS along the φ -ring can be calculated as , which only relates to the metric 𝒢 φ φ on the Bloch band and the concrete expression of the ground state is not required. On the other hand, the topological Euler number index χ of the Bloch band is (see Ref.  [48]) given by , where is the Ricci scalar curvature associate to the Bloch state | un(k)〉 of the n-th Bloch band. In the present model, the expression of Euler number χ can be simplified as

The relation between the cyclic ground-state distance lGS and the Euler number χ of the Bloch band illustrates that the cyclic ground-state distance is an intrinsic property of the Bloch bands. Meanwhile, this relation also tells us that the cyclic ground-state distance lGS in the ferromagnetic phase of a uniform Ising chain is equivalent to the topological Euler number of the band in the quasi-momentum space.

4. Conclusions

In summary, we study the ground-state geometric properties of a two-sublattice quantum Ising chain in a transverse field. In particular, we introduce a gauge transformation to the spin system by a twist operator, which endows the Hamiltonian of the system with a topology of a ring S1 without changing its energy spectrum. The advantages of introducing the twist parameter is that the Hamiltonian describing the system is periodic of the parameter φ , i.e., . Hence, we can naturally discuss the Berry phase and the cyclic quantum distance on the φ -ring S1. On the other hand, our approach can be easily extended to two dimensions because the fermionized system has a natural S1 topology (the periodic gauge | un (k)〉 = eiG· r | un(k + G)〉 , G denotes the reciprocal lattice vector) on the first Brillouin zone. Finally, the spin Hamiltonian can be naturally regarded as a fermionic Hamiltonian on a 2-torus T2 = S1 × S1 in the parameter space (k, φ ). However, it is hard to find a natural S1 or T2 topological structure for the Hamiltonian if we take the Hamiltonian parameters Ja, Jb or h as the parameters manifold. Note that a nontrivial topological number for the quantum states requires a nontrivial topological structure of the Hamiltonian. This is the reason why we introduce a twist parameter to the Hamiltonian instead of using the Hamiltonian parameters Ja, Jb and h as the parameter manifold.

On the φ parameter manifold S1, we derive the exact expression of the ground-state Riemannian metric and define a cyclic quantum distance of the ground state. We study extensively the ground-state Riemannian metric and the cyclic quantum distance in different parameter regions of the two-sublattice Ising chain. Furthermore, we show that the ferromagnetic phase of a uniform Ising chain can be characterized by an invariant quantum metric Gφ φ = 1/2 and this metric will decay to zero with h− 2 in the paramagnetic phase. While the metric Gφ φ is continuous in the critical point, the singularity will appear when we take the derivative of Gφ φ with respect to the external field h. This approach provides an intrinsic description of the geometric properties of the ground state. We hope that the current work will raise renewed interest in the understanding of the geometric nature of the ground state in quantum condensed-matter systems.

Reference
1 Sachdev S 1999 Quantum Phase Transitions Cambridge Cambridge University Press [Cited within:1]
2 Sondhi S L, Girvin S M, Carini J P and Shahar D 1997 Rev. Mod. Phys. 69 315 DOI:10.1103/RevModPhys.69.315 [Cited within:1]
3 Vojta M 2003 Rep. Prog. Phys. 66 2069 DOI:10.1088/0034-4885/66/12/R01 [Cited within:1]
4 Bengtsson I and Zyczkowski K 2008 Geometry of Quantum States: An Introduction to Quantum Entanglement Cambridge Cambridge University Press [Cited within:1]
5 Ortiz G 2010 Understand ing in Quantum Phase TransitionsCarr Led. Boca Raton Taylor & Francis [Cited within:1]
6 Osterloh A, Amico L, Falci G and Fazio R 2002 Nature 416 608 DOI:10.1038/416608a [Cited within:1]
7 Osborne T J and Nielsen M A 2002 Phys. Rev. A 66 032110 DOI:10.1103/PhysRevA.66.032110 [Cited within:1]
8 Vidal G, Latorre J I, Rico E and Kitaev A 2003 Phys. Rev. Lett. 90 227902 DOI:10.1103/PhysRevLett.90.227902 [Cited within:1]
9 Gu S J, Deng S S, Li Y Q and Lin H Q 2004 Phys. Rev. Lett. 93 086402 DOI:10.1103/PhysRevLett.93.086402 [Cited within:1]
10 Kitaev A and Preskill J 2006 Phys. Rev. Lett. 96 110404 DOI:10.1103/PhysRevLett.96.110404 [Cited within:1]
11 Levin M and Wen X G 2006 Phys. Rev. Lett. 96 110405 DOI:10.1103/PhysRevLett.96.110405 [Cited within:1]
12 Ollivier H and Zurek W H 2001 Phys. Rev. Lett. 88 017901 DOI:10.1103/PhysRevLett.88.017901 [Cited within:1]
13 Datta A, Shaji A and Caves C M 2008 Phys. Rev. Lett. 100 050502 DOI:10.1103/PhysRevLett.100.050502 [Cited within:1]
14 Campbell S, Mazzola L, de Chiara G, Apollaro T J G, Plastina F, Busch T and Paternostro M 2013 New J. Phys. 15 043033 DOI:10.1088/1367-2630/15/4/043033 [Cited within:1]
15 Zanardi P, Giorda P and Cozzini M 2007 Phys. Rev. Lett. 99 100603 DOI:10.1103/PhysRevLett.99.100603 [Cited within:1]
16 Venuti L C and Zanardi P 2007 Phys. Rev. Lett. 99 095701 DOI:10.1103/PhysRevLett.99.095701 [Cited within:1]
17 You W L, Li Y W and Gu S J 2007 Phys. Rev. E 76 022101 [Cited within:1]
18 Gu S J 2010 Int. J. Mod. Phys. B 24 4371 DOI:10.1142/S0217979210056335 [Cited within:1]
19 Chen S, Wang L, Gu S J and Wang Y 2007 Phys. Rev. E 76 061108 DOI:10.1103/PhysRevE.76.061108 [Cited within:1]
20 Yang S, Gu S J, Sun C P and Lin H Q 2008 Phys. Rev. A 78 012304 DOI:10.1103/PhysRevA.78.012304 [Cited within:1]
21 Zhao J H and Zhou H Q 2009 Phys. Rev. B 80 014403 DOI:10.1103/PhysRevB.80.014403 [Cited within:1]
22 Berry M V 1984 Proc. R. Soc. London A 392 45 DOI:10.1098/rspa.1984.0023 [Cited within:1]
23 Simon B 1983 Phys. Rev. Lett. 51 2167 DOI:10.1103/PhysRevLett.51.2167 [Cited within:1]
24 Carollo A C M and Pachos J K 2005 Phys. Rev. Lett. 95 157203 DOI:10.1103/PhysRevLett.95.157203 [Cited within:1]
25 Hamma A 2006arXiv: quant-ph/0602091 [Cited within:1]
26 Zhu S L 2006 Phys. Rev. Lett. 96 077206 DOI:10.1103/PhysRevLett.96.077206 [Cited within:1]
27 Ma Y Q and Chen S 2009 Phys. Rev. A 79 022116 DOI:10.1103/PhysRevA.79.022116 [Cited within:1]
28 Hirano T, Katsura H and Hatsugai Y 2008 Phys. Rev. B 77 094431 DOI:10.1103/PhysRevB.77.094431 [Cited within:1]
29 Hatsugai Y 2010 New J. Phys. 12 065004 DOI:10.1088/1367-2630/12/6/065004 [Cited within:1]
30 Fukui T and Fujiwara T 2009 J. Phys. Soc. Jpn. 78 093001 DOI:10.1143/JPSJ.78.093001 [Cited within:1]
31 Ma Y Q, Yu Z X, Wang D S, Xie B H and Li X G 2012 EPL 100 60001 DOI:10.1209/0295-5075/100/60001 [Cited within:1]
32 Ma Y Q, Yu Z X, Wang D S and Li X G 2013 Phys. Lett. A 377 1250 DOI:10.1016/j.physleta.2013.03.021 [Cited within:1]
33 Provost J P and Vallee G 1980 Commun. Math. Phys. 76 289 DOI:10.1007/BF02193559 [Cited within:1]
34 Berry M V 1989 Geometric Phases in PhysicsShapere A Wilczek F eds. Singapore World Scientific [Cited within:1]
35 Resta R 2005 Phys. Rev. Lett. 95 196805 DOI:10.1103/PhysRevLett.95.196805 [Cited within:1]
36 Haldane F D M 2011 Phys. Rev. Lett. 107 116801 DOI:10.1103/PhysRevLett.107.116801 [Cited within:1]
37 Ma Y Q, Chen S, Fan H and Liu W M 2010 Phys. Rev. B 81 245129 DOI:10.1103/PhysRevB.81.245129 [Cited within:1]
38 Rezakhani A T, Abasto D F, Lidar D A and Zanardi P 2010 Phys. Rev. A 82 012321 DOI:10.1103/PhysRevA.82.012321 [Cited within:1]
39 Matsuura S and Ryu S 2010 Phys. Rev. B 82 245113 DOI:10.1103/PhysRevB.82.245113 [Cited within:1]
40 Neupert T, Chamon C and Mudry C 2013 Phys. Rev. B 87 245103 DOI:10.1103/PhysRevB.87.245103 [Cited within:1]
41 Legner M and Neupert T 2013 Phys. Rev. B 88 115114 DOI:10.1103/PhysRevB.88.115114 [Cited within:2]
42 Thouless D J, Kohmoto M, Nightingale M P and Nijs M 1982 Phys. Rev. Lett. 49 405 DOI:10.1103/PhysRevLett.49.405 [Cited within:1]
43 Niu Q and Thouless D J 1984 J. Phys. A 17 2453 DOI:10.1088/0305-4470/17/12/016 [Cited within:1]
44 Kane C L and Mele E J 2005 Phys. Rev. Lett. 95 146802 DOI:10.1103/PhysRevLett.95.146802 [Cited within:1]
45 Fu L, Kane C L and Mele E J 2007 Phys. Rev. Lett. 98 106803 DOI:10.1103/PhysRevLett.98.106803 [Cited within:1]
46 Hasan M Z and Kane C L 2010 Rev. Mod. Phys. 82 3045 DOI:10.1103/RevModPhys.82.3045 [Cited within:1]
47 Qi X L and Zhang S C 2011 Rev. Mod. Phys. 83 1057 DOI:10.1103/RevModPhys.83.1057 [Cited within:1]
48 Ma Y Q, Gu S J, Chen S, Fan H and Liu W M 2012 EPL 103 10008 DOI:10.1209/0295-5075/103/10008 [Cited within:2]
49 Kolodrubetz M, Gritsev V and Polkovnikov A 2013 Phys. Rev. B 88 064304 DOI:10.1103/PhysRevB.88.064304 [Cited within:1]
50 Ma Y Q 2014 Phys. Rev. E 90 042133 DOI:10.1103/PhysRevE.90.042133 [Cited within:1]
51 Yang L, Ma Y Q and Li X G 2015 Physica B 456 359 DOI:10.1016/j.physb.2014.09.022 [Cited within:1]
52 Wootters W K 1981 Phys. Rev. D 23 357 DOI:10.1103/PhysRevD.23.357 [Cited within:1]
53 Braunstein S L and Caves C M 1994 Phys. Rev. Lett. 72 3439 DOI:10.1103/PhysRevLett.72.3439 [Cited within:1]
54 de Lima J P, Goncalves L L and Alves T F A 2007 Phys. Rev. B 75 214406 DOI:10.1103/PhysRevB.75.214406 [Cited within:1]
55 Wang L C, Yan J Y and Yi X X 2010 Chin. Phys. B 19 040512 DOI:10.1088/1674-1056/19/4/040512 [Cited within:1]
56 Marzari N and Vand erbilt D 1997 Phys. Rev. B 56 12847 DOI:10.1103/PhysRevB.56.12847 [Cited within:1]