Search

Fermionic Gaussian states are eigenstates of quadratic Hamiltonians and widely used in quantum many-body problems. We propose a highly efficient algorithm that converts fermionic Gaussian states to matrix product states. It can be formulated for finite-size systems without translation invariance, but becomes particularly appealing when applied to infinite systems with translation invariance. If the ground states of a topologically ordered system on infinite cylinders are expressed as matrix product states, then the fixed points of the transfer matrix can be harnessed to filter out the anyon eigenbasis, also known as minimally entangled states. This allows for efficient computation of universal properties such as entanglement spectrum and modular matrices. The potential of our method is demonstrated by numerical calculations in two chiral spin liquids that have the same topological orders as the bosonic Laughlin and Moore-Read states, respectively. The anyon eigenbasis for the first one has been worked out before and serves as a useful benchmark. The anyon eigenbasis of the second one is, however, not transparent and its successful construction provides a nontrivial corroboration of our method., Comment: 13 pages, 7 figures

The low-energy physics of Tomonaga-Luttinger liquids (TLLs) is controlled by the Luttinger parameter. We demonstrate that this parameter can be extracted from a single wave function for one-component TLLs with periodic boundary condition. This method relies on the fact that TLLs are described by conformal field theory in which crosscap states can be constructed. The overlaps between the crosscap states and the ground state as well as some excited states are proved to be universal numbers that directly reveal the Luttinger parameter. In microscopic lattice models, crosscap states are formed by putting each pair of antipodal sites into a maximally entangled state. Analytical and numerical calculations are performed in a few representative models to substantiate the conformal field theory prediction. The extracted Luttinger parameters are generally quite accurate in finite-size systems with moderate lengths, so there is no need to perform data fitting and/or finite-size scaling., Comment: 6+9 pages, 2 figures

We propose a novel mechanism of impurity screening in (1+1)$d$ quantum critical states described by conformal field theories (CFTs). An impurity can be screened if it has the same quantum numbers as some gapless degrees of freedom of the CFT. The common source of these degrees of freedom is the chiral primary fields of the CFT, but we uncover that topological defect lines of the CFT may also take this role. Theoretical analysis relies on the insight that the impurities can be interpreted as edge modes of certain symmetry-protected topological (SPT) states. By stacking a SPT state with a CFT, one or two interfaces on which the SPT edge modes reside are created. If screening occurs due to topological defect lines, a symmetry-enriched CFT with exotic boundary states are obtained. The boundary conditions that appear in these cases are difficult to achieve using previously known methods. Numerical simulations are performed in two quantum spin chains whose bulks are described by the SU(3)$_{1}$ and Spin(5)$_{1}$ CFTs and edges are coupled to spin-1/2 impurities. In both cases, low-energy eigenvalues are consistent with analytical predictions., Comment: 10 pages, 3 figures

Approaching the long-time dynamics of non-Markovian open quantum systems presents a challenging task if the bath is strongly coupled. Recent proposals address this problem through a representation of the so-called process tensor in terms of a tensor network. We show that for Gaussian environments highly efficient contraction to matrix product operator (MPO) form can be achieved with infinite MPO evolution methods, leading to significant computational speed-up over existing proposals. The result structurally resembles open system evolution with carefully designed auxiliary degrees of freedom, as in hierarchical or pseudomode methods. Here, however, these degrees of freedom are generated automatically by the MPO evolution algorithm. Moreover, the semi-group form of the resulting propagator enables us to explore steady-state physics, such as phase transitions.

We introduce a family of SO($n$)-symmetric spin chains which generalize the transverse-field Ising chain for $n=1$. These spin chains are defined with Gamma matrices and can be exactly solved by mapping to $n$ species of itinerant Majorana fermions coupled to a static $\mathbb{Z}_2$ gauge field. Their phase diagrams include a critical point described by the $\mathrm{Spin}(n)_{1}$ conformal field theory as well as two distinct gapped phases. We show that one of the gapped phases is a trivial phase and the other realizes a symmetry-protected topological phase when $n \geq 2$. These two gapped phases are proved to be related to each other by a Kramers-Wannier duality. Furthermore, other elegant structures in the transverse-field Ising chain, such as the infinite-dimensional Onsager algebra, also carry over to our models., Comment: 12 pages, 3 figures, published version

We study quantum phase transitions in Bose-Fermi mixtures driven by interspecies interaction in the quantum Hall regime. In the absence of such an interaction, the bosons and fermions form their respective fractional quantum Hall (FQH) states at certain filling factors. A symmetry-protected topological (SPT) state is identified as the ground state for strong interspecies interaction. The phase transitions between them are proposed to be described by Chern-Simons-Higgs field theories. For a simple microscopic Hamiltonian, we present numerical evidence for the existence of the SPT state and a continuous transition to the FQH state. It is also found that the entanglement entropy between the bosons and fermions exhibits scaling behavior in the vicinity of this transition., Comment: 9 pages, 6 figures

We show that the Klein bottle entropy [Phys. Rev. Lett. 119, 261603 (2017)] for conformal field theories (CFTs) perturbed by a relevant operator is a universal function of the dimensionless coupling constant. The universal scaling of the Klein bottle entropy near criticality provides an efficient approach to extract the scaling dimension of lattice operators via data collapse. As paradigmatic examples, we validate the universal scaling of the Klein bottle entropy for Ising and Z3 parafermion CFTs with various perturbations using numerical simulation with continuous matrix product operator approach., Comment: 5+5 pages, 3 figures, published version

We develop a formalism for constructing particle-number-conserving Gaussian fermionic projected entangled pair states [U(1)-GfPEPS] and show that these states can describe ground states of band insulators and gapless fermions with band touching points. When using them as variational Ans\"{a}tze for two Dirac fermion systems ($\pi$-flux model on the square lattice and $[0,\pi]$-flux model on the kagome lattice), we find that the U(1)-GfPEPS, even with a relatively small bond dimension, can accurately approximate the Dirac Fermi sea ground states. By applying Gutzwiller projectors on top of these U(1)-GfPEPS, we obtain PEPS representation of U(1)-Dirac spin liquid states for spin-1/2 systems. With state-of-the-art tensor network numerics, the critical exponent in the spin-spin correlation function of the Gutzwiller-projected $\pi$-flux state is estimated to be $\eta \approx 1.7$.

We investigate the physics of projected d-wave pairing states using their fermionic projected entangled pair state (fPEPS) representation. First, we approximate a d-wave Bardeen-Cooper-Schrieffer state using the Gaussian fPEPS. Next, we translate the resulting state into fPEPS tensors and implement the Gutzwiller projection which removes double occupancy by modifying the local tensor elements. The tensor network representation of the projected d-wave pairing state allows us to evaluate physical quantities in the thermodynamic limit without employing the Gutzwiller approximation. Despite having very few variational parameters, such physically motivated tensor network states are shown to exhibit competitive energies for the doped t-J model. We expect that such construction offers useful initial states and guidance for variational tensor network calculations., Comment: 9 pages, 7 figures

We consider a chiral spin liquid constructed using the parton theory. This state supports non-Abelian anyons and neutral fermions, which share some similarities with the celebrated Moore-Read state. The edge physics of the parton state and the Moore-Read state is very different. Based on conformal field theory (CFT) analysis, it is proposed that the edge states exhibit an emergent SO(5) symmetry. The counting of edge states in the low-lying SO(5)$_{1}$ CFT towers is computed. The chiral spin liquid is generated using tensor network methods, which allows us to confirm the SO(5)$_{1}$ counting using entanglement spectrum. An additional feature of multiple branches in the entanglement spectrum is observed and analyzed.

We propose a series of models with inverse-square interactions, for which the ground states have an exact closed form, that describe one localized spin-$1/2$ Kondo impurity coupled to Luttinger liquids of itinerant bosons or fermions in the continuum or on a lattice. The model Hamiltonians contain two parts: a two-component Calogero-Sutherland model with open boundary condition in the continuum or a long-range $t \textendash J$-type model on the lattice, and long-range Kondo couplings between the impurity and the particles. The wave functions of the ground states in a certain basis have Jastrow product forms with powers depending on the particle statistics. It is shown that the wave functions can be expressed as chiral correlators of certain conformal field theories, which reveals that the particles are fractionalized to fully screen the Kondo impurity. This insight would be very useful in future studies on the very difficult problem of Kondo physics with interacting baths., Comment: (5.5 + 12.5) pages, (2 + 2) figures, published version

We introduce a many-body topological invariant, called the topological disorder parameter (TDP), to characterize gapped quantum phases with global internal symmetry in (2+1)d. TDP is defined as the constant correction that appears in the ground state expectation value of a partial symmetry transformation applied to a connected spatial region $M$, the absolute value of which scales generically as $\exp(-\alpha l+\gamma)$ where $l$ is the perimeter of $M$ and $\gamma$ is the TDP. Motivated by a topological quantum field theory interpretation of the operator, we show that $e^\gamma$ can be related to the quantum dimension of the symmetry defect, and provide a general formula for $\gamma$ when the entanglement Hamiltonian of the topological phase can be described by a (1+1)d conformal field theory (CFT). A special case of TDP is equivalent to the topological R\'enyi entanglement entropy when the symmetry is the cyclic permutation of the replica of the gapped phase. We then investigate several examples of lattice models of topological phases, both analytically and numerically, in particular when the assumption of having a CFT edge theory is not satisfied. We also consider an example of partial translation symmetry in Wen's plaquette model and show that the result can be understood using the edge CFT. Our results establish a new tool to detect quantum topological order., Comment: 21 pages, 6 figures; v2 published version

The nature of the ground state for the $S = 1/2$ kagome Heisenberg antiferromagnet (KHAF) has been elusive. We revisit this challenging problem and provide numerical evidence that its ground state might be a chiral spin liquid. Combining the density matrix renormalization group method and analytical analyses, we demonstrate that the previously observed chiral spin liquid phase in the KHAF with longer-range couplings is stable in a broader region. We characterize the nature of the ground state by computing energy derivatives, revealing ground-state degeneracy arising from spontaneous breaking of time-reversal symmetry, and targeting the semion sector. We further investigate the phase diagram in the vicinity of the KHAF and observe a $\sqrt{3}\times\sqrt{3}$ magnetically ordered phase and two valence-bond crystal phases.

We provide an efficient and accurate method for converting Hartree-Fock-Bogoliubov wave functions into matrix product states (MPSs). These wave functions, also known as Bogoliubov vacua, exhibit a peculiar entanglement structure that the eigenvectors of the reduced density matrix are also Bogoliubov vacua. We exploit this important feature to obtain their optimal MPS approximation and derive an explicit formula for corresponding MPS matrices. The performance of our method is benchmarked with the Kitaev chain and the Majorana-Hubbard model on the honeycomb lattice. The approach facilitates the applications of Hartree-Fock-Bogoliubov wave functions and is ideally suited for combining with the density-matrix renormalization group method., Comment: 5 pages, 3 figures, published version

We perform a tensor network simulation of the (1+1)-dimensional $O(3)$ nonlinear $\sigma$-model with $\theta=\pi$ term. Within the Hamiltonian formulation, this field theory emerges as the finite-temperature partition function of a modified quantum rotor model decorated with magnetic monopoles. Using the monopole harmonics basis, we derive the matrix representation for this modified quantum rotor model, which enables tensor network simulations. We employ our recently developed continuous matrix product operator method [Tang et al., Phys. Rev. Lett. 125, 170604 (2020)] to study the finite-temperature properties of this model and reveal its massless nature. The central charge as a function of the coupling constant is directly extracted in our calculations and compared with field theory predictions., Comment: Published version. 13 pages, 7 figures, code: https://github.com/TensorBFS/U1cMPO

We explore the topological defects of the critical three-state Potts spin system on the torus, Klein bottle and cylinder. A complete characterization is obtained by breaking down the Fuchs-Runkel-Schweigert construction of 2d rational CFT to the lattice setting. This is done by applying the strange correlator prescription to the recently obtained tensor network descriptions of string-net ground states in terms of bimodule categories [Lootens, Fuchs, Haegeman, Schweigert, Verstraete, SciPost Phys. 10, 053 (2021)]. The symmetries are represented by matrix product operators (MPO), as well as intertwiners between the diagonal tetracritical Ising model and the non-diagonal three-state Potts model. Our categorical construction lifts the global transfer matrix symmetries and intertwiners, previously obtained by solving Yang-Baxter equations, to MPO symmetries and intertwiners that can be locally deformed, fused and split. This enables the extraction of conformal characters from partition functions and yields a comprehensive picture of all boundary conditions., Comment: 71 pages, many figures, includes supplementary material

Model wave functions constructed from (1+1)D conformal field theory (CFT) have played a vital role in studying chiral topologically ordered systems. There usually exist multiple degenerate ground states when such states are placed on the torus. The common practice for dealing with this degeneracy within the CFT framework is to take a full correlator on the torus, which includes both holomorphic and antiholomorphic sectors, and decompose it into several conformal blocks. In this paper, we propose a pure chiral approach for the torus wave function construction. By utilizing the operator formalism, the wave functions are written as chiral correlators of holormorphic fields restricted to each individual topological sector. This method is not only conceptually much simpler, but also automatically provides us the anyon eigenbasis of the degenerate ground states (also known as the "minimally entangled states"). As concrete examples, we construct the full set of degenerate ground states for SO($n$)$_1$ and SU($n$)$_1$ chiral spin liquids on the torus, the former of which provide a complete wave function realization of Kitaev's sixteenfold way of anyon theories. We further characterize their topological orders by analytically computing the associated modular $S$ and $T$ matrices., Comment: 50 pages, 5 figures; published version, with note added and minor modifications

The simplest spin-orbital model can host a nematic spin-orbital liquid state on the triangular lattice. We provide clear evidence that the ground state of the SU(4) Kugel-Khomskii model on the triangular lattice can be well described by a "single" Gutzwiller projected wave function with an emergent parton Fermi surface, despite it exhibits strong finite-size effect in quasi-one-dimensional cylinders. The finite-size effect can be resolved by the fact that the parton Fermi surface consists of open orbits in the reciprocal space. Thereby, a stripy liquid state is expected in the two-dimensional limit, which preserves the SU(4) symmetry while breaks the translational symmetry by doubling the unit cell along one of the lattice vector directions. It is indicative that these stripes are critical and the central charge is $c=3$, in agreement with the SU(4)$_1$ Wess-Zumino-Witten conformal field theory. All these results are consistent with the Lieb-Schultz-Mattis-Oshikawa-Hastings theorem., Comment: 4.5 page main text + supplementary materials. Science Bulletin (2022)

In the physics of the Fractional Quantum Hall (FQH) effect, a zoo of Abelian topological phases can be obtained by varying the magnetic field. Aiming to reach the same phenomenology in spin-like systems, we propose a family of SU($N$)-symmetric models in the fundamental representation, on the square lattice with short-range interactions restricted to triangular units, a natural generalization for arbitrary $N$ of an SU($3$) model studied previously where time-reversal symmetry is broken explicitly. Guided by the recent discovery of SU($2$)$_1$ and SU($3$)$_1$ chiral spin liquids (CSL) on similar models we search for topological SU($N$)$_1$ CSL in some range of the Hamiltonian parameters via a combination of complementary numerical methods such as exact diagonalizations (ED), infinite density matrix renormalization group (iDMRG) and infinite Projected Entangled Pair State (iPEPS). Extensive ED on small (periodic and open) clusters up to $N=10$ and an innovative SU($N$)-symmetric version of iDMRG to compute entanglement spectra on (infinitely-long) cylinders in all topological sectors provide unambiguous signatures of the SU($N$)$_1$ character of the chiral liquids. An SU($4$)-symmetric chiral PEPS, constructed in a manner similar to its SU($2$) and SU($3$) analogs, is shown to give a good variational ansatz of the $N=4$ ground state, with chiral edge modes originating from the PEPS holographic bulk-edge correspondence. Finally, we discuss the possible observation of such Abelian CSL in ultracold atom setups where the possibility of varying $N$ provides a tuning parameter similar to the magnetic field in the physics of the FQH effect., Comment: 41 pages, 21 figures, references added, figure 2(j) modified, abstract and introduction revised

Spin-orbital liquids are quantum disordered states in systems with entangled spin and orbital degrees of freedom. We study exactly solvable spin-orbital models in two dimensions with selected Heisenberg-, Kitaev-, and $\Gamma$-type interactions, as well as external magnetic fields. These models realize a variety of spin-orbital-liquid phases featuring dispersing Majorana fermions with Fermi surfaces, nodal Dirac or quadratic band touching points, or full gaps. In particular, we show that Zeeman magnetic fields can stabilize nontrivial flux patterns and induce metamagnetic transitions between states with different topological character. Solvable nearest-neighbor biquadratic spin-orbital perturbations can be tuned to stabilize zero-energy flat bands. We discuss in detail the examples of $\mathrm{SO}(2)$- and $\mathrm{SO}(3)$-symmetric spin-orbital models on the square and honeycomb lattices, and use group-theoretical arguments to generalize to $\mathrm{SO}(\nu)$-symmetric models with arbitrary integer $\nu > 1$. These results extend the list of exactly solvable models with spin-orbital-liquid ground states and highlight the intriguing general features of such exotic phases. Our models are thus excellent starting points for more realistic modellings of candidate materials., Comment: 18 pages, 11 figures; v3 (expanded conclusion and appendix, additional minor changes)

We propose resonating valence bond wave functions for a spin-1 system on the torus that realize a non-Abelian chiral spin liquid. The wave functions take the form of infinite dimensional matrix product states constructed from conformal blocks of the $\mathrm{SO}(3)_{1}$ Wess-Zumino-Witten model. This means that they are lattice analogues of the bosonic Moore-Read state introduced in fractional quantum Hall systems. The topological order of this system is revealed by explicit construction of three-fold degenerate ground states and analytical computation of the modular S and T matrices., Comment: 11 pages, 2 figures

We study transitions between topological phases featuring emergent fractionalized excitations in two-dimensional models for Mott insulators with spin and orbital degrees of freedom. The models realize fermionic quantum critical points in fractionalized Gross-Neveu$^\ast$ universality classes in (2+1) dimensions. They are characterized by the same set of critical exponents as their ordinary Gross-Neveu counterparts, but feature a different energy spectrum, reflecting the nontrivial topology of the adjacent phases. We exemplify this in a square-lattice model, for which an exact mapping to a $t$-$V$ model of spinless fermions allows us to make use of large-scale numerical results, as well as in a honeycomb-lattice model, for which we employ $\epsilon$-expansion and large-$N$ methods to estimate the critical behavior. Our results are potentially relevant for Mott insulators with $d^1$ electronic configurations and strong spin-orbit coupling, or for twisted bilayer structures of Kitaev materials., Comment: 6+6 pages, 2+3 figures; v3 (minor changes, discussion on strong-coupling limit)

We propose to boost the performance of the density matrix renormalization group (DMRG) in two dimensions by using Gutzwiller projected states as the initialization ansatz. When the Gutzwiller projected state is properly chosen, the notorious "local minimum" issue in DMRG can be circumvented and the precision of DMRG can be improved by orders of magnitude without extra computational cost. Moreover, this method allows to quantify the closeness of the initial Gutzwiller projected state and the final converged state after DMRG sweeps, thereby sheds light on whether the Gutzwiller ansatz captures the essential entanglement features of the actual ground state for a given Hamiltonian. The Kitaev honeycomb model has been exploited to demonstrate and benchmark this new method., Comment: 4.5 pages of main text + 9 pages of supplementary material; Fig. 3 in main text is updated

In two dimensions, the topological order described by $\mathbb{Z}_2$ gauge theory coupled to free or weakly interacting fermions with a nonzero spectral Chern number $\nu$ is classified by $\nu \; \mathrm{mod}\; 16$ as predicted by Kitaev [Ann. Phys. 321, 2 (2006)]. Here we provide a systematic and complete construction of microscopic models realizing this so-called sixteenfold way of anyon theories. These models are defined by $\Gamma$ matrices satisfying the Clifford algebra, enjoy a global $\mathrm{SO}(\nu)$ symmetry, and live on either square or honeycomb lattices depending on the parity of $\nu$. We show that all these models are exactly solvable by using a Majorana representation and characterize the topological order by calculating the topological spin of an anyonic quasiparticle and the ground-state degeneracy. The possible relevance of the $\nu=2$ and $\nu=3$ models to materials with Kugel-Khomskii-type spin-orbital interactions is discussed., Comment: 6+9 pages, 2+1 figures, published version

We present an algorithm for studying quantum systems at finite temperature using continuous matrix product operator representation. The approach handles both short-range and long-range interactions in the thermodynamic limit without incurring any time discretization error. Moreover, the approach provides direct access to physical observables including the specific heat, local susceptibility, and local spectral functions. After verifying the method using the prototypical quantum XXZ chains, we apply it to quantum Ising models with power-law decaying interactions and on the infinite cylinder, respectively. The approach offers predictions that are relevant to experiments in quantum simulators and the nuclear magnetic resonance spin-lattice relaxation rate., Comment: 4.5+6 pages, 3+6 figures, published version, code: https://github.com/TensorBFS/cMPO

Recent work by Wu {\em et al.} [arXiv:1910.11011] proposed a numerical method, so-called matrix product operator-matrix product state (MPO-MPS) method, by which several types of quantum many-body wave functions, in particular, the projected Fermi sea state, can be efficiently represented as a tensor network. In this paper, we generalize the MPO-MPS method to study Gutzwiller projected paired states of fermions, where the maximally localized Wannier orbitals for Bogoliubov quasiparticles/quasiholes have been adapted to improve the computational performance. The study of $SO(3)$-symmetric spin-1 chains reveals that this new method has better performance than variational Monte Carlo for gapped states and similar performance for gapless states. Moreover, we demonstrate that dynamic correlation functions can be easily evaluated by this method cooperating with other MPS-based accurate approaches, such as the Chebyshev MPS method., Comment: 12 pages, 7 figures, 2 tables, published version

We perform the state-of-the-art tensor network simulations directly in the thermodynamic limit to clarify the critical properties of the $q$-state clock model on the square lattice. We determine accurately the two phase transition temperatures through the singularity of the classical analog of the entanglement entropy, and provide extensive numerical evidences to show that both transitions are of the Berezinskii-Kosterlitz-Thouless (BKT) type for $q\ge 5$ and that the low-energy physics of this model is well described by the $\mathbb{Z}_q$-deformed sine-Gordon theory. We also determine the characteristic conformal parameters, especially the compactification radius, that govern the critical properties of the intermediate BKT phase., Comment: 5 pages, 5 figures plus supplemental material; accepted for publication in Phys. Rev. E (Rapid Com.), minor typo corrections

We study two-component bosons on the Harper-Hofstadter model with two legs. The synthetic magnetic fields for the two types of bosons point to either the same direction or opposite directions. The bosons have hardcore intra-species interaction such that there can be no more than one boson of the same type on each lattice site. For certain filling factors in the absence of inter-species interaction, each component realizes a vortex Mott insulator with rung current or a Meissner superfluid without rung current. The system undergoes phase transitions to other phases as inter-species interaction is turned on, which are characterized numerically using the density matrix renormalization group method and supplemented with analytical studies when possible. The vortex Mott insulator transits to a gapped Meissner phase without rung current and the Meissner superfluid transits to a gapped vortex phase with rung current. In both cases, we observe gapped spin density wave states that break certain ${\mathbb Z}_{2}$ symmetries., Comment: 11 pages, 10 figures

Tensor network states and parton wave functions are two pivotal methods for studying quantum many-body systems. This work connects these two subjects as we demonstrate that a variety of parton wave functions, such as projected Fermi sea and projected fermionic or bosonic paired states, can be represented exactly as tensor networks. The results can be compressed into matrix product states with moderate bond dimensions so various physical quantities can be computed efficiently. For the projected Fermi sea, we develop an excellent compression scheme with high fidelity using maximally localized Wannier orbitals. Numerical calculations on two parton wave functions demonstrate that our method exceeds commonly adopted Monte Carlo methods in some aspects. It produces energy and correlation function with very high accuracy that is difficult to achieve using Monte Carlo method. The entanglement measures that were almost impossible to compute before can also be obtained easily using our method., Comment: 6+4 pages, 3+1 figures, published version

The coupling between localized magnetic moments and itinerant electrons presents a plethora of interesting physics. The low-energy physics of some quantum impurity systems can be described using conformal field theory (CFT). In this paper, the connection between quantum impurity models and CFT is further strengthened as we construct a class of exactly solvable models with ground states given by CFT correlators. The method developed here is completely analytical and can be applied to fermions with an arbitrary number of colors and multiple impurities. Numerical calculations are performed to characterize certain aspects of our models for which we do not have analytical results., Comment: 25 pages

The quantized thermal Hall effect is an important probe for detecting chiral topological order and revealing the nature of chiral gapless edge states. The standard Kubo formula approach for the thermal Hall conductance $\kappa_{xy}$ based on the linear-response theory faces difficulties in practical application due to the lack of a reliable numerical method for calculating dynamical quantities in microscopic models at finite temperature. In this work, we propose an approach for calculating $\kappa_{xy}$ in two-dimensional lattice models displaying chiral topological order. Our approach targets at the edge current localized at the boundary which involves only thermal averages of local operators in equilibrium, thus drastically lowering the barrier for the calculation of $\kappa_{xy}$. We use the chiral $p$-wave superconductor (with and without disorder) and the Hofstadter model as benchmark examples to illustrate several sources of finite-size effects, and we suggest the infinite (or sufficiently long) strip as the best geometry for carrying out numerical simulations., Comment: 12 pages, 9 figures, published version

We propose and study an SU(3) trimer resonating-valence-bond (tRVB) state with $C_{4v}$ point-group symmetry on the square lattice. By devising a projected entangled-pair state representation, we show that all (connected) correlation functions between local operators in this SU(3) tRVB state decay exponentially, indicating its gapped nature. We further calculate the modular $S$ and $T$ matrices by constructing all nine topological sectors on a torus and establish the existence of $\mathbb{Z}_3$ topological order in this SU(3) tRVB state., Comment: 6 pages, 6 figures

We extend the scope of the Klein bottle entropy, originally introduced by [Tu, Phys. Rev. Lett. 119, 261603 (2017)] in the rational conformal field theory (CFT), to the compactified boson CFT, which are relevant to the studies of Luttinger liquids. We first review the Klein bottle entropy in rational CFT and discuss details of how to extract the Klein bottle entropy from lattice models using the example of the transverse field Ising model. We then go beyond the scope of rational CFT and study the Klein bottle entropy $\ln g$ in the compactified boson CFT, which turns out to have a straightforward relation to the compactification radius $R$, $\ln g = \ln R$. This relation indicates a convenient and efficient method to extract the Luttinger parameter from lattice model calculations. Our numerical results on the Klein bottle entropy in the spin-$1/2$ XXZ chain show excellent agreement with the CFT predictions, up to some small deviations near the isotropic point, which we attribute to the marginally irrelevant terms. For the $S = 1$ XXZ chain that cannot be exactly solved, our numerical results provide an accurate numerical determination of the Luttinger parameter in this model., Comment: 18 pages, 6 figures, published version

We construct an exactly solvable quantum impurity model which consists of spin-1/2 conduction fermions and the spin-1/2 magnetic moment. The ground state is a Gutzwiller projected Fermi sea with non-orthonormal modes and its wave function in the site-occupation basis is a Jastrow-type homogeneous polynomial. The parent Hamiltonian has all-to-all inverse-square hopping terms between the conduction fermions and inverse-square spin-exchange terms between the conduction fermions and the magnetic moments. The low-lying energy levels, spin-spin correlation function, and von Neumann entanglement entropy of our model demonstrate that it exhibits the essential aspects of spin-1/2 Kondo physics. The machinery developed in this work can generate many other exactly solvable quantum impurity models., Comment: 6+9 pages, 4 figures, published version

Recently, entropy corrections on nonorientable manifolds such as the Klein bottle are proposed as a universal characterization of critical systems with an emergent conformal field theory (CFT). We show that entropy correction on the Klein bottle can be interpreted as a boundary effect via transforming the Klein bottle into an orientable manifold with nonlocal boundary interactions. The interpretation reveals the conceptual connection of the Klein bottle entropy with the celebrated Affleck-Ludwig entropy in boundary CFT. We propose a generic scheme to extract these universal boundary entropies from quantum Monte Carlo calculation of partition function ratios in lattice models. Our numerical results on the Affleck-Ludwig entropy and Klein bottle entropy for the $q$-state quantum Potts chains with $q=2,3$ show excellent agreement with the CFT predictions. For the quantum Potts chain with $q=4$, the Klein bottle entropy slightly deviates from the CFT prediction, which is possibly due to marginally irrelevant terms in the low-energy effective theory., Comment: 10 pages, 4 figures. Published version

We show that rational conformal field theories in 1+1 dimensions on a Klein bottle, with length $L$ and width $\beta$, satisfying $L \gg \beta$, have a universal entropy. This universal entropy is a topological invariant depending on the quantum dimensions of the primary fields and can be accurately extracted by taking a proper ratio between the Klein bottle and torus partition functions, enabling a characterization of conformal critical theories. The result is checked against exact calculations in quantum spin-1/2 XY and Ising chains., Comment: 5 pages, 2 figures, published version

We introduce a family of many-body quantum states that describe interacting spin one-half hard-core particles with bosonic or fermionic statistics on arbitrary one- and two-dimensional lattices. The wave functions at lattice filling fraction $\nu=2/(2m+1)$ are derived from deformations of the Wess-Zumino-Witten model $\mathfrak{su}(3)_1$ and are related to the $(m+1,m+1,m)$ Halperin fractional quantum Hall states. We derive long-range SU(2) invariant parent Hamiltonians for these states which in two dimensions are chiral $t$-$J$-$V$ models with additional three-body interaction terms. In one dimension we obtain a generalisation to open chains of a periodic inverse-square $t$-$J$-$V$ model proposed in [Z. N. C. Ha and F. D. M. Haldane, Phys. Rev. B $\mathbf{46}$, 9359 (1992)]. We observe that the gapless low-energy spectrum of this model and its open-boundary generalisation can be described by rapidity sets with the same generalised Pauli exclusion principle. A two-component compactified free boson conformal field theory is identified that has the same central charge and scaling dimensions as the periodic bosonic inverse-square $t$-$J$-$V$ model., Comment: 19 pages, 2 figures. v2: minor corrections and partial rewriting of section IV B 2

Motivated by a recent experiment which synthesizes Landau levels for photons on cones [Schine {\em et al.}, Nature 534, 671 (2016)], and more generally the interest in understanding gravitational responses of quantum Hall states, we study fractional quantum Hall states of bosons on cones. A variety of trial wave functions for conical systems are constructed and compared with exact diagonalization results. The tip of a cone is a localized geometrical defect with singular curvature which can modify the density profiles of quantum Hall states. The density profiles on cones can be used to extract some universal information about quantum Hall states. The values of certain quantities are computed numerically using the density profiles of some quantum Hall states and they agree with analytical predictions., Comment: 9.5 pages, 8 figures

We introduce a general framework to construct many-body lattice wavefunctions starting from the conformal blocks (CBs) of rational conformal field theories (RCFTs). We discuss the different ways of encoding the physical degrees of freedom of the lattice system using both the internal symmetries of the theory and the fusion channels of the CBs. We illustrate this construction both by revisiting the known Haldane-Shastry model and by providing a novel implementation for the Ising RCFT. In the latter case, we find a connection to the Ising transverse field (ITF) spin chain via the Kramers-Wannier duality and the Temperley-Lieb-Jones algebra. We also find evidence that the ground state of the finite-size critical ITF Hamiltonian corresponds exactly to the wavefunction obtained from CBs of spin fields.

We propose an SU(3) symmetric Hamiltonian with short-range interactions on the Kagome lattice and show that it hosts an Abelian chiral spin liquid (CSL) state. We provide numerical evidence based on exact diagonalization to show that this CSL state is stabilized in an extended region of the parameter space and can be viewed as a lattice version of the Halperin 221 fractional quantum Hall (FQH) state of two-component bosons. We also construct a parton wave function for this CSL state and demonstrate that its variational energies are in good agreement with exact diagonalization results. The parton description further supports that the CSL is characterized by a chiral edge conformal field theory (CFT) of the SU(3)_1 Wess-Zumino-Witten type., Comment: 5 pages, 4 figures; v4: typo in the caption of Fig. 3 in the published version corrected

Symmetry-protected trivial (SPt) phases of matter are the product-state analogue of symmetry-protected topological (SPT) phases. This means, SPt phases can be adiabatically connected to a product state by some path that preserves the protecting symmetry. Moreover, SPt and SPT phases can be adiabatically connected to each other when interaction terms that break the symmetries protecting the SPT order are added in the Hamiltonian. It is also known that spin-1 SPT phases in quantum spin chains can emerge as effective intermediate phases of spin-2 Hamiltonians. In this paper we show that a similar scenario is also valid for SPt phases. More precisely, we show that for a given spin-2 quantum chain, effective intermediate spin-1 SPt phases emerge in some regions of the phase diagram, these also being adiabatically connected to non-trivial intermediate SPT phases. We characterize the phase diagram of our model by studying quantities such as the entanglement entropy, symmetry-related order parameters, and 1-site fidelities. Our numerical analysis uses Matrix Product States (MPS) and the infinite Time-Evolving Block Decimation (iTEBD) method to approximate ground states of the system in the thermodynamic limit. Moreover, we provide a field theory description of the possible quantum phase transitions between the SPt phases. Together with the numerical results, such a description shows that the transitions may be described by Conformal Field Theories (CFT) with central charge c=1. Our results are in agreement, and further generalize, those in [Y. Fuji, F. Pollmann, M. Oshikawa, Phys. Rev. Lett. 114, 177204 (2015)]., Comment: 7 pages, 5 figures, 1 table, revised version. Accepted in PRB

We study lattice wave functions obtained from the SU(2)$_1$ Wess-Zumino-Witten conformal field theory. Following Moore and Read's construction, the Kalmeyer-Laughlin fractional quantum Hall state is defined as a correlation function of primary fields. By an additional insertion of Kac-Moody currents, we associate a wave function to each state of the conformal field theory. These wave functions span the complete Hilbert space of the lattice system. On the cylinder, we study global properties of the lattice states analytically and correlation functions numerically using a Metropolis Monte Carlo method. By comparing short-range bulk correlations, numerical evidence is provided that the states with one current operator represent edge states in the thermodynamic limit. We show that the edge states with one Kac-Moody current of lowest order have a good overlap with low-energy excited states of a local Hamiltonian, for which the Kalmeyer-Laughlin state approximates the ground state. For some states, exact parent Hamiltonians are derived on the cylinder. These Hamiltonians are SU(2) invariant and nonlocal with up to four-body interactions., Comment: 17 pages, 8 figures, 3 tables; v4: error bars and fig. 6 corrected + minor corrections

We study $S=1$ spin liquid states on the kagome lattice constructed by Gutzwiller-projected $p_x+ip_y$ superconductors. We show that the obtained spin liquids are either non-Abelian or Abelian topological phases, depending on the topology of the fermionic mean-field state. By calculating the modular matrices $S$ and $T$, we confirm that projected topological superconductors are non-Abelian chiral spin liquid (NACSL). The chiral central charge and the spin Hall conductance we obtained agree very well with the $SO(3)_1$ (or, equivalently, $SU(2)_2$) field theory predictions. We propose a local Hamiltonian which may stabilize the NACSL. From a variational study we observe a topological phase transition from the NACSL to the $Z_2$ Abelian spin liquid., Comment: 12 pages, 7 figures, 1 table

We consider a spin-1/2 fermionic ladder with spin-orbit coupling and a perpendicular magnetic field, which shares important similarities with topological superconducting wires. We fully characterize the symmetry-protected topological phase of this ladder through the identification of fractionalized edge modes and non-trivial spin winding numbers. We propose an experimental scheme to engineer such a ladder system with cold atoms in optical lattices, and we present two protocols that can be used to extract the topological signatures from density and momentum-distribution measurements. We then consider the presence of interactions and discuss the effects of a contact on-site repulsion on the topological phase. We find that such interactions could enhance the extension of the topological phase in certain parameters regimes., Comment: 27 pages, 10 figures, 1 appendix

Unlike the random radial orientation distribution of field elliptical galaxies, galaxies in a cluster are expected to point preferentially towards the center of the cluster, as a result of the cluster's tidal force on its member galaxies. In this work an analytic model is formulated to simulate this effect. The deformation time scale of a galaxy in a cluster is usually much shorter than the time scale of change of the tidal force; the dynamical process of the tidal interaction within the galaxy can thus be ignored. An equilibrium shape of a galaxy is then assumed to be the surface of equipotential, which is the sum of the self-gravitational potential of the galaxy and the tidal potential of the cluster at this location. We use a Monte-Carlo method to calculate the radial orientation distribution of these galaxies, by assuming the NFW mass profile of the cluster and the initial ellipticity of field galaxies. The radial angles show a single peak distribution centered at zero. The Monte-Carlo simulations also show that a shift of the reference center from the real cluster center weakens the anisotropy of the radial angle distribution. Therefore, the expected radial alignment cannot be revealed if the distribution of spatial position angle is used instead of that of radial angle. The observed radial orientations of elliptical galaxies in cluster Abell~2744 are consistent with the simulated distribution., Comment: 8 pages, 6 figures, 2 tables. MNRAS in press

$f(R)$ gravity is one of the simplest theories of modified gravity to explain the accelerated cosmic expansion. Although it is usually assumed that the quasi-Newtonian approach (a combination of the quasi-static approximation and sub-Hubble limit) for cosmic perturbations is good enough to describe the evolution of large scale structure in $f(R)$ models, some studies have suggested that this method is not valid for all $f(R)$ models. Here, we show that in the matter-dominated era, the pressure and shear equations alone, which can be recast into four first-order equations to solve for cosmological perturbations exactly, are sufficient to solve for the Newtonian potential, $\Psi$, and the curvature potential, $\Phi$. Based on these two equations, we are able to clarify how the exact linear perturbations fit into different limits. We find that the Compton length controls the quasi-static behaviours in $f(R)$ gravity. In addition, regardless the validity of quasi-static approximation, a strong version of the sub-Hubble limit alone is sufficient to reduce the exact linear perturbations in any viable $f(R)$ gravity to second order. Our findings disagree with some previous studies where we find little difference between our exact and quasi-Newtonian solutions even up to $k=10 c^{-1} \mathcal{H}_0$.

We show that infinite Matrix Product States (MPS) constructed from conformal field theories can describe ground states of one-dimensional critical systems with open boundary conditions. To illustrate this, we consider a simple infinite MPS for a spin-1/2 chain and derive an inhomogeneous open Haldane-Shastry model. For the spin-1/2 open Haldane-Shastry model, we derive an exact expression for the two-point spin correlation function. We also provide an SU($n$) generalization of the open Haldane-Shastry model and determine its twisted Yangian generators responsible for the highly degenerate multiplets in the energy spectrum., Comment: 5+7 pages, 4 figures, published version, a typo in the twisted Yangian generators corrected (thanks to the authors of arXiv:1511.08613 for pointing out this typo)

We develop a method of constructing excited states in one dimensional spin chains which are derived from the $SU(2)_1$ Wess-Zumino-Witten Conformal Field Theory (CFT) using a parent Hamiltonian approach. The resulting systems are equivalent to the Haldane-Shastry model. In our ansatz, correlation functions between primary fields correspond to the ground state of the spin system, whereas excited states are obtained by insertion of descendant fields. Our construction is based on the current algebra of the CFT and emphasizes the close relation between the spectrum of the spin system and the underlying CFT. This general structure might imply that the method could be applied to a wider range of model systems., Comment: 15 pages, 4 figures, v2: significantly improved, v3: accepted version

We perform a systematic investigation on the hexagon-singlet solid (HSS) states, which are a class of spin liquid candidates for the spin-1 kagome antiferromagnet. With the Schwinger boson representation, we show that all HSS states have exponentially decaying correlations and can be interpreted as a (special) subset of the resonating Affleck-Kennedy-Lieb-Tasaki (AKLT) loop states. We provide a compact tensor network representation of the HSS states, with which we are able to calculate physical observables efficiently. We find that the HSS states have vanishing topological entanglement entropy, suggesting the absence of intrinsic topological order. We also employ the HSS states to perform a variational study of the spin-1 kagome Heisenberg antiferromagnetic model. When we use a restricted HSS ansatz preserving lattice symmetry, the best variational energy per site is found to be $e_0 = -1.3600$. In contrast, when allowing lattice symmetry breaking, we find a trimerized simplex valence bond crystal with a lower energy, $e_0=-1.3871$., Comment: 14 pages, 12 figures, published version

Here we study the emergence of different Symmetry-Protected Topological (SPT) phases in a spin-2 quantum chain. We consider a Heisenberg-like model with bilinear, biquadratic, bicubic, and biquartic nearest-neighbor interactions, as well as uniaxial anisotropy. We show that this model contains four different effective spin-1 SPT phases, corresponding to different representations of the $(\mathbb{Z}_2 \times \mathbb{Z}_2) + T$ symmetry group, where $\mathbb{Z}_2$ is some $\pi$-rotation in the spin internal space and $T$ is time-reversal. One of these phases is equivalent to the usual spin-1 Haldane phase, while the other three are different but also typical of spin-1 systems. The model also exhibits an $SO(5)$-Haldane phase. Moreover, we also find that the transitions between the different effective spin-1 SPT phases are continuous, and can be described by a $c=2$ conformal field theory. At such transitions, indirect evidence suggests a possible effective field theory of four massless Majorana fermions. The results are obtained by approximating the ground state of the system in the thermodynamic limit using Matrix Product States via the infinite Time Evolving Block Decimation method, as well as by effective field theory considerations. Our findings show, for the first time, that different large effective spin-1 SPT phases separated by continuous quantum phase transitions can be stabilized in a simple quantum spin chain., Comment: 7 pages, 6 figures, revised version. To appear in PRB