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:0903.2519 · Physical Review B

Quantum Monte Carlo study of the two-dimensional fermion Hubbard model

Christopher N. Varney, C.-R. Lee, Zhaojun Bai, Simone Chiesa, Mark Jarrell, Richard T. Scalettar

Published Aug 18, 2009• 164 citations
arXivHTMLPDF

Abstract

We report large scale determinant quantum Monte Carlo calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian at half-filling. The sharp Fermi surface of the noninteracting limit is significantly broadened by the electronic correlations but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite-size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak-coupling to the strong-coupling Heisenberg limit. Our lattices provide improved resolution of the Green's function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.

DOI
10.1103/physrevb.80.075116

Related papers

Determinant quantum Monte Carlo study of the two-dimensional single-band Hubbard-Holstein model

Steven Johnston, Elizabeth Nowadnick, Y. F. Kung +3 more

We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here, we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic orders at half-filling.

arXiv:1306.2968•Physical Review B•Jun 24, 2013•82 citations

Phase transitions in strongly interacting quantum systems

Christopher N. Varney

Strongly correlated materials exhibit some of the most complex and technologically useful phenomena in condensed matter. We use the world line and determinant quantum Monte Carlo techniques to explore the phases of the Hubbard Hamiltonian, a model describing strong interactions, in a number of different contexts. We first consider the spin, charge, and bond order correlations of the one-dimensional extended fermion Hubbard model in the presence of a coupling to the lattice. A static alternating lattice distortion leads to enhanced charge density wave correlations at the expense of antiferromagnetic order. When the lattice degrees of freedom are dynamic, we show that a similar effect occurs even though the charge asymmetry must arise spontaneously. Although the evolution of the total energy with lattice coupling is smooth, the individual components exhibit sharp crossovers at the phase boundaries. Finally, we observe a tendency for bond order in the region between the charge and spin density wave phases. Second, we examine mixtures of bosons and fermions in one-dimensional optical lattices. We evaluate the density profiles and bosonic visibility [special characters omitted] resolving the discrepancy between theory and experiment by identifying parameter regimes where [special characters omitted] is reduced and increased. We present a simple qualitative picture of the different response to the fermion admixture in terms of the superfluid and Mott-insulating domains before and after the fermions are included. Finally, we show that [special characters omitted] exhibits kinks which are tied to the domain evolution present in the pure bosonic case, and also additional structure arising from the formation of boson-fermion molecules, a prediction for future experiments. Third, we report large scale calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak-coupling to the strong-coupling Heisenberg limit. Our lattices provide improved resolution of the momentum distribution, allowing a more quantitative comparison with time-of-flight optical lattice experiments.

•Jan 1, 2009•0 citations

Quantum Monte Carlo study of the half-filled Hubbard model on the honeycomb lattice

Yuichi Otsuka, Seiji Yunoki, Sandro Sorella

The half-filled Hubbard model on the honeycomb lattice is investigated by numerically exact large-scale quantum Monte Carlo simulations for lattice sizes up to 2592 sites. By performing careful finite-size scaling for the spin-spin correlation functions, calculated with a high degree of accuracy, we find that the ground state is antiferromagnetically long-range ordered at U / t = 4, where U is the on-site Hubbard interaction and t is the nearest neighbor hopping. Our result is in sharp contrast to a recent report [Meng et al. , Nature 464, 847 (2010)], where instead strong evidence of spin liquid behavior is found.

Journal of Physics Conference Series•Aug 12, 2013•2 citations

Quantum Monte Carlo study of the one-dimensional Holstein model of spinless fermions

Ross H. McKenzie, Chris J. Hamer, D. W. Murray

The Holstein model of spinless fermions interacting with dispersionless phonons in one dimension is studied by a Green's function Monte Carlo technique. The ground-state energy, first fermionic excited state, density wave correlations, and mean lattice displacements are calculated for lattices of up to 16 sites, for one fermion per two sites, i.e., a half-filled band. Results are obtained for values of the fermion hopping parameter of t=0.1\ensuremath{\omega}, \ensuremath{\omega}, and 10\ensuremath{\omega}, where \ensuremath{\omega} is the phonon frequency. At a finite fermion-phonon coupling g there is a transition from a metallic phase to an insulating phase in which there is charge-density-wave order. Finite-size scaling is found to hold in the metallic phase and is used to extract the coupling dependence of the Luttinger liquid parameters, ${\mathit{u}}_{\mathrm{\ensuremath{\rho}}}$ and ${\mathit{K}}_{\mathrm{\ensuremath{\rho}}}$, the velocity of charge excitations, and the correlation exponent, respectively. For free fermions (g=0) and for strong coupling (${\mathit{g}}^{2}$\ensuremath{\gg}t\ensuremath{\omega}) our results agree well with known analytic results. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9512046•Physical review. B, Condensed matter•Apr 15, 1996•62 citations

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Mingpu Qin, Hao Shi, Shiwei Zhang

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of $U/t$ from 2 to 8, for a range of densities including half-filling and $n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75$, and $0.875$. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.

arXiv:1605.09421•Physical review. B./Physical review. B•Aug 4, 2016•100 citations

Two-Dimensional Hubbard Model –Metal Insulator Transition Studied by Monte Carlo Calculation–

Nobuo Furukawa, Masatoshi Imada

The Hubbard model on a square lattice in the ground state is investigated. Various aspects of Mott transition at half-filling are clarified from the quantum Monte Carlo study. Critical exponents of the transition for charge and spin correlations are estimated. For the doping concentration δ, the charge susceptibility is proportional to δ -1 indicating the divergence of charge mass for δ→0, while the spin susceptibility is always finite and not singular near δ=0. Several theoretical consequences of the above results are discussed. Incommensurate peak amplitude of spin correlation scales as (δ-δ c ) -1 with δ c <0.01 in disagreement with RPA results ln (δ-δ c ). The antiferromagnetic order at the half-filling also shows strong correlation character and is not sensitive to the shape of the fermi surface, while the incommensurate peak observed away from the half-filling is sensitively suppressed by the loss of partial nesting consistently with the weak coupling picture. Comparisons with experimental indications of high- T c oxides are made.

Journal of the Physical Society of Japan•Sep 15, 1992•218 citations
We report large scale determinant quantum Monte Carlo calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian at half-filling. The sharp Fermi surface of the noninteracting limit is significantly broadened by the electronic correlations but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite-size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak-coupling to the strong-coupling Heisenberg limit. Our lattices provide improved resolution of the Green's function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.
We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here, we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic orders at half-filling.
Strongly correlated materials exhibit some of the most complex and technologically useful phenomena in condensed matter. We use the world line and determinant quantum Monte Carlo techniques to explore the phases of the Hubbard Hamiltonian, a model describing strong interactions, in a number of different contexts. We first consider the spin, charge, and bond order correlations of the one-dimensional extended fermion Hubbard model in the presence of a coupling to the lattice. A static alternating lattice distortion leads to enhanced charge density wave correlations at the expense of antiferromagnetic order. When the lattice degrees of freedom are dynamic, we show that a similar effect occurs even though the charge asymmetry must arise spontaneously. Although the evolution of the total energy with lattice coupling is smooth, the individual components exhibit sharp crossovers at the phase boundaries. Finally, we observe a tendency for bond order in the region between the charge and spin density wave phases. Second, we examine mixtures of bosons and fermions in one-dimensional optical lattices. We evaluate the density profiles and bosonic visibility [special characters omitted] resolving the discrepancy between theory and experiment by identifying parameter regimes where [special characters omitted] is reduced and increased. We present a simple qualitative picture of the different response to the fermion admixture in terms of the superfluid and Mott-insulating domains before and after the fermions are included. Finally, we show that [special characters omitted] exhibits kinks which are tied to the domain evolution present in the pure bosonic case, and also additional structure arising from the formation of boson-fermion molecules, a prediction for future experiments. Third, we report large scale calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak-coupling to the strong-coupling Heisenberg limit. Our lattices provide improved resolution of the momentum distribution, allowing a more quantitative comparison with time-of-flight optical lattice experiments.
The half-filled Hubbard model on the honeycomb lattice is investigated by numerically exact large-scale quantum Monte Carlo simulations for lattice sizes up to 2592 sites. By performing careful finite-size scaling for the spin-spin correlation functions, calculated with a high degree of accuracy, we find that the ground state is antiferromagnetically long-range ordered at U / t = 4, where U is the on-site Hubbard interaction and t is the nearest neighbor hopping. Our result is in sharp contrast to a recent report [Meng et al. , Nature 464, 847 (2010)], where instead strong evidence of spin liquid behavior is found.
The Holstein model of spinless fermions interacting with dispersionless phonons in one dimension is studied by a Green's function Monte Carlo technique. The ground-state energy, first fermionic excited state, density wave correlations, and mean lattice displacements are calculated for lattices of up to 16 sites, for one fermion per two sites, i.e., a half-filled band. Results are obtained for values of the fermion hopping parameter of t=0.1\ensuremath{\omega}, \ensuremath{\omega}, and 10\ensuremath{\omega}, where \ensuremath{\omega} is the phonon frequency. At a finite fermion-phonon coupling g there is a transition from a metallic phase to an insulating phase in which there is charge-density-wave order. Finite-size scaling is found to hold in the metallic phase and is used to extract the coupling dependence of the Luttinger liquid parameters, ${\mathit{u}}_{\mathrm{\ensuremath{\rho}}}$ and ${\mathit{K}}_{\mathrm{\ensuremath{\rho}}}$, the velocity of charge excitations, and the correlation exponent, respectively. For free fermions (g=0) and for strong coupling (g2{\mathit{g}}^{2}g2\ensuremath{\gg}t\ensuremath{\omega}) our results agree well with known analytic results. \textcopyright{} 1996 The American Physical Society.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/tU/tU/t from 2 to 8, for a range of densities including half-filling and n=0.3,0.5,0.6,0.75n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75n=0.3,0.5,0.6,0.75, and 0.8750.8750.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
The Hubbard model on a square lattice in the ground state is investigated. Various aspects of Mott transition at half-filling are clarified from the quantum Monte Carlo study. Critical exponents of the transition for charge and spin correlations are estimated. For the doping concentration δ, the charge susceptibility is proportional to δ -1 indicating the divergence of charge mass for δ→0, while the spin susceptibility is always finite and not singular near δ=0. Several theoretical consequences of the above results are discussed. Incommensurate peak amplitude of spin correlation scales as (δ-δ c ) -1 with δ c <0.01 in disagreement with RPA results ln (δ-δ c ). The antiferromagnetic order at the half-filling also shows strong correlation character and is not sensitive to the shape of the fermi surface, while the incommensurate peak observed away from the half-filling is sensitively suppressed by the loss of partial nesting consistently with the weak coupling picture. Comparisons with experimental indications of high- T c oxides are made.