Phases of the hard-plate lattice gas on a three-dimensional cubic lattice.

We study the phase diagram of a lattice gas of 2×2×1 hard plates on the three-dimensional cubic lattice. Each plate covers an elementary plaquette of the cubic lattice, with the constraint that a site can belong to utmost one plate. We focus on the isotropic system, with equal fugacities for the three orientations of plates. We show, using grand canonical Monte Carlo simulations, that the system undergoes two phase transitions when the density of plates is increased: the first from a disordered fluid phase to a layered phase, and the second from the layered phase to a sublattice-ordered phase. In the layered phase, the system breaks up into disjoint slabs of thickness two along one spontaneously chosen Cartesian direction, corresponding to a twofold (Z_{2}) symmetry breaking of translation symmetry along the layering direction. Plates with normals perpendicular to this layering direction are preferentially contained entirely within these slabs, while plates straddling two adjacent slabs have a lower density, thus breaking the symmetry between the three types of plates. We show that the slabs exhibit two-dimensional power-law columnar order even in the presence of a nonzero density of vacancies. In contrast, interslab correlations of the two-dimensional columnar order parameter decay exponentially with the separation between the slabs. In the sublattice-ordered phase, there is twofold symmetry breaking of lattice translation symmetry along all three Cartesian directions. We present numerical evidence that the disordered to layered transition is continuous and consistent with universality class of the three-dimensional O(3) model with cubic anisotropy, while the layered to sublattice transition is first-order in nature.

Spontaneous layering and power-law order in the three-dimensional fully packed hard-plate lattice gas.

We obtain the phase diagram of fully packed hard plates on the cubic lattice. Each plate covers an elementary plaquette of the cubic lattice and occupies its four vertices, with each vertex of the cubic lattice occupied by exactly one such plate. We consider the general case with fugacities s_{µ} for "µ plates," whose normal is the µ direction (µ=x,y,z). At and close to the isotropic point, we find, consistent with previous work, a phase with long-range sublattice order. When two of the fugacities s_{µ_{1}} and s_{µ_{2}} are comparable, and the third fugacity s_{µ_{3}} is much smaller, we find a spontaneously layered phase. In this phase, the system breaks up into disjoint slabs of width two stacked along the µ_{3} axis. µ_{1} and µ_{2} plates are preferentially contained entirely within these slabs, while plates straddling two successive slabs have a lower density. This corresponds to a twofold breaking of translation symmetry along the µ_{3} axis. In the opposite limit, with µ_{3}≫µ_{1}∼µ_{2}, we find a phase with long-range columnar order, corresponding to simultaneous twofold symmetry breaking of lattice translation symmetry in directions µ_{1} and µ_{2}. The spontaneously layered phases display critical behavior, with power-law decay of correlations in the µ_{1} and µ_{2} directions when the slabs are stacked in the µ_{3} direction, and represent examples of "floating phases" discussed earlier in the context of coupled Luttinger liquids and quasi-two-dimensional classical systems. We ascribe this remarkable behavior to the constrained motion of defects in this phase, and we sketch a coarse-grained effective field theoretical understanding of the stability of power-law order in this unusual three-dimensional floating phase.

Emergent Moments and Random Singlet Physics in a Majorana Spin Liquid.

We exhibit an exactly solvable example of a SU(2) symmetric Majorana spin liquid phase, in which quenched disorder leads to random-singlet phenomenology of emergent magnetic moments. More precisely, we argue that a strong-disorder fixed point controls the low temperature susceptibility χ(T) of an exactly solvable S=1/2 model on the decorated honeycomb lattice with vacancy and/or bond disorder, leading to χ(T)=C/T+DT^{α(T)-1}, where α(T)→0 slowly as the temperature T→0. The first term is a Curie tail that represents the emergent response of vacancy-induced spin textures spread over many unit cells: it is an intrinsic feature of the site-diluted system, rather than an extraneous effect arising from isolated free spins. The second term, common to both vacancy and bond disorder [with different α(T) in the two cases] is the response of a random singlet phase, familiar from random antiferromagnetic spin chains and the analogous regime in phosphorus-doped silicon (Si:P).

Fractional Brownian motion of worms in worm algorithms for frustrated Ising magnets.

We study the distribution of lengths and other statistical properties of worms constructed by Monte Carlo worm algorithms in the power-law three-sublattice ordered phase of frustrated triangular and kagome lattice Ising antiferromagnets. Viewing each step of the worm construction as a position increment (step) of a random walker, we demonstrate that the persistence exponent θ and the dynamical exponent z of this random walk depend only on the universal power-law exponents of the underlying critical phase and not on the details of the worm algorithm or the microscopic Hamiltonian. Further, we argue that the detailed balance condition obeyed by such worm algorithms and the power-law correlations of the underlying equilibrium system together give rise to two related properties of this random walk: First, the steps of the walk are expected to be power-law correlated in time. Second, the position distribution of the walker relative to its starting point is given by the equilibrium position distribution of a particle in an attractive logarithmic central potential of strength η_{m}, where η_{m} is the universal power-law exponent of the equilibrium defect-antidefect correlation function of the underlying spin system. We derive a scaling relation, z=(2-η_{m})/(1-θ), that allows us to express the dynamical exponent z(η_{m}) of this process in terms of its persistence exponent θ(η_{m}). Our measurements of z(η_{m}) and θ(η_{m}) are consistent with this relation over a range of values of the universal equilibrium exponent η_{m} and yield subdiffusive (z>2) values of z in the entire range. Thus, we demonstrate that the worms represent a discrete-time realization of a fractional Brownian motion characterized by these properties.

Bilayer Coulomb phase of two-dimensional dimer models: Absence of power-law columnar order.

Using renormalization group (RG) analyses and Monte Carlo (MC) simulations, we study the fully packed dimer model on the bilayer square lattice with fugacity equal to z (1) for interlayer (intralayer) dimers, and intralayer interaction V between neighboring parallel dimers on any elementary plaquette in either layer. For a range of not-too-large z>0 and repulsive interactions 00 destroys the power-law correlations of the z=0 decoupled layers, and leads immediately to a short-range correlated state, albeit with a slow crossover for small |V|. For V_{c}

Phase diagram of a system of hard cubes on the cubic lattice.

We study the phase diagram of a system of 2×2×2 hard cubes on a three-dimensional cubic lattice. Using Monte Carlo simulations, we show that the system exhibits four different phases as the density of cubes is increased: disordered, layered, sublattice ordered, and columnar ordered. In the layered phase, the system spontaneously breaks up into parallel slabs of size 2×L×L where only a very small fraction cubes do not lie wholly within a slab. Within each slab, the cubes are disordered; translation symmetry is thus broken along exactly one principal axis. In the solidlike sublattice-ordered phase, the hard cubes preferentially occupy one of eight sublattices of the cubic lattice, breaking translational symmetry along all three principal directions. In the columnar phase, the system spontaneously breaks up into weakly interacting parallel columns of size 2×2×L, where only a very small fraction cubes do not lie wholly within a column. Within each column, the system is disordered, and thus translational symmetry is broken only along two principal directions. Using finite-size scaling, we show that the disordered-layered phase transition is continuous, while the layered-sublattice and sublattice-columnar transitions are discontinuous. We construct a Landau theory written in terms of the layering and columnar order parameters which is able to describe the different phases that are observed in the simulations and the order of the transitions. Additionally, our results near the disordered-layered transition are consistent with the O(3) universality class perturbed by cubic anisotropy as predicted by the Landau theory.

Cluster algorithms for frustrated two-dimensional Ising antiferromagnets via dual worm constructions.

We report on the development of two dual worm constructions that lead to cluster algorithms for efficient and ergodic Monte Carlo simulations of frustrated Ising models with arbitrary two-spin interactions that extend up to third-neighbors on the triangular lattice. One of these algorithms generalizes readily to other frustrated systems, such as Ising antiferromagnets on the Kagome lattice with further neighbor couplings. We characterize the performance of both these algorithms in a challenging regime with power-law correlations at finite wave vector.

Sign-Problem-Free Monte Carlo Simulation of Certain Frustrated Quantum Magnets.

We introduce a quantum Monte Carlo (QMC) method for efficient sign-problem-free simulations of a broad class of frustrated S=1/2 antiferromagnets using the basis of spin eigenstates of clusters to avoid the severe sign problem faced by other QMC methods. We demonstrate the utility of the method in several cases with competing exchange interactions and flag important limitations as well as possible extensions of the method.

Classical Spin Liquid on the Maximally Frustrated Honeycomb Lattice.

We show that the honeycomb Heisenberg antiferromagnet with J_{1}/2=J_{2}=J_{3}, where J_{1}, J_{2}, and J_{3} are first-, second-, and third-neighbor couplings, respectively, forms a classical spin liquid with pinch-point singularities in the structure factor at the Brillouin zone corners. Upon dilution with nonmagnetic ions, fractionalized degrees of freedom carrying 1/3 of the free moment emerge. Their effective description in the limit of low temperature is that of spins randomly located on a triangular lattice, with a frustrated sublattice-sensitive interaction of long-ranged logarithmic form. The XY version of this magnet exhibits nematic thermal order by disorder. This comes with a clear experimental diagnostic in neutron scattering, which turns out to apply also to the case of the celebrated planar order by disorder of the kagome Heisenberg antiferromagnet.

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-neighbor 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, but not-too-low, |ε|/t in this compensated case, we show that the density of states ρ(ε) exhibits a strong divergence of the form ρ_{Dyson}(ε)∼|ε|^{-1}/[log(t/|ε|)]^{(y+1)}, which crosses over to the universal low-energy asymptotic form (modified Gade-Wegner scaling) expected on symmetry grounds ρ_{GW}(ε)∼|ε|^{-1}e^{-b[log(t/|ε|)]^{2/3}} below a crossover scale ε_{c}≪t. ε_{c} is found to decrease rapidly with decreasing n_{v}, while y decreases much more slowly.

Melting of Three-Sublattice Order in Easy-Axis Antiferromagnets on Triangular and Kagome Lattices.

When the constituent spins have an energetic preference to lie along an easy axis, triangular and kagome lattice antiferromagnets often develop long-range order that distinguishes the three sublattices of the underlying triangular Bravais lattice. In zero magnetic field, this three-sublattice order melts either in a two-step manner, i.e., via an intermediate phase with power-law three-sublattice order controlled by a temperature-dependent exponent η(T)∈(1/9,1/4), or via a transition in the three-state Potts universality class. Here, I predict that the uniform susceptibility to a small easy-axis field B diverges as χ(B)∼|B|^{-[(4-18η)/(4-9η)]} in a large part of the intermediate power-law ordered phase [corresponding to η(T)∈(1/9,2/9)], providing an easy-to-measure thermodynamic signature of two-step melting. I also show that these two melting scenarios can be generically connected via an intervening multicritical point and obtain numerical estimates of multicritical exponents.

Columnar order and Ashkin-Teller criticality in mixtures of hard squares and dimers.

We show that critical exponents of the transition to columnar order in a mixture of 2×1 dimers and 2×2 hard squares on the square lattice depends on the composition of the mixture in exactly the manner predicted by the theory of Ashkin-Teller criticality, including in the hard-square limit. This result settles the question regarding the nature of the transition in the hard-square lattice gas. It also provides the first example of a polydisperse system whose critical properties depend on composition. Our ideas also lead to some interesting predictions for a class of frustrated quantum magnets that exhibit columnar ordering of the bond energies at low temperature.

Néel-state to valence-bond-solid transition on the honeycomb lattice: evidence for deconfined criticality.

We study a spin-1/2 SU(2) model on the honeycomb lattice with nearest-neighbor antiferromagnetic exchange J that favors Néel order and competing six-spin interactions Q that favor a valence-bond-solid (VBS) state in which the bond energies order at the "columnar" wave vector K=(2π/3,-2π/3). We present quantum Monte Carlo evidence for a direct continuous quantum phase transition between Néel and VBS states, with exponents and logarithmic violations of scaling consistent with those at analogous deconfined critical points on the square lattice. Although this strongly suggests a description in terms of deconfined criticality, the measured threefold anisotropy of the phase of the VBS order parameter shows unusual near-marginal behavior at the critical point.

Resonating valence bond wave functions and classical interacting dimer models.

We relate properties of nearest-neighbor resonating valence-bond (NNRVB) wave functions for SU(g) spin systems on two-dimensional bipartite lattices to those of fully packed interacting classical dimer models on the same lattice. The interaction energy can be expressed as a sum of n-body potentials V(n), which are recursively determined from the NNRVB wave function on finite subgraphs of the original lattice. The magnitude of the n-body interaction V(n) (n>1) is of order O(g(-(n-1))) for small g(-1). The leading term is a two-body nearest-neighbor interaction V2(g) favoring two parallel dimers on elementary plaquettes. For SU(2) spins, using our calculated value of V2(g=2), we find that the long-distance behavior of the bond-energy correlation function is dominated by an oscillatory term that decays as 1/|r|α with α≈1.22. This result is in remarkable quantitative agreement with earlier direct numerical studies of the corresponding wave function, which give α≈1.20.

Phase transitions in the distribution of the Andreev conductance of superconductor-metal junctions with multiple transverse modes.

We compute analytically the full distribution of Andreev conductance G(NS) of a metal-superconductor interface with a large number N(c) of transverse modes, using a random matrix approach. The probability distribution P(G(NS),N(c) in the limit of large N(c) displays a Gaussian behavior near the average value =(2-√2)N(c) and asymmetric power-law tails in the two limits of very small and very large G(NS). In addition, we find a novel third regime sandwiched between the central Gaussian peak and the power-law tail for large G(NS). Weakly nonanalytic points separate these four regimes-these are shown to be consequences of three phase transitions in an associated Coulomb gas problem.

Fractional spin textures in the frustrated magnet SrCr(9p)Ga(12-9p)O19.

We consider the archetypal frustrated antiferromagnet SrCr(9p)Ga(12-9p)O19 in its well-known spin-liquid state, and demonstrate that a Cr(3+) spin S=3/2 ion in direct proximity to a pair of vacancies (in disordered p<1 samples) is cloaked by a spatially extended spin texture that encodes the correlations of the parent spin liquid. In this spin-liquid regime, our analytic theory predicts that the combined object has a magnetic response identical to a classical spin of length S/2=3/4, which dominates over the small intrinsic susceptibility of the pure system. This fractional-spin texture leaves an unmistakable imprint on the measured 7¹Ga nuclear magnetic resonance line shapes, which we compute using Monte Carlo simulations and compare with experimental data.

Triangular and kagome antiferromagnets with a strong easy-axis anisotropy.

We consider S>3/2 kagome and triangular lattice magnets with strong easy-axis single-ion anisotropy D and antiferromagnetic exchange J. When D>>J, the low energy states selected by the anisotropy map onto configurations of the corresponding classical Ising antiferromagnet. Subleading O(J;{3}S/D;{2}) multispin interaction arising from the transverse quantum dynamics makes the low-temperature behavior very different from the well-known classical case: The kagome magnet goes into a semiclassical spin-liquid state with distinctive and unusual short-range correlations below a crossover temperature T;{*} approximately 0.08J;{3}S/D;{2}, while the triangular magnet undergoes a first-order transition at T_{c} approximately 0.1J;{3}S/D;{2} to an orientationally ordered collinear state that gives rise to a novel zero-magnetization plateau for small magnetic fields along the easy axis. Possible experimental implications are also discussed.

Variational wave-function study of the triangular lattice supersolid.

We present a variational wave function which explains the behavior of the supersolid state formed by hard-core bosons on the triangular lattice. The wave function is a linear superposition of only and all configurations minimizing the repulsion between the bosons (which it thus implements as a hard constraint). Its properties can be evaluated exactly--in particular, the variational minimization of the energy yields (i) the surprising and initially controversial spontaneous density deviation from half-filling (ii) a quantitatively accurate estimate of the corresponding density wave (solid) order parameter.

Unusual liquid state of hard-core bosons on the pyrochlore lattice.

We study the physics of hard-core bosons with unfrustrated hopping (t) and nearest-neighbor repulsion (V) on the three dimensional pyrochlore lattice. At half-filling, we demonstrate that the small V/t superfluid state eventually becomes unstable at large enough V/t to an unusual insulating state which displays no broken lattice translation symmetry. Equal time and static density correlators in this insulator are well described by a mapping to electric field correlators in the Coulomb phase of a U(1) lattice gauge theory, allowing us to identify this insulator with a U(1) fractionalized Mott-insulating state. The possibility of observing this phase in suitably designed atom-trap experiments with ultracold atoms is also discussed, as are specific experimental signatures.

Magnetization plateaus and sublattice ordering in easy-axis kagome lattice antiferromagnets.

We study kagome lattice antiferromagnets where the effects of easy-axis single-ion anisotropy (D) dominates over the Heisenberg exchange J. For S> or =3/2, virtual quantum fluctuations help lift the extensive classical degeneracy. We demonstrate the presence of a one-third magnetization plateau for a broad range of magnetic fields J3/D2 < or = B < or = JS along the easy axis. The fully equilibrated system at low temperature on this plateau develops an unusual nematic order that breaks sublattice rotation symmetry but not translation symmetry; however, extremely slow dynamics associated with this ordering is expected to lead to glassy freezing of the system on intermediate time scales.