Source author record

Olexei I. Motrunich

Olexei I. Motrunich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

Infinite randomness with continuously varying critical exponents in the random XYZ spin chain

We study the antiferromagnetic XYZ spin chain with quenched bond randomness, focusing on a critical line between localized Ising magnetic phases. A previous calculation using the spectrum-bifurcation renormalization group, and assuming marginal many-body localization, proposed that critical indices vary continuously. In this work we solve the low-energy physics using an unbiased numerically exact tensor network method named the "rigorous renormalization group." We find a line of fixed points consistent with infinite-randomness phenomenology, with indeed continuously varying critical exponents for average spin correlations. A self-consistent Hartree-Fock-type treatment of the $z$ couplings as interactions added to the free-fermion random XY model captures much of the important physics including the varying exponents; we provide an understanding of this as a result of local correlation induced between the mean-field couplings. We solve the problem of the locally-correlated XY spin chain with arbitrary degree of correlation and provide analytical strong-disorder renormalization group proofs of continuously varying exponents based on an associated classical random walk problem. This is also an example of a line of fixed points with continuously varying exponents in the equivalent disordered free-fermion chain. We argue that this line of fixed points also controls an extended region of the critical interacting XYZ spin chain.

preprint2021arXiv

Hilbert Space Fragmentation and Commutant Algebras

We study the phenomenon of Hilbert space fragmentation in isolated Hamiltonian and Floquet quantum systems using the language of commutant algebras, the algebra of all operators that commute with each term of the Hamiltonian or each gate of the circuit. We provide a precise definition of Hilbert space fragmentation in this formalism as the case where the dimension of the commutant algebra grows exponentially with the system size. Fragmentation can hence be distinguished from systems with conventional symmetries such as $U(1)$ or $SU(2)$, where the dimension of the commutant algebra grows polynomially with the system size. Further, the commutant algebra language also helps distinguish between "classical" and "quantum" Hilbert space fragmentation, where the former refers to fragmentation in the product state basis. We explicitly construct the commutant algebra in several systems exhibiting classical fragmentation, including the $t-J_z$ model and the spin-1 dipole-conserving model, and we illustrate the connection to previously-studied "Statistically Localized Integrals of Motion" (SLIOMs). We also revisit the Temperley-Lieb spin chains, including the spin-1 biquadratic chain widely studied in the literature, and show that they exhibit quantum Hilbert space fragmentation. Finally, we study the contribution of the full commutant algebra to the Mazur bounds in various cases. In fragmented systems, we use expressions for the commutant to analytically obtain new or improved Mazur bounds for autocorrelation functions of local operators that agree with previous numerical results. In addition, we are able to rigorously show the localization of the on-site spin operator in the spin-1 dipole-conserving model.

preprint2020arXiv

Eta-pairing states as true scars in an extended Hubbard Model

The eta-pairing states are a set of exactly known eigenstates of the Hubbard model on hypercubic lattices, first discovered by Yang [Phys. Rev. Lett. 63, 2144 (1989)]. These states are not many-body scar states in the Hubbard model because they occupy unique symmetry sectors defined by the so-called "eta-pairing SU(2)" symmetry. We study an extended Hubbard model with bond-charge interactions, popularized by Hirsch [Physica C 158, 326 (1989)], where the eta-pairing states survive without the eta-pairing symmetry and become true scar states. We also discuss similarities between the eta-pairing states and exact scar towers in the spin-1 XY model found by Schecter and Iadecola [Phys. Rev. Lett. 123, 147201 (2019)], and systematically arrive at all nearest-neighbor terms that preserve such scar towers in 1D. We also generalize these terms to arbitrary bipartite lattices. Our study of the spin-1 XY model also leads us to several new scarred models, including a spin-1/2 $J_1-J_2$ model with Dzyaloshinkskii-Moriya interaction, in realistic quantum magnet settings in 1D and 2D.

preprint2020arXiv

New exact eigenstates in the Lesanovsky model, proximity to integrability and the PXP model, and approximate scar states

We study a model of Rydberg atoms in a nearest-neighbor Rydberg blockaded regime, introduced by Lesanovsky in Phys. Rev. Lett. 108, 105301 (2012). This many-body model (which has one parameter $z$) has an exactly known gapped liquid ground state, and two exactly known low-lying excitations. We discover two new exact low-lying eigenstates. We also discuss behavior of the model at small parameter $z$ and its proximity to an integrable model. Lastly, we discuss connections between the Lesanovsky model at intermediate $z$ and so-called PXP model. The PXP model describes a recent experiment that observed unusual revivals from a charge density wave initial state, which are attributed to a set of many-body "scar states" which do not obey the eigenstate thermalization hypothesis. We discuss the possibility of scar states in the Lesanovsky model and present two approximations for them.

preprint2020arXiv

Slow Thermalization of Exact Quantum Many-Body Scar States Under Perturbations

Quantum many-body scar states are exceptional finite energy density eigenstates in an otherwise thermalizing system that do not satisfy the eigenstate thermalization hypothesis. We investigate the fate of exact many-body scar states under perturbations. At small system sizes, deformed scar states described by perturbation theory survive. However, we argue for their eventual thermalization in the thermodynamic limit from the finite-size scaling of the off-diagonal matrix elements. Nevertheless, we show numerically and analytically that the nonthermal properties of the scars survive for a parametrically long time in quench experiments. We present a rigorous argument that lower-bounds the thermalization time for any scar state as $t^{*} \sim O(λ^{-1/(1+d)})$, where $d$ is the spatial dimension of the system and $λ$ is the perturbation strength.

preprint2020arXiv

Unified structure for exact towers of scar states in the AKLT and other models

Quantum many-body scar states are many-body states with finite energy density in non-integrable models that do not obey the eigenstate thermalization hypothesis. Recent works have revealed "towers" of scar states that are exactly known and are equally spaced in energy, specifically in the AKLT model, the spin-1 XY model, and a spin-1/2 model that conserves number of domain walls. We provide a common framework to understand and prove known exact towers of scars in these systems, by evaluating the commutator of the Hamiltonian and a ladder operator. In particular we provide a simple proof of the scar towers in the integer-spin 1d AKLT models by studying two-site spin projectors. Through this picture we deduce a family of Hamiltonians that share the scar tower with the AKLT model, and also find common parent Hamiltonians for the AKLT and XY model scars. We also introduce new towers of exact states, organized in a "pyramid" structure, in the spin-1/2 model through successive application of a non-local ladder operator.

preprint2016arXiv

Bosonic Analogue of Dirac Composite Fermi Liquid

We introduce a particle-hole-symmetric metallic state of bosons in a magnetic field at odd-integer filling. This state hosts composite fermions whose energy dispersion features a quadratic band touching and corresponding $2π$ Berry flux protected by particle-hole and discrete rotation symmetries. We also construct an alternative particle-hole symmetric state---distinct in the presence of inversion symmetry---without Berry flux. As in the Dirac composite Fermi liquid introduced by Son, breaking particle-hole symmetry recovers the familiar Chern-Simons theory. We discuss realizations of this phase both in 2D and on bosonic topological insulator surfaces, as well as signatures in experiments and simulations.

preprint2016arXiv

Entanglement entropy of composite Fermi liquid states on the lattice: In support of the Widom formula

Quantum phases characterized by surfaces of gapless excitations are known to violate the otherwise ubiquitous boundary law of entanglement entropy in the form of a multiplicative log correction: $S\sim L^{d-1} \log L$. Using variational Monte Carlo, we calculate the second Rényi entropy for a model wavefunction of the $ν=1/2$ composite Fermi liquid (CFL) state defined on the two-dimensional triangular lattice. By carefully studying the scaling of the total Rényi entropy and, crucially, its contributions from the modulus and sign of the wavefunction on various finite-size geometries, we argue that the prefactor of the leading $L \log L$ term is equivalent to that in the analogous free fermion wavefunction. In contrast to the recent results of Shao et al. [PRL 114, 206402 (2015)], we thus conclude that the "Widom formula" holds even in this non-Fermi liquid CFL state. More generally, our results further elucidate---and place on a more quantitative footing---the relationship between nontrivial wavefunction sign structure and $S\sim L \log L$ entanglement scaling in such highly entangled gapless phases.

preprint2016arXiv

Explicit derivation of duality between a free Dirac cone and quantum electrodynamics in (2+1) dimensions

We explicitly derive the duality between a free electronic Dirac cone and quantum electrodynamics in $(2+1)$ dimensions (QED$_3$) with $N = 1$ fermion flavors. The duality proceeds via an exact, non-local mapping from electrons to dual fermions with long-range interactions encoded by an emergent gauge field. This mapping allows us to construct parent Hamiltonians for exotic topological-insulator surface phases, derive the particle-hole-symmetric field theory of a half-filled Landau level, and nontrivially constrain QED$_3$ scaling dimensions. We similarly establish duality between bosonic topological insulator surfaces and $N = 2$ QED$_3$.

preprint2016arXiv

Vacancy-induced low-energy states in undoped graphene

We demonstrate that a nonzero concentration $n_v$ of static, randomly-placed vacancies in graphene leads to a density $w$ of zero-energy quasiparticle states at the band-center $ε=0$ within a tight-binding description with nearest-neighbour hopping $t$ on the honeycomb lattice. We show that $w$ remains generically nonzero in the compensated case (exactly equal number of vacancies on the two sublattices) even in the presence of hopping disorder, and depends sensitively on $n_v$ and correlations between vacancy positions. For low, {\em but not-too-low} $|ε|/t$ in this compensated case, we show that the density of states (DOS) $ρ(ε)$ exhibits a strong divergence of the form $ρ_{\rm 1D}(ε) \sim |ε|^{-1}/ [\log(t/|ε|)]^{(y+1)} $, which crosses over to the universal low-energy asymptotic form expected on symmetry grounds $ρ_{\rm GW}(ε) \sim |ε|^{-1}e^{-b[\log(t/|ε|)]^{2/3} }$ below a crossover scale $ε_c \ll t$. $ε_c$ is found to decrease rapidly with decreasing $n_v$, while $y$ decreases much more slowly.

preprint2015arXiv

A continuous Mott transition between a metal and a quantum spin liquid

More than half a century after first being proposed by Sir Nevill Mott, the deceptively simple question of whether the interaction-driven electronic metal-insulator transition may be continuous remains enigmatic. Recent experiments on two-dimensional materials suggest that when the insulator is a quantum spin liquid, lack of magnetic long-range order on the insulating side may cause the transition to be continuous, or only very weakly first order. Motivated by this, we study a half-filled extended Hubbard model on a triangular lattice strip geometry. We argue, through use of large-scale numerical simulations and analytical bosonization, that this model harbors a continuous (Kosterlitz-Thouless-like) quantum phase transition between a metal and a gapless spin liquid characterized by a spinon Fermi surface, i.e., a "spinon metal." These results may provide a rare insight into the development of Mott criticality in strongly interacting two-dimensional materials and represent one of the first numerical demonstrations of a Mott insulating quantum spin liquid phase in a genuinely electronic microscopic model.

preprint2015arXiv

Monte Carlo study of phase transitions out of Symmetry-Enriched Topological phases of bosons in two dimensions

We study a statistical mechanics model of two species of bosons with mutual statistics $θ=2π/n$ in (2+1) dimensions. This model realizes a fractionalized topological phase of bosons, which is a fractionalized version of the boson integer quantum Hall effect. The model can be studied with sign-free Monte Carlo simulations. We study the phase transitions between the fractionalized topological phase and a trivial insulator, and between different topological phases. We find that these transitions are continuous, and we measure their critical exponents.

preprint2015arXiv

The half-filled Landau level: the case for Dirac composite fermions

One of the most spectacular experimental findings in the fractional quantum Hall effect is evidence for an emergent Fermi surface when the electron density is nearly half the density of magnetic flux quanta ($ν= 1/2$). The seminal work of Halperin, Lee, and Read (HLR) first predicted that at $ν= 1/2$ composite fermions--bound states of an electron and a pair of vortices--experience zero net magnetic field and can form a "composite Fermi liquid" with an emergent Fermi surface. In this paper we use infinite cylinder DMRG to provide compelling numerical evidence for the existence of a Fermi sea of composite fermions for realistic interactions between electrons at $ν= 1/2$. Moreover, we show that the state is particle-hole symmetric, in contrast to the construction of HLR. Instead, our findings are consistent if the composite fermions are massless Dirac particles, at finite density, similar to the surface state of a 3D topological insulator. Exploiting this analogy we devise a numerical test and successfully observe the suppression of $2k_F$ backscattering characteristic of Dirac particles.

preprint2014arXiv

Exact Models for Symmetry-Protected Topological Phases in One Dimension

We present an exactly solvable model for one-dimensional symmetry-protected topological phases with $\mathbb{Z}_N\times\mathbb{Z}_N$ symmetry. The model works by binding point topological defects (domain walls) of one symmetry to charges of the other and condensing these bound states. Binding single topological defects to charges leads to symmetry-protected topological phases, while binding multiple topological defects to charges leads to phases with a combination of symmetry-breaking and topological properties.

preprint2014arXiv

Model of Fractionalization of Faraday Lines in Compact Electrodynamics

Motivated by ideas of fractionalization and intrinsic topological order in bosonic models with short-range interactions, we consider similar phenomena in formal lattice gauge theory models. Specifically, we show that a compact quantum electrodynamics (CQED) can have, besides the familiar Coulomb and confined phases, additional unusual confined phases where excitations are quantum lines carrying fractions of the elementary unit of electric field strength. We construct a model that has $N$-tupled monopole condensation and realizes $1/N$ fractionalization of the quantum Faraday lines. This phase has another excitation which is a $Z_N$ quantum surface in spatial dimensions five and higher, but can be viewed as a quantum line or a quantum particle in four or three spatial dimensions respectively. These excitation have statistical interactions with the fractionalized Faraday lines; for example, in three spatial dimensions, the particle excitation picks up a Berry phase of $e^{i2π/N}$ when going around the fractionalized Faraday line excitation. We demonstrate the existence of this phase by Monte Carlo simulations in (3+1) space-time dimensions.

preprint2014arXiv

Plaquette Ordered Phase and Quantum Phase Diagram in the Spin-1/2 J1-J2 Square Heisenberg Model

We study the spin-1/2 Heisenberg model on the square lattice with first- and second-neighbor antiferromagnetic interactions J1 and J2, which possesses a nonmagnetic region that has been debated for many years and might realize the interesting Z2 spin liquid. We use the density matrix renormalization group approach with explicit implementation of SU(2) spin rotation symmetry and study the model accurately on open cylinders with different boundary conditions. With increasing J2, we find a Neel phase, a plaquette valence-bond (PVB) phase with a finite spin gap, and a possible spin liquid in a small region of J2 between these two phases. From the finite-size scaling of the magnetic order parameter, we estimate that the Neel order vanishes at J2/J1~0.44. For 0.5<J2/J1<0.61, we find dimer correlations and PVB textures whose decay lengths grow strongly with increasing system width, consistent with a long-range PVB order in the two-dimensional limit. The dimer-dimer correlations reveal the s-wave character of the PVB order. For 0.44<J2/J1<0.5, spin order, dimer order, and spin gap are small on finite-size systems and appear to scale to zero with increasing system width, which is consistent with a possible gapless SL or a near-critical behavior. We compare and contrast our results with earlier numerical studies.

preprint2013arXiv

Exact realization of Integer and Fractional Quantum Hall Phases in U(1)xU(1) models in (2+1)d

In this work we present a set of microscopic U(1)xU(1) models which realize insulating phases with a quantized Hall conductivity σ_{xy}. The models are defined in terms of physical degrees of freedom, and can be realized by local Hamiltonians. For one set of these models, we find that σ_{xy} is quantized to be an even integer. The origin of this effect is a condensation of objects made up of bosons of one species bound to a single vortex of the other species. For other models, the Hall conductivity can be quantized as a rational number times two. For these systems, the condensed objects contain bosons of one species bound to multiple vortices of the other species. These systems have excitations carrying fractional charges and non-trivial mutual statistics. We present sign-free reformulations of these models which can be studied in Monte Carlo, and we use such reformulations to numerically detect a gapless boundary between the quantum Hall and trivial insulator states. We also present the broader phase diagrams of the models.

preprint2013arXiv

Non-Fermi-liquid d-wave metal phase of strongly interacting electrons

Developing a theoretical framework for conducting electronic fluids qualitatively distinct from those described by Landau's Fermi-liquid theory is of central importance to many outstanding problems in condensed matter physics. One such problem is that, above the transition temperature and near optimal doping, high-transition-temperature copper-oxide superconductors exhibit `strange metal' behaviour that is inconsistent with being a traditional Landau Fermi liquid. Indeed, a microscopic theory of a strange-metal quantum phase could shed new light on the interesting low-temperature behaviour in the pseudogap regime and on the d-wave superconductor itself. Here we present a theory for a specific example of a strange metal---the 'd-wave metal'. Using variational wavefunctions, gauge theoretic arguments, and ultimately large-scale density matrix renormalization group calculations, we show that this remarkable quantum phase is the ground state of a reasonable microscopic Hamiltonian---the usual t-J model with electron kinetic energy $t$ and two-spin exchange $J$ supplemented with a frustrated electron `ring-exchange' term, which we here examine extensively on the square lattice two-leg ladder. These findings constitute an explicit theoretical example of a genuine non-Fermi-liquid metal existing as the ground state of a realistic model.

preprint2013arXiv

Phase diagram of spin-1/2 J1-J2 Heisenberg model on honeycomb lattice

We use density matrix renormalization group (DMRG) algorithm to study the phase diagram of the spin-1/2 Heisenberg model on honeycomb lattice with first (J1) and second (J2) neighbor antiferromagnetic interactions, where a Z2 spin liquid region has been proposed. By implementing SU(2) symmetry in the DMRG code, we are able to obtain accurate results for long cylinders with width slightly over 15 lattice spacings and torus up to the size N=2*6*6. With increasing J2, we find a Neel phase with vanishing spin gap and a plaquette valence-bond (PVB) phase with non-zero spin gap. By extrapolating the square of the staggered magnetic moment ms^2 on finite-size cylinders to thermodynamic limit, we find the Neel order vanishing at J2/J1~0.22. For 0.25<J2/J1<0.35, we find a possible PVB order, which shows a fast growing PVB decay length with increasing system width. For 0.22<J2/J1<0.25, both spin and dimer orders are weak for all systems we have studied, which is consistent with a possible spin liquid phase or a deconfined quantum critical point. We present calculations of the topological entanglement entropy, compare the DMRG results with the variational Monte Carlo, and discuss possible scenarios in the thermodynamic limit for this region.

preprint2012arXiv

Bose Metals and Insulators on Multi-Leg Ladders with Ring Exchange

We establish compelling evidence for the existence of new quasi-one-dimensional descendants of the d-wave Bose liquid (DBL), an exotic two-dimensional quantum phase of uncondensed itinerant bosons characterized by surfaces of gapless excitations in momentum space [O. I. Motrunich and M. P. A. Fisher, Phys. Rev. B {\bf 75}, 235116 (2007)]. In particular, motivated by a strong-coupling analysis of the gauge theory for the DBL, we study a model of hard-core bosons moving on the $N$-leg square ladder with frustrating four-site ring exchange. Here, we focus on four- and three-leg systems where we have identified two novel phases: a compressible gapless Bose metal on the four-leg ladder and an incompressible gapless Mott insulator on the three-leg ladder. The former is conducting along the ladder and has five gapless modes, one more than the number of legs. This represents a significant step forward in establishing the potential stability of the DBL in two dimensions. The latter, on the other hand, is a fundamentally quasi-one-dimensional phase that is insulating along the ladder but has two gapless modes and incommensurate power law transverse density-density correlations. In both cases, we can understand the nature of the phase using slave-particle-inspired variational wave functions consisting of a product of two distinct Slater determinants, the properties of which compare impressively well to a density matrix renormalization group solution of the model Hamiltonian. Stability arguments are made in favor of both quantum phases by accessing the universal low-energy physics with a bosonization analysis of the appropriate quasi-1D gauge theory. We will briefly discuss the potential relevance of these findings to high-temperature superconductors, cold atomic gases, and frustrated quantum magnets.

preprint2012arXiv

Line of continuous phase transitions in a three dimensional U(1) model with 1/r^2 current-current interactions

We study a lattice model of interacting loops in three dimensions with a $1/r^2$ interaction. Using Monte Carlo, we find that the phase diagram contains a line of second-order phase transitions between a phase where the loops are gapped and a phase where they proliferate. The correlation length exponent and critical conductivity vary continuously along this line. Our model is exactly self-dual at a special point on the critical line, which allows us to calculate the critical conductivity exactly at this point.

preprint2012arXiv

Majorana Spin Liquids on a two-leg ladder

We realize a gapless Majorana Orbital Liquid (MOL) using orbital degrees of freedom and also an SU(2)-invariant Majorana Spin Liquid (MSL) using both spin and orbital degrees of freedom in Kitaev-type models on a 2-leg ladder. The models are exactly solvable by Kitaev's parton approach, and we obtain long-wavelength descriptions for both Majorana liquids. The MOL has one gapless mode and power law correlations in energy at incommensuare wavevectors, while the SU(2) MSL has three gapless modes and power law correlations in spin, spin-nematic, and local energy observables. We study the stability of such states to perturbations away from the exactly solvable points. We find that both MOL and MSL can be stable against allowed short-range parton interactions. We also argue that both states persist upon allowing $Z_2$ gauge field fluctuations, in that the number of gapless modes is retained, although with an expanded set of contributions to observables compared to the free parton mean field.

preprint2012arXiv

Monte Carlo Study of a $U(1)\times U(1)$ Loop Model with Modular Invariance

We study a $U(1)\times U(1)$ system in (2+1)-dimensions with long-range interactions and mutual statistics. The model has the same form after the application of operations from the modular group, a property which we call modular invariance. Using the modular invariance of the model, we propose a possible phase diagram. We obtain a sign-free reformulation of the model and study it in Monte Carlo. This study confirms our proposed phase diagram. We use the modular invariance to analytically determine the current-current correlation functions and conductivities in all the phases in the diagram, as well as at special "fixed" points which are unchanged by an operation from the modular group. We numerically determine the order of the phase transitions, and find segments of second-order transitions. For the statistical interaction parameter $θ=π$, these second-order transitions are evidence of a critical loop phase obtained when both loops are trying to condense simulataneously. We also measure the critical exponents of the second-order transitions.

preprint2012arXiv

Phases and phase transitions in a $U(1)\times U(1)$ system with $θ=2π/3$ mutual statistics

We study a $U(1)\times U(1)$ system with short-range interactions and mutual $θ=2π/3$ statistics in (2+1) dimensions. We are able to reformulate the model to eliminate the sign problem, and perform a Monte Carlo study. We find a phase diagram containing a phase with only small loops and two phases with one species of proliferated loop. We also find a phase where both species of loop condense, but without any gapless modes. Lastly, when the energy cost of loops becomes small we find a phase which is a condensate of bound states, each made up of three particles of one species and a vortex of the other. We define several exact reformulations of the model, which allow us to precisely describe each phase in terms of gapped excitations. We propose field-theoretic descriptions of the phases and phase transitions, which are particularly interesting on the "self-dual" line where both species have identical interactions. We also define irreducible responses useful for describing the phases.

preprint2011arXiv

Exotic Gapless Mott Insulators of Bosons on Multi-Leg Ladders

We present evidence for an exotic gapless insulating phase of hard-core bosons on multi-leg ladders with a density commensurate with the number of legs. In particular, we study in detail a model of bosons moving with direct hopping and frustrating ring exchange on a 3-leg ladder at $ν=1/3$ filling. For sufficiently large ring exchange, the system is insulating along the ladder but has two gapless modes and power law transverse density correlations at incommensurate wave vectors. We propose a determinantal wave function for this phase and find excellent comparison between variational Monte Carlo and density matrix renormalization group calculations on the model Hamiltonian, thus providing strong evidence for the existence of this exotic phase. Finally, we discuss extensions of our results to other $N$-leg systems and to $N$-layer two-dimensional structures.

preprint2011arXiv

Monte Carlo Study of a U(1)xU(1) system with π-statistical Interaction

We study a $U(1)\times U(1)$ system with two species of loops with mutual $π$-statistics in (2+1) dimensions. We are able to reformulate the model in a way that can be studied by Monte Carlo and we determine the phase diagram. In addition to a phase with no loops, we find two phases with only one species of loop proliferated. The model has a self-dual line, a segment of which separates these two phases. Everywhere on the segment, we find the transition to be first-order, signifying that the two loop systems behave as immiscible fluids when they are both trying to condense. Moving further along the self-dual line, we find a phase where both loops proliferate, but they are only of even strength, and therefore avoid the statistical interactions. We study another model which does not have this phase, and also find first-order behavior on the self-dual segment.

preprint2011arXiv

Power-Law Behavior of Bond Energy Correlators in a Kitaev-type Model with a Stable Parton Fermi Surface

We study bond energy correlation functions in an exactly solvable quantum spin model of Kitaev type on the kagome lattice with stable Fermi surface of partons proposed recently by Chua et al, Ref.\[arXiv:1010.1035]. Even though any spin correlations are ultra-short ranged, we find that the bond energy correlations have power law behavior with a $1/|{\bm r}|^3$ envelope and oscillations at incommensurate wavevectors. We determine the corresponding singular surfaces in momentum space, which provide a gauge-invariant characterization of this gapless spin liquid.

preprint2011arXiv

Sign structures for short-range RVB states on small kagome clusters

We extended the Schwinger boson construction to obtain wave functions that are resonating valence bond (RVB) counterparts of the degenerate coplanar classical states on the Kagome lattice. We examined all 84 of them on the 36-site cluster and found that they form a narrow energy band. On the 12-site cluster, there are only four such states and their superpositions accurately account for the second through fourth exact diagonalization (ED) states, while the ED ground state is accurately reproduced by allowing a particular two-vison insertion on top of the q=0 RVB state. Thus, we have established the RVB sign structures for the low energy states on this cluster.

preprint2011arXiv

Spin Bose-Metal and Valence Bond Solid phases in a spin-1/2 model with ring exchanges on a four-leg triangular ladder

We study a spin-1/2 system with Heisenberg plus ring exchanges on a four-leg triangular ladder using the density matrix renormalization group and Gutzwiller variational wave functions. Near an isotropic lattice regime, for moderate to large ring exchanges we find a spin Bose-metal phase with a spinon Fermi sea consisting of three partially filled bands. Going away from the triangular towards the square lattice regime, we find a staggered dimer phase with dimers in the transverse direction, while for small ring exchanges the system is in a featureless rung phase. We also discuss parent states and a possible phase diagram in two dimensions.

preprint2011arXiv

SU(2)-invariant Majorana spin liquid with stable parton Fermi surfaces in an exactly solvable model

We construct an exactly solvable spin-orbital model on a decorated square lattice that realizes an SU(2)-invariant Majorana spin liquid with parton Fermi surfaces, of the kind discussed recently by Biswas et al. [ Phys. Rev. B. {\bf 83}, 245131(2011)]. We find power-law spin correlations as well as power-law spin-nematic correlations with the same dominant $1/|{\bf r}|^3$ envelope. The model is solvable also in the presence of Zeeman magnetic field. One fermion species carries $S^z=0$ quantum number and its Fermi surface is not altered in the field, while the Fermi surfaces of the other species evolve and can disappear. In particular, we find an interesting half-magnetization plateau phase in which spin excitations are gapful while there remain spinless gapless excitations that still produce metal-like thermal properties. In the fully-magnetized phase, the model reduces to the one proposed by Baskaran et al., [arXiv:0908.1614v3] in terms of the orbital degrees of freedom.

preprint2011arXiv

Variational study of J1-J2 Heisenberg model on Kagome lattice using projected Schwinger boson wave functions

Motivated by the unabating interest in the spin-1/2 Heisenberg antiferromagnetic model on the Kagome lattice, we investigate the energetics of projected Schwinger boson (SB) wave functions in the $J_1$--$J_2$ model with antiferromagnetic $J_2$ coupling. Our variational Monte Carlo results show that Sachdev's $Q_1=Q_2$ SB ansatz has a lower energy than the Dirac spin liquid for $J_2\gtrsim 0.08 J_1$ and the ${\bf q=0}$ Jastrow type magnetically ordered state. This work demonstrates that the projected SB wave functions can be tested on the same footing as their fermionic counterparts.

preprint2010arXiv

Effects of impurities in Spin Bose-Metal phase on a two-leg triangular strip

We study effects of nonmagnetic impurities in a Spin Bose-Metal (SBM) phase discovered in a two-leg triangular strip spin-1/2 model with ring exchanges (D. N. Shenget al, arXiv:0902.4210). This phase is a quasi-1D descendant of a 2D spin liquid with spinon Fermi sea, and the present study aims at interpolating between the 1D and 2D cases. Different types of defects can be treated as local energy perturbations, which we find are always relevant. As a result, a nonmagnetic impurity generically cuts the system into two decoupled parts. We calculate bond energy and local spin susceptibility near the defect, both of which can be measured in experiments. The Spin Bose-Metal has dominant correlations at characteristic incommensurate wavevectors that are revealed near the defect. Thus, the bond energy shows a static texture oscillating as a function of distance from the defect and decaying as a slow power law. The local spin susceptibility also oscillates and actually {\it increases} as a function of distance from the defect, similar to the effect found in the 1D chain [S. Eggert and I. Affleck, Phys. Rev. Lett. {\bf 75}, 934 (1995)]. We calculate the corresponding power law exponents for the textures as a function of one Luttinger parameter of the SBM theory.

preprint2010arXiv

Effects of Zeeman field on a Spin Bose-Metal phase

We consider Zeeman field effects on a Spin Bose-Metal (SBM) phase on a two-leg triangular ladder. This phase was found in a spin-1/2 model with ring exchanges [D. N. Sheng et. al., Phys. Rev. B {\bf 79}, 205112 (2009)], and was also proposed to appear in an interacting electronic model with longer-ranged repulsion [Lai et. al., Phys. Rev. B {\bf 81}, 045105 (2010)]. Using bosonization of a spinon-gauge theory, we study the stability of the SBM phase and its properties under the field. We also explore phases arising from potential instabilities of the SBM; in all cases, we find a gap to spin-1 excitations while spin-nematic correlations are power law. We discuss two-dimensional analogues of these phases where spinons can pair with their own species.

preprint2010arXiv

Failure of Gutzwiller-type wave function to capture gauge fluctuations: Case study in the Exciton Bose Liquid context

Slave particle approaches are widely used in studies of exotic quantum phases. A complete description beyond mean field also contains dynamical gauge fields, while a simplified procedure considers Gutzwiller-projected trial states. We apply this in the context of bosonic models with ring exchanges realizing so-called Exciton Bose Liquid (EBL) phase and compare a Gutzwiller wave function against an accurate EBL wave function. We solve the parton-gauge theory and show that dynamical fluctuations of the spatial gauge fields are necessary for obtaining qualitatively accurate EBL description. On the contrary, just the Gutzwiller projection leads to a state with subtle differences in the long-wavelength properties, thus suggesting that Gutzwiller wave functions may generally fail to capture long-wavelength physics.

preprint2010arXiv

Insulating phases of electrons on a zigzag strip in the orbital magnetic field

We consider electrons on a two-leg triangular ladder at half-filling and in an orbital magnetic field. In a two-band regime in the absence of the field, the electronic system remains conducting for weak interactions since there is no four-fermion Umklapp term. We find that in the presence of the orbital field there is a four-fermion Umklapp and it is always relevant for repulsive interactions. Thus in this special ladder, the combination of the orbital magnetic field and interactions provides a mechanism to drive metal-insulator transition already at weak coupling. We discuss properties of the possible resulting phases C0S2 and various C0S1 and C0S0.

preprint2010arXiv

Possible realization of the Exciton Bose Liquid phase in a hard-core boson model with ring-only exchange interactions

We investigate a hard-core boson model with ring-only exchanges on a square lattice, where a $K_1$ term acts on 1$\times$1 plaquettes and a $K_2$ term acts on 1$\times$2 and 2$\times$1 plaquettes, with a goal of realizing a novel Exciton Bose Liquid (EBL) phase first proposed by Paramekanti, et al [Phys. Rev. B {\bf 66}, 054526 (2002)]. We construct Jastrow-type variational wave functions for the EBL, study their formal properties, and then use them as seeds for a projective Quantum Monte Carlo study. Using Green's Function Monte Carlo, we obtain an unbiased phase diagram which at half-filling reveals CDW for small $K_2$, valence bond solid for intermediate $K_2$, and possibly for large $K_2$ the EBL phase. Away from half-filling, the EBL phase is present for intermediate $K_2$ and remains stable for a range of densities below 1/2 before phase separation occurs at lower densities.

preprint2010arXiv

Study of a hard-core boson model with ring-only interactions

We present a Quantum Monte Carlo study of a hardcore boson model with ring-only exchanges on a square lattice, where a $K_1$ term acts on 1$\times$1 plaquettes and a $K_2$ term acts on 1$\times$2 and 2$\times$1 plaquettes. At half-filling, the phase diagram reveals charge density wave for small $K_2$, valence bond solid for intermediate $K_2$, and possibly for large $K_2$ the novel Exciton Bose Liquid (EBL) phase first proposed by Paramekanti, et al[Phys. Rev. B {\bf 66}, 054526 (2002)]. Away from half-filling, the EBL phase is present already for intermediate $K_2$ and remains stable for a range of densities below 1/2 before phase separation sets in at lower densities

preprint2010arXiv

Variational Studies of Triangular Heisenberg Antiferromagnet in Magnetic Field

We present a variational study of the Heisenberg antiferromagnet on the spatially anisotropic triangular lattice in magnetic field. First we construct a simple yet accurate wavefunction for the 1/3-magnetization plateau uud phase on the isotropic lattice. Beginning with this state, we obtain natural extensions to nearby commensurate coplanar phases on either side of the plateau. The latter occur also for low lattice anisotropy, while the uud state extends to much larger anisotropy. Far away from the 1/3 plateau and for significant anisotropy, incommensurate states have better energetics, and we address competition between coplanar and non-coplanar states in the high field regime. For very strong anisotropy, our study is dominated by quasi-1d physics. The variational study is supplemented by exact diagonalization calculations which provide a reference for testing the energetics of our trial wavefunctions as well as helping to identify candidate phases.

preprint2009arXiv

Interlayer coherent composite Fermi liquid phase in quantum Hall bilayers

Composite fermions have played a seminal role in understanding the quantum Hall effect, particularly the formation of a compressible `composite Fermi liquid' (CFL) at filling factor nu = 1/2. Here we suggest that in multi-layer systems interlayer Coulomb repulsion can similarly generate `metallic' behavior of composite fermions between layers, even if the electrons remain insulating. Specifically, we propose that a quantum Hall bilayer with nu = 1/2 per layer at intermediate layer separation may host such an interlayer coherent CFL, driven by exciton condensation of composite fermions. This phase has a number of remarkable properties: the presence of `bonding' and `antibonding' composite Fermi seas, compressible behavior with respect to symmetric currents, and fractional quantum Hall behavior in the counterflow channel. Quantum oscillations associated with the Fermi seas give rise to a new series of incompressible states at fillings nu = p/[2(p \pm 1)] per layer (p an integer), which is a bilayer analogue of the Jain sequence.

preprint2009arXiv

Spin Bose-Metal phase in a spin-1/2 model with ring exchange on a two-leg triangular strip

Recent experiments on triangular lattice organic Mott insulators have found evidence for a 2D spin liquid in proximity to the metal-insulator transition. A Gutzwiller wavefunction study of the triangular lattice Heisenberg model with appropriate four-spin ring exchanges has found that the projected spinon Fermi sea state has a low variational energy. This wavefunction, together with a slave particle gauge theory, suggests that such spin liquid possesses spin correlations that are singular along surfaces in momentum space ("Bose surfaces"). Signatures of this state, which we refer to as a "Spin Bose-Metal" (SBM), are expected to be manifest in quasi-1D ladder systems: The discrete transverse momenta cut through the 2D Bose surface leading to a distinct pattern of 1D gapless modes. Here we search for a quasi-1D descendant of the triangular lattice SBM state by exploring the Heisenberg plus ring model on a two-leg strip (zigzag chain). Using DMRG, variational wavefunctions, and a Bosonization analysis, we map out the full phase diagram. Without ring exchange the model is equivalent to the J_1 - J_2 Heisenberg chain, and we find the expected Bethe-chain and dimerized phases. Remarkably, moderate ring exchange reveals a new gapless phase over a large swath of the phase diagram. Spin and dimer correlations possess particular singular wavevectors and allow us to identify this phase as the hoped for quasi-1D descendant SBM state. We derive a low energy theory and find three gapless modes and one Luttinger parameter controlling all power laws. Potential instabilities out of the zigzag SBM give rise to other interesting phases such as a period-3 VBS or a period-4 Chirality order, which we discover in the DMRG; we also find an interesting SBM state with partial ferromagnetism.

preprint2009arXiv

Two-band electronic metal and neighboring spin liquid (spin Bose-metal) on a zigzag strip with longer-ranged repulsion

We consider an electronic model for realizing the Spin Bose-metal (SBM) phase on a 2-leg triangular strip --a spin liquid phase found by D. N. Sheng et al [Phys. Rev. B {\bf 79}, 205112 (2009)] in a spin-1/2 model with ring exchanges. The SBM can be viewed as a "C1S2" Mott insulator of electrons where the overall charge transporting mode is gapped out. We start from a two-band "C2S2" metal and consider extended repulsion motivated by recent ab initio derivation of electronic model for $κ$-ET spin liquid material [K. Nakamura et al, J. Phys. Soc. Jpn. {\bf 78}, 083710(2009)]. Using weak coupling renormalization group analysis, we find that the extended interactions allow much wider C2S2 metallic phase than in the Hubbard model with on-site repulsion only. An eight-fermion Umklapp term plays a crucial role in producing a Mott insulator but can not be treated in weak coupling. We use Bosonization to extend the analysis to intermediate coupling and study phases obtained out of the C2S2 metal upon increasing overall repulsion strength, finding that the SBM phase is a natural outcome for extended interactions.

preprint2008arXiv

Impurities in triangular lattice spin 1/2 antiferromagnet

We study effects of nonmagnetic impurities in a spin-1/2 frustrated triangular antiferromagnet with the aim of understanding the observed broadening of $^{13}$C NMR lines in the organic spin liquid material $κ$-(ET)$_2$Cu$_2$(CN)$_3$. For high temperatures down to $J/3$, we calculate local susceptibility near a nonmagnetic impurity and near a grain boundary for the nearest neighbor Heisenberg model in high temperature series expansion. We find that the local susceptibility decays to the uniform one in few lattice spacings, and for a low density of impurities we would not be able to explain the line broadening present in the experiments already at elevated temperatures. At low temperatures, we assume a gapless spin liquid with a Fermi surface of spinons. We calculate the local susceptibility in the mean field and also go beyond the mean field by Gutzwiller projection. The zero temperature local susceptibility decays as a power law and oscillates at $2 k_F$. As in the high temperature analysis we find that a low density of impurities is not able to explain the observed broadening of the lines. We are thus led to conclude that there is more disorder in the system. We find that a large density of point-like disorder gives broadening that is consistent with the experiment down to about 5K, but that below this temperature additional mechanism is likely needed.