Source author record

E. M. Stoudenmire

E. M. Stoudenmire 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

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

15 published item(s)

preprint2020arXiv

DMRG Approach to Optimizing Two-Dimensional Tensor Networks

Tensor network algorithms have been remarkably successful solving a variety of problems in quantum many-body physics. However, algorithms to optimize two-dimensional tensor networks known as PEPS lack many of the aspects that make the seminal density matrix renormalization group (DMRG) algorithm so powerful for optimizing one-dimensional tensor networks known as matrix product states. We implement a framework for optimizing two-dimensional PEPS tensor networks which includes all of steps that make DMRG so successful for optimizing one-dimension tensor networks. We present results for several 2D spin models and discuss possible extensions and applications.

preprint2015arXiv

Assembling Fibonacci Anyons From a $\mathbb{Z}_3$ Parafermion Lattice Model

Recent concrete proposals suggest it is possible to engineer a two-dimensional bulk phase supporting non-Abelian Fibonacci anyons out of Abelian fractional quantum Hall systems. The low-energy degrees of freedom of such setups can be modeled as $\mathbb{Z}_3$ parafermions "hopping" on a two-dimensional lattice. We use the density matrix renormalization group to study a model of this type interpolating between the decoupled-chain, triangular-lattice, and square-lattice limits. The results show clear evidence of the Fibonacci phase over a wide region of the phase diagram, most notably including the isotropic triangular lattice point. We also study the broader phase diagram of this model and show that elsewhere it supports an Abelian state with semionic excitations.

preprint2015arXiv

Many-body localization in disorder-free systems: the importance of finite-size constraints

Recently it has been suggested that many-body localization (MBL) can occur in translation-invariant systems, and candidate 1D models have been proposed. We find that such models, in contrast to MBL systems with quenched disorder, typically exhibit much more severe finite-size effects due to the presence of two or more vastly different energy scales. In a finite system, this can easily create an artificial splitting of the density of states (DOS) into bands separated by large energy gaps. We argue that in order for such models to faithfully represent the physics of the thermodynamic limit, the ratio of the relevant coupling parameters must be larger than a certain cutoff that depends on system size, and should be chosen in such a way that various bands in the DOS of a given model overlap with one another. By setting the parameters in this way to minimize the finite-size effects, we then perform exact diagonalization studies of several translation-invariant MBL candidate models. Based on the variety of diagnostics, including entanglement properties and the behaviour of local observables, we find the systems exhibit thermal (ergodic), rather than MBL-like behaviour. Our results suggest that MBL in translation-invariant systems with two or more very different energy scales is less robust than perturbative arguments suggest, possibly pointing to the importance of non-perturbative effects which induce delocalization in the thermodynamic limit.

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

Kohn-Sham calculations with the exact functional

As a proof of principle, self-consistent Kohn--Sham calculations are performed with the exact exchange-correlation functional. Finding the exact functional for even one trial density requires solving the interacting Schrödinger equation many times. The density matrix renormalization group method makes this possible for one-dimensional, real-space systems of more than two interacting electrons. We illustrate and explore the convergence properties of the exact KS scheme for both weakly and strongly correlated systems. We also explore the spin-dependent generalization and densities for which the functional is ill defined.

preprint2013arXiv

Guaranteed convergence of the Kohn-Sham equations

A sufficiently damped iteration of the Kohn-Sham equations with the exact functional is proven to always converge to the true ground-state density, regardless of the initial density or the strength of electron correlation, for finite Coulomb systems. We numerically implement the exact functional for one-dimensional continuum systems and demonstrate convergence of the damped KS algorithm. More strongly correlated systems converge more slowly.

preprint2013arXiv

Real-Space Parallel Density Matrix Renormalization Group

We demonstrate how to parallelize the density matrix renormalization group (DMRG) algorithm in real space through a straightforward modification of serial DMRG. This makes it possible to apply at least an order of magnitude more computational power to challenging simulations, greatly accelerating investigations of two-dimensional systems and large parameter spaces. We discuss details of the algorithm and present benchmark results including a study of valence-bond-solid order within the square-lattice Q2 model and Neel order within the triangular lattice Heisenberg model. The parallel DMRG algorithm also motivates an alternative canonical form for matrix product states.

preprint2012arXiv

One-dimensional Continuum Electronic Structure with the Density Matrix Renormalization Group and Its Implications For Density Functional Theory

We extend the density matrix renormalization group to compute exact ground states of continuum many-electron systems in one dimension with long-range interactions. We find the exact ground state of a chain of 100 strongly correlated artificial hydrogen atoms. The method can be used to simulate 1d cold atom systems and to study density functional theory in an exact setting. To illustrate, we find an interacting, extended system which is an insulator but whose Kohn-Sham system is metallic.

preprint2012arXiv

Reference electronic structure calculations in one dimension

Large strongly correlated systems provide a challenge to modern electronic structure methods, because standard density functionals usually fail and traditional quantum chemical approaches are too demanding. The density-matrix renormalization group method, an extremely powerful tool for solving such systems, has recently been extended to handle long-range interactions on real-space grids, but is most efficient in one dimension where it can provide essentially arbitrary accuracy. Such 1d systems therefore provide a theoretical laboratory for studying strong correlation and developing density functional approximations to handle strong correlation, {\em if} they mimic three-dimensional reality sufficiently closely. We demonstrate that this is the case, and provide reference data for exact and standard approximate methods, for future use in this area.

preprint2012arXiv

Topological phases in ultracold polar-molecule quantum magnets

We show how to use polar molecules in an optical lattice to engineer quantum spin models with arbitrary spin S >= 1/2 and with interactions featuring a direction-dependent spin anisotropy. This is achieved by encoding the effective spin degrees of freedom in microwave-dressed rotational states of the molecules and by coupling the spins through dipolar interactions. We demonstrate how one of the experimentally most accessible anisotropies stabilizes symmetry protected topological phases in spin ladders. Using the numerically exact density matrix renormalization group method, we find that these interacting phases -- previously studied only in the nearest-neighbor case -- survive in the presence of long-range dipolar interactions. We also show how to use our approach to realize the bilinear-biquadratic spin-1 and the Kitaev honeycomb models. Experimental detection schemes and imperfections are discussed.

preprint2011arXiv

Interaction Effects in Topological Superconducting Wires Supporting Majorana Fermions

Among the broad spectrum of systems predicted to exhibit topological superconductivity and Majorana fermions, one-dimensional wires with strong spin-orbit coupling provide one of the most promising experimental candidates. Here we investigate the fate of the topological superconducting phase in such wires when repulsive interactions are present. Using a combination of Density Matrix Renormalization Group, bosonization, and Hartree-Fock techniques, we demonstrate that while interactions degrade the bulk gap -consistent with recent results of Gangadharaiah et al.- they also greatly expand the parameter range over which the topological phase arises. In particular, we show that with interactions this phase can be accessed over a broader chemical potential window, thereby leading to greater immunity against disorder-induced chemical potential fluctuations in the wire. We also suggest that in certain wires strong interactions may allow Majorana fermions to be generated without requiring a magnetic field.

preprint2011arXiv

Studying Two Dimensional Systems With the Density Matrix Renormalization Group

The Density Matrix Renormalization Group (DMRG) method scales exponentially in the system width for models in two dimensions, but remains one of the most powerful methods for studying 2D systems with a sign problem. Reviewing past applications of DMRG in 2D demonstrates its success in treating a wide variety of problems, although it remains underutilized in this setting. We present techniques for performing cutting edge 2D DMRG studies including methods for ensuring convergence, extrapolating finite-size data and extracting gaps and excited states. Finally, we compare the current performance of a recently developed tensor network method to 2D DMRG.

preprint2009arXiv

Tuning magnetic frustration on the diamond lattice of the A-site magnetic spinels CoAl$_{2-x}$Ga$_x$O$_4$: Lattice expansion and site disorder

The spinels CoB$_2$O$_4$ with magnetic Co$^{2+}$ ions on the diamond lattice A site can be frustrated because of competing near-neighbor ($J_1$) and next-near neighbor ($J_2$) interactions. Here we describe attempts to tune the relative strengths of these interactions by substitution on the non-magnetic B-site. The system we employ is CoAl$_{2-x}$Ga$_x$O$_4$, where Al is systematically replaced by the larger Ga, ostensibly on the B site. As expected, Ga substitution expands the lattice, resulting in Co atoms on the A-site being pushed further from one other and thereby weakening magnetic interactions. In addition, Ga distributes between the B and the A site in a concentration dependent manner displacing an increasing amount of Co from the A site with increasing $x$. This increased inversion, which is confirmed by neutron diffraction studies carried out at room temperature, affects magnetic ordering very significantly, and changes the nature of the ground state. Modeling of the magnetic coupling illustrates the complexity that arises from the cation site disorder.

preprint2004arXiv

Magnetoresistive Effects in Ferromagnet-Superconductor Multilayers

We consider a nanoscale system consisting of Manganite-ferromagnet and Cuprate-superconductor multilayers in a spin valve configuration. The magnetization of the bottom Manganite-ferromagnet is pinned by a Manganite-antiferromagnet. The magnetization of the top Manganite-ferromagnet is coupled to the bottom one via indirect exchange through the superconducting layers. We study the behavior of the critical temperature and the magnetoresistance as a function of an externally applied parallel magnetic field, when the number of Cuprate-superconductor layers are changed. There are two typical behaviors in the case of a few monolayers of the Cuprates: a) For small magnetic fields, the critical temperature and the magnetoresistance change abruptly when the flipping field of the top Manganite-ferromagnet is reached. b) For large magnetic fields, the multilayered system re-enters the zero-resistance (superconducting) state after having become resistive (normal).