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:1407.2784 · Physical Review C

Fermion-dimer scattering using an impurity lattice Monte Carlo approach and the adiabatic projection method

Serdar Elhatisari, Dean Lee

Published Dec 1, 2014• 23 citations
arXivHTMLPDF

Abstract

We present lattice Monte Carlo calculations of fermion-dimer scattering in the limit of zero-range interactions using the adiabatic projection method. The adiabatic projection method uses a set of initial cluster states and Euclidean time projection to give a systematically improvable description of the low-lying scattering cluster states in a finite volume. We use L\"uscher's finite-volume relations to determine the $s$-wave, $p$-wave, and $d$-wave phase shifts. For comparison, we also compute exact lattice results using Lanczos iteration and continuum results using the Skorniakov-Ter-Martirosian equation. For our Monte Carlo calculations we use a new lattice algorithm called impurity lattice Monte Carlo. This algorithm can be viewed as a hybrid technique which incorporates elements of both worldline and auxiliary-field Monte Carlo simulations.

DOI
10.1103/physrevc.90.064001

Related papers

Low Energy Continuum and Lattice Effective Field Theories

Serdar Elhatisari

In the first part of the thesis we consider the constraints of causality and unitarity for particles interacting via strictly finite-range interactions. We generalize Wigner's causality bound to the case of non-vanishing partial-wave mixing. Specifically we analyze the system of the low-energy interactions between protons and neutrons. We also analyze low-energy scattering for systems with arbitrary short-range interactions plus an attractive $1/r^α$ tail for $α\geq2$. In particular, we focus on the case of $α=6$ and we derive the constraints of causality and unitarity also for these systems and find that the van der Waals length scale dominates over parameters characterizing the short-distance physics of the interaction. This separation of scales suggests a separate universality class for physics characterizing interactions with an attractive $1/r^{6}$ tail. We argue that a similar universality class exists for any attractive potential $1/r^α$ for $α\geq2$. In the second part of the thesis we present lattice Monte Carlo calculations of fermion-dimer scattering in the limit of zero-range interactions using the adiabatic projection method. The adiabatic projection method uses a set of initial cluster states and Euclidean time projection to give a systematically improvable description of the low-lying scattering cluster states in a finite volume. We use Lüscher's finite-volume relations to determine the $s$-wave, $p$-wave, and $d$-wave phase shifts. For comparison, we also compute exact lattice results using Lanczos iteration and continuum results using the Skorniakov-Ter-Martirosian equation. For our Monte Carlo calculations we use a new lattice algorithm called impurity lattice Monte Carlo. This algorithm can be viewed as a hybrid technique which incorporates elements of both worldline and auxiliary-field Monte Carlo simulations.

arXiv:1409.4048•Sep 14, 2014•0 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

Benchmark calculations for elastic fermion-dimer scattering

Shahin Bour, H.‐W. Hammer, Dean Lee +1 more

We present continuum and lattice calculations for elastic scattering between a fermion and a bound dimer in the shallow binding limit. For the lattice calculation we use the finite-volume method of L\"uscher. We take into account topological finite-volume corrections to the dimer binding energy which depend on the momentum of the dimer. After subtracting these effects, we find from the lattice calculation $\ensuremath{\kappa}{a}_{fd}=1.174(9)$ and $\ensuremath{\kappa}{r}_{fd}=\ensuremath{-}0.029(13)$. These results agree well with the continuum values $\ensuremath{\kappa}{a}_{fd}=1.17907(1)$ and $\ensuremath{\kappa}{r}_{fd}=\ensuremath{-}0.0383(3)$ obtained from the STM equation. We discuss applications to cold atomic Fermi gases, deuteron-neutron scattering in the spin-quartet channel, and lattice calculations of scattering for nuclei and hadronic molecules at finite volume.

arXiv:1206.1765•Physical Review C•Sep 14, 2012•44 citations

Worldline Monte Carlo method for few-body nuclear physics

Shailesh Chandrasekharan A, Son T. Nguyen, Thomas R. Richardson

In this work we introduce a worldline-based fermion Monte Carlo algorithm for studying few body quantum mechanics of self-interacting fermions in the Hamiltonian lattice formulation. Our motivation to construct the method comes from our interest in studying renormalization of chiral nuclear effective-field theory with lattice regularization. In particular we wish to apply our method to compute the lattice spacing dependence of local lattice interactions as we take the continuum limit of the lattice theory. Our algorithm can compute matrix elements of the operator $exp(\ensuremath{-}\ensuremath{\beta}H)$, where $H$ is the lattice Hamiltonian and $\ensuremath{\beta}$ is a free real parameter. These elements help us compute deep bound states that are well separated from scattering states even at values of $\ensuremath{\beta}$ which are not very large. Computing these bound-state energies accurately can help us study renormalization of the lattice theory. In addition to developing the algorithm, in this work we also introduce a finite-volume renormalization scheme for the lattice Hamiltonian of the leading pionless effective-field theory and show how it would work in the one- and two-body sectors.

arXiv:2402.15377•Physical Review C•Aug 14, 2024•4 citations

Ab initio Lattice Results for Fermi Polarons in Two Dimensions

Shahin Bour, Dean Lee, H.‐W. Hammer +1 more

We investigate the attractive Fermi polaron problem in two dimensions using nonperturbative Monte Carlo simulations. We introduce a new Monte Carlo algorithm called the impurity lattice Monte Carlo method. This algorithm samples the path integral in a computationally efficient manner and has only small sign oscillations for systems with a single impurity. As a benchmark of the method, we calculate the universal polaron energy in three dimensions in the scale-invariant unitarity limit and find agreement with published results. We then present the first fully nonperturbative calculations of the polaron energy in two dimensions and density correlations between the impurity and majority particles in the limit of zero-range interactions. We find evidence for a smooth crossover transition from fermionic quasiparticle to molecular state as a function of the interaction strength.

arXiv:1412.8175•Physical Review Letters•Oct 27, 2015•32 citations

Quantum Monte Carlo impurity solvers for multi-orbital problems and frequency-dependent interactions

Hiroshi Shinaoka, Fakher F. Assaad, Nils Blümer +1 more

The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.

The European Physical Journal Special Topics•Jul 1, 2017•11 citations
We present lattice Monte Carlo calculations of fermion-dimer scattering in the limit of zero-range interactions using the adiabatic projection method. The adiabatic projection method uses a set of initial cluster states and Euclidean time projection to give a systematically improvable description of the low-lying scattering cluster states in a finite volume. We use L"uscher's finite-volume relations to determine the sss-wave, ppp-wave, and ddd-wave phase shifts. For comparison, we also compute exact lattice results using Lanczos iteration and continuum results using the Skorniakov-Ter-Martirosian equation. For our Monte Carlo calculations we use a new lattice algorithm called impurity lattice Monte Carlo. This algorithm can be viewed as a hybrid technique which incorporates elements of both worldline and auxiliary-field Monte Carlo simulations.
In the first part of the thesis we consider the constraints of causality and unitarity for particles interacting via strictly finite-range interactions. We generalize Wigner's causality bound to the case of non-vanishing partial-wave mixing. Specifically we analyze the system of the low-energy interactions between protons and neutrons. We also analyze low-energy scattering for systems with arbitrary short-range interactions plus an attractive 1/rα1/r^α1/rα tail for α≥2α\geq2α≥2. In particular, we focus on the case of α=6α=6α=6 and we derive the constraints of causality and unitarity also for these systems and find that the van der Waals length scale dominates over parameters characterizing the short-distance physics of the interaction. This separation of scales suggests a separate universality class for physics characterizing interactions with an attractive 1/r61/r^{6}1/r6 tail. We argue that a similar universality class exists for any attractive potential 1/rα1/r^α1/rα for α≥2α\geq2α≥2. In the second part of the thesis we present lattice Monte Carlo calculations of fermion-dimer scattering in the limit of zero-range interactions using the adiabatic projection method. The adiabatic projection method uses a set of initial cluster states and Euclidean time projection to give a systematically improvable description of the low-lying scattering cluster states in a finite volume. We use Lüscher's finite-volume relations to determine the sss-wave, ppp-wave, and ddd-wave phase shifts. For comparison, we also compute exact lattice results using Lanczos iteration and continuum results using the Skorniakov-Ter-Martirosian equation. For our Monte Carlo calculations we use a new lattice algorithm called impurity lattice Monte Carlo. This algorithm can be viewed as a hybrid technique which incorporates elements of both worldline and auxiliary-field Monte Carlo simulations.
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).
We present continuum and lattice calculations for elastic scattering between a fermion and a bound dimer in the shallow binding limit. For the lattice calculation we use the finite-volume method of L\"uscher. We take into account topological finite-volume corrections to the dimer binding energy which depend on the momentum of the dimer. After subtracting these effects, we find from the lattice calculation $\ensuremath{\kappa}{a}_{fd}=1.174(9)$ and $\ensuremath{\kappa}{r}_{fd}=\ensuremath{-}0.029(13)$. These results agree well with the continuum values $\ensuremath{\kappa}{a}_{fd}=1.17907(1)$ and $\ensuremath{\kappa}{r}_{fd}=\ensuremath{-}0.0383(3)$ obtained from the STM equation. We discuss applications to cold atomic Fermi gases, deuteron-neutron scattering in the spin-quartet channel, and lattice calculations of scattering for nuclei and hadronic molecules at finite volume.
In this work we introduce a worldline-based fermion Monte Carlo algorithm for studying few body quantum mechanics of self-interacting fermions in the Hamiltonian lattice formulation. Our motivation to construct the method comes from our interest in studying renormalization of chiral nuclear effective-field theory with lattice regularization. In particular we wish to apply our method to compute the lattice spacing dependence of local lattice interactions as we take the continuum limit of the lattice theory. Our algorithm can compute matrix elements of the operator $exp(\ensuremath{-}\ensuremath{\beta}H)$, where HHH is the lattice Hamiltonian and $\ensuremath{\beta}$ is a free real parameter. These elements help us compute deep bound states that are well separated from scattering states even at values of $\ensuremath{\beta}$ which are not very large. Computing these bound-state energies accurately can help us study renormalization of the lattice theory. In addition to developing the algorithm, in this work we also introduce a finite-volume renormalization scheme for the lattice Hamiltonian of the leading pionless effective-field theory and show how it would work in the one- and two-body sectors.
We investigate the attractive Fermi polaron problem in two dimensions using nonperturbative Monte Carlo simulations. We introduce a new Monte Carlo algorithm called the impurity lattice Monte Carlo method. This algorithm samples the path integral in a computationally efficient manner and has only small sign oscillations for systems with a single impurity. As a benchmark of the method, we calculate the universal polaron energy in three dimensions in the scale-invariant unitarity limit and find agreement with published results. We then present the first fully nonperturbative calculations of the polaron energy in two dimensions and density correlations between the impurity and majority particles in the limit of zero-range interactions. We find evidence for a smooth crossover transition from fermionic quasiparticle to molecular state as a function of the interaction strength.
The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.