We now analyze the physics of a single impurity immersed in the BEC and the formation of the polaron. The spectral properties of an impurity with quasi-momentum
k can be described by the impurity Green’s function
where
is the self-energy of the impurity. Due to the great complexity of many-body systems, an exact calculation of the self-energy is not feasible, for this reason, one has to make certain approximations. To calculate the self-energy, we employ the T-matrix approximation which has been successfully used to describe Bose polaron experiments [
10,
12,
56,
57,
58]. In this approximation, the self-energy of the impurity is simply the product of the matrix
and the equilibrium density of the bosons
, that is
. In contrast to the two-body problem, the scattering matrix
has to include the presence of the BEC. To do so, we assume that the BEC can be accurately described by the Bogoliubov theory. This procedure gives the chemical potential
, the excitation spectrum
, and the modified two-particle propagator
where
is the usual Bogoliubov coherence factor
Substitution of
into Equation (
4) gives the scattering matrix
. As expected, for
, the scattering matrix
is identical to the two-body scattering matrix
. In
Figure 3, we plot the spectral function of the polaron
as a function of the impurity-boson interaction
for vanishing quasi-momentum
and several values of the power
. In each panel, the yellow curve is associated with the energy of the dimer states, that is the poles of the
matrix shown in
Figure 2, the white dashed lines enclose the Bogoliubov continuum with energies
, which is essentially indistinguishable from the two-particle continuum. For the numerical calculations, we take the value
and
. The reason for considering a small boson-boson interaction is to ensure that the bosonic system is far from the transition to the Mott phase. In all numerical calculations, it was ensured that the spectral function is properly normalized
As one can notice from
Figure 3, there are well-defined quasi-particle branches for both, positive and negative interactions. The poles of
are associated with the energies of the polaron
. This energy can be found by solving the following self-consistent equation
where
denotes the real part. For
, the quasiparticle has a higher energy than the repulsive dimer and, like the latter, it is dynamically stable since there are no available states in which it can decay by converting the interaction energy into kinetic energy. As the range of the hopping increases, the energy of the impurity approaches the energy of the dimer for small
, that is, the polaron becomes a dimer into a BEC. For
, the quasiparticle branch has a smaller energy than the attractively bound pair. Remarkably, the attractive polaron is still well-defined even when the attractively bound pair is already absorbed into the continuum (see
Figure 3(d)). Furthermore, as
decreases, the broadening of the attractive polaron inside the Bogoliubov continuum is reduced. This can be physically understood from the asymmetry of the density of states (see
Figure 1(d)). As the range of the hopping increases, the density of available states into which the polaron can decay decreases, and therefore the broadening of the impurity inside the Bogoliubov continuum is reduced.
Within the quasiparticle picture and in the vicinity of a pole, the impurity Green’s function in Equation (
7) can be approximated as follows [
59]
where
is the quasiparticle residue and
is the damping rate. These quantities are given by
In
Figure 4, we show the residue and damping rate of a zero-quasi-momentum polaron as a function of the impurity-boson interaction for several values of
. The residue of the repulsive branch decreases as the range of the hopping increases, this behavior is in agreement with the previous observation that the polaron branch is pushed towards the bound state branch. Since
, it is not possible to discern this branch in current experiments. In stark contrast, the residue of the attractive branch increases when the hopping range increases, making the polaron picture more robust and feasible to be observed. As shown in
Figure 4b, the damping rate for the attractive polaron branch decreases as
decreases, that is, the polaron becomes more long-lived. The damping of the repulsive branch does not change significantly with the tunneling range.