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

Investigation of Spin-Polaron Breakdown in the Single-Hole Doped Two-Dimensional t–t′–Jz Model

Matthias Bernhart

Published Apr 1, 2026• 0 citations
Publisher

Abstract

Solid state physics, a very broad field in natural science, aims to explain physical properties of materials modeled as a many-body ensemble of the order of ∼ 10^23 interacting particles. For interacting metallic systems, one of the most powerful tools for unraveling the consequences of the collective interaction is the Fermi liquid theory. While still being an approximation, it excellently models many materials especially in the low-temperature and low-energy regime by treating themany-body problem as being equivalent to effectively non-interacting ”particles”, which is done by introducing the concept of quasiparticles. During the last few decades, several exotic phases of matter have been discovered, which originate from strong Coulomb interactions between the electrons, whose description is beyond the Fermi liquid theory. This thesis addresses the question of the validity of Fermi liquid theory once a single hole is introduced into an antiferromagnet with completely localized spins due to the strong interactions. Here, the most-studied case concerns the so-called 2D t − J model, which hosts a well-defined quasiparticle, the spin polaron. Focusing on a more tractable case, the t − Jz model without quantum spin fluctuations, we study the stability of the spin polaron solution once the nextnearest neighbor hopping t′ is introduced into the system. We focus on observations at the Γ-point and the X-point of systems with fixed Jz = 0.4t and varying t′ ∈ [−0.5t, 0.5t]. The model is numerically solved using exact diagonalization on finite-size lattices of up to 32 sites. A study of the quasiparticle weight under finite-size scaling reveals an “anomalous” regime at the Γ-point fort′ ≤ −0.3t, indicated by a vanishing quasiparticle weight. In order to characterize the low-energy dynamics within this “anomalous” regime, we investigated the spin correlation centered around the hole, as well as the magnon number distribution of a single-doped 20-site system. Based on the results, we present a preliminary picture of the antiferromagnetic ground state doped with one hole within that “anomalous” regime (t′ ≤ −0.3t) that deviates from the standard spin polaron as the cloud of magnetic fluctuations, typically localized around the hole, delocalizes over the whole finite-size lattice. Consequently, such systems show a weakly suppressed antiferromagnetic configuration under single-hole doping.

DOI
10.34726/hss.2026.137502

Related papers

Numerical Studies of Quantum Spin Systems

C. Brünger

The understanding of matter does not only require a precise analysis of individual characteristics of the fundamental constituents, but also the study of their collective behavior.Solid state physics, and in particular research in strongly correlated electron systems, deals with such collective phenomena.The properties of the electrons, charge and spin quantum numbers, do not on their own explain the rich world of states of matter, however mutual interactions lead to various states or phases of matter.This work focuses on phases in quantum spin systems.On the basis of the Heisenberg model one (1D) and two dimensional (2D) spin lattices are investigated.In a first part the bilayer Heisenberg Model and the 2D Kondo necklace model are studied.Both models exhibit a quantum phase transition between an ordered and disordered phase.For small interlayer couplings the models have antiferromagnetic long-range order, which is destroyed by increasing the coupling.However, even in the disordered phase the models show global features in terms of collective excitations called magnons.Coming from the disordered phase the gapped magnons condense at a critical coupling and form the order of the ordered phase.Here, another interesting question arises, namely that of the coupling of a single doped hole to such critical fluctuations.A self-consistent Born approximation predicts that the doped hole couples to the magnons such that the quasiparticle residue vanishes at the quantum critical point.In this work the delicate question about the fate of the quasiparticle residue across the quantum phase transition is also tackled by means of large scale quantum Monte Carlo simulations.Furthermore the dynamics of a single hole doped in the magnetic background is investigated.

•Jan 1, 2007•0 citations

Spin polarons in the t - J model

Gerardo Martı́nez, Peter Horsch

The motion of a single hole in a two-dimensional Heisenberg antiferromagnet (hAF) is studied in a representation where holes are described as spinless fermions (holons) and spins as normal bosons. Assuming long-range AF order the spin dynamics is treated in linear spin-wave theory. The formulation highlights the close relation with the conventional polaron problem. The holon Green's function is calculated self-consistently within the Born approximation using finite-cluster geometries for the numerical solution. As a remarkable result we find close agreement with the spectral function A(k,\ensuremath{\omega}) of a hole calculated by exact diagonalization methods. A(k,\ensuremath{\omega}) is characterized by a narrow quasiparticle (QP) peak at the low-energy side of the spectrum, which is well separated from the incoherent part for large enough J values. A complete characterization of our solution is given, including the spectral weight, the dispersion relation, and effective masses of the QP state. A finite-size-scaling study gives a nonvanishing spectral weight of the QP in the thermodynamic limit for values J/t typical for copper oxide superconductors. Our calculations indicate that the self-consistent Born approximation is a valuable scheme for characterizing the dynamics of a hole in the t-J model, even in the strong-coupling regime.

Physical review. B, Condensed matter•Jul 1, 1991•498 citations

A model study of momentum-selective Mott physics

Sebastian Albert Erwin Huber

The Fermi liquid theory is a central concept in modern condensed matter physics used to describe conventional metals. A state of this universality class consists of well-defined quasiparticles, which occupy a finite number of states within a Fermi volume defined in momentum space. According to the Luttinger’s theorem, such a state encloses a Fermi volume proportional to the electron density modulo the filled bands if no symmetry is broken. In the past decades strong deviations from Fermi liquid theory have been observed in the scaling behavior of thermodynamic observables in different materials, among which are for instance cuprates and iron pnictides. This thesis is about two prototypical systems that can not be described by Fermi liquid theory. The first part of this work investigates two-dimensional effective models, which are in a Mott insulating state at half-filling and have a finite conductivity when doped with holes. In such a Mott insulating state, the charge carriers are strongly localized due to the Coulomb repulsion. Hence, strong correlations are assumed to have a strong impact on the formation of the ground state, also in the regime of small hole concentration. In a first project we have used a so-called spinon dopon mean-field theory to represent the two-dimensional Fermi-Hubbard model with strong on-site repulsive interaction in effective degrees of freedom, in which the holes can be embedded into a quantum spin liquid. The corresponding SU(2) invariant ground state belongs to the class of fractionalized Fermi liquids. In a second project we investigate a quantum dimer model, an effective model based on a Hilbert space spanned by short range singlets and bound states of holes and spins. The focus here is on the calculation of the hole-part of the electron spectral function by using exact diagonalization and its comparison with two analytic methods, a diagrammatic computation based on the Bethe-Salpeter equation and a so-called two-mode approximation. The electron spectral function shows a similar analytic form in momentum space between nodal and antinodal point when compared to results from photoemission spectroscopy experiments on cuprates. Furthermore, in a subsequent work we calculate the exact ground state of the quantum dimer model along a certain parameter line. In order to analyze the behavior of the elctron spectral function when increasing the density of holes, we investigate the Fermi-Hubbard model in a current project using a dynamical mean-field approach. To solve the 4-site cluster impurity problem, we use a numerical renormalization group approach. However, numerical limitations force us to restrict the analysis to spin-polarized baths. According to the spectral data, the system is similar to the SU(2) invariant case at half-filling in a Mott insulating state and posseses a momentum-selective energy gap at finite doping. The topology of the Fermi surface shows a Lifshitz transition when increasing the hole concentration. Here, the curvature of the Fermi surface changes from electron- to hole-like. For comparison we apply the dynamical mean-field theory also to the quantum dimer model and observe that the electron spectral functions at finite doping are qualitatively similar to that of the two-dimensional Fermi-Hubbard model. The second half of the work is about Tomonaga-Luttinger liquid theory, which is used to describe the low-energy effective degrees of freedom of one-dimensional systems. Here, we first discuss a conceptional extension of the operator-based bosonization theory for one-dimensional systems. This extension is especially suited for inhomogeneous one-dimensional systems. First, we investigate a one-dimensional system with a local interation potential and compute an exact solution of the single particle propagator at T = 0. The critical exponent of the single particle propagator has an unconventional form as a function of the microscopic Tomonaga-Luttinger parameters, which is not covered by the original Luttinger paradigm postulated by F. Duncan M. Haldane. In a second project on one-dimensional systems, we study the impact of scattering processes among bosonic low-energy excitations on the thermalization process. Such scattering processes are irrelevant on large length scales, however strongly affect the dynamics. In our analytic analysis we focus on a experimental setup, where a one-dimensional Bose gas is instantaneously splitted in two identical, however strongly correlated, halves of one-dimensional electronic systems. In the following, the corresponding non-equilibrium state runs through multiple regimes in time, such as a metastable prethermalization regime. However, above a certain threshold in time such scattering processes cause an effective thermalization of the system. In order to demonstrate this, we compute the kinetic equation in the Keldysh field integral formalism from a diagrammatic expansion based on a self-consistent Born approximation.

•Jan 1, 2020•0 citations

Sub-dimensional magnetic polarons in the one-hole doped SU(3) $t$-$J$ model

Henning Schlömer, Fabian Grusdt, Ulrich Schollwöck +2 more

The physics of doped Mott insulators is at the heart of strongly correlated materials and is believed to constitute an essential ingredient for high-temperature superconductivity. In systems with higher SU(N) spin symmetries, even richer magnetic ground states appear at a filling of one particle per site compared to the case of SU(2) spins, but their fate upon doping remains largely unexplored. Here we address this question by studying a single hole in the SU(3) $t$-$J$ model, whose undoped ground state features long-range, diagonal spin stripes. By analyzing both ground state and dynamical properties utilizing the density matrix renormalization group, we establish the appearence of magnetic polarons consisting of chargons and flavor defects, whose dynamics is constrained to a single effective dimension along the ordered diagonal. We semi-analytically describe the system using geometric string theory, where paths of hole motion are the fundamental degrees of freedom. With recent advances in the realization and control of SU(N) Fermi-Hubbard models with ultracold atoms in optical lattices, our results can directly be observed in quantum gas microscopes with single-site resolution. Our work suggests the appearance of intricate ground states at finite doping constituted by emergent, coupled Luttinger liquids along diagonals, and is a first step towards exploring a wealth of physics in doped SU(N) Fermi-Hubbard models on various geometries.

arXiv:2312.14137•Dec 21, 2023•1 citation

Hole propagation in the Kitaev-Heisenberg model: From quasiparticles in quantum Néel states to non-Fermi liquid in the Kitaev phase

Fabien Trousselet, Peter Horsch, Andrzej M. Oleś +1 more

We explore with exact diagonalization the propagation of a single hole in four magnetic phases of the $t\text{\ensuremath{-}}J$-like Kitaev-Heisenberg model on a honeycomb lattice: the N\'eel antiferromagnetic, stripe, zigzag, and Kitaev spin-liquid phases. We find coherent propagation of spin-polaron quasiparticles in the antiferromagnetic phase by a mechanism similar to that in the $t\text{\ensuremath{-}}J$ model for high-${T}_{c}$ cuprates. In the stripe and zigzag phases clear quasiparticles features appear in spectral functions of those propagators where holes are created and annihilated on one sublattice, while they remain largely hidden in those spectral functions that correspond to photoemission experiments. As the most surprising result, we find a totally incoherent spectral weight distribution for the spectral function of a hole moving in the Kitaev spin-liquid phase in the strong-coupling regime $t\ensuremath{\gg}J$ relevant for iridates. At intermediate coupling the finite systems calculation reveals a well-defined quasiparticle at the $\ensuremath{\Gamma}$ point; however, we find that the gapless spin excitations wipe out quasiparticles at finite momenta. Also for this more subtle case we conclude that in the thermodynamic limit the lightly doped Kitaev liquid phase does not support quasiparticle states in the neighborhood of $\ensuremath{\Gamma}$, and therefore yields a non-Fermi liquid, contrary to earlier suggestions based on slave-boson studies. These observations are supported by the presented study of the dynamic spin-structure factor in the Kitaev spin liquid regime.

arXiv:1308.3373•Physical Review B•Jul 10, 2014•23 citations

Polarons and dimerons in the two-dimensional attractive Hubbard model

Gerard Pascual, Jordi Boronat, Kris Van Houcke

A two-dimensional (2D) spin-up ideal Fermi gas interacting attractively with a spin-down impurity in the continuum undergoes, at zero temperature, a first-order phase transition from a polaron to a dimeron state. Here, we study a similar system on a square lattice, by considering the attractive 2D Fermi-Hubbard model with a single spin-down and a finite filling fraction of spin-up fermions. We study polaron and dimeron quasiparticle properties via variational ansatz up to one particle-hole excitation. Moreover, we develop a determinant diagrammatic Monte Carlo algorithm for this problem based on expansion in bare on-site coupling U . This algorithm turns out to be sign-problem-free at any filling of spin-up fermions, allowing one to sample very high diagram order (larger than 200 in our study) and to do simulations for large U / t (we go up to U / t = − 20 with t the hopping strength). Both methods give qualitatively consistent results. With variational ansatz we go to even larger on-site attraction. At very low spin-up filling fraction, we observe the polaron-to-dimeron transition, in agreement with the continuum case. Upon increasing the filling fraction, however, the transition shifts to higher values of | U | / t and the transition disappears beyond a filling fraction of about 20 % . In this region, the polaron state always gives a lower energy and has a finite quasiparticle residue. Our findings are directly relevant to cold atom experiments with 2D optical lattices: A small and finite density of spin-down impurities in the ground state will form a superfluid at strong coupling at low spin-up filling fraction. Above some critical spin-up filling fraction, the system is expected to remain a normal Fermi liquid, even in the strong-coupling limit.

arXiv:2411.19725•Physical Review Research•Sep 11, 2025•2 citations
Solid state physics, a very broad field in natural science, aims to explain physical properties of materials modeled as a many-body ensemble of the order of ∼ 10^23 interacting particles. For interacting metallic systems, one of the most powerful tools for unraveling the consequences of the collective interaction is the Fermi liquid theory. While still being an approximation, it excellently models many materials especially in the low-temperature and low-energy regime by treating themany-body problem as being equivalent to effectively non-interacting ”particles”, which is done by introducing the concept of quasiparticles. During the last few decades, several exotic phases of matter have been discovered, which originate from strong Coulomb interactions between the electrons, whose description is beyond the Fermi liquid theory. This thesis addresses the question of the validity of Fermi liquid theory once a single hole is introduced into an antiferromagnet with completely localized spins due to the strong interactions. Here, the most-studied case concerns the so-called 2D t − J model, which hosts a well-defined quasiparticle, the spin polaron. Focusing on a more tractable case, the t − Jz model without quantum spin fluctuations, we study the stability of the spin polaron solution once the nextnearest neighbor hopping t′ is introduced into the system. We focus on observations at the Γ-point and the X-point of systems with fixed Jz = 0.4t and varying t′ ∈ [−0.5t, 0.5t]. The model is numerically solved using exact diagonalization on finite-size lattices of up to 32 sites. A study of the quasiparticle weight under finite-size scaling reveals an “anomalous” regime at the Γ-point fort′ ≤ −0.3t, indicated by a vanishing quasiparticle weight. In order to characterize the low-energy dynamics within this “anomalous” regime, we investigated the spin correlation centered around the hole, as well as the magnon number distribution of a single-doped 20-site system. Based on the results, we present a preliminary picture of the antiferromagnetic ground state doped with one hole within that “anomalous” regime (t′ ≤ −0.3t) that deviates from the standard spin polaron as the cloud of magnetic fluctuations, typically localized around the hole, delocalizes over the whole finite-size lattice. Consequently, such systems show a weakly suppressed antiferromagnetic configuration under single-hole doping.
The understanding of matter does not only require a precise analysis of individual characteristics of the fundamental constituents, but also the study of their collective behavior.Solid state physics, and in particular research in strongly correlated electron systems, deals with such collective phenomena.The properties of the electrons, charge and spin quantum numbers, do not on their own explain the rich world of states of matter, however mutual interactions lead to various states or phases of matter.This work focuses on phases in quantum spin systems.On the basis of the Heisenberg model one (1D) and two dimensional (2D) spin lattices are investigated.In a first part the bilayer Heisenberg Model and the 2D Kondo necklace model are studied.Both models exhibit a quantum phase transition between an ordered and disordered phase.For small interlayer couplings the models have antiferromagnetic long-range order, which is destroyed by increasing the coupling.However, even in the disordered phase the models show global features in terms of collective excitations called magnons.Coming from the disordered phase the gapped magnons condense at a critical coupling and form the order of the ordered phase.Here, another interesting question arises, namely that of the coupling of a single doped hole to such critical fluctuations.A self-consistent Born approximation predicts that the doped hole couples to the magnons such that the quasiparticle residue vanishes at the quantum critical point.In this work the delicate question about the fate of the quasiparticle residue across the quantum phase transition is also tackled by means of large scale quantum Monte Carlo simulations.Furthermore the dynamics of a single hole doped in the magnetic background is investigated.
The motion of a single hole in a two-dimensional Heisenberg antiferromagnet (hAF) is studied in a representation where holes are described as spinless fermions (holons) and spins as normal bosons. Assuming long-range AF order the spin dynamics is treated in linear spin-wave theory. The formulation highlights the close relation with the conventional polaron problem. The holon Green's function is calculated self-consistently within the Born approximation using finite-cluster geometries for the numerical solution. As a remarkable result we find close agreement with the spectral function A(k,\ensuremath{\omega}) of a hole calculated by exact diagonalization methods. A(k,\ensuremath{\omega}) is characterized by a narrow quasiparticle (QP) peak at the low-energy side of the spectrum, which is well separated from the incoherent part for large enough J values. A complete characterization of our solution is given, including the spectral weight, the dispersion relation, and effective masses of the QP state. A finite-size-scaling study gives a nonvanishing spectral weight of the QP in the thermodynamic limit for values J/t typical for copper oxide superconductors. Our calculations indicate that the self-consistent Born approximation is a valuable scheme for characterizing the dynamics of a hole in the t-J model, even in the strong-coupling regime.
The Fermi liquid theory is a central concept in modern condensed matter physics used to describe conventional metals. A state of this universality class consists of well-defined quasiparticles, which occupy a finite number of states within a Fermi volume defined in momentum space. According to the Luttinger’s theorem, such a state encloses a Fermi volume proportional to the electron density modulo the filled bands if no symmetry is broken. In the past decades strong deviations from Fermi liquid theory have been observed in the scaling behavior of thermodynamic observables in different materials, among which are for instance cuprates and iron pnictides. This thesis is about two prototypical systems that can not be described by Fermi liquid theory. The first part of this work investigates two-dimensional effective models, which are in a Mott insulating state at half-filling and have a finite conductivity when doped with holes. In such a Mott insulating state, the charge carriers are strongly localized due to the Coulomb repulsion. Hence, strong correlations are assumed to have a strong impact on the formation of the ground state, also in the regime of small hole concentration. In a first project we have used a so-called spinon dopon mean-field theory to represent the two-dimensional Fermi-Hubbard model with strong on-site repulsive interaction in effective degrees of freedom, in which the holes can be embedded into a quantum spin liquid. The corresponding SU(2) invariant ground state belongs to the class of fractionalized Fermi liquids. In a second project we investigate a quantum dimer model, an effective model based on a Hilbert space spanned by short range singlets and bound states of holes and spins. The focus here is on the calculation of the hole-part of the electron spectral function by using exact diagonalization and its comparison with two analytic methods, a diagrammatic computation based on the Bethe-Salpeter equation and a so-called two-mode approximation. The electron spectral function shows a similar analytic form in momentum space between nodal and antinodal point when compared to results from photoemission spectroscopy experiments on cuprates. Furthermore, in a subsequent work we calculate the exact ground state of the quantum dimer model along a certain parameter line. In order to analyze the behavior of the elctron spectral function when increasing the density of holes, we investigate the Fermi-Hubbard model in a current project using a dynamical mean-field approach. To solve the 4-site cluster impurity problem, we use a numerical renormalization group approach. However, numerical limitations force us to restrict the analysis to spin-polarized baths. According to the spectral data, the system is similar to the SU(2) invariant case at half-filling in a Mott insulating state and posseses a momentum-selective energy gap at finite doping. The topology of the Fermi surface shows a Lifshitz transition when increasing the hole concentration. Here, the curvature of the Fermi surface changes from electron- to hole-like. For comparison we apply the dynamical mean-field theory also to the quantum dimer model and observe that the electron spectral functions at finite doping are qualitatively similar to that of the two-dimensional Fermi-Hubbard model. The second half of the work is about Tomonaga-Luttinger liquid theory, which is used to describe the low-energy effective degrees of freedom of one-dimensional systems. Here, we first discuss a conceptional extension of the operator-based bosonization theory for one-dimensional systems. This extension is especially suited for inhomogeneous one-dimensional systems. First, we investigate a one-dimensional system with a local interation potential and compute an exact solution of the single particle propagator at T = 0. The critical exponent of the single particle propagator has an unconventional form as a function of the microscopic Tomonaga-Luttinger parameters, which is not covered by the original Luttinger paradigm postulated by F. Duncan M. Haldane. In a second project on one-dimensional systems, we study the impact of scattering processes among bosonic low-energy excitations on the thermalization process. Such scattering processes are irrelevant on large length scales, however strongly affect the dynamics. In our analytic analysis we focus on a experimental setup, where a one-dimensional Bose gas is instantaneously splitted in two identical, however strongly correlated, halves of one-dimensional electronic systems. In the following, the corresponding non-equilibrium state runs through multiple regimes in time, such as a metastable prethermalization regime. However, above a certain threshold in time such scattering processes cause an effective thermalization of the system. In order to demonstrate this, we compute the kinetic equation in the Keldysh field integral formalism from a diagrammatic expansion based on a self-consistent Born approximation.
Sub-dimensional magnetic polarons in the one-hole doped SU(3) ttt-JJJ model
The physics of doped Mott insulators is at the heart of strongly correlated materials and is believed to constitute an essential ingredient for high-temperature superconductivity. In systems with higher SU(N) spin symmetries, even richer magnetic ground states appear at a filling of one particle per site compared to the case of SU(2) spins, but their fate upon doping remains largely unexplored. Here we address this question by studying a single hole in the SU(3) ttt-JJJ model, whose undoped ground state features long-range, diagonal spin stripes. By analyzing both ground state and dynamical properties utilizing the density matrix renormalization group, we establish the appearence of magnetic polarons consisting of chargons and flavor defects, whose dynamics is constrained to a single effective dimension along the ordered diagonal. We semi-analytically describe the system using geometric string theory, where paths of hole motion are the fundamental degrees of freedom. With recent advances in the realization and control of SU(N) Fermi-Hubbard models with ultracold atoms in optical lattices, our results can directly be observed in quantum gas microscopes with single-site resolution. Our work suggests the appearance of intricate ground states at finite doping constituted by emergent, coupled Luttinger liquids along diagonals, and is a first step towards exploring a wealth of physics in doped SU(N) Fermi-Hubbard models on various geometries.
We explore with exact diagonalization the propagation of a single hole in four magnetic phases of the $t\text{\ensuremath{-}}J$-like Kitaev-Heisenberg model on a honeycomb lattice: the N\'eel antiferromagnetic, stripe, zigzag, and Kitaev spin-liquid phases. We find coherent propagation of spin-polaron quasiparticles in the antiferromagnetic phase by a mechanism similar to that in the $t\text{\ensuremath{-}}J$ model for high-Tc{T}_{c}Tc​ cuprates. In the stripe and zigzag phases clear quasiparticles features appear in spectral functions of those propagators where holes are created and annihilated on one sublattice, while they remain largely hidden in those spectral functions that correspond to photoemission experiments. As the most surprising result, we find a totally incoherent spectral weight distribution for the spectral function of a hole moving in the Kitaev spin-liquid phase in the strong-coupling regime $t\ensuremath{\gg}J$ relevant for iridates. At intermediate coupling the finite systems calculation reveals a well-defined quasiparticle at the $\ensuremath{\Gamma}$ point; however, we find that the gapless spin excitations wipe out quasiparticles at finite momenta. Also for this more subtle case we conclude that in the thermodynamic limit the lightly doped Kitaev liquid phase does not support quasiparticle states in the neighborhood of $\ensuremath{\Gamma}$, and therefore yields a non-Fermi liquid, contrary to earlier suggestions based on slave-boson studies. These observations are supported by the presented study of the dynamic spin-structure factor in the Kitaev spin liquid regime.
A two-dimensional (2D) spin-up ideal Fermi gas interacting attractively with a spin-down impurity in the continuum undergoes, at zero temperature, a first-order phase transition from a polaron to a dimeron state. Here, we study a similar system on a square lattice, by considering the attractive 2D Fermi-Hubbard model with a single spin-down and a finite filling fraction of spin-up fermions. We study polaron and dimeron quasiparticle properties via variational ansatz up to one particle-hole excitation. Moreover, we develop a determinant diagrammatic Monte Carlo algorithm for this problem based on expansion in bare on-site coupling U . This algorithm turns out to be sign-problem-free at any filling of spin-up fermions, allowing one to sample very high diagram order (larger than 200 in our study) and to do simulations for large U / t (we go up to U / t = − 20 with t the hopping strength). Both methods give qualitatively consistent results. With variational ansatz we go to even larger on-site attraction. At very low spin-up filling fraction, we observe the polaron-to-dimeron transition, in agreement with the continuum case. Upon increasing the filling fraction, however, the transition shifts to higher values of | U | / t and the transition disappears beyond a filling fraction of about 20 % . In this region, the polaron state always gives a lower energy and has a finite quasiparticle residue. Our findings are directly relevant to cold atom experiments with 2D optical lattices: A small and finite density of spin-down impurities in the ground state will form a superfluid at strong coupling at low spin-up filling fraction. Above some critical spin-up filling fraction, the system is expected to remain a normal Fermi liquid, even in the strong-coupling limit.