Source author record

Roger G. Melko

Roger G. Melko 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

51works
15topics
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

51 published item(s)

preprint2025arXiv

Probing Defects with Quantum Simulator Snapshots

Snapshots, i.e. projective measurements of local degrees of freedom, are the most standard data taken in experiments on quantum simulators. Snapshots are usually used to probe local physics. In this work we propose a simple protocol to experimentally probe physics of defects with these snapshots. Our protocol relies only on snapshots from the bulk system, without introducing the defect explicitly; as such, the physics of different kinds of defects can be probed using the same dataset. In particular, we demonstrate that with snapshots of local spin configurations of, for example, the $1d$ Rydberg atom realization of the quantum Ising criticality, we can (1) extract the ``defect entropy", and (2) access the continuous line of fixed points of effective defect conformal field theory, which was recently discussed in the context of the ``weak-measurement altered criticality".

preprint2022arXiv

Data-Enhanced Variational Monte Carlo Simulations for Rydberg Atom Arrays

Rydberg atom arrays are programmable quantum simulators capable of preparing interacting qubit systems in a variety of quantum states. Due to long experimental preparation times, obtaining projective measurement data can be relatively slow for large arrays, which poses a challenge for state reconstruction methods such as tomography. Today, novel groundstate wavefunction ansätze like recurrent neural networks (RNNs) can be efficiently trained not only from projective measurement data, but also through Hamiltonian-guided variational Monte Carlo (VMC). In this paper, we demonstrate how pretraining modern RNNs on even small amounts of data significantly reduces the convergence time for a subsequent variational optimization of the wavefunction. This suggests that essentially any amount of measurements obtained from a state prepared in an experimental quantum simulator could provide significant value for neural-network-based VMC strategies.

preprint2022arXiv

Neural annealing and visualization of autoregressive neural networks in the Newman-Moore model

Artificial neural networks have been widely adopted as ansatzes to study classical and quantum systems. However, some notably hard systems such as those exhibiting glassiness and frustration have mainly achieved unsatisfactory results despite their representational power and entanglement content, thus, suggesting a potential conservation of computational complexity in the learning process. We explore this possibility by implementing the neural annealing method with autoregressive neural networks on a model that exhibits glassy and fractal dynamics: the two-dimensional Newman-Moore model on a triangular lattice. We find that the annealing dynamics is globally unstable because of highly chaotic loss landscapes. Furthermore, even when the correct ground state energy is found, the neural network generally cannot find degenerate ground-state configurations due to mode collapse. These findings indicate that the glassy dynamics exhibited by the Newman-Moore model caused by the presence of fracton excitations in the configurational space likely manifests itself through trainability issues and mode collapse in the optimization landscape.

preprint2022arXiv

Persistent Homology of $\mathbb{Z}_2$ Gauge Theories

Topologically ordered phases of matter display a number of unique characteristics, including ground states that can be interpreted as patterns of closed strings. In this paper, we consider the problem of detecting and distinguishing closed strings in Ising spin configurations sampled from the classical $\mathbb{Z}_2$ gauge theory. We address this using the framework of persistent homology, which computes the size and frequency of general loop structures in spin configurations via the formation of geometric complexes. Implemented numerically on finite-size lattices, we show that the first Betti number of the Vietoris-Rips complexes achieves a high density at low temperatures in the $\mathbb{Z}_2$ gauge theory. In addition, it displays a clear signal at the finite-temperature deconfinement transition of the three-dimensional theory. We argue that persistent homology should be capable of interpreting prominent loop structures that occur in a variety of systems, making it an useful tool in theoretical and experimental searches for topological order.

preprint2022arXiv

Toward Orbital-Free Density Functional Theory with Small Data Sets and Deep Learning

We use voxel deep neural networks to predict energy densities and functional derivatives of electron kinetic energies for the Thomas-Fermi model and Kohn-Sham density functional theory calculations. We show that the ground-state electron density can be found via direct minimization for a graphene lattice without any projection scheme using a voxel deep neural network trained with the Thomas-Fermi model. Additionally, we predict the kinetic energy of a graphene lattice within chemical accuracy after training from only 2 Kohn-Sham density functional theory (DFT) calculations. We identify an important sampling issue inherent in Kohn-Sham DFT calculations and propose future work to rectify this problem. Furthermore, we demonstrate an alternative, functional derivative-free, Monte Carlo based orbital free density functional theory algorithm to calculate an accurate 2-electron density in a double inverted Gaussian potential with a machine-learned kinetic energy functional.

preprint2021arXiv

Bulk and Boundary Quantum Phase Transitions in a Square Rydberg Atom Array

Motivated by recent experimental realizations of exotic phases of matter on programmable quantum simulators, we carry out a comprehensive theoretical study of quantum phase transitions in a Rydberg atom array on a square lattice, with both open and periodic boundary conditions. In the bulk, we identify several first-order and continuous phase transitions by performing large-scale quantum Monte Carlo simulations and develop an analytical understanding of the nature of these transitions using the framework of Landau-Ginzburg-Wilson theory. Remarkably, we find that under open boundary conditions, the boundary itself undergoes a second-order quantum phase transition, independent of the bulk. These results explain recent experimental observations and provide important new insights into both the adiabatic state preparation of novel quantum phases and quantum optimization using Rydberg atom array platforms.

preprint2021arXiv

Pruning a restricted Boltzmann machine for quantum state reconstruction

Restricted Boltzmann machines (RBMs) have proven to be a powerful tool for learning quantum wavefunction representations from qubit projective measurement data. Since the number of classical parameters needed to encode a quantum wavefunction scales rapidly with the number of qubits, the ability to learn efficient representations is of critical importance. In this paper we study magnitude-based pruning as a way to compress the wavefunction representation in an RBM, focusing on RBMs trained on data from the transverse-field Ising model in one dimension. We find that pruning can reduce the total number of RBM weights, but the threshold at which the reconstruction accuracy starts to degrade varies significantly depending on the phase of the model. In a gapped region of the phase diagram, the RBM admits pruning over half of the weights while still accurately reproducing relevant physical observables. At the quantum critical point however, even a small amount of pruning can lead to significant loss of accuracy in the physical properties of the reconstructed quantum state. Our results highlight the importance of tracking all relevant observables as their sensitivity varies strongly with pruning. Finally, we find that sparse RBMs are trainable and discuss how a successful sparsity pattern can be created without pruning.

preprint2020arXiv

Reconstructing quantum molecular rotor ground states

Nanomolecular assemblies of C$_{60}$ can be synthesized to enclose dipolar molecules. The low-temperature states of such endofullerenes are described by quantum mechanical rotors, which are candidates for quantum information devices with higher-dimensional local Hilbert spaces. The experimental exploration of endofullerene arrays comes at a time when machine learning techniques are rapidly being adopted to characterize, verify, and reconstruct quantum states from measurement data. In this paper, we develop a strategy for reconstructing the ground state of chains of dipolar rotors using restricted Boltzmann machines (RBMs) adapted to train on data from higher-dimensional Hilbert spaces. We demonstrate accurate generation of energy expectation values from an RBM trained on data in the free-rotor eigenstate basis, and explore the learning resources required for various chain lengths and dipolar interaction strengths. Finally, we show evidence for fundamental limitations in the accuracy achievable by RBMs due to the difficulty in imposing symmetries in the sampling procedure. We discuss possible avenues to overcome this limitation in the future, including the further development of autoregressive models such as recurrent neural networks for the purposes of quantum state reconstruction.

preprint2020arXiv

Recurrent Neural Network Wave Functions

A core technology that has emerged from the artificial intelligence revolution is the recurrent neural network (RNN). Its unique sequence-based architecture provides a tractable likelihood estimate with stable training paradigms, a combination that has precipitated many spectacular advances in natural language processing and neural machine translation. This architecture also makes a good candidate for a variational wave function, where the RNN parameters are tuned to learn the approximate ground state of a quantum Hamiltonian. In this paper, we demonstrate the ability of RNNs to represent several many-body wave functions, optimizing the variational parameters using a stochastic approach. Among other attractive features of these variational wave functions, their autoregressive nature allows for the efficient calculation of physical estimators by providing independent samples. We demonstrate the effectiveness of RNN wave functions by calculating ground state energies, correlation functions, and entanglement entropies for several quantum spin models of interest to condensed matter physicists in one and two spatial dimensions.

preprint2020arXiv

Wavefunction positivization via automatic differentiation

We introduce a procedure to systematically search for a local unitary transformation that maps a wavefunction with a non-trivial sign structure into a positive-real form. The transformation is parametrized as a quantum circuit compiled into a set of one and two qubit gates. We design a cost function that maximizes the average sign of the output state and removes its complex phases. The optimization of the gates is performed through automatic differentiation algorithms, widely used in the machine learning community. We provide numerical evidence for significant improvements in the average sign for a two-leg triangular Heisenberg ladder with next-to-nearest neighbour and ring-exchange interactions. This model exhibits phases where the sign structure can be removed by simple local one-qubit unitaries, but also an exotic Bose-metal phase whose sign structure induces "Bose surfaces" with a fermionic character and a higher entanglement that requires deeper circuits.

preprint2019arXiv

Machine learning quantum states in the NISQ era

We review the development of generative modeling techniques in machine learning for the purpose of reconstructing real, noisy, many-qubit quantum states. Motivated by its interpretability and utility, we discuss in detail the theory of the restricted Boltzmann machine. We demonstrate its practical use for state reconstruction, starting from a classical thermal distribution of Ising spins, then moving systematically through increasingly complex pure and mixed quantum states. Intended for use on experimental noisy intermediate-scale quantum (NISQ) devices, we review recent efforts in reconstruction of a cold atom wavefunction. Finally, we discuss the outlook for future experimental state reconstruction using machine learning, in the NISQ era and beyond.

preprint2016arXiv

A Large-$S$ Study of Quantum Kagome Ice

We present a large-$S$ study of a quantum spin ice Hamiltonian, introduced by Huang, et al., [PRL 112, 167203 (2014)], on the kagome lattice. This model involves a competition between the frustrating Ising term of classical kagome ice, a Zeeman magnetic field $h$,and a nearest-neighbor transverse spin-flip term $S^x_iS^x_j - S^y_iS^y_j$. Recent Quantum Monte Carlo (QMC) simulations by Carrasquilla, et al., [Nature Communications, 6, 7421, (2015)], uncovered lobes of a disordered phase for large Ising interaction and $h\neq 0$ -- a putative quantum spin liquid phase. Here, we examine the nature of this model using large-$S$ expansion. We show that the ground state properties generally have the same trends with those observed in QMC simulations. In particular, the large-$S$ ground state phase diagram captures the existence of the disordered lobes.

preprint2016arXiv

Geometrical mutual information at the tricritical point of the two-dimensional Blume-Capel model

The spin-1 classical Blume-Capel model on a square lattice is known to exhibit a finite-temperature phase transition described by the tricritical Ising CFT in 1+1 space-time dimensions. This phase transition can be accessed with classical Monte Carlo simulations, which, via a replica-trick calculation, can be used to study the shape-dependence of the classical Rényi entropies for a torus divided into two cylinders. From the second Rényi entropy, we calculate the Geometrical Mutual Information (GMI) introduced by Stéphan et. al. [Phys. Rev. Lett. 112, 127204 (2014)] and use it to extract a numerical estimate for the value of the central charge near the tricritical point. By comparing to the known CFT result, $c=7/10$, we demonstrate how this type of GMI calculation can be used to estimate the position of the tricritical point in the phase diagram.

preprint2016arXiv

Learning Thermodynamics with Boltzmann Machines

A Boltzmann machine is a stochastic neural network that has been extensively used in the layers of deep architectures for modern machine learning applications. In this paper, we develop a Boltzmann machine that is capable of modelling thermodynamic observables for physical systems in thermal equilibrium. Through unsupervised learning, we train the Boltzmann machine on data sets constructed with spin configurations importance-sampled from the partition function of an Ising Hamiltonian at different temperatures using Monte Carlo (MC) methods. The trained Boltzmann machine is then used to generate spin states, for which we compare thermodynamic observables to those computed by direct MC sampling. We demonstrate that the Boltzmann machine can faithfully reproduce the observables of the physical system. Further, we observe that the number of neurons required to obtain accurate results increases as the system is brought close to criticality.

preprint2016arXiv

Shape dependence of two-cylinder Renyi entropies for free bosons on a lattice

Universal scaling terms occurring in Renyi entanglement entropies have the potential to bring new understanding to quantum critical points in free and interacting systems. Quantitative comparisons between analytical continuum theories and numerical calculations on lattice models play a crucial role in advancing such studies. In this paper, we exactly calculate the universal two-cylinder shape dependence of entanglement entropies for free bosons on finite-size square lattices, and compare to approximate functions derived in the continuum using several different ansatzes. Although none of these ansatzes are exact in the thermodynamic limit, we find that numerical fits are in good agreement with continuum functions derived using the AdS/CFT correspondence, an extensive mutual information model, and a quantum Lifshitz model. We use fits of our lattice data to these functions to calculate universal scalars defined in the thin-cylinder limit, and compare to values previously obtained for the free boson field theory in the continuum.

preprint2016arXiv

Spinon walk in quantum spin ice

We study a minimal model for the dynamics of spinons in quantum spin ice. The model captures the essential strong coupling between the spinon and the disordered background spins. We demonstrate that the spinon motion can be mapped to a random walk with an entropy-induced memory in imaginary time. Our numerical simulation of the spinon walk indicates that the spinon propagates as a massive quasiparticle at low energy despite its strong coupling to the spin background at the microscopic energy scale. We discuss the experimental implications of our findings.

preprint2016arXiv

Universal corner entanglement of Dirac fermions and gapless bosons from the continuum to the lattice

A quantum critical (QC) fluid exhibits universal subleading corrections to the area law of its entanglement entropies. In two dimensions when the partition involves a corner of angle $θ$, the subleading term is logarithmic with coefficient $a_α(θ)$ for the $α$-Rényi entropy. In the smooth limit $θ\!\to\!π$, $a_1(θ)$ yields the central charge of the stress tensor when the QC point is described by a conformal field theory (CFT). For general Rényi indices and angles, $a_α(θ)$ is richer and few general results exist. We study $a_α(θ)$ focusing on two benchmark CFTs, the free Dirac fermion and boson. We perform numerical lattice calculations to obtain high precision results in $θ,α$ regimes hitherto unexplored. We derive field theory estimates for $a_α(θ)$, including new exact results, and demonstrate an excellent quantitative match with our numerical calculations. We also develop and test strong lower bounds, which apply to both free and interacting QC systems. Finally, we comment on the near collapse of $a_α(θ)$ for various theories, including interacting $O(N)$ models.

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

Detecting Goldstone Modes with Entanglement Entropy

In the face of mounting numerical evidence, Metlitski and Grover [arXiv:1112.5166] have given compelling analytical arguments that systems with spontaneous broken continuous symmetry contain a sub-leading contribution to the entanglement entropy that diverges logarithmically with system size. They predict that the coefficient of this log is a universal quantity that depends on the number of Goldstone modes. In this paper, we confirm the presence of this log term through quantum Monte Carlo calculations of the second Rényi entropy on the spin 1/2 XY model. Devising an algorithm to facilitate convergence of entropy data at extremely low temperatures, we demonstrate that the single Goldstone mode in the ground state can be identified through the coefficient of the log term. Furthermore, our simulation accuracy allows us to obtain an additional geometric constant additive to the Rényi entropy, that matches a predicted fully-universal form obtained from a free bosonic field theory with no adjustable parameters.

preprint2015arXiv

Fluctuating orders and quenched randomness in the cuprates

We study a quasi-2D classical Landau-Ginzburg-Wilson effective field theory in the presence of quenched disorder in which incommensurate charge-density wave and superconducting orders are intertwined. The disorder precludes long-range charge-density wave order, but not superconducting or nematic order. We select three representative sets of input parameters and compute the corresponding charge-density wave structure factors using both large-$N$ techniques and classical Monte Carlo simulations. Where nematicity and superconductivity coexist at low temperature, the peak height of the charge-density wave structure factor decreases monotonically as a function of increasing temperature, unlike what is seen in X-ray experiments on YBa$_2$Cu$_3$O$_{6+x}$. Conversely, where the thermal evolution of the charge-density wave structure factor qualitatively agrees with experiments, the nematic correlation length, computed to one-loop order, is shorter than the charge-density wave correlation length.

preprint2015arXiv

Unusual Corrections to Scaling and Convergence of Universal Renyi Properties at Quantum Critical Points

At a quantum critical point, bipartite entanglement entropies have universal quantities which are subleading to the ubiquitous area law. For Renyi entropies, these terms are known to be similar to the von Neumann entropy, while being much more amenable to numerical and even experimental measurement. We show here that when calculating universal properties of Renyi entropies, it is important to account for unusual corrections to scaling that arise from relevant local operators present at the conical singularity in the multi-sheeted Riemann surface. These corrections grow in importance with increasing Renyi index. We present studies of Renyi correlation functions in the 1+1 transverse-field Ising model (TFIM) using conformal field theory, mapping to free fermions, and series expansions, and the logarithmic entropy singularity at a corner in 2+1 for both free bosonic field theory and the TFIM, using numerical linked cluster expansions. In all numerical studies, accurate results are only obtained when unusual corrections to scaling are taken into account. In the worst case, an analysis ignoring these corrections can get qualitatively incorrect answers, such as predicting a decrease in critical exponents with the Renyi index, when they are actually increasing. We discuss a two-step extrapolation procedure that can be used to account for the unusual corrections to scaling.

preprint2014arXiv

Corner contribution to the entanglement entropy of an O(3) quantum critical point in 2+1 dimensions

The entanglement entropy for a quantum critical system across a boundary with a corner exhibits a sub-leading logarithmic scaling term with a scale-invariant coefficient. Using a Numerical Linked Cluster Expansion, we calculate this universal quantity for a square-lattice bilayer Heisenberg model at its quantum critical point. We find, for this 2+1 dimensional O(3) universality class, that it is thrice the value calculated previously for the Ising universality class. This relation gives substantial evidence that this coefficient provides a measure of the number of degrees of freedom of the theory, analogous to the central charge in a 1+1 dimensional conformal field theory.

preprint2014arXiv

Corner contribution to the entanglement entropy of strongly-interacting O(2) quantum critical systems in 2+1 dimensions

In a D=2+1 quantum critical system, the entanglement entropy across a boundary with a corner contains a subleading logarithmic scaling term with a universal coefficient. It has been conjectured that this coefficient is, to leading order, proportional to the number of field components N in the associated O(N) continuum $ϕ^4$ field theory. Using density matrix renormalization group calculations combined with the powerful numerical linked cluster expansion technique, we confirm this scenario for the O(2) Wilson-Fisher fixed point in a striking way, through direct calculation at the quantum critical points of two very different microscopic models. The value of this corner coefficient is, to within our numerical precision, twice the coefficient of the Ising fixed point. Our results add to the growing body of evidence that this universal term in the Rényi entanglement entropy reflects the number of low-energy degrees of freedom in a system, even for strongly interacting theories.

preprint2014arXiv

Diamagnetism and density wave order in the pseudogap regime of YBa$_2$Cu$_3$O$_{6+x}$

Clear experimental evidence of charge density wave correlations competing with superconducting order in YBCO have thrust their relationship with the pseudogap regime into the spotlight. To aid in characterizing the pseudogap regime, we propose a dimensionless ratio of the diamagnetic susceptibility to the correlation length of the charge density wave correlations. Using Monte Carlo simulations, we compute this ratio on the classical model of Hayward et. al. (Science 343, 1336 (2014)), which describes angular fluctuations of a multicomponent order, capturing both superconducting and density wave correlations. We compare our results with available data on YBa$_2$Cu$_3$O$_{6+x}$, and propose experiments to clarify the value of this dimensionless ratio using existing samples and techniques.

preprint2014arXiv

Geometric mutual information at classical critical points

A practical use of the entanglement entropy in a 1d quantum system is to identify the conformal field theory describing its critical behavior. It is exactly $(c/3)\ln \ell$ for an interval of length $\ell$ in an infinite system, where $c$ is the central charge of the conformal field theory. Here we define the geometric mutual information, an analogous quantity for classical critical points. We compute this for 2d conformal field theories in an arbitrary geometry, and show in particular that for a rectangle cut into two rectangles, it is proportional to $c$. This makes it possible to extract $c$ in classical simulations, which we demonstrate for the critical Ising and 3-state Potts models.

preprint2014arXiv

Quantum Kagome Ice

Actively shought since the turn of the century, two-dimensional quantum spin liquids (QSLs) are exotic phases of matter where magnetic moments remain disordered even at extremely low temperatures. Despite ongoing searches, QSLs remain elusive, due to a lack of concrete knowledge of the microscopic mechanisms that inhibit magnetic order in real materials. Here, we study a theoretical model for a broad class of frustrated magnetic rare-earth pyrochlore materials called "quantum spin ices". When subject to an external magnetic field along the [111] crystallographic direction, the resulting spin interactions contain a mix of geometric frustration and quantum fluctuations in decoupled two-dimensional kagome planes. Using large-scale quantum Monte Carlo simulations, we identify a simple set of interactions sufficient to promote a groundstate with no magnetic long-range order, and a gap to excitations, consistent with a $Z_2$ spin liquid phase. This suggests a systematic experimental procedure to search for two-dimensional QSLs within the broader class of three-dimensional pyrochlore quantum spin ice materials.

preprint2013arXiv

Angular fluctuations of a multi-component order describe the pseudogap regime of the cuprate superconductors

The hole-doped cuprate high temperature superconductors enter the pseudogap regime as their superconducting critical temperature, $T_c$, falls with decreasing hole density. Experiments have probed this regime for over two decades, but we argue that decisive new information has emerged from recent X-ray scattering experiments (arXiv:1207.0915, arXiv:1206.4333, arXiv:1207.3667). The experiments observe incommensurate charge density wave fluctuations whose strength rises gradually over a wide temperature range above $T_c$, but then decreases as the temperature is lowered below $T_c$. We propose a theory in which the superconducting and charge-density wave orders exhibit angular fluctuations in a 6-dimensional space. The theory provides a natural quantitative fit to the X-ray data, and can be a basis for understanding other characteristics of the pseudogap.

preprint2013arXiv

Anomalies in the Entanglement Properties of the Square Lattice Heisenberg Model

We compute the bipartite entanglement properties of the spin-half square-lattice Heisenberg model by a variety of numerical techniques that include valence bond quantum Monte Carlo (QMC), stochastic series expansion QMC, high temperature series expansions and zero temperature coupling constant expansions around the Ising limit. We find that the area law is always satisfied, but in addition to the entanglement entropy per unit boundary length, there are other terms that depend logarithmically on the subregion size, arising from broken symmetry in the bulk and from the existence of corners at the boundary. We find that the numerical results are anomalous in several ways. First, the bulk term arising from broken symmetry deviates from an exact calculation that can be done for a mean-field Neel state. Second, the corner logs do not agree with the known results for non-interacting Boson modes. And, third, even the finite temperature mutual information shows an anomalous behavior as T goes to zero, suggesting that T->0 and L->infinity limits do not commute. These calculations show that entanglement entropy demonstrates a very rich behavior in d>1, which deserves further attention.

preprint2013arXiv

Detecting Classical Phase Transitions with Renyi Mutual Information

By developing a method to represent the Renyi entropies via a replica-trick on classical statistical mechanical systems, we introduce a procedure to calculate the Renyi Mutual Information in any Monte Carlo simulation. Through simulations on several classical models, we demonstrate that the Renyi Mutual Information can detect finite-temperature critical points, and even identify their universality class, without knowledge of an order parameter or other thermodynamic estimators. Remarkably, in addition to critical points mediated by symmetry breaking, the Renyi Mutual Information is able to detect topological vortex-unbinding transitions, as we explicitly demonstrate on simulations of the XY model.

preprint2013arXiv

Entanglement at a Two-Dimensional Quantum Critical Point: a Numerical Linked Cluster Expansion Study

We develop a method to calculate the bipartite entanglement entropy of quantum models, in the thermodynamic limit, using a Numerical Linked Cluster Expansion (NLCE) involving only rectangular clusters. It is based on exact diagonalization of all n x m rectangular clusters at the interface between entangled subsystems A and B. We use it to obtain the Renyi entanglement entropy of the two-dimensional transverse field Ising model, for arbitrary real Renyi index alpha. Extrapolating these results as a function of the order of the calculation, we obtain universal pieces of the entanglement entropy associated with lines and corners at the quantum critical point. They show NLCE to be one of the few methods capable of accurately calculating universal properties of arbitrary Renyi entropies at higher dimensional critical points.

preprint2013arXiv

Entanglement at a Two-Dimensional Quantum Critical Point: a T=0 Projector Quantum Monte Carlo Study

Although the leading-order scaling of entanglement entropy is non-universal at a quantum critical point (QCP), sub-leading scaling can contain universal behaviour. Such universal quantities are commonly studied in non-interacting field theories, however it typically requires numerical calculation to access them in interacting theories. In this paper, we use large-scale T=0 quantum Monte Carlo simulations to examine in detail the second Rényi entropy of entangled regions at the QCP in the transverse-field Ising model in 2+1 space-time dimensions -- a fixed point for which there is no exact result for the scaling of entanglement entropy. We calculate a universal coefficient of a vertex-induced logarithmic scaling for a polygonal entangled subregion, and compare the result to interacting and non-interacting theories. We also examine the shape-dependence of the Rényi entropy for finite-size toroidal lattices divided into two entangled cylinders by smooth boundaries. Remarkably, we find that the dependence on cylinder length follows a shape-dependent function calculated previously by Stephan {\it et al.} [New J. Phys., 15, 015004, (2013)] at the QCP corresponding to the 2+1 dimensional quantum Lifshitz free scalar field theory. The quality of the fit of our data to this scaling function, as well as the apparent cutoff-independent coefficient that results, presents tantalizing evidence that this function may reflect universal behaviour across these and other very disparate QCPs in 2+1 dimensional systems.

preprint2013arXiv

Entanglement in gapless resonating valence bond states

We study resonating-valence-bond (RVB) states on the square lattice of spins and of dimers, as well as SU(N)-invariant states that interpolate between the two. These states are ground states of gapless models, although the SU(2)-invariant spin RVB state is also believed to be a gapped liquid in its spinful sector. We show that the gapless behavior in spin and dimer RVB states is qualitatively similar by studying the Rényi entropy for splitting a torus into two cylinders, We compute this exactly for dimers, showing it behaves similarly to the familiar one-dimensional log term, although not identically. We extend the exact computation to an effective theory believed to interpolate among these states. By numerical calculations for the SU(2) RVB state and its SU(N)-invariant generalizations, we provide further support for this belief. We also show how the entanglement entropy behaves qualitatively differently for different values of the Rényi index $n$, with large values of $n$ proving a more sensitive probe here, by virtue of exhibiting a striking even/odd effect.

preprint2013arXiv

Fate of CP(N-1) fixed points with q-monopoles

We present an extensive quantum Monte Carlo study of the Néel-valence bond solid (VBS) phase transition on rectangular and honeycomb lattice SU($N$) antiferromagnets in sign problem free models. We find that in contrast to the honeycomb lattice and previously studied square lattice systems, on the rectangular lattice for small $N$ a first order Néel-VBS transition is realized. On increasing $N\geq 4$, we observe that the transition becomes continuous and with the {\em same} universal exponents as found on the honeycomb and square lattices (studied here for $N=5,7,10$), providing strong support for a deconfined quantum critical point. Combining our new results with previous numerical and analytical studies we present a general phase diagram of the stability of $\mathbb{CP}^{N-1}$ fixed points with $q$-monopoles.

preprint2013arXiv

Self-Correcting Quantum Memories Beyond the Percolation Threshold

We analyze several high dimensional generalizations of the toric code at nonzero temperature. We find that in large enough dimension, there can be a distinct separation between the critical temperature $T_c$, given by thermodynamic singularities, and the percolation temperature $T_p$, given by the percolation of defects. We argue that the regime $T_p<T<T_c$ is a range of temperatures where a self-correcting quantum memory can operate despite having percolating defects. We present analytic arguments and numerical evidence in support of this scenario, including a mean-field treatment and Monte Carlo simulations. Near $T_c$, simulations observe a large hysteretic behavior, which may have applications by allowing the self-correcting phase to survive in a "superheated" regime.

preprint2012arXiv

A Wang-Landau method for calculating Renyi entropies in finite-temperature quantum Monte Carlo simulations

We implement a Wang-Landau sampling technique in quantum Monte Carlo (QMC) for the purpose of calculating the Renyi entanglement entropies and associated mutual information. The algorithm converges an estimate for an analogue to the density of states for Stochastic Series Expansion QMC allowing a direct calculation of Renyi entropies without explicit thermodynamic integration. We benchmark results for the mutual information on two-dimensional (2D) isotropic and anisotropic Heisenberg models, 2D transverse field Ising model, and 3D Heisenberg model, confirming a critical scaling of the mutual information in cases with a finite-temperature transition. We discuss the benefits and limitations of broad sampling techniques compared to standard importance sampling methods.

preprint2012arXiv

Bridging lattice-scale physics and continuum field theory with quantum Monte Carlo simulations

We discuss designer Hamiltonians---lattice models tailored to be free from sign problems ("de-signed") when simulated with quantum Monte Carlo methods but which still host complex many-body states and quantum phase transitions of interest in condensed matter physics. We focus on quantum spin systems in which competing interactions lead to non-magnetic ground states. These states and the associated quantum phase transitions can be studied in great detail, enabling direct access to universal properties and connections with low-energy effective quantum field theories. As specific examples, we discuss the transition from a Neel antiferromagnet to either a uniform quantum paramagnet or a spontaneously symmetry-broken valence-bond solid in SU(2) and SU(N) invariant spin models. We also discuss anisotropic (XXZ) systems harboring topological Z2 spin liquids and the XY* transition. We briefly review recent progress on quantum Monte Carlo algorithms, including ground state projection in the valence-bond basis and direct computation of the Renyi variants of the entanglement entropy.

preprint2012arXiv

Entanglement scaling in two-dimensional gapless systems

We numerically determine subleading scaling terms in the ground-state entanglement entropy of several two-dimensional (2D) gapless systems, including a Heisenberg model with Néel order, a free Dirac fermion in the π-flux phase, and the nearest-neighbor resonating-valence-bond wavefunction. For these models, we show that the entanglement entropy between cylindrical regions of length x and L - x, extending around a torus of length L, depends upon the dimensionless ratio x/L. This can be well-approximated on finite-size lattices by a function ln(sin(πx/L)), akin to the familiar chord-length dependence in one dimension. We provide evidence, however, that the precise form of this bulk-dependent contribution is a more general function in the 2D thermodynamic limit.

preprint2012arXiv

Thermodynamic singularities in the entanglement entropy at a 2D quantum critical point

We study the bipartite entanglement entropy of the two-dimensional (2D) transverse-field Ising model in the thermodynamic limit using series expansion methods. Expansions are developed for the Renyi entropy around both the small-field and large-field limits, allowing the separate calculation of the entanglement associated with lines and corners at the boundary between sub-systems. Series extrapolations are used to extract subleading power laws and logarithmic singularities as the quantum critical point is approached. In 1D, we find excellent agreement with exact results as well as quantum Monte Carlo simulations. In 2D, we find compelling evidence that the entanglement at a corner is significantly different from a free boson field theory. These results demonstrate the power of the series expansion method for calculating entanglement entropy in interacting systems, a fact that will be particularly useful in future searches for exotic quantum criticality in models with and without the sign problem.

preprint2012arXiv

Topological Sector Fluctuations and Curie Law Crossover in Spin Ice

At low temperatures, a spin ice enters a Coulomb phase - a state with algebraic correlations and topologically constrained spin configurations. In Ho2Ti2O7, we have observed experimentally that this process is accompanied by a non-standard temperature evolution of the wave vector dependent magnetic susceptibility, as measured by neutron scattering. Analytical and numerical approaches reveal signatures of a crossover between two Curie laws, one characterizing the high temperature paramagnetic regime, and the other the low temperature topologically constrained regime, which we call the spin liquid Curie law. The theory is shown to be in excellent agreement with neutron scattering experiments. On a more general footing, i) the existence of two Curie laws appears to be a general property of the emergent gauge field for a classical spin liquid, and ii) sheds light on the experimental difficulty of measuring a precise Curie-Weiss temperature in frustrated materials; iii) the mapping between gauge and spin degrees of freedom means that the susceptibility at finite wave vector can be used as a local probe of fluctuations among topological sectors.

preprint2011arXiv

Finite Temperature Critical Behavior of Mutual Information

We study mutual information for Renyi entropy of arbitrary index n, in interacting quantum systems at finite-temperature critical points, using high-temperature expansion, quantum Monte Carlo simulations and scaling theory. We find that for n>1, the critical behavior is manifest at two temperatures T_c and n*T_c. For the XXZ model with Ising anisotropy, the coefficient of the area-law has a t*ln(t) singularity, whereas the subleading correction from corners has a logarithmic divergence, with a coefficient related to the exact results of Cardy and Peschel. For T<n*T_c there is a constant term associated with broken symmetries that jumps at both T_c and n*T_c, which can be understood in terms of a scaling function analogous to the boundary entropy of Affleck and Ludwig.

preprint2011arXiv

Quantum spin liquid in a kagome lattice spin-1/2 XY model with four-site exchange

We study the ground state phase diagram of a two-dimensional kagome lattice spin-1/2 XY model (J) with a four-site ring exchange interaction (K) using quantum Monte Carlo simulations. We find that the superfluid phase, existing in the regime of small ring exchange, undergoes a direct transition to a Z_2 quantum spin liquid phase at (K/J)_c ~ 22, which is related to the phase proposed by Balents, Girvin and Fisher [Phys. Rev. B, 65 224412 (2002)]. The quantum phase transition between the superfluid and the spin liquid phase has exponents z and νfalling in the 3D XY universality class, making it a candidate for an exotic XY* quantum critical point, mediated by the condensation of bosonic spinons.

preprint2011arXiv

Topological Entanglement Entropy of a Bose-Hubbard Spin Liquid

The Landau paradigm of classifying phases by broken symmetries was demonstrated to be incomplete when it was realized that different quantum Hall states could only be distinguished by more subtle, topological properties. Today, the role of topology as an underlying description of order has branched out to include topological band insulators, and certain featureless gapped Mott insulators with a topological degeneracy in the groundstate wavefunction. Despite intense focus, very few candidates for these topologically ordered "spin liquids" exist. The main difficulty in finding systems that harbour spin liquid states is the very fact that they violate the Landau paradigm, making conventional order parameters non-existent. Here, we uncover a spin liquid phase in a Bose-Hubbard model on the kagome lattice, and measure its topological order directly via the topological entanglement entropy. This is the first smoking-gun demonstration of a non-trivial spin liquid, identified through its entanglement entropy as a gapped groundstate with emergent Z2 gauge symmetry.

preprint2011arXiv

Universal Signatures of Fractionalized Quantum Critical Points

Groundstates of certain materials can support exotic excitations with a charge that's a fraction of the fundamental electron charge. The condensation of these fractionalized particles has been predicted to drive novel quantum phase transitions, which haven't yet been observed in realistic systems. Through numerical and theoretical analysis of a physical model of interacting lattice bosons, we establish the existence of such an exotic critical point, called XY*. We measure a highly non-classical critical exponent eta = 1.49(2), and construct a universal scaling function of winding number distributions that directly demonstrates the distinct topological sectors of an emergent Z_2 gauge field. The universal quantities used to establish this exotic transition can be used to detect other fractionalized quantum critical points in future model and material systems.

preprint2010arXiv

A generalized Monte Carlo loop algorithm for frustrated Ising models

We introduce a Generalized Loop Move (GLM) update for Monte Carlo simulations of frustrated Ising models on two-dimensional lattices with bond-sharing plaquettes. The GLM updates are designed to enhance Monte Carlo sampling efficiency when the system's low-energy states consist of an extensive number of degenerate or near-degenerate spin configurations, separated by large energy barriers to single spin flips. Through implementation on several frustrated Ising models, we demonstrate the effectiveness of the GLM updates in cases where both degenerate and near-degenerate sets of configurations are favored at low temperatures. The GLM update's potential to be straightforwardly extended to different lattices and spin interactions allow it to be readily adopted on many other frustrated Ising models of physical relevance.

preprint2010arXiv

Classical Topological Order in Kagome Ice

We examine the onset of classical topological order in a nearest-neighbor kagome ice model. Using Monte Carlo simulations, we characterize the topological sectors of the groundstate using a non-local cut measure which circumscribes the toroidal geometry of the simulation cell. We demonstrate that simulations which employ global loop updates that are allowed to wind around the periodic boundaries cause the topological sector to fluctuate, while restricted local loop updates freeze the simulation into one topological sector. The freezing into one topological sector can also be observed in the susceptibility of the real magnetic spin vectors projected onto the kagome plane. The ability of the susceptibility to distinguish between fluctuating and non-fluctuating topological sectors should motivate its use as a local probe of topological order in a variety of related model and experimental systems.

preprint2010arXiv

Finite Size Scaling of Mutual Information: A Scalable Simulation

We develop a quantum Monte Carlo procedure to compute the Renyi mutual information of an interacting quantum many-body system at non-zero temperature. Performing simulations on a spin-1/2 XXZ model, we observe that for a subregion of fixed size embedded in a system of size L, the mutual information converges at large L to a limiting function which displays non-monotonic temperature behavior corresponding to the onset of correlations. For a region of size L/2 embedded in a system of size L, the mutual information divided by L converges to a limiting function of temperature, with apparently nontrivial corrections near critical points.

preprint2010arXiv

Measuring Renyi Entanglement Entropy with Quantum Monte Carlo

We develop a quantum Monte Carlo procedure, in the valence bond basis, to measure the Renyi entanglement entropy of a many-body ground state as the expectation value of a unitary {\it Swap} operator acting on two copies of the system. An improved estimator involving the ratio of {\it Swap} operators for different subregions enables simulations to converge the entropy in a time polynomial in the system size. We demonstrate convergence of the Renyi entropy to exact results for a Heisenberg chain. Finally, we calculate the scaling of the Renyi entropy in the two-dimensional Heisenberg model and confirm that the Néel groundstate obeys the expected area law for systems up to linear size L=28.

preprint2010arXiv

Mott Insulators in a Fully-Frustrated Bose Hubbard Model on the Honeycomb Lattice

We examine the effects of quantum fluctuations on a classical spin liquid state in the fully-frustrated honeycomb lattice Bose Hubbard model using quantum Monte Carlo simulations. Frustration is induced explicitly in the model by modulating the sign of the interaction spatially around each lattice hexagon. A superfluid to Mott insulating quantum phase transition can be induced by varying the relative strength of the classical interaction and quantum hopping. In the cases where the interaction has a regular spatial modulation, hopping promotes a phase transition to a symmetry-broken valence-bond solid state. When the interaction is forced to have no regular pattern, the Mott insulating phase is found to be featureless and gapped, making it an interesting candidate state for a quantum spin liquid arising in a Hamiltonian with only nearest-neighbor interactions.

preprint2008arXiv

Superfluid and supersolid phases of lattice bosons with ring-exchange interaction

We examine the superfluid phase of a hard-core boson model with nearest-neighbor exchange J and four-particle ring-exchange K at half-filling on the square lattice. At zero temperature we find that the superfluid in the pure-J model is quickly destroyed by the inclusion of negative-K ring-exchange interactions, favoring a state with a (pi,pi) ordering wavevector. Minimization of the mean-field energy suggests that a supersolid state with coexisting superfluidity, charge-density wave, and valence-bond-like order is formed. We also study the behavior of the finite-T Kosterlitz-Thouless phase transition in the superfluid phase, by forcing the Nelson-Kosterlitz universal jump condition on the finite-T spin wave superfluid density. Away from the pure J point, T_{KT} decreases rapidly for negative K, while for positive K, T_{KT} reaches a maximum at some K \neq 0 in agreement with recent quantum Monte Carlo simulations.

preprint2007arXiv

A short-loop algorithm for quantum Monte Carlo simulations

We present an algorithmic framework for a variant of the quantum Monte Carlo operator-loop algorithm, where non-local cluster updates are constructed in a way that makes each individual loop smaller. The algorithm is designed to increase simulation efficiency in cases where conventional loops become very large, do not close altogether, or otherwise behave poorly. We demonstrate and characterize some aspects of the short-loop on a square lattice spin-1/2 XXZ model where, remarkably, a significant increase in simulation efficiency is observed in some parameter regimes. The simplicity of the model provides a prototype for the use of short-loops on more complicated quantum systems.

preprint2003arXiv

Thermodynamic Properties of the Dipolar Spin Ice Model

We present a detailed theoretical overview of the thermodynamic properties of the dipolar spin ice model, which has been shown to be an excellent quantitative descriptor of the Ising pyrochlore materials Dy_2Ti_2O_7 and Ho_2Ti_2O_7. We show that the dipolar spin ice model can reproduce an effective quasi macroscopically degenerate ground state and spin-ice behavior of these materials when the long-range nature of dipole-dipole interaction is handled carefully using Ewald summation techniques. This degeneracy is, however, ultimately lifted at low temperature. The long-range ordered state is identified via mean field theory and Monte Carlo simulation techniques. Finally, we investigate the behavior of the dipolar spin ice model in an applied magnetic field, and compare our predictions with experimental results. We find that a number of different long-range ordered states are favored by the model depending on field direction.