Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/9705069 · Journal of the Physical Society of Japan

A Quantum Monte Carlo Method and Its Applications to Multi-Orbital Hubbard Models

Yukitoshi Motome, Masatoshi Imada

Published Jul 15, 1997• 38 citations
arXivHTMLPDF

Abstract

We present a framework of an auxiliary field quantum Monte Carlo (QMC) method for multi-orbital Hubbard models. Our formulation can be applied to a Hamiltonian which includes terms for on-site Coulomb interaction for both intra- and inter-orbitals, intra-site exchange interaction and energy differences between orbitals. Based on our framework, we point out possible ways to investigate various phase transitions such as metal-insulator, magnetic and orbital order-disorder transitions without the minus sign problem. As an application, a two-band model is investigated by the projection QMC method and the ground state properties of this model are presented.

DOI
10.1143/jpsj.66.1872

Related papers

Auxiliary field quantum Monte Carlo for multiband Hubbard models: Controlling the sign and phase problems to capture Hund's physics

Hongxia Hao, Brenda M. Rubenstein, Hao Shi

In the study of strongly correlated, many-electron systems, the Hubbard-Kanamori (HK) model has emerged as one of the prototypes for transition metal oxide physics. The model is multiband in nature and contains Hund's coupling terms, which have pronounced effects on metal-insulator transitions, high-temperature superconductivity, and other physical properties. In the following, we present a complete theoretical framework for treating the HK model using the ground-state auxiliary field quantum Monte Carlo (AFQMC) method and analyze its performance on few-band models whose parameters approximate those observed in ruthenates, rhodates, and other materials exhibiting Hund's physics. Unlike previous studies, the constrained path and phaseless approximations are used to respectively control the sign and phase problems, which enables high-accuracy modeling of the HK model's ground-state properties within parameter regimes of experimental interest. We demonstrate that, after careful consideration of the Hubbard-Stratonovich transformations and trial wave functions employed, relative errors in the energy of less than $1%$ can routinely be achieved for moderate to large values of the Hund's coupling constant. Crucially, our methodology also accurately predicts magnetic ordering and phase transitions. The results presented open the door to more predictive modeling of Hund's physics within a wide range of strongly correlated materials using AFQMC.

arXiv:1902.01463•Physical review. B./Physical review. B•Jun 20, 2019•20 citations

Numerical Study for the Ground State of Multi-Orbital Hubbard Models

Yukitoshi Motome, Masatoshi Imada

Ground state properties of multi-orbital Hubbard models are investigated by the auxiliary field quantum Monte Carlo method. A Monte Carlo technique generalized to the multi-orbital systems is introduced and examined in detail. The algorithm contains non-trivial cases where the negative sign problem does not exist. We investigate one-dimensional systems with doubly degenerate orbitals by this new technique. Properties of the Mott insulating state are quantitatively clarified as the strongly correlated insulator, where the charge gap amplitude is much larger than the spin gap. The insulator-metal transitions driven by the chemical potential shows a universality class with the correlation length exponent $ν=1/2$, which is consistent with the scaling arguments. Increasing level split between two orbitals drives crossover from the Mott insulator with high spin state to the band insulator with low spin state, where the spin gap amplitude increases and becomes closer to the charge gap. Experimental relevance of our results especially to Haldane materials is discussed.

arXiv:cond-mat/9806279•Journal of the Physical Society of Japan•Sep 15, 1998•10 citations

Computational Studies Of Strongly Correlated Quantum Matter

Hao Shi

The study of strongly correlated quantum many-body systems is an outstanding challenge. Highly accurate results are needed for the understanding of practical and fundamental problems in condensed-matter physics, high energy physics, material science, quantum chemistry and so on. Our familiar mean-field or perturbative methods tend to be ineffective. Numerical simulations provide a promising approach for studying such systems. The fundamental difficulty of numerical simulation is that the dimension of the Hilbert space needed to describe interacting systems increases exponentially with the system size. Quantum Monte Carlo (QMC) methods are one of the best approaches to tackle the problem of enormous Hilbert space. They have been highly successful for boson systems and unfrustrated spin models. For systems with fermions, the exchange symmetry in general causes the infamous sign problem, making the statistical noise in the computed results grow exponentially with the system size. This hinders our understanding of interesting physics such as high-temperature superconductivity, metal-insulator phase transition. In this thesis, we present a variety of new developments in the auxiliary-field quantum Monte Carlo (AFQMC) methods, including the incorporation of symmetry in both the trial wave function and the projector, developing the constraint release method, using the force-bias to drastically improve the efficiency in Metropolis framework, identifying and solving the infinite variance problem, and sampling Hartree-Fock-Bogoliubov wave function. With these developments, some of the most challenging many-electron problems are now under control. We obtain an exact numerical solution of two-dimensional strongly interacting Fermi atomic gas, determine the ground state properties of the 2D Fermi gas with Rashba spin-orbit coupling, provide benchmark results for the ground state of the two-dimensional Hubbard model, and establish that the Hubbard model has a stripe order in the underdoped region.

•Jan 1, 2017•0 citations

Auxiliary-Field Quantum Monte Carlo Calculations of Excited States and Strongly Correlated Systems

S. Zhang

to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.

•Jan 1, 2012•0 citations

Continuous Time Quantum Monte Carlo in Combination with Machine Learning on the Hubbard Model

Hunpyo Lee

The acceleration of exact continuous time quantum Monte Carlo (CTQMC) approaches in multi-site or multi-orbital systems is extremely interesting work, because these approaches are very time-consuming in terms of numerical computation and might account for the nature of exotic behaviors such as high-temperature superconductivity and Mott insulator behavior observed in the strongly correlated materials. We extend the recently developed interaction-expansion CTQMC method in combination with a machine learning (CTQMC+ML) approach for the single-site and single-orbital systems to multi-site and multi-orbital ones. This method can be applied to explore the nonlocal correlation effects in lattice models and to study the electronic structure of real materials via an ab-initio density functional theory plus dynamical mean field theory approach. We find that our CTQMC+ML method for multi-site (and multi-orbital) systems accurately predicts the impurity Green’s function with less computational time than the CTQMC approaches, as in the case of the single-site and single-orbital version.

Journal of the Korean Physical Society•Nov 1, 2019•1 citation

Efficient local energy evaluation for multi-Slater wave functions in orbital space quantum Monte Carlo

Ankit Mahajan, Sandeep Sharma

We present an algorithm for calculating the local energy of a multi-Slater wave function in orbital space quantum Monte Carlo (QMC). Recent developments in selected configuration interaction methods have led to increased interest in using multi-Slater trial wave functions in various QMC methods. For an ab initio Hamiltonian, our algorithm has a cost scaling of O(n5 + nc), as opposed to the O(n4nc) scaling of existing orbital space algorithms, where n is the system size and nc is the number of configurations in the wave function. We present our method using variational Monte Carlo calculations with the Jastrow multi-Slater wave function, although the formalism should be applicable for auxiliary field QMC. We apply it to polyacetylene and demonstrate the possibility of using a much larger number of configurations than possible using existing methods.

arXiv:2008.06477•The Journal of Chemical Physics•Nov 18, 2020•22 citations
We present a framework of an auxiliary field quantum Monte Carlo (QMC) method for multi-orbital Hubbard models. Our formulation can be applied to a Hamiltonian which includes terms for on-site Coulomb interaction for both intra- and inter-orbitals, intra-site exchange interaction and energy differences between orbitals. Based on our framework, we point out possible ways to investigate various phase transitions such as metal-insulator, magnetic and orbital order-disorder transitions without the minus sign problem. As an application, a two-band model is investigated by the projection QMC method and the ground state properties of this model are presented.
In the study of strongly correlated, many-electron systems, the Hubbard-Kanamori (HK) model has emerged as one of the prototypes for transition metal oxide physics. The model is multiband in nature and contains Hund's coupling terms, which have pronounced effects on metal-insulator transitions, high-temperature superconductivity, and other physical properties. In the following, we present a complete theoretical framework for treating the HK model using the ground-state auxiliary field quantum Monte Carlo (AFQMC) method and analyze its performance on few-band models whose parameters approximate those observed in ruthenates, rhodates, and other materials exhibiting Hund's physics. Unlike previous studies, the constrained path and phaseless approximations are used to respectively control the sign and phase problems, which enables high-accuracy modeling of the HK model's ground-state properties within parameter regimes of experimental interest. We demonstrate that, after careful consideration of the Hubbard-Stratonovich transformations and trial wave functions employed, relative errors in the energy of less than 11%1 can routinely be achieved for moderate to large values of the Hund's coupling constant. Crucially, our methodology also accurately predicts magnetic ordering and phase transitions. The results presented open the door to more predictive modeling of Hund's physics within a wide range of strongly correlated materials using AFQMC.
Ground state properties of multi-orbital Hubbard models are investigated by the auxiliary field quantum Monte Carlo method. A Monte Carlo technique generalized to the multi-orbital systems is introduced and examined in detail. The algorithm contains non-trivial cases where the negative sign problem does not exist. We investigate one-dimensional systems with doubly degenerate orbitals by this new technique. Properties of the Mott insulating state are quantitatively clarified as the strongly correlated insulator, where the charge gap amplitude is much larger than the spin gap. The insulator-metal transitions driven by the chemical potential shows a universality class with the correlation length exponent ν=1/2ν=1/2ν=1/2, which is consistent with the scaling arguments. Increasing level split between two orbitals drives crossover from the Mott insulator with high spin state to the band insulator with low spin state, where the spin gap amplitude increases and becomes closer to the charge gap. Experimental relevance of our results especially to Haldane materials is discussed.
The study of strongly correlated quantum many-body systems is an outstanding challenge. Highly accurate results are needed for the understanding of practical and fundamental problems in condensed-matter physics, high energy physics, material science, quantum chemistry and so on. Our familiar mean-field or perturbative methods tend to be ineffective. Numerical simulations provide a promising approach for studying such systems. The fundamental difficulty of numerical simulation is that the dimension of the Hilbert space needed to describe interacting systems increases exponentially with the system size. Quantum Monte Carlo (QMC) methods are one of the best approaches to tackle the problem of enormous Hilbert space. They have been highly successful for boson systems and unfrustrated spin models. For systems with fermions, the exchange symmetry in general causes the infamous sign problem, making the statistical noise in the computed results grow exponentially with the system size. This hinders our understanding of interesting physics such as high-temperature superconductivity, metal-insulator phase transition. In this thesis, we present a variety of new developments in the auxiliary-field quantum Monte Carlo (AFQMC) methods, including the incorporation of symmetry in both the trial wave function and the projector, developing the constraint release method, using the force-bias to drastically improve the efficiency in Metropolis framework, identifying and solving the infinite variance problem, and sampling Hartree-Fock-Bogoliubov wave function. With these developments, some of the most challenging many-electron problems are now under control. We obtain an exact numerical solution of two-dimensional strongly interacting Fermi atomic gas, determine the ground state properties of the 2D Fermi gas with Rashba spin-orbit coupling, provide benchmark results for the ground state of the two-dimensional Hubbard model, and establish that the Hubbard model has a stripe order in the underdoped region.
to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.
The acceleration of exact continuous time quantum Monte Carlo (CTQMC) approaches in multi-site or multi-orbital systems is extremely interesting work, because these approaches are very time-consuming in terms of numerical computation and might account for the nature of exotic behaviors such as high-temperature superconductivity and Mott insulator behavior observed in the strongly correlated materials. We extend the recently developed interaction-expansion CTQMC method in combination with a machine learning (CTQMC+ML) approach for the single-site and single-orbital systems to multi-site and multi-orbital ones. This method can be applied to explore the nonlocal correlation effects in lattice models and to study the electronic structure of real materials via an ab-initio density functional theory plus dynamical mean field theory approach. We find that our CTQMC+ML method for multi-site (and multi-orbital) systems accurately predicts the impurity Green’s function with less computational time than the CTQMC approaches, as in the case of the single-site and single-orbital version.
We present an algorithm for calculating the local energy of a multi-Slater wave function in orbital space quantum Monte Carlo (QMC). Recent developments in selected configuration interaction methods have led to increased interest in using multi-Slater trial wave functions in various QMC methods. For an ab initio Hamiltonian, our algorithm has a cost scaling of O(n5 + nc), as opposed to the O(n4nc) scaling of existing orbital space algorithms, where n is the system size and nc is the number of configurations in the wave function. We present our method using variational Monte Carlo calculations with the Jastrow multi-Slater wave function, although the formalism should be applicable for auxiliary field QMC. We apply it to polyacetylene and demonstrate the possibility of using a much larger number of configurations than possible using existing methods.