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:1708.04383

Constrained Path Monte Carlo method for spin $ 1/2 $ fermions at unitarity limit

Faisal Etminan, Mohammad Mehdi Firoozabadi

Published Aug 15, 2017• 0 citations
arXivHTMLPDF

Abstract

We present calculations for spin $ 1/2 $ fermions at unitarity limit, where the effective range of the interaction is zero and the scattering length is infinite. We compute the ground-state energy for a system of 6, 10,14,18 and 20 particles, with equal numbers of up and down spins in a periodic cube in the full ground-state constrained-path Monte Carlo (CPMC) method using the extended, attractive Hubbard model. Our results in a careful extrapolation to the thermodynamic limit may suggest that the ratio of the ground-state energy to that of a free Fermi gas is $ ξ= 0.43(4) $, which can compare with recent experimental results and consistent with the fixed node Green's function Monte Carlo and novel lattice approaches results. We also obtain results for interactions with different effective ranges and find that the energy is consistent with a universal linear dependence on the product of the Fermi momentum and the effective range.

DOI
10.48550/arxiv.1708.04383

Related papers

Ground-state energy of spin- 12 fermions in the unitary limit

Dean Lee

We present lattice results for the ground-state energy of a spin-$\frac{1}{2}$ fermion system in the unitary limit, where the effective range of the interaction is zero and the scattering length is infinite. We compute the ground-state energy for a system of 6, 10, 14, 18, and 22 particles, with equal numbers of up and down spins in a periodic cube. We estimate that in the limit of large number of particles, the ground-state energy is 0.25(3) times the ground-state energy of the free Fermi system.

arXiv:cond-mat/0511332•Physical Review B•Mar 17, 2006•66 citations

Ground state energy at unitarity

Dean Lee

We consider two-component fermions on the lattice in the unitarity limit. This is an idealized limit of attractive fermions where the range of the interaction is zero and the scattering length is infinite. Using Euclidean time projection, we compute the ground state energy using four computationally different but physically identical auxiliary-field methods. The best performance is obtained using a bounded continuous auxiliary field and a nonlocal updating algorithm called the hybrid Monte Carlo. With this method, we calculate results for 10 and 14 fermions at lattice volumes ${4}^{3},{5}^{3},{6}^{3},{7}^{3},{8}^{3}$ and extrapolate to the continuum limit. For 10 fermions in a periodic cube, the ground state energy is $0.292(12)$ times the ground state energy for noninteracting fermions. For 14 fermions, the ratio is $0.329(5)$.

arXiv:0803.1280•Physical Review C•Aug 13, 2008•67 citations

Precision benchmark calculations for four particles at unitarity

Shahin Bour, Xin Li, Dean Lee +2 more

The unitarity limit describes interacting particles where the range of the interaction is zero and the scattering length is infinite. We present precision benchmark calculations for two-component fermions at unitarity using three different ab initio methods: Hamiltonian lattice formalism using iterated eigenvector methods, Euclidean lattice formalism with auxiliary-field projection Monte Carlo methods, and continuum diffusion Monte Carlo methods with fixed and released nodes. We have calculated the ground-state energy of the unpolarized four-particle system in a periodic cube as a dimensionless fraction of the ground-state energy for the noninteracting system. We obtain values of $0.211(2)$ and $0.210(2)$ using two different Hamiltonian lattice representations, $0.206(9)$ using Euclidean lattice formalism, and an upper bound of $0.212(2)$ from fixed-node diffusion Monte Carlo methods. Released-node calculations starting from the fixed-node result yield a decrease of less than $0.002$ over a propagation of $0.4{E}_{F}^{\ensuremath{-}1}$ in Euclidean time, where ${E}_{F}$ is the Fermi energy. We find good agreement among all three ab initio methods.

arXiv:1104.2102•Physical Review A•Jun 15, 2011•22 citations

Atomic Fermi gas in the unitary limit by quantum Monte Carlo methods: Effects of the interaction range

Xin Li, Jindřich Kolorenč, Luboš Mitáš

We calculate the ground-state properties of an unpolarized two-component Fermi gas with the aid of the diffusion quantum Monte Carlo (DMC) methods. Using an extrapolation to the zero effective range of the attractive two-particle interaction, we find $E/{E}_{\mathrm{free}}$ in the unitary limit to be 0.212(2), 0.407(2), 0.409(3), and 0.398(3) for 4, 14, 38, and 66 atoms, respectively. Our calculations indicate that the dependence of the total energy on the effective range of the interaction ${R}_{\mathrm{eff}}$ is sizable and the extrapolation to ${R}_{\mathrm{eff}}=0$ is therefore important for reaching the true unitary limit. To test the quality of nodal surfaces and to estimate the impact of the fixed-node approximation, we perform released-node DMC calculations for 4 and 14 atoms. Analysis of the released-node and the fixed-node results suggests that the main sources of the fixed-node errors are long-range correlations, which are difficult to sample in the released-node approaches due to the fast growth of the bosonic noise. Besides energies, we evaluate the two-body density matrix and the condensate fraction. We find that the condensate fraction for the 66-atom system converges to 0.56(1) after the extrapolation to the zero interaction range.

arXiv:1105.1748•Physical Review A•Aug 8, 2011•20 citations

Configuration-interaction Monte Carlo method and its application to the trapped unitary Fermi gas

Abhishek Mukherjee, Y. Alhassid

We develop a quantum Monte Carlo method to estimate the ground-state energy of a fermionic many-particle system in the configuration-interaction shell model approach. The fermionic sign problem is circumvented by using a guiding wave function in Fock space. The method provides an upper bound on the ground-state energy whose tightness depends on the choice of the guiding wave function. We argue that the antisymmetric geminal product class of wave functions is a good choice for guiding wave functions. We demonstrate our method for the trapped two-species fermionic cold atom system in the unitary regime of infinite scattering length using the particle-number projected Hartree-Fock-Bogoliubov wave function as the guiding wave function. We estimate the ground-state energy and energy-staggering pairing gap as a function of the number of particles. We compare our results with exact numerical diagonalization results and with previous fixed-node coordinate-space Monte Carlo calculations.

arXiv:1304.1645•Physical Review A•Nov 18, 2013•16 citations

The Ground State Of Two-Dimensional Hubbard-Like Models

Allman, Eric

We present results from a constrained path Monte Carlo (CPMC) study of a modified two-dimensional Hubbard model. We include more general forms of the band structure and electron interaction in order to examine their effects on ground-state properties, such as electron pairing correlations. Both next-nearest neighbor hopping, t', and third-nearest neighbor hopping, t″, are introduced in the Hamiltonian. A nearest neighbor interaction of strength V is also included. We carry out CPMC calculations on system sizes up to 16 x 16, at various electron fillings, to investigate the ground state of the model for different values of these parameters. For realistic systems, these calculations indicate that Hubbard-like models are not capable of showing enhanced pairing. The modified Hamiltonian also presents an opportunity to more closely examine the accuracy and robustness of the CPMC algorithm. Results of further benchmark calculations involving CPMC are presented. These benchmarks confirm that CPMC is able to show enhanced pairing in systems where such behavior is known to exist.

•Jan 1, 2002•0 citations
We present calculations for spin 1/21/21/2 fermions at unitarity limit, where the effective range of the interaction is zero and the scattering length is infinite. We compute the ground-state energy for a system of 6, 10,14,18 and 20 particles, with equal numbers of up and down spins in a periodic cube in the full ground-state constrained-path Monte Carlo (CPMC) method using the extended, attractive Hubbard model. Our results in a careful extrapolation to the thermodynamic limit may suggest that the ratio of the ground-state energy to that of a free Fermi gas is ξ=0.43(4)ξ= 0.43(4)ξ=0.43(4), which can compare with recent experimental results and consistent with the fixed node Green's function Monte Carlo and novel lattice approaches results. We also obtain results for interactions with different effective ranges and find that the energy is consistent with a universal linear dependence on the product of the Fermi momentum and the effective range.
We present lattice results for the ground-state energy of a spin-12\frac{1}{2}21​ fermion system in the unitary limit, where the effective range of the interaction is zero and the scattering length is infinite. We compute the ground-state energy for a system of 6, 10, 14, 18, and 22 particles, with equal numbers of up and down spins in a periodic cube. We estimate that in the limit of large number of particles, the ground-state energy is 0.25(3) times the ground-state energy of the free Fermi system.
We consider two-component fermions on the lattice in the unitarity limit. This is an idealized limit of attractive fermions where the range of the interaction is zero and the scattering length is infinite. Using Euclidean time projection, we compute the ground state energy using four computationally different but physically identical auxiliary-field methods. The best performance is obtained using a bounded continuous auxiliary field and a nonlocal updating algorithm called the hybrid Monte Carlo. With this method, we calculate results for 10 and 14 fermions at lattice volumes 43,53,63,73,83{4}^{3},{5}^{3},{6}^{3},{7}^{3},{8}^{3}43,53,63,73,83 and extrapolate to the continuum limit. For 10 fermions in a periodic cube, the ground state energy is 0.292(12)0.292(12)0.292(12) times the ground state energy for noninteracting fermions. For 14 fermions, the ratio is 0.329(5)0.329(5)0.329(5).
The unitarity limit describes interacting particles where the range of the interaction is zero and the scattering length is infinite. We present precision benchmark calculations for two-component fermions at unitarity using three different ab initio methods: Hamiltonian lattice formalism using iterated eigenvector methods, Euclidean lattice formalism with auxiliary-field projection Monte Carlo methods, and continuum diffusion Monte Carlo methods with fixed and released nodes. We have calculated the ground-state energy of the unpolarized four-particle system in a periodic cube as a dimensionless fraction of the ground-state energy for the noninteracting system. We obtain values of 0.211(2)0.211(2)0.211(2) and 0.210(2)0.210(2)0.210(2) using two different Hamiltonian lattice representations, 0.206(9)0.206(9)0.206(9) using Euclidean lattice formalism, and an upper bound of 0.212(2)0.212(2)0.212(2) from fixed-node diffusion Monte Carlo methods. Released-node calculations starting from the fixed-node result yield a decrease of less than 0.0020.0020.002 over a propagation of $0.4{E}_{F}^{\ensuremath{-}1}$ in Euclidean time, where EF{E}_{F}EF​ is the Fermi energy. We find good agreement among all three ab initio methods.
We calculate the ground-state properties of an unpolarized two-component Fermi gas with the aid of the diffusion quantum Monte Carlo (DMC) methods. Using an extrapolation to the zero effective range of the attractive two-particle interaction, we find E/EfreeE/{E}_{\mathrm{free}}E/Efree​ in the unitary limit to be 0.212(2), 0.407(2), 0.409(3), and 0.398(3) for 4, 14, 38, and 66 atoms, respectively. Our calculations indicate that the dependence of the total energy on the effective range of the interaction Reff{R}_{\mathrm{eff}}Reff​ is sizable and the extrapolation to Reff=0{R}_{\mathrm{eff}}=0Reff​=0 is therefore important for reaching the true unitary limit. To test the quality of nodal surfaces and to estimate the impact of the fixed-node approximation, we perform released-node DMC calculations for 4 and 14 atoms. Analysis of the released-node and the fixed-node results suggests that the main sources of the fixed-node errors are long-range correlations, which are difficult to sample in the released-node approaches due to the fast growth of the bosonic noise. Besides energies, we evaluate the two-body density matrix and the condensate fraction. We find that the condensate fraction for the 66-atom system converges to 0.56(1) after the extrapolation to the zero interaction range.
We develop a quantum Monte Carlo method to estimate the ground-state energy of a fermionic many-particle system in the configuration-interaction shell model approach. The fermionic sign problem is circumvented by using a guiding wave function in Fock space. The method provides an upper bound on the ground-state energy whose tightness depends on the choice of the guiding wave function. We argue that the antisymmetric geminal product class of wave functions is a good choice for guiding wave functions. We demonstrate our method for the trapped two-species fermionic cold atom system in the unitary regime of infinite scattering length using the particle-number projected Hartree-Fock-Bogoliubov wave function as the guiding wave function. We estimate the ground-state energy and energy-staggering pairing gap as a function of the number of particles. We compare our results with exact numerical diagonalization results and with previous fixed-node coordinate-space Monte Carlo calculations.
We present results from a constrained path Monte Carlo (CPMC) study of a modified two-dimensional Hubbard model. We include more general forms of the band structure and electron interaction in order to examine their effects on ground-state properties, such as electron pairing correlations. Both next-nearest neighbor hopping, t', and third-nearest neighbor hopping, t″, are introduced in the Hamiltonian. A nearest neighbor interaction of strength V is also included. We carry out CPMC calculations on system sizes up to 16 x 16, at various electron fillings, to investigate the ground state of the model for different values of these parameters. For realistic systems, these calculations indicate that Hubbard-like models are not capable of showing enhanced pairing. The modified Hamiltonian also presents an opportunity to more closely examine the accuracy and robustness of the CPMC algorithm. Results of further benchmark calculations involving CPMC are presented. These benchmarks confirm that CPMC is able to show enhanced pairing in systems where such behavior is known to exist.