BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Steven R. White

Steven R. White appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elphysics.chem-phcond-mat.supr-conquant-phcond-mat.mtrl-scicond-mat.stat-mechphysics.comp-phcond-matcond-mat.otherphysics.atm-clus

Catalog footprint

What is connected

34works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

34 published item(s)

preprint2022arXiv

Conditional probability density functional theory

We present conditional probability (CP) density functional theory (DFT) as a formally exact theory. In essence, CP-DFT determines the ground-state energy of a system by finding the CP density from a series of independent Kohn-Sham (KS) DFT calculations. By directly calculating CP densities, we bypass the need for an approximate XC energy functional. In this work we discuss and derive several key properties of the CP density and corresponding CP-KS potential. Illustrative examples are used throughout to help guide the reader through the various concepts and theory presented. We explore a suitable CP-DFT approximation and discuss exact conditions, limitations, and results for selected examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Absence of superconductivity in the pure two-dimensional Hubbard model

We study the superconducting pairing correlations in the ground state of the doped Hubbard model -- in its original form without hopping beyond nearest neighbor or other perturbing parameters -- in two dimensions at intermediate to strong coupling and near optimal doping. The nature of such correlations has been a central question ever since the discovery of cuprate high-temperature superconductors. Despite unprecedented effort and tremendous progress in understanding the properties of this fundamental model, a definitive answer to whether the ground state is superconducting in the parameter regime most relevant to cuprates has proved exceedingly difficult to establish. In this work, we employ two complementary, state-of-the-art many-body computational methods, constrained path (CP) auxiliary-field quantum Monte Carlo (AFQMC) and density matrix renormalization group (DMRG) methods, deploying the most recent algorithmic advances in each. Systematic and detailed comparisons between the two methods are performed. The DMRG is extremely reliable on small width cylinders, where we use it to validate the AFQMC. The AFQMC is then used to study wide systems as well as fully periodic systems, to establish that we have reached the thermodynamic limit. The ground state is found to be non-superconducting in the moderate to strong coupling regime in the vicinity of optimal hole doping.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bypassing the energy functional in density functional theory: Direct calculation of electronic energies from conditional probability densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP-DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP-DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP-DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ground-state properties of the hydrogen chain: insulator-to-metal transition, dimerization, and magnetic phases

Accurate and predictive computations of the quantum-mechanical behavior of many interacting electrons in realistic atomic environments are critical for the theoretical design of materials with desired properties, and require solving the grand-challenge problem of the many-electron Schrodinger equation. An infinite chain of equispaced hydrogen atoms is perhaps the simplest realistic model for a bulk material, embodying several central themes of modern condensed matter physics and chemistry, while retaining a connection to the paradigmatic Hubbard model. Here we report a combined application of cutting-edge computational methods to determine the properties of the hydrogen chain in its quantum-mechanical ground state. Varying the separation between the nuclei leads to a rich phase diagram, including a Mott phase with quasi long-range antiferromagnetic order, electron density dimerization with power-law correlations, an insulator-to-metal transition and an intricate set of intertwined magnetic orders.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Plaquette versus ordinary $d$-wave pairing in the $t'$-Hubbard model on a width 4 cylinder

The Hubbard model and its extensions are important microscopic models for understanding high- $T_c$ superconductivity in cuprates. In the model with next-nearest-neighbor hopping $t'$ (the $t'$- Hubbard model), pairing is strongly influenced by $t'$ . In particular, a recent study on a width-4 cylinder observed quasi-long-rage superconducting order, associated with a negative $t'$ , which was taken to imply superconductivity in the two-dimensional (2D) limit. In this work we study more carefully pairing in the width-4 $t'$-Hubbard model. We show that in this specific system, the pairing symmetry with $t'<0$ is not the ordinary $d$-wave one would expect in the 2D limit. Instead we observe a so-called plaquette d-wave pairing. The plaquette d-wave exists only on a width-4 cylinder, and so is not representative of the 2D limit. We find that a negative $t'$ suppresses the conventional d-wave, leading to plaquette pairing. In contrast, a different $t''$ coupling acting diagonally on the plaquettes suppresses plaquette pairing, leading to conventional $d$-wave pairing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electronic quasiparticles in the quantum dimer model: density matrix renormalization group results

We study a recently proposed quantum dimer model for the pseudogap metal state of the cuprates. The model contains bosonic dimers, representing a spin-singlet valence bond between a pair of electrons, and fermionic dimers, representing a quasiparticle with spin-$1/2$ and charge $+e$. By density matrix renormalization group calculations on a long but finite cylinder, we obtain the ground-state density distribution of the fermionic dimers for a number of different total densities. From the Friedel oscillations at open boundaries, we deduce that the Fermi surface consists of small hole pockets near $(π/2, π/2)$, and this feature persists up to a doping density of $1/16$. We also compute the entanglement entropy and find that it closely matches the sum of the entanglement entropies of a critical boson and a low density of free fermions. Our results support the existence of a fractionalized Fermi liquid in this model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Entanglement renormalization and wavelets

We establish a precise connection between discrete wavelet transforms (WTs) and entanglement renormalization (ER), a real-space renormalization group transformation for quantum systems on the lattice, in the context of free particle systems. Specifically, we employ Daubechies wavelets to build approximations to the ground state of the critical Ising model, then demonstrate that these states correspond to instances of the multi-scale entanglement renormalization ansatz (MERA), producing the first known analytic MERA for critical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Matrix Product Operators, Matrix Product States, and ab initio Density Matrix Renormalization Group algorithms

Current descriptions of the ab initio DMRG algorithm use two superficially different languages: an older language of the renormalization group and renormalized operators, and a more recent language of matrix product states and matrix product operators. The same algorithm can appear dramatically different when written in the two different vocabularies. In this work, we carefully describe the translation between the two languages in several contexts. First, we describe how to efficiently implement the ab-initio DMRG sweep using a matrix product operator based code, and the equivalence to the original renormalized operator implementation. Next we describe how to implement the general matrix product operator/matrix product state algebra within a pure renormalized operator-based DMRG code. Finally, we discuss two improvements of the ab initio DMRG sweep algorithm motivated by matrix product operator language: Hamiltonian compression, and a sum over operators representation that allows for perfect computational parallelism. The connections and correspondences described here serve to link the future developments with the past, and are important in the efficient implementation of continuing advances in ab initio DMRG and related algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

One Dimensional Mimicking of Electronic Structure: The Case for Exponentials

An exponential interaction is constructed so that one-dimensional atoms and chains of atoms mimic the general behavior of their three-dimensional counterparts. Relative to the more commonly used soft-Coulomb interaction, the exponential greatly diminishes the computational time needed for calculating highly accurate quantities with the density matrix renormalization group. This is due to the use of a small matrix product operator and to exponentially vanishing tails. Furthermore, its more rapid decay closely mimics the screened Coulomb interaction in three dimensions. Choosing parameters to best match earlier calculations, we report results for the one dimensional hydrogen atom, uniform gas, and small atoms and molecules both exactly and in the local density approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pure density functional for strong correlations and the thermodynamic limit from machine learning

We use density-matrix renormalization group, applied to a one-dimensional model of continuum Hamiltonians, to accurately solve chains of hydrogen atoms of various separations and numbers of atoms. We train and test a machine-learned approximation to $F[n]$, the universal part of the electronic density functional, to within quantum chemical accuracy. Our calculation (a) bypasses the standard Kohn-Sham approach, avoiding the need to find orbitals, (b) includes the strong correlation of highly-stretched bonds without any specific difficulty (unlike all standard DFT approximations) and (c) is so accurate that it can be used to find the energy in the thermodynamic limit to quantum chemical accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Compression of Correlation Matrices and an Efficient Method for Forming Matrix Product States of Fermionic Gaussian States

Here we present an efficient and numerically stable procedure for compressing a correlation matrix into a set of local unitary single-particle gates, which leads to a very efficient way of forming the matrix product state (MPS) approximation of a pure fermionic Gaussian state, such as the ground state of a quadratic Hamiltonian. The procedure involves successively diagonalizing subblocks of the correlation matrix to isolate local states which are purely occupied or unoccupied. A small number of nearest neighbor unitary gates isolates each local state. The MPS of this state is formed by applying the many-body version of these gates to a product state. We treat the simple case of compressing the correlation matrix of spinless free fermions with definite particle number in detail, though the procedure is easily extended to fermions with spin and more general BCS states (utilizing the formalism of Majorana modes). We also present a DMRG-like algorithm to obtain the compressed correlation matrix directly from a hopping Hamiltonian. In addition, we discuss a slight variation of the procedure which leads to a simple construction of the multiscale entanglement renormalization ansatz (MERA) of a fermionic Gaussian state, and present a simple picture of orthogonal wavelet transforms in terms of the gate structure we present in this paper. As a simple demonstration we analyze the Su-Schrieffer-Heeger model (free fermions on a 1D lattice with staggered hopping amplitudes).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Doping Asymmetry of a 3-orbital CuO$_2$ Hubbard Model

While both the hole and electron doped cuprates can exhibit $d_{x^2-y^2}$-wave superconductivity, the local distribution of the doped carriers is known to be significantly different with the doped holes going primarily on the O sites while the doped electrons go on the Cu sites. Here we report the results of density matrix renormalization calculations for a three-orbital model of a CuO$_2$ lattice. In addition to the asymmetric dependence of the intra-unit-cell occupation of the Cu and O for hole and electron doping, we find important differences in the longer range spin and charge correlations. As expected, the pair-field response has a $d_{x^2-y^2}$-like structure for both the hole and electron doped systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Space group symmetry fractionalization in a chiral kagome Heisenberg antiferromagnet

The anyonic excitations of a spin-liquid can feature fractional quantum numbers under space-group symmetries. Detecting these fractional quantum numbers, which are analogs of the fractional charge of Laughlin quasiparticles, may prove easier than the direct observation of anyonic braiding and statistics. Motivated by the recent numerical discovery of spin-liquid phases in the kagome Heisenberg antiferromagnet, we theoretically predict the pattern of space group symmetry fractionalization in the kagome lattice chiral spin liquid. We provide a method to detect these fractional quantum numbers in finite-size numerics which is simple to implement in DMRG. Applying these developments to the chiral spin liquid phase of a kagome Heisenberg model, we find perfect agreement between our theoretical prediction and numerical observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spin liquid phase of the $S=\frac{1}{2}$ $J_1-J_2$ Heisenberg model on the triangular lattice

We study the $S=1/2$ Heisenberg model on the triangular lattice with nearest neighbor interaction $J_1$ and next nearest neighbor interaction $J_2$ with the density matrix renormalization group. We are able to study long open cylinders with widths up to 9 lattice spacings. At an intermediate $J_2$ region $0.06 \lesssim J_2/J_1 \lesssim 0.17$, we find evidence for a spin liquid (SL) state with short range spin-spin, bond-bond and chiral correlation lengths, bordered by a classical $120^\circ$ Néel ordered state at small $J_2$ and by a two sub-lattice collinear magnetically ordered state at larger $J_2$. Focusing on $J_2/J_1 = 0.1$, we find a number of signatures of a gapped SL phase: two quasi-degenerate ground states on even cylinders, with an energy gap that decreases exponentially with the cylinder width; a dimerization effect on odd cylinders; and large spin triplet and singlet bulk gaps.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum phases of the frustrated XY models on the honeycomb lattice

Searching for spin liquid states has long been attracting both experimentalists and theorists. In this article, we review recent density matrix renormalization group studies of the spin-1/2 XY model on the honeycomb lattice, with first-neighbor ($J_1 = 1$) and frustrating second-neighbor ($J_2>0$) interactions. For the intermediate frustration regime $0.22\lesssim J_2\lesssim0.36$, there exists a surprising antiferromagnetic Ising phase, with ordered moments pointing along the $z$ axis, despite the absence of any $S_z S_z$ interactions in the Hamiltonian. Surrounding this phase as a function of $J_2$ are antiferromagnetic phases with the moments pointing in the $\textit{x-y}$ plane for small $J_2$ and a close competition between an $\textit{x-y}$ plane magnetic collinear phase and a dimer phase for large values of $J_2$. No spin-liquid phases was found in the XY model even with the third neighbor ($J_3>0$) interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unexpected z-Direction Ising Antiferromagnetic Order in a frustrated Spin-1/2 $J_1-J_2$ XY Model on the Honeycomb Lattice

Using the density matrix renormalization group (DMRG) on wide cylinders, we study the phase diagram of the spin-1/2 XY model on the honeycomb lattice, with first-neighbor ($J_1 = 1$) and frustrating second-neighbor ($J_2>0$) interactions. For the intermediate frustration regime $0.22\lesssim J_2\lesssim0.36$, we find a surprising antiferromagnetic Ising phase, with ordered moments pointing along the z axis, despite the absence of any S_z_z interactions in the Hamiltonian. Surrounding this phase as a function of $J_2$ are antiferromagnetic phases with the moments pointing in the $x-y$ plane for small $J_2$ and a close competition between an $x-y$ plane magnetic collinear phase and a dimer phase for large values of $J_2$. We do not find any spin liquid phases in this model.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Weak plaquette valence bond order in the $S=1/2$ honeycomb $J_1-J_2$ Heisenberg model

Using the density matrix renormalization group, we investigate the $S=1/2$ Heisenberg model on the honeycomb lattice with first- ($J_1$) and second-neighbor ($J_2$) interactions. We are able to study long open cylinders with widths up to 12 lattice spacings. For $J_2/J_1$ near 0.3, we find an apparently paramagnetic phase, bordered by an antiferromagnetic phase for $J_2\lesssim 0.26$ and by a valence bond crystal for $J_2\gtrsim 0.36$. The longest correlation length that we find in this intermediate phase is for plaquette valence bond (PVB) order. This correlation length grows strongly with cylinder circumference, indicating either quantum criticality or weak PVB order.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Enhanced Pairing in the "Checkerboard" Hubbard Ladder

We study signatures of superconductivity in a 2--leg "checkerboard" Hubbard ladder model, defined as a one--dimensional (period 2) array of square plaquettes with an intra-plaquette hopping $t$ and inter-plaquette hopping $t'$, using the density matrix renormalization group method. The highest pairing scale (characterized by the spin gap or the pair binding energy, extrapolated to the thermodynamic limit) is found for doping levels close to half filling, $U\approx 6t$ and $t'/t \approx 0.6$. Other forms of modulated hopping parameters, with periods of either 1 or 3 lattice constants, are also found to enhance pairing relative to the uniform two--leg ladder, although to a lesser degree. A calculation of the phase stiffness of the ladder reveals that in the regime with the strongest pairing, the energy scale associated with phase ordering is comparable to the pairing scale.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Incommensurate correlations in the anisotropic triangular Heisenberg lattice

We study the anisotropic spin-half antiferromagnetic triangular Heisenberg lattice in two dimensions, seen as a set of chains with couplings J (J') along (in between) chains, respectively. Our focus is on the incommensurate correlation that emerges in this system in a wide parameter range due to the intrinsic frustration of the spins. We study this system with traditional DMRG using cylindrical boundary conditions to least constrain possible incommensurate order. Despite that the limit of essentially decoupled chains J'/J < 0.5 is not very accessible numerically, it appears that the spin-spin correlations remain incommensurate for any finite 0 < J' < Jc', where Jc'/J > 1. The incommensurate wave vector q_J, however, approaches the commensurate value corresponding to the antiferromagnetic correlation of a single chain very rapidly with decreasing J'/J, roughly as q_J ~ pi - c_1 (J'/J)^n exp(-c_2 J/J').

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient Matrix Product State Method for periodic boundary conditions

We introduce an efficient method to calculate the ground state of one-dimensional lattice models with periodic boundary conditions. The method works in the representation of Matrix Product States (MPS), related to the Density Matrix Renormalization Group (DMRG) method. It improves on a previous approach by Verstraete et al. We introduce a factorization procedure for long products of MPS matrices, which reduces the computational effort from m^5 to m^3, where m is the matrix dimension, and m ~ 100 - 1000 in typical cases. We test the method on the S=1/2 and S=1 Heisenberg chains. It is also applicable to non-translationally invariant cases. The new method makes ground state calculations with periodic boundary conditions about as efficient as traditional DMRG calculations for systems with open boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spin Liquid Ground State of the $S=1/2$ Kagome Heisenberg Model

Condensed matter physicists have long sought a realistic two-dimensional (2D) magnetic system whose ground state is a {\it spin liquid}---a zero temperature state in which quantum fluctuations have melted away any form of magnetic order. The nearest-neighbor $S=1/2$ Heisenberg model on the kagome lattice has seemed an ideal candidate, but in recent years some approximate numerical approaches to it have yielded instead a valence bond crystal. We have used the density matrix renormalization group to perform very accurate simulations on numerous cylinders with circumferences up to 12 lattice spacings, finding instead of the valence bond crystal a singlet-gapped spin liquid with substantially lower energy that appears to have $Z_2$ topological order. Our results, through a combination of very low energy, short correlation lengths and corresponding small finite size effects, a new rigorous energy bound, and consistent behavior on many cylinders, provide strong evidence that the 2D ground state of this model is a gapped spin liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Minimally entangled typical quantum states at finite temperature

We introduce a class of states, called minimally entangled typical thermal states (METTS), designed to resemble a typical state of a quantum system at finite temperature with a bias towards classical (minimally entangled) properties. These states reveal in an intuitive way properties such as short-range order which may be hidden in correlation functions. An algorithm is presented which, when used with the density matrix renormalization group (DMRG), is faster by a factor of $10^3 - 10^{10}$ than previous heat-bath approaches for thermally averaged quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Pairing on striped $t$-$t'$-$J$ lattices

Results are given from a density matrix renormalization group study of pairing on a striped $t$-$t'$-$J$ lattice in the presence of boundary magnetic and pair fields. We find that pairing on a stripe depends sensitively on both $J/t$ and $t'/t$. In the strong-pairing model-parameter regime the stripes are easily coupled by the pair field, and have a uniform phase. There is a small but measurable energy cost to create anti-phase superconducting domain walls.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Neel order in square and triangular lattice Heisenberg models

Using examples of the square- and triangular-lattice Heisenberg models we demonstrate that the density matrix renormalization group method (DMRG) can be effectively used to study magnetic ordering in two-dimensional lattice spin models. We show that local quantities in DMRG calculations, such as the on-site magnetization M, should be extrapolated with the truncation error, not with its square root, as previously assumed. We also introduce convenient sequences of clusters, using cylindrical boundary conditions and pinning magnetic fields, which provide for rapidly converging finite-size scaling. This scaling behavior on our clusters is clarified using finite-size analysis of the effective sigma-model and finite-size spin-wave theory. The resulting greatly improved extrapolations allow us to determine the thermodynamic limit of M for the square lattice with an error comparable to quantum Monte Carlo. For the triangular lattice, we verify the existence of three-sublattice magnetic order, and estimate the order parameter to be M = 0.205(15).

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Density matrix renormalization group algorithms with a single center site

We develop a correction to the density matrix used in density matrix renormalization group calculations to take into account the incompleteness of the environment block. The correction allows successful calculations using only a single site in the center of the system, rather than the standard two sites, improving typical computation times by a factor of two to four. In addition, in many cases where ordinary DMRG can get stuck in metastable configurations, the correction eliminates the sticking. We test the new method on the Heisenberg S=1 chain.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

On the Nagaoka polaron in the t-J model

It is widely believed that a single hole in the two (or three) dimensional t-J model, for sufficiently small exchange coupling J, creates a ferromagnetic bubble around itself, a finite J remnant of the ferromagnetic groundstate at J=0 (the infinite U Hubbard model), first established by Nagaoka. We investigate this phenomenon in two dimensions using the density matrix renormalization group, for system sizes up to 9x9. We find that the polaron forms for J/t<0.02-0.03 (a somewhat larger value than estimated previously). Although finite-size effects appear large, our data seems consistent with the expected 1.1(J/t)^{-1/4} variation of polarion radius. We also test the Brinkman-Rice model of non-retracing paths in a Neel background, showing that it is quite accurate, at larger J. Results are also presented in the case where the Heisenberg interaction is dropped (the t-J^z model). Finally we discuss a "dressed polaron" picture in which the hole propagates freely inside a finite region but makes only self-retracing excursions outside this region.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Equivalence of the antiferromagnetic Heisenberg ladder to a single $S=1$ chain

I introduce two continuous transformations between the $S=1$ Heisenberg chain and the antiferromagnetic $S=1/2$ Heisenberg ladder. Both transformations couple diagonally situated {\it next nearest neighbor} $S=1/2$'s to form each $S=1$. Using the density matrix renormalization group, I demonstrate that the two systems are in the same phase. Furthermore, I find that the hidden topological long-range order characterizing the $S=1$ system is even stronger in the isotropic two-chain system.

0 citations0 reviewsNo code linkedOpen access