Researcher profile

Stephen Inglis

Stephen Inglis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2016arXiv

Accessing many-body localized states through the Generalized Gibbs Ensemble

We show how the thermodynamic properties of large many-body localized systems can be studied using quantum Monte Carlo simulations. To this end we devise a heuristic way of constructing local integrals of motion of very high quality, which are added to the Hamiltonian in conjunction with Lagrange multipliers. The ground state simulation of the shifted Hamiltonian corresponds to a high-energy state of the original Hamiltonian in case of exactly known local integrals of motion. We can show that the inevitable mixing between eigenstates as a consequence of non-perfect integrals of motion is weak enough such that the characteristics of many-body localized systems are not averaged out in our approach, unlike the standard ensembles of statistical mechanics. Our method paves the way to study higher dimensions and indicates that a full many-body localized phase in 2d, where (nearly) all eigenstates are localized, is likely to exist.

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

The classical mutual information in mean-field spin glass models

We investigate the classical Renyi entropy S_n and the associated mutual information I_n in the Sherrington-Kirkpatrick (S-K) model, which is the paradigm model of mean-field spin glasses. Using classical Monte Carlo simulations and analytical tools we investigate the S-K model on the n-sheets booklet. This is obtained by gluing together n independent copies of the model, and it is the main ingredient to construct the Renyi entanglement-related quantities. We find a glassy phase at low temperature, whereas at high temperature the model exhibits paramagnetic behavior, consistent with the regular S-K model. The temperature of the paramagnetic-glassy transition depends non-trivially on the geometry of the booklet. At high-temperatures we provide the exact solution of the model by exploiting the replica symmetry. This is the permutation symmetry among the fictitious replicas that are used to perform disorder averages (via the replica trick). In the glassy phase the replica symmetry has to be broken. Using a generalization of the Parisi solution, we provide analytical results for S_n and I_n, and for standard thermodynamic quantities. Both S_n and I_n exhibit a volume law in the whole phase diagram. We characterize the behavior of the corresponding densities S_n/N, I_n/N, in the thermodynamic limit. Interestingly, at the critical point the mutual information does not exhibit any crossing for different system sizes, in contrast with local spin models.

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.

preprint2014arXiv

A path integral Monte Carlo method for Rényi entanglement entropies

We introduce a quantum Monte Carlo algorithm to measure the Rényi entanglement entropies in systems of interacting bosons in the continuum. This approach is based on a path integral ground state method that can be applied to interacting itinerant bosons in any spatial dimension with direct relevance to experimental systems of quantum fluids. We demonstrate how it may be used to compute spatial mode entanglement, particle partitioned entanglement, and the entanglement of particles, providing insights into quantum correlations generated by fluctuations, indistinguishability and interactions. We present proof-of-principle calculations, and benchmark against an exactly soluble model of interacting bosons in one spatial dimension. As this algorithm retains the fundamental polynomial scaling of quantum Monte Carlo when applied to sign-problem-free models, future applications should allow for the study of entanglement entropy in large scale many-body systems of interacting bosons.

preprint2014arXiv

Destroying a topological quantum bit by condensing Ising vortices

The imminent realization of topologically-protected qubits in fabricated systems will provide not only an elementary implementation of fault-tolerant quantum computing architecture, but also an experimental vehicle for the general study of topological order. The simplest topological qubit harbors what is known as a Z$_2$ liquid phase, which encodes information via a degeneracy depending on the system's topology. Elementary excitations of the phase are fractionally charged objects called {\it spinons}, or Ising flux vortices called {\it visons}. At zero temperature a Z$_2$ liquid is stable under deformations of the Hamiltonian until spinon or vison condensation induces a quantum phase transition destroying the topological order. In this paper, we use quantum Monte Carlo to study a vison-induced transition from a Z$_2$ liquid to a valence-bond solid in a quantum dimer model on the kagome lattice. Our results indicate that this critical point is controlled by a new universality class beyond the standard Landau paradigm.

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.

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 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.

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.

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.

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.