-
Unusual energy spectra of matrix product states
Authors:
J. Maxwell Silvester,
Giuseppe Carleo,
Steven R. White
Abstract:
In the simulation of ground states of strongly-correlated quantum systems, the decomposition of an approximate solution into the exact eigenstates of the Hamiltonian -- the energy spectrum of the state -- determines crucial aspects of the simulation's performance. For example, in approaches based on imaginary-time evolution, the spectrum falls off exponentially with the energy, ensuring rapid conv…
▽ More
In the simulation of ground states of strongly-correlated quantum systems, the decomposition of an approximate solution into the exact eigenstates of the Hamiltonian -- the energy spectrum of the state -- determines crucial aspects of the simulation's performance. For example, in approaches based on imaginary-time evolution, the spectrum falls off exponentially with the energy, ensuring rapid convergence. Here we consider the energy spectra of approximate matrix product state ground states, such as those obtained with the density matrix renormalization group. Despite the high accuracy of these states, contributions to the spectra are roughly constant out to surprisingly high energy, with an increase in bond dimension reducing the amplitude but not the extent of these high-energy tails. The unusual spectra, which appear to be a general feature of compressed wavefunctions, have a strong effect on sampling-based methods, yielding large fluctuations. For example, estimating the energy variance using sampling performs much more poorly than one might expect. Bounding the most extreme samples makes the variance estimate much less noisy but introduces a strong bias. However, we find that this biased variance estimator is an excellent surrogate for the variance when extrapolating the ground-state energy, and this approach outperforms competing extrapolation methods in both accuracy and computational cost.
△ Less
Submitted 24 August, 2024;
originally announced August 2024.
-
Nested Gausslet Basis Sets
Authors:
Steven R. White,
Michael J. Lindsey
Abstract:
We introduce nested gausslet (NG) bases, an improvement on previous gausslet bases which can treat systems containing atoms with much larger atomic number. We also introduce pure Gaussian distorted gausslet bases, which allow the Hamiltonian integrals to be performed analytically, as well as hybrid bases in which the gausslets are combined with standard Gaussian-type bases. All these bases feature…
▽ More
We introduce nested gausslet (NG) bases, an improvement on previous gausslet bases which can treat systems containing atoms with much larger atomic number. We also introduce pure Gaussian distorted gausslet bases, which allow the Hamiltonian integrals to be performed analytically, as well as hybrid bases in which the gausslets are combined with standard Gaussian-type bases. All these bases feature the diagonal approximation for the electron-electron interactions, so that the Hamiltonian is completely defined by two $N_b\times N_b$ matrices, where $N_b \approx 10^4$ is small enough to permit fast calculations at the Hartree-Fock level. In constructing these bases we have gained new mathematical insight into the construction of one-dimensional diagonal bases. In particular we have proved an important theorem relating four key basis set properties: completeness, orthogonality, zero-moment conditions, and diagonalization of the coordinate operator matrix. We test our basis sets on small systems with a focus on high accuracy, obtaining, for example, an accuracy of $2\times10^{-5}$ Ha for the total Hartree-Fock energy of the neon atom in the complete basis set limit.
△ Less
Submitted 19 September, 2023;
originally announced September 2023.
-
A Novel Method of Function Extrapolation Inspired by Techniques in Low-entangled Many-body Physics
Authors:
Lambert Lin,
Steven R White
Abstract:
We introduce a novel extrapolation algorithm inspired by quantum mechanics and evaluate its performance against linear prediction. Our method involves mapping function values onto a quantum state and estimating future function values by minimizing entanglement entropy. We demonstrate the effectiveness of our approach on various simple functions, both with and without noise, comparing it to linear…
▽ More
We introduce a novel extrapolation algorithm inspired by quantum mechanics and evaluate its performance against linear prediction. Our method involves mapping function values onto a quantum state and estimating future function values by minimizing entanglement entropy. We demonstrate the effectiveness of our approach on various simple functions, both with and without noise, comparing it to linear prediction. Our results show that the proposed algorithm produces extrapolations comparable to linear prediction, while exhibiting improved performance for functions with sharp features.
△ Less
Submitted 12 October, 2023; v1 submitted 17 August, 2023;
originally announced August 2023.
-
Validating phase-space methods with tensor networks in two-dimensional spin models with power-law interactions
Authors:
Sean R. Muleady,
Mingru Yang,
Steven R. White,
Ana Maria Rey
Abstract:
Using a recently developed extension of the time-dependent variational principle for matrix product states, we evaluate the dynamics of 2D power-law interacting XXZ models, implementable in a variety of state-of-the-art experimental platforms. We compute the spin squeezing as a measure of correlations in the system, and compare to semiclassical phase-space calculations utilizing the discrete trunc…
▽ More
Using a recently developed extension of the time-dependent variational principle for matrix product states, we evaluate the dynamics of 2D power-law interacting XXZ models, implementable in a variety of state-of-the-art experimental platforms. We compute the spin squeezing as a measure of correlations in the system, and compare to semiclassical phase-space calculations utilizing the discrete truncated Wigner approximation (DTWA). We find the latter efficiently and accurately captures the scaling of entanglement with system size in these systems, despite the comparatively resource-intensive tensor network representation of the dynamics. We also compare the steady-state behavior of DTWA to thermal ensemble calculations with tensor networks. Our results open a way to benchmark dynamical calculations for two-dimensional quantum systems, and allow us to rigorously validate recent predictions for the generation of scalable entangled resources for metrology in these systems.
△ Less
Submitted 30 May, 2023; v1 submitted 26 May, 2023;
originally announced May 2023.
-
Quantum Phases in the Honeycomb-Lattice $J_1$--$J_3$ Ferro-Antiferromagnetic Model
Authors:
Shengtao Jiang,
Steven R. White,
A. L. Chernyshev
Abstract:
Using large-scale density-matrix renormalzation group calculations and minimally augmented spin-wave theory, we demonstrate that the phase diagram of the quantum $S\!=\!\frac12$ $J_1$--$J_3$ ferro-antiferromagnetic model on the honeycomb lattice differs dramatically from the classical one. It hosts the double-zigzag and Ising-z phases as unexpected intermediaries between ferromagnetic and zigzag s…
▽ More
Using large-scale density-matrix renormalzation group calculations and minimally augmented spin-wave theory, we demonstrate that the phase diagram of the quantum $S\!=\!\frac12$ $J_1$--$J_3$ ferro-antiferromagnetic model on the honeycomb lattice differs dramatically from the classical one. It hosts the double-zigzag and Ising-z phases as unexpected intermediaries between ferromagnetic and zigzag states that are also extended beyond their classical regions of stability. In broad agreement with quantum order-by-disorder arguments, these collinear phases replace the classical spiral state.
△ Less
Submitted 29 March, 2024; v1 submitted 12 April, 2023;
originally announced April 2023.
-
Coexistence of superconductivity with partially filled stripes in the Hubbard model
Authors:
Hao Xu,
Chia-Min Chung,
Mingpu Qin,
Ulrich Schollwöck,
Steven R. White,
Shiwei Zhang
Abstract:
Combining the complementary capabilities of two of the most powerful modern computational methods, we find superconductivity in both the electron- and hole-doped regimes of the two-dimensional Hubbard model (with next nearest neighbor hopping). In the electron-doped regime, superconductivity is weaker and is accompanied by antiferromagnetic Néel correlations at low doping. The strong superconducti…
▽ More
Combining the complementary capabilities of two of the most powerful modern computational methods, we find superconductivity in both the electron- and hole-doped regimes of the two-dimensional Hubbard model (with next nearest neighbor hopping). In the electron-doped regime, superconductivity is weaker and is accompanied by antiferromagnetic Néel correlations at low doping. The strong superconductivity on the hole-doped side coexists with stripe order, which persists into the overdoped region with weaker hole density modulation. These stripe orders, neither filled as in the pure Hubbard model (no next nearest neighbor hopping) nor half-filled as seen in previous state-of-the-art calculations, vary in fillings between 0.6 and 0.8. The resolution of the tiny energy scales separating competing orders requires exceedingly high accuracy combined with averaging and extrapolating with a wide range of system sizes and boundary conditions. These results validate the applicability of this iconic model for describing cuprate high-$T_c$ superconductivity.
△ Less
Submitted 15 March, 2023;
originally announced March 2023.
-
Density-matrix-renormalization-group-based downfolding of the three-band Hubbard model: the importance of density-assisted hopping
Authors:
Shengtao Jiang,
Douglas J. Scalapino,
Steven R. White
Abstract:
Typical Wannier-function downfolding starts with a mean-field or density functional set of bands to construct the Wannier functions. Here we carry out a controlled approach, using DMRG-computed natural orbital bands, to downfold the three-band Hubbard model to an effective single band model. A sharp drop-off in the natural orbital occupancy at the edge of the first band provides a clear justificat…
▽ More
Typical Wannier-function downfolding starts with a mean-field or density functional set of bands to construct the Wannier functions. Here we carry out a controlled approach, using DMRG-computed natural orbital bands, to downfold the three-band Hubbard model to an effective single band model. A sharp drop-off in the natural orbital occupancy at the edge of the first band provides a clear justification for a single-band model. Constructing Wannier functions from the first band, we compute all possible two-particle terms and retain those with significant magnitude. The resulting single-band model includes two-site density-assisted hopping terms with $t_n \sim 0.6 t$. These terms lead to a reduction of the ratio $U/t_{\rm eff}$, and are important in capturing the doping-asymmetric carrier mobility, as well as in enhancing the pairing in a single-band model for the hole-doped cuprates.
△ Less
Submitted 18 July, 2023; v1 submitted 1 March, 2023;
originally announced March 2023.
-
Variational Benchmarks for Quantum Many-Body Problems
Authors:
Dian Wu,
Riccardo Rossi,
Filippo Vicentini,
Nikita Astrakhantsev,
Federico Becca,
Xiaodong Cao,
Juan Carrasquilla,
Francesco Ferrari,
Antoine Georges,
Mohamed Hibat-Allah,
Masatoshi Imada,
Andreas M. Läuchli,
Guglielmo Mazzola,
Antonio Mezzacapo,
Andrew Millis,
Javier Robledo Moreno,
Titus Neupert,
Yusuke Nomura,
Jannes Nys,
Olivier Parcollet,
Rico Pohle,
Imelda Romero,
Michael Schmid,
J. Maxwell Silvester,
Sandro Sorella
, et al. (8 additional authors not shown)
Abstract:
The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, ide…
▽ More
The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, identifying cases where state-of-the-art numerical approaches show limited accuracy, and novel algorithms or computational platforms, such as quantum computing, could provide improved accuracy. The V-score can be used as a metric to assess the progress of quantum variational methods towards quantum advantage for ground-state problems, especially in regimes where classical verifiability is impossible.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
Quantum Fourier Transform Has Small Entanglement
Authors:
Jielun Chen,
E. M. Stoudenmire,
Steven R. White
Abstract:
The Quantum Fourier Transform (QFT) is a key component of many important quantum algorithms, most famously as being the essential ingredient in Shor's algorithm for factoring products of primes. Given its remarkable capability, one would think it can introduce large entanglement to qubit systems and would be difficult to simulate classically. While early results showed QFT indeed has maximal opera…
▽ More
The Quantum Fourier Transform (QFT) is a key component of many important quantum algorithms, most famously as being the essential ingredient in Shor's algorithm for factoring products of primes. Given its remarkable capability, one would think it can introduce large entanglement to qubit systems and would be difficult to simulate classically. While early results showed QFT indeed has maximal operator entanglement, we show that this is entirely due to the bit reversal in the QFT. The core part of the QFT has Schmidt coefficients decaying exponentially quickly, and thus it can only generate a constant amount of entanglement regardless of the number of qubits. In addition, we show the entangling power of the QFT is the same as the time evolution of a Hamiltonian with exponentially decaying interactions, and thus a variant of the area law for dynamics can be used to understand the low entanglement intuitively. Using the low entanglement property of the QFT, we show that classical simulations of the QFT on a matrix product state with low bond dimension only take time linear in the number of qubits, providing a potential speedup over the classical fast Fourier transform (FFT) on many classes of functions. We demonstrate this speedup in test calculations on some simple functions. For data vectors of length $10^6$ to $10^8$, the speedup can be a few orders of magnitude.
△ Less
Submitted 27 October, 2023; v1 submitted 16 October, 2022;
originally announced October 2022.
-
Where is the Quantum Spin Nematic?
Authors:
Shengtao Jiang,
Judit Romhányi,
Steven R. White,
M. E. Zhitomirsky,
A. L. Chernyshev
Abstract:
We provide strong evidence of the spin-nematic state in a paradigmatic ferro-antiferromagnetic $J_1$--$J_2$ model using analytical and density-matrix renormalization group methods. In zero field, the attraction of spin-flip pairs leads to a first-order transition and no nematic state, while pair-repulsion at larger $J_2$ stabilizes the nematic phase in a narrow region near the pair-condensation fi…
▽ More
We provide strong evidence of the spin-nematic state in a paradigmatic ferro-antiferromagnetic $J_1$--$J_2$ model using analytical and density-matrix renormalization group methods. In zero field, the attraction of spin-flip pairs leads to a first-order transition and no nematic state, while pair-repulsion at larger $J_2$ stabilizes the nematic phase in a narrow region near the pair-condensation field. A devil's staircase of multi-pair condensates is conjectured for weak pair-attraction. A suppression of the spin-flip gap by many-body effects leads to an order-of-magnitude contraction of the nematic phase compared to naïve expectations. The proposed phase diagram should be broadly valid.
△ Less
Submitted 15 March, 2023; v1 submitted 31 August, 2022;
originally announced September 2022.
-
Pairing Properties of the $t$-$t'$-$t''$-$J$ model
Authors:
Shengtao Jiang,
Douglas J. Scalapino,
Steven R. White
Abstract:
We study the pairing properties of the two-dimensional $t$-$t'$-$t''$-$J$ model, where $t'$ and $t''$ are second and third neighbor hoppings, at a doping level $x\approx0.1$. Recent studies of the $t$-$t'$-$J$ model find strong pairing for $t'>0$, associated with electron doping, but an absence of pairing for $t'<0$ associated with hole doping. This is in contrast to the cuprates, where the highes…
▽ More
We study the pairing properties of the two-dimensional $t$-$t'$-$t''$-$J$ model, where $t'$ and $t''$ are second and third neighbor hoppings, at a doping level $x\approx0.1$. Recent studies of the $t$-$t'$-$J$ model find strong pairing for $t'>0$, associated with electron doping, but an absence of pairing for $t'<0$ associated with hole doping. This is in contrast to the cuprates, where the highest transition temperatures appear for hole doping. Model parameterizations for the cuprates estimate a $t''$ comparable to $t'$, which, in principle, might fix this discrepancy. However, we find that it does not; we observe a suppression of pairing for the hole-doped system ($t'<0, t''>0$) while for the electron-doped system ($t'>0, t''<0$) $d$-wave pairing is robust. Extended hoppings appear to be insufficient to make the one-band $t$-$t'$-$J$ model capable of describing the pairing in the hole doped system.
△ Less
Submitted 15 June, 2022;
originally announced June 2022.
-
Conditional probability density functional theory
Authors:
Ryan Pederson,
Jielun Chen,
Steven R. White,
Kieron Burke
Abstract:
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 p…
▽ More
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.
△ Less
Submitted 24 June, 2022; v1 submitted 17 March, 2022;
originally announced March 2022.
-
Bayesian profile regression for clustering analysis involving a longitudinal response and explanatory variables
Authors:
Anaïs Rouanet,
Rob Johnson,
Magdalena E Strauss,
Sylvia Richardson,
Brian D Tom,
Simon R White,
Paul D W Kirk
Abstract:
The identification of sets of co-regulated genes that share a common function is a key question of modern genomics. Bayesian profile regression is a semi-supervised mixture modelling approach that makes use of a response to guide inference toward relevant clusterings. Previous applications of profile regression have considered univariate continuous, categorical, and count outcomes. In this work, w…
▽ More
The identification of sets of co-regulated genes that share a common function is a key question of modern genomics. Bayesian profile regression is a semi-supervised mixture modelling approach that makes use of a response to guide inference toward relevant clusterings. Previous applications of profile regression have considered univariate continuous, categorical, and count outcomes. In this work, we extend Bayesian profile regression to cases where the outcome is longitudinal (or multivariate continuous) and provide PReMiuMlongi, an updated version of PReMiuM, the R package for profile regression. We consider multivariate normal and Gaussian process regression response models and provide proof of principle applications to four simulation studies. The model is applied on budding yeast data to identify groups of genes co-regulated during the Saccharomyces cerevisiae cell cycle. We identify 4 distinct groups of genes associated with specific patterns of gene expression trajectories, along with the bound transcriptional factors, likely involved in their co-regulation process.
△ Less
Submitted 8 November, 2021;
originally announced November 2021.
-
Constructing Hubbard Models for the Hydrogen Chain using Sliced Basis DMRG
Authors:
Randy C. Sawaya,
Steven R. White
Abstract:
Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like mod…
▽ More
Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like model, whose parameters are determined by minimizing the expectation value of the Wannier Hamiltonian in the ground state of the Hubbard Hamiltonian. This indirect variational procedure not only yields compact and simple models for the hydrogen chain, but also allows us to explore the importance of constraints in the effective Hamiltonian, such as the restricting the range of the single-particle hopping and two-particle interactions, and to assess the reliability of more conventional downfolding. The entanglement entropy for a model's ground state, cut in the middle, is an important property determining the ability of DMRG and tensor networks to simulate the model, and we study its variation with the range of the interactions. Counterintuitively, we find that shorter ranged interactions often have larger entanglement.
△ Less
Submitted 10 September, 2021;
originally announced September 2021.
-
Ground State Phase Diagram of the $t$-$t'$-$J$ model
Authors:
Shengtao Jiang,
Douglas J. Scalapino,
Steven R. White
Abstract:
We report results of large scale ground state density matrix renormalization group(DMRG) calculations on $t$-$t'$-$J$ cylinders with circumferences 6 and 8. We determine a rough phase diagram which appears to approximate the 2D system. While for many properties, positive and negative $t'$ values ($t'/t = \pm 0.2$) appear to correspond to electron and hole doped cuprate systems, respectively, the b…
▽ More
We report results of large scale ground state density matrix renormalization group(DMRG) calculations on $t$-$t'$-$J$ cylinders with circumferences 6 and 8. We determine a rough phase diagram which appears to approximate the 2D system. While for many properties, positive and negative $t'$ values ($t'/t = \pm 0.2$) appear to correspond to electron and hole doped cuprate systems, respectively, the behavior of superconductivity itself shows an inconsistency between the model and the materials. The $t'<0$ (hole doped) region shows antiferromagnetism limited to very low doping, stripes more generally, and the familiar Fermi-surface of the hole doped cuprates. However, we find $t'<0$ strongly suppresses superconductivity. The $t'>0$ (electron doped) region shows the expected circular Fermi pocket of holes around the $(π,π)$ point and a broad low-doped region of coexisting antiferromagnetism and $d$-wave pairing with a strong triplet $p$ component at wave-vector $(π,π)$ induced by the antiferromagnetism and $d$-wave pairing. The pairing for the electron low-doped system with $t'>0$ is strong and unambiguous in the DMRG simulations. At larger doping another broad region with stripes in addition to weaker $d$-wave pairing and striped $p$-wave pairing appears. In a small doping region near $x=0.08$ for $t'\sim-0.2$, we find a new type of stripe involving unpaired holes located predominantly on chains spaced three lattice spacings apart. The undoped two-leg ladder regions in between mimic the short-ranged spin correlations seen in two-leg Heisenberg ladders.
△ Less
Submitted 19 May, 2021; v1 submitted 20 April, 2021;
originally announced April 2021.
-
Hybrid gausslet/Gaussian basis sets
Authors:
Yiheng Qiu,
Steven R. White
Abstract:
We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part…
▽ More
We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian, so that it continues to scale as the second power of the number of basis functions, rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction which compensates for the incompleteness of the basis stemming from the electron-electron cusps, based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems, and Hartree Fock on a ten-atom hydrogen chain, to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results, for the two electron cases, at the level of several microHartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.
△ Less
Submitted 3 March, 2021;
originally announced March 2021.
-
Matrix product state recursion methods for strongly correlated quantum systems
Authors:
Yifan Tian,
Steven R. White
Abstract:
We present a method for extrapolation of real-time dynamical correlation functions which can improve the capability of matrix product state methods to compute spectral functions. Unlike the widely used linear prediction method, which ignores the origin of the data being extrapolated, our recursion methods utilize a representation of the wavefunction in terms of an expansion of the same wavefunctio…
▽ More
We present a method for extrapolation of real-time dynamical correlation functions which can improve the capability of matrix product state methods to compute spectral functions. Unlike the widely used linear prediction method, which ignores the origin of the data being extrapolated, our recursion methods utilize a representation of the wavefunction in terms of an expansion of the same wavefunction and its translations at earlier times. This recursion method is exact for a noninteracting Fermi system. Surprisingly, the recursion method is also more robust than linear prediction at large interaction strength. We test this method on the Hubbard two-leg ladder and present more accurate results for the spectral function than previous studies.
△ Less
Submitted 3 December, 2020; v1 submitted 1 October, 2020;
originally announced October 2020.
-
Stripes, Antiferromagnetism, and the Pseudogap in the Doped Hubbard Model at Finite Temperature
Authors:
Alexander Wietek,
Yuan-Yao He,
Steven R. White,
Antoine Georges,
E. Miles Stoudenmire
Abstract:
The interplay between thermal and quantum fluctuations controls the competition between phases of matter in strongly correlated electron systems. We study finite-temperature properties of the strongly coupled two-dimensional doped Hubbard model using the minimally-entangled typical thermal states (METTS) method on width $4$ cylinders. We discover that a phase characterized by commensurate short-ra…
▽ More
The interplay between thermal and quantum fluctuations controls the competition between phases of matter in strongly correlated electron systems. We study finite-temperature properties of the strongly coupled two-dimensional doped Hubbard model using the minimally-entangled typical thermal states (METTS) method on width $4$ cylinders. We discover that a phase characterized by commensurate short-range antiferromagnetic correlations and no charge ordering occurs at temperatures above the half-filled stripe phase extending to zero temperature. The transition from the antiferromagnetic phase to the stripe phase takes place at temperature $T/t \approx 0.05$ and is accompanied by a step-like feature of the specific heat. We find the single-particle gap to be smallest close to the nodal point at $\mathbf{k}=(π/2, π/2)$ and detect a maximum in the magnetic susceptibility. These features bear a strong resemblance to the pseudogap phase of high-temperature cuprate superconductors. The simulations are verified using a variety of different unbiased numerical methods in the three limiting cases of zero temperature, small lattice sizes, and half-filling. Moreover, we compare to and confirm previous determinantal quantum Monte Carlo results on incommensurate spin-density waves at finite doping and temperature.
△ Less
Submitted 16 July, 2021; v1 submitted 22 September, 2020;
originally announced September 2020.
-
The ITensor Software Library for Tensor Network Calculations
Authors:
Matthew Fishman,
Steven R. White,
E. Miles Stoudenmire
Abstract:
ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagram notation, which allows users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of tensor network algorithms. After discussing the philosophy behind the ITensor…
▽ More
ITensor is a system for programming tensor network calculations with an interface modeled on tensor diagram notation, which allows users to focus on the connectivity of a tensor network without manually bookkeeping tensor indices. The ITensor interface rules out common programming errors and enables rapid prototyping of tensor network algorithms. After discussing the philosophy behind the ITensor approach, we show examples of each part of the interface including Index objects, the ITensor product operator, tensor factorizations, tensor storage types, algorithms for matrix product state (MPS) and matrix product operator (MPO) tensor networks, quantum number conserving block-sparse tensors, and the NDTensors library. We also review publications that have used ITensor for quantum many-body physics and for other areas where tensor networks are increasingly applied. To conclude we discuss promising features and optimizations to be added in the future.
△ Less
Submitted 20 December, 2021; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Bypassing the energy functional in density functional theory: Direct calculation of electronic energies from conditional probability densities
Authors:
Ryan J. McCarty,
Dennis Perchak,
Ryan Pederson,
Robert Evans,
Yiheng Qiu,
Steven R. White,
Kieron Burke
Abstract:
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 e…
▽ More
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.
△ Less
Submitted 21 October, 2020; v1 submitted 3 July, 2020;
originally announced July 2020.
-
Time Dependent Variational Principle with Ancillary Krylov Subspace
Authors:
Mingru Yang,
Steven R. White
Abstract:
We propose an improved scheme to do the time dependent variational principle (TDVP) in finite matrix product states (MPS) for two-dimensional systems or one-dimensional systems with long range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state-averaging with global Krylov vectors. We show that the projection error is significantly reduc…
▽ More
We propose an improved scheme to do the time dependent variational principle (TDVP) in finite matrix product states (MPS) for two-dimensional systems or one-dimensional systems with long range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state-averaging with global Krylov vectors. We show that the projection error is significantly reduced so that precise time evolution can still be obtained even if a larger time step is used. Combined with the one-site TDVP, our approach provides a way to dynamically increase the bond dimension while still preserving unitarity for real time evolution. Our method can be more accurate and exhibit slower bond dimension growth than the conventional two-site TDVP.
△ Less
Submitted 29 September, 2020; v1 submitted 12 May, 2020;
originally announced May 2020.
-
Plaquette versus ordinary $d$-wave pairing in the $t'$-Hubbard model on a width 4 cylinder
Authors:
Chia-Min Chung,
Mingpu Qin,
Shiwei Zhang,
Ulrich Schollwöck,
Steven R. White
Abstract:
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 take…
▽ More
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.
△ Less
Submitted 13 April, 2020; v1 submitted 6 April, 2020;
originally announced April 2020.
-
Critical properties of a comb lattice
Authors:
Natalia Chepiga,
Steven R. White
Abstract:
In this paper we study the critical properties of the Heisenberg spin-1/2 model on a comb lattice -- a 1D backbone decorated with finite 1D chains -- the teeth. We address the problem numerically by a comb tensor network that duplicates the geometry of a lattice. We observe a fundamental difference between the states on a comb with even and odd number of sites per tooth, which resembles an even-od…
▽ More
In this paper we study the critical properties of the Heisenberg spin-1/2 model on a comb lattice -- a 1D backbone decorated with finite 1D chains -- the teeth. We address the problem numerically by a comb tensor network that duplicates the geometry of a lattice. We observe a fundamental difference between the states on a comb with even and odd number of sites per tooth, which resembles an even-odd effect in spin-1/2 ladders. The comb with odd teeth is always critical, not only along the teeth, but also along the backbone, which leads to a competition between two critical regimes in orthogonal directions. In addition, we show that in a weak-backbone limit the excitation energy scales as $1/(NL)$, and not as $1/N$ or $1/L$ typical for 1D systems. For even teeth in the weak backbone limit the system corresponds to a collection of decoupled critical chains of length $L$, while in the strong backbone limit, one spin from each tooth forms the backbone, so the effective length of a critical tooth is one site shorter, $L-1$. Surprisingly, these two regimes are connected via a state where a critical chain spans over two nearest neighbor teeth, with an effective length $2L$.
△ Less
Submitted 8 June, 2020; v1 submitted 26 February, 2020;
originally announced February 2020.
-
The Alzheimer's Disease Prediction Of Longitudinal Evolution (TADPOLE) Challenge: Results after 1 Year Follow-up
Authors:
Razvan V. Marinescu,
Neil P. Oxtoby,
Alexandra L. Young,
Esther E. Bron,
Arthur W. Toga,
Michael W. Weiner,
Frederik Barkhof,
Nick C. Fox,
Arman Eshaghi,
Tina Toni,
Marcin Salaterski,
Veronika Lunina,
Manon Ansart,
Stanley Durrleman,
Pascal Lu,
Samuel Iddi,
Dan Li,
Wesley K. Thompson,
Michael C. Donohue,
Aviv Nahon,
Yarden Levy,
Dan Halbersberg,
Mariya Cohen,
Huiling Liao,
Tengfei Li
, et al. (71 additional authors not shown)
Abstract:
We present the findings of "The Alzheimer's Disease Prediction Of Longitudinal Evolution" (TADPOLE) Challenge, which compared the performance of 92 algorithms from 33 international teams at predicting the future trajectory of 219 individuals at risk of Alzheimer's disease. Challenge participants were required to make a prediction, for each month of a 5-year future time period, of three key outcome…
▽ More
We present the findings of "The Alzheimer's Disease Prediction Of Longitudinal Evolution" (TADPOLE) Challenge, which compared the performance of 92 algorithms from 33 international teams at predicting the future trajectory of 219 individuals at risk of Alzheimer's disease. Challenge participants were required to make a prediction, for each month of a 5-year future time period, of three key outcomes: clinical diagnosis, Alzheimer's Disease Assessment Scale Cognitive Subdomain (ADAS-Cog13), and total volume of the ventricles. The methods used by challenge participants included multivariate linear regression, machine learning methods such as support vector machines and deep neural networks, as well as disease progression models. No single submission was best at predicting all three outcomes. For clinical diagnosis and ventricle volume prediction, the best algorithms strongly outperform simple baselines in predictive ability. However, for ADAS-Cog13 no single submitted prediction method was significantly better than random guesswork. Two ensemble methods based on taking the mean and median over all predictions, obtained top scores on almost all tasks. Better than average performance at diagnosis prediction was generally associated with the additional inclusion of features from cerebrospinal fluid (CSF) samples and diffusion tensor imaging (DTI). On the other hand, better performance at ventricle volume prediction was associated with inclusion of summary statistics, such as the slope or maxima/minima of biomarkers. TADPOLE's unique results suggest that current prediction algorithms provide sufficient accuracy to exploit biomarkers related to clinical diagnosis and ventricle volume, for cohort refinement in clinical trials for Alzheimer's disease. However, results call into question the usage of cognitive test scores for patient selection and as a primary endpoint in clinical trials.
△ Less
Submitted 27 December, 2021; v1 submitted 9 February, 2020;
originally announced February 2020.
-
Ground-state properties of the hydrogen chain: insulator-to-metal transition, dimerization, and magnetic phases
Authors:
Mario Motta,
Claudio Genovese,
Fengjie Ma,
Zhi-Hao Cui,
Randy Sawaya,
Garnet Kin-Lic Chan,
Natalia Chepiga,
Phillip Helms,
Carlos Jimenez-Hoyos,
Andrew J. Millis,
Ushnish Ray,
Enrico Ronca,
Hao Shi,
Sandro Sorella,
Edwin M. Stoudenmire,
Steven R. White,
Shiwei Zhang
Abstract:
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 bul…
▽ More
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.
△ Less
Submitted 13 July, 2020; v1 submitted 4 November, 2019;
originally announced November 2019.
-
Absence of superconductivity in the pure two-dimensional Hubbard model
Authors:
Mingpu Qin,
Chia-Min Chung,
Hao Shi,
Ettore Vitali,
Claudius Hubig,
Ulrich Schollwöck,
Steven R. White,
Shiwei Zhang
Abstract:
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 superconductor…
▽ More
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.
△ Less
Submitted 22 July, 2020; v1 submitted 20 October, 2019;
originally announced October 2019.
-
Discontinuous Galerkin discretization for quantum simulation of chemistry
Authors:
Jarrod R. McClean,
Fabian M. Faulstich,
Qinyi Zhu,
Bryan O'Gorman,
Yiheng Qiu,
Steven R. White,
Ryan Babbush,
Lin Lin
Abstract:
Methods for electronic structure based on Gaussian and molecular orbital discretizations offer a well established, compact representation that forms much of the foundation of correlated quantum chemistry calculations on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growt…
▽ More
Methods for electronic structure based on Gaussian and molecular orbital discretizations offer a well established, compact representation that forms much of the foundation of correlated quantum chemistry calculations on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from $O(N^{4.5})$ in molecular orbital representations to $O(N^{2.6})$ in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1-2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.
△ Less
Submitted 30 August, 2019;
originally announced September 2019.
-
Comb tensor networks
Authors:
Natalia Chepiga,
Steven R. White
Abstract:
In this paper we propose a special type of a tree tensor network that has the geometry of a comb---a 1D backbone with finite 1D teeth projecting out from it. This tensor network is designed to provide an effective description of higher dimensional objects with special limited interactions, or, alternatively, one-dimensional systems composed of complicated zero-dimensional objects. We provide detai…
▽ More
In this paper we propose a special type of a tree tensor network that has the geometry of a comb---a 1D backbone with finite 1D teeth projecting out from it. This tensor network is designed to provide an effective description of higher dimensional objects with special limited interactions, or, alternatively, one-dimensional systems composed of complicated zero-dimensional objects. We provide details on the best numerical procedures for the proposed network, including an algorithm for variational optimization of the wave-function as a comb tensor network, and the transformation of the comb into a matrix product state. We compare the complexity of using a comb versus alternative matrix product state representations using density matrix renormalization group (DMRG) algorithms. As an application, we study a spin-1 Heisenberg model system which has a comb geometry. In the case where the ends of the teeth are terminated by spin-1/2 spins, we find that Haldane edge states of the teeth along the backbone form a critical spin-1/2 chain, whose properties can be tuned by the coupling constant along the backbone. By adding next-nearest-neighbor interactions along the backbone, the comb can be brought into a gapped phase with a long-range dimerization along the backbone. The critical and dimerized phases are separated by a Kosterlitz-Thouless phase transition, the presence of which we confirm numerically. Finally, we show that when the teeth contain an odd number of spins and are not terminated by spin-1/2's, a special type of comb edge states emerge.
△ Less
Submitted 1 March, 2019;
originally announced March 2019.
-
Anisotropic-Exchange Magnets on a Triangular Lattice: Spin Waves, Accidental Degeneracies, and Dual Spin Liquids
Authors:
P. A. Maksimov,
Zhenyue Zhu,
Steven R. White,
A. L. Chernyshev
Abstract:
We present an extensive overview of the phase diagram, spin-wave excitations, and finite-temperature transitions of the anisotropic-exchange magnets on an ideal nearest-neighbor triangular lattice. We investigate transitions between five principal classical phases of the corresponding model: ferromagnetic, Néel, its dual, and the two stripe phases. Transitions are identified by the spin-wave insta…
▽ More
We present an extensive overview of the phase diagram, spin-wave excitations, and finite-temperature transitions of the anisotropic-exchange magnets on an ideal nearest-neighbor triangular lattice. We investigate transitions between five principal classical phases of the corresponding model: ferromagnetic, Néel, its dual, and the two stripe phases. Transitions are identified by the spin-wave instabilities and by the Luttinger-Tisza approach. Some of the transitions are direct and others occur via intermediate phases with more complicated forms of ordering. In a portion of the Néel phase, we find spin-wave instabilities to a long-range spiral-like state. In the stripe phases, quantum fluctuations are mostly negligible, leaving the ordered moment nearly saturated even for the $S=1/2$ case. However, for a two-dimensional surface of the full 3D parameter space, the spin-wave spectrum in one of the stripe phases exhibits an enigmatic accidental degeneracy manifested by pseudo-Goldstone modes. As a result, despite the nearly classical ground state, the ordering transition temperature in a wide region of the phase diagram is significantly suppressed from the mean-field expectation. We identify this accidental degeneracy as due to an exact correspondence to an extended Kitaev-Heisenberg model with emergent symmetries that naturally lead to the pseudo-Goldstone modes. There are previously studied dualities within the Kitaev-Heisenberg model on the triangular lattice that are exposed here in a wider parameter space. One important implication of this correspondence for the $S=1/2$ case is the existence of a region of the spin-liquid phase that is dual to the spin-liquid phase discovered recently by us. We complement our studies by the density-matrix renormalization group of the $S=1/2$ model to confirm some of the duality relations and to verify the existence of the dual spin-liquid phase.
△ Less
Submitted 14 December, 2021; v1 submitted 14 November, 2018;
originally announced November 2018.
-
Multi-sliced Gausslet Basis Sets for Electronic Structure
Authors:
Steven R. White,
E. Miles Stoudenmire
Abstract:
We introduce highly local basis sets for electronic structure which are very efficient for correlation calculations near the complete basis set limit. Our approach is based on gausslets, recently introduced wavelet-like smooth orthogonal functions. We adapt the gausslets to particular systems using one dimensional coordinate transformations, putting more basis functions near nuclei, while maintain…
▽ More
We introduce highly local basis sets for electronic structure which are very efficient for correlation calculations near the complete basis set limit. Our approach is based on gausslets, recently introduced wavelet-like smooth orthogonal functions. We adapt the gausslets to particular systems using one dimensional coordinate transformations, putting more basis functions near nuclei, while maintaining orthogonality. Three dimensional basis functions are composed out of products of the 1D functions in an efficient way called multislicing. We demonstrate the new bases with both Hartree Fock and density matrix renormalization group (DMRG) calculations on hydrogen chain systems. With both methods, we can go to higher accuracy in the complete basis set limit than is practical for conventional Gaussian basis sets, with errors near 0.1 mH per atom.
△ Less
Submitted 26 September, 2018;
originally announced September 2018.
-
Density Matrix Renormalization Group Study of a One Dimensional Diatomic Molecule beyond the Born-Oppenheimer Approximation
Authors:
Mingru Yang,
Steven R. White
Abstract:
We study one dimensional models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate DMRG efficiently with this system, a three-site algor…
▽ More
We study one dimensional models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate DMRG efficiently with this system, a three-site algorithm has been implemented. We also use a compression method to treat the long-range interactions between charged particles. We find that 1D diatomic molecules with spin-1/2 nuclei in the spin-triplet state will unbind when the mass of the nuclei reduces to only a few times larger than the electron mass, while the molecule with nuclei in the singlet state always binds, given the two electrons in their singlet state in both cases.
△ Less
Submitted 30 September, 2020; v1 submitted 16 September, 2018;
originally announced September 2018.
-
Topography of Spin Liquids on a Triangular Lattice
Authors:
Zhenyue Zhu,
P. A. Maksimov,
Steven R. White,
A. L. Chernyshev
Abstract:
Spin systems with frustrated anisotropic interactions are of significant interest due to possible exotic ground states. We have explored their phase diagram on a nearest-neighbor triangular lattice using the density-matrix renormalization group and mapped out the topography of the region that can harbor a spin liquid. We find that this spin-liquid phase is continuously connected to a previously di…
▽ More
Spin systems with frustrated anisotropic interactions are of significant interest due to possible exotic ground states. We have explored their phase diagram on a nearest-neighbor triangular lattice using the density-matrix renormalization group and mapped out the topography of the region that can harbor a spin liquid. We find that this spin-liquid phase is continuously connected to a previously discovered spin-liquid phase of the isotropic $J_1\!-\!J_2$ model. The two limits show nearly identical spin correlations, making the case that their respective spin liquids are isomorphic to each other.
△ Less
Submitted 9 May, 2018; v1 submitted 3 January, 2018;
originally announced January 2018.
-
The t-t'-J model in one dimension using extremely correlated Fermi liquid theory and time dependent density matrix renormalization group
Authors:
Peizhi Mai,
Steven R. White,
B. Sriram Shastry
Abstract:
We study the one dimensional t-t'-J model for generic couplings using two complementary theories, the extremely correlated Fermi liquid theory and time-dependent density matrix renormalization group over a broad energy scale. The two methods provide a unique insight into the strong momentum dependence of the self-energy of this prototypical non-Fermi liquid, described at low energies as a Tomonaga…
▽ More
We study the one dimensional t-t'-J model for generic couplings using two complementary theories, the extremely correlated Fermi liquid theory and time-dependent density matrix renormalization group over a broad energy scale. The two methods provide a unique insight into the strong momentum dependence of the self-energy of this prototypical non-Fermi liquid, described at low energies as a Tomonaga-Luttinger liquid. We also demonstrate its intimate relationship to spin-charge separation, i.e. the splitting of Landau quasiparticles of higher dimensions into two constituents, driven by strong quantum fluctuations inherent in one dimension. The momentum distribution function, the spectral function, and the excitation dispersion of these two methods also compare well.
△ Less
Submitted 23 June, 2018; v1 submitted 14 December, 2017;
originally announced December 2017.
-
Hybrid grid/basis set discretizations of the Schrödinger equation
Authors:
Steven R. White
Abstract:
We present a new kind of basis function for discretizing the Schrödinger equation in electronic structure calculations, called a gausslet, which has wavelet-like features but is composed of a sum of Gaussians. Gausslets are placed on a grid and combine advantages of both grid and basis set approaches. They are orthogonal, infinitely smooth, symmetric, polynomially complete, and with a high degree…
▽ More
We present a new kind of basis function for discretizing the Schrödinger equation in electronic structure calculations, called a gausslet, which has wavelet-like features but is composed of a sum of Gaussians. Gausslets are placed on a grid and combine advantages of both grid and basis set approaches. They are orthogonal, infinitely smooth, symmetric, polynomially complete, and with a high degree of locality. Because they are formed from Gaussians, they are easily combined with traditional atom-centered Gaussian bases. We also introduce diagonal approximations which dramatically reduce the computational scaling of two-electron Coulomb terms in the Hamiltonian.
△ Less
Submitted 23 September, 2017;
originally announced September 2017.
-
Chemical accuracy from small, system-adapted basis functions
Authors:
Thomas E. Baker,
Kieron Burke,
Steven R. White
Abstract:
We propose a general method for constructing system-dependent basis functions for correlated quantum chemical calculations. Our construction combines features from several traditional approaches: plane waves, localized basis functions, and wavelets. In a one-dimensional mimic of Coulomb systems, it requires only 2-3 basis functions per electron to achieve chemical accuracy, and reproduces the natu…
▽ More
We propose a general method for constructing system-dependent basis functions for correlated quantum chemical calculations. Our construction combines features from several traditional approaches: plane waves, localized basis functions, and wavelets. In a one-dimensional mimic of Coulomb systems, it requires only 2-3 basis functions per electron to achieve chemical accuracy, and reproduces the natural orbitals. We illustrate its effectiveness for molecular energy curves and chains of many atoms. We discuss the promise and challenges for realistic quantum chemical calculations.
△ Less
Submitted 21 February, 2018; v1 submitted 11 September, 2017;
originally announced September 2017.
-
Matrix product state techniques for two-dimensional systems at finite temperature
Authors:
Benedikt Bruognolo,
Zhenyue Zhu,
Steven R. White,
E. Miles Stoudenmire
Abstract:
The density matrix renormalization group is one of the most powerful numerical methods for computing ground-state properties of two-dimensional (2D) quantum lattice systems. Here we show its finite-temperature extensions are also viable for 2D, using the following strategy: At high temperatures, we combine density-matrix purification and numerical linked-cluster expansions to extract static observ…
▽ More
The density matrix renormalization group is one of the most powerful numerical methods for computing ground-state properties of two-dimensional (2D) quantum lattice systems. Here we show its finite-temperature extensions are also viable for 2D, using the following strategy: At high temperatures, we combine density-matrix purification and numerical linked-cluster expansions to extract static observables directly in the thermodynamic limit. At low temperatures inaccessible to purification, we use the minimally entangled typical thermal state (METTS) algorithm on cylinders. We consider the triangular Heisenberg antiferromagnet as a first application, finding excellent agreement with other state of the art methods. In addition, we present a METTS-based approach that successfully extracts critical temperatures, and apply it to a frustrated lattice model. On a technical level, we compare two different schemes for performing imaginary-time evolution of 2D clusters, finding that a Suzuki-Trotter decomposition with swap gates is currently the most accurate and efficient.
△ Less
Submitted 22 August, 2017; v1 submitted 16 May, 2017;
originally announced May 2017.
-
Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods
Authors:
Mario Motta,
David M. Ceperley,
Garnet Kin-Lic Chan,
John A. Gomez,
Emanuel Gull,
Sheng Guo,
Carlos Jimenez-Hoyos,
Tran Nguyen Lan,
Jia Li,
Fengjie Ma,
Andrew J. Millis,
Nikolay V. Prokof'ev,
Ushnish Ray,
Gustavo E. Scuseria,
Sandro Sorella,
Edwin M. Stoudenmire,
Qiming Sun,
Igor S. Tupitsyn,
Steven R. White,
Dominika Zgid,
Shiwei Zhang
Abstract:
We present numerical results for the equation of state of an infinite chain of hydrogen atoms. A variety of modern many-body methods are employed, with exhaustive cross-checks and validation. Approaches for reaching the continuous space limit and the thermodynamic limit are investigated, proposed, and tested. The detailed comparisons provide a benchmark for assessing the current state of the art i…
▽ More
We present numerical results for the equation of state of an infinite chain of hydrogen atoms. A variety of modern many-body methods are employed, with exhaustive cross-checks and validation. Approaches for reaching the continuous space limit and the thermodynamic limit are investigated, proposed, and tested. The detailed comparisons provide a benchmark for assessing the current state of the art in many-body computation, and for the development of new methods. The ground-state energy per atom in the linear chain is accurately determined versus bondlength, with a confidence bound given on all uncertainties.
△ Less
Submitted 6 November, 2017; v1 submitted 1 May, 2017;
originally announced May 2017.
-
Infinite Sparse Structured Factor Analysis
Authors:
Matthew C. Pearce,
Simon R. White
Abstract:
Matrix factorisation methods decompose multivariate observations as linear combinations of latent feature vectors. The Indian Buffet Process (IBP) provides a way to model the number of latent features required for a good approximation in terms of regularised reconstruction error. Previous work has focussed on latent feature vectors with independent entries. We extend the model to include nondiagon…
▽ More
Matrix factorisation methods decompose multivariate observations as linear combinations of latent feature vectors. The Indian Buffet Process (IBP) provides a way to model the number of latent features required for a good approximation in terms of regularised reconstruction error. Previous work has focussed on latent feature vectors with independent entries. We extend the model to include nondiagonal latent covariance structures representing characteristics such as smoothness. This is done by . Using simulations we demonstrate that under appropriate conditions a smoothness prior helps to recover the true latent features, while denoising more accurately. We demonstrate our method on a real neuroimaging dataset, where computational tractability is a sufficient challenge that the efficient strategy presented here is essential.
△ Less
Submitted 13 April, 2017;
originally announced April 2017.
-
Disorder-Induced Mimicry of a Spin Liquid in YbMgGaO$_4$
Authors:
Zhenyue Zhu,
P. A. Maksimov,
Steven R. White,
A. L. Chernyshev
Abstract:
We suggest that a randomization of the pseudo-dipolar interaction in the spin-orbit-generated low-energy Hamiltonian of YbMgGaO$_4$ due to an inhomogeneous charge environment from a natural mixing of Mg$^{2+}$ and Ga$^{3+}$ can give rise to orientational spin disorder and mimic a spin-liquid-like state. In the absence of such quenched disorder, $1/S$ and density matrix renormalization group calcul…
▽ More
We suggest that a randomization of the pseudo-dipolar interaction in the spin-orbit-generated low-energy Hamiltonian of YbMgGaO$_4$ due to an inhomogeneous charge environment from a natural mixing of Mg$^{2+}$ and Ga$^{3+}$ can give rise to orientational spin disorder and mimic a spin-liquid-like state. In the absence of such quenched disorder, $1/S$ and density matrix renormalization group calculations both show robust ordered states for the physically relevant phases of the model. Our scenario is consistent with the available experimental data and further experiments are proposed to support it.
△ Less
Submitted 14 September, 2017; v1 submitted 8 March, 2017;
originally announced March 2017.
-
Sliced Basis Density Matrix Renormalization Group for Electronic Structure
Authors:
E. Miles Stoudenmire,
Steven R. White
Abstract:
We introduce a hybrid approach to applying the density matrix renormalization group (DMRG) to continuous systems, combining a grid approximation along one direction with a finite Gaussian basis set along the remaining two directions. This approach is especially useful for chain-like molecules, where the grid is used in the long direction, and we demonstrate the approach with results for hydrogen c…
▽ More
We introduce a hybrid approach to applying the density matrix renormalization group (DMRG) to continuous systems, combining a grid approximation along one direction with a finite Gaussian basis set along the remaining two directions. This approach is especially useful for chain-like molecules, where the grid is used in the long direction, and we demonstrate the approach with results for hydrogen chains. The computational time for this system scales approximately linearly with the length of the chain, as we demonstrate with minimal basis set calculations with up to 1000 atoms, which are near-exact within the basis. The linear scaling comes from the combination of localization of the basis and a compression method with controlled accuracy for the long-ranged Coulomb terms in the Hamiltonian.
△ Less
Submitted 13 February, 2017;
originally announced February 2017.
-
Hybrid-space density matrix renormalization group study of the doped two-dimensional Hubbard model
Authors:
G. Ehlers,
S. R. White,
R. M. Noack
Abstract:
The performance of the density matrix renormalization group (DMRG) is strongly influenced by the choice of the local basis of the underlying physical lattice. We demonstrate that, for the two-dimensional Hubbard model, the hybrid real-momentum space formulation of the DMRG is computationally more efficient than the standard real-space formulation. In particular, we show that the computational cost…
▽ More
The performance of the density matrix renormalization group (DMRG) is strongly influenced by the choice of the local basis of the underlying physical lattice. We demonstrate that, for the two-dimensional Hubbard model, the hybrid real-momentum space formulation of the DMRG is computationally more efficient than the standard real-space formulation. In particular, we show that the computational cost for fixed bond dimension of the hybrid-space DMRG is approximately independent of the width of the lattice, in contrast to the real-space DMRG, for which it is proportional to the width squared. We apply the hybrid-space algorithm to calculate the ground state of the doped two-dimensional Hubbard model on cylinders of width four and six sites; at $n=0.875$ filling, the ground state exhibits a striped charge-density distribution with a wavelength of eight sites for both $U/t=4.0$ and $U/t=8.0$. We find that the strength of the charge ordering depends on $U/t$ and on the boundary conditions.Furthermore, we investigate the magnetic ordering as well as the decay of the static spin, charge, and pair-field correlation functions.
△ Less
Submitted 21 March, 2017; v1 submitted 13 January, 2017;
originally announced January 2017.
-
Stripe order in the underdoped region of the two-dimensional Hubbard model
Authors:
Bo-Xiao Zheng,
Chia-Min Chung,
Philippe Corboz,
Georg Ehlers,
Ming-Pu Qin,
Reinhard M. Noack,
Hao Shi,
Steven R. White,
Shiwei Zhang,
Garnet Kin-Lic Chan
Abstract:
Competing inhomogeneous orders are a central feature of correlated electron materials including the high-temperature superconductors. The two- dimensional Hubbard model serves as the canonical microscopic physical model for such systems. Multiple orders have been proposed in the underdoped part of the phase diagram, which corresponds to a regime of maximum numerical difficulty. By combining the la…
▽ More
Competing inhomogeneous orders are a central feature of correlated electron materials including the high-temperature superconductors. The two- dimensional Hubbard model serves as the canonical microscopic physical model for such systems. Multiple orders have been proposed in the underdoped part of the phase diagram, which corresponds to a regime of maximum numerical difficulty. By combining the latest numerical methods in exhaustive simulations, we uncover the ordering in the underdoped ground state. We find a stripe order that has a highly compressible wavelength on an energy scale of a few Kelvin, with wavelength fluctuations coupled to pairing order. The favored filled stripe order is different from that seen in real materials. Our results demonstrate the power of modern numerical methods to solve microscopic models even in challenging settings.
△ Less
Submitted 7 December, 2017; v1 submitted 30 December, 2016;
originally announced January 2017.
-
Pure density functional for strong correlations and the thermodynamic limit from machine learning
Authors:
Li Li,
Thomas E. Baker,
Steven R. White,
Kieron Burke
Abstract:
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-S…
▽ More
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.
△ Less
Submitted 13 September, 2016;
originally announced September 2016.
-
Electronic quasiparticles in the quantum dimer model: density matrix renormalization group results
Authors:
Junhyun Lee,
Subir Sachdev,
Steven R. White
Abstract:
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 dens…
▽ More
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.
△ Less
Submitted 8 September, 2016; v1 submitted 13 June, 2016;
originally announced June 2016.
-
Representation and design of wavelets using unitary circuits
Authors:
Glen Evenbly,
Steven R. White
Abstract:
The representation of discrete, compact wavelet transformations (WTs) as circuits of local unitary gates is discussed. We employ a similar formalism as used in the multi-scale representation of quantum many-body wavefunctions using unitary circuits, further cementing the relation established in [Phys. Rev. Lett. 116, 140403 (2016)] between classical and quantum multi-scale methods. An algorithm fo…
▽ More
The representation of discrete, compact wavelet transformations (WTs) as circuits of local unitary gates is discussed. We employ a similar formalism as used in the multi-scale representation of quantum many-body wavefunctions using unitary circuits, further cementing the relation established in [Phys. Rev. Lett. 116, 140403 (2016)] between classical and quantum multi-scale methods. An algorithm for constructing the circuit representation of known orthogonal, dyadic, discrete WTs is presented, and the explicit representation for Daubechies wavelets, coiflets, and symlets is provided. Furthermore, we demonstrate the usefulness of the circuit formalism in designing novel WTs, including various classes of symmetric wavelets and multi-wavelets, boundary wavelets and biorthogonal wavelets.
△ Less
Submitted 24 May, 2016;
originally announced May 2016.
-
Matrix Product Operators, Matrix Product States, and ab initio Density Matrix Renormalization Group algorithms
Authors:
Garnet Kin-Lic Chan,
Anna Keselman,
Naoki Nakatani,
Zhendong Li,
Steven R. White
Abstract:
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 transl…
▽ More
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.
△ Less
Submitted 18 June, 2016; v1 submitted 9 May, 2016;
originally announced May 2016.
-
Entanglement renormalization and wavelets
Authors:
Glen Evenbly,
Steven R. White
Abstract:
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 corres…
▽ More
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.
△ Less
Submitted 16 February, 2016; v1 submitted 2 February, 2016;
originally announced February 2016.
-
Space group symmetry fractionalization in a chiral kagome Heisenberg antiferromagnet
Authors:
Michael P. Zaletel,
Zhenyue Zhu,
Yuan-Ming Lu,
Ashvin Vishwanath,
Steven R. White
Abstract:
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 ant…
▽ More
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.
△ Less
Submitted 6 November, 2015; v1 submitted 4 November, 2015;
originally announced November 2015.
-
Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms
Authors:
J. P. F. LeBlanc,
Andrey E. Antipov,
Federico Becca,
Ireneusz W. Bulik,
Garnet Kin-Lic Chan,
Chia-Min Chung,
Youjin Deng,
Michel Ferrero,
Thomas M. Henderson,
Carlos A. Jiménez-Hoyos,
E. Kozik,
Xuan-Wen Liu,
Andrew J. Millis,
N. V. Prokof'ev,
Mingpu Qin,
Gustavo E. Scuseria,
Hao Shi,
B. V. Svistunov,
Luca F. Tocchio,
I. S. Tupitsyn,
Steven R. White,
Shiwei Zhang,
Bo-Xiao Zheng,
Zhenyue Zhu,
Emanuel Gull
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 December, 2015; v1 submitted 9 May, 2015;
originally announced May 2015.
-
Compression of Correlation Matrices and an Efficient Method for Forming Matrix Product States of Fermionic Gaussian States
Authors:
Matthew T. Fishman,
Steven R. White
Abstract:
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 co…
▽ More
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).
△ Less
Submitted 28 April, 2015;
originally announced April 2015.