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

ISTA Thesis

Krzysztof Myśliwy

Published Jan 1, 2022• 0 citations
OpenAlex

Abstract

The polaron model is a basic model of quantum field theory describing a single particle interacting with a bosonic field. It arises in many physical contexts. We are mostly concerned with models applicable in the context of an impurity atom in a Bose-Einstein condensate as well as the problem of electrons moving in polar crystals. The model has a simple structure in which the interaction of the particle with the field is given by a term linear in the field’s creation and annihilation operators. In this work, we investigate the properties of this model by providing rigorous estimates on various energies relevant to the problem. The estimates are obtained, for the most part, by suitable operator techniques which constitute the principal mathematical substance of the thesis. The first application of these techniques is to derive the polaron model rigorously from first principles, i.e., from a full microscopic quantum-mechanical many-body problem involving an impurity in an otherwise homogeneous system. We accomplish this for the N + 1 Bose gas in the mean-field regime by showing that a suitable polaron-type Hamiltonian arises at weak interactions as a low-energy effective theory for this problem. In the second part, we investigate rigorously the ground state of the model at fixed momentum and for large values of the coupling constant. Qualitatively, the system is expected to display a transition from the quasi-particle behavior at small momenta, where the dispersion relation is parabolic and the particle moves through the medium dragging along a cloud of phonons, to the radiative behavior at larger momenta where the polaron decelerates and emits free phonons. At the same time, in the strong coupling regime, the bosonic field is expected to behave purely classically. Accordingly, the effective mass of the polaron at strong coupling is conjectured to be asymptotically equal to the one obtained from the semiclassical counterpart of the problem, first studied by Landau and Pekar in the 1940s. For polaron models with regularized form factors and phonon dispersion relations of superfluid type, i.e., bounded below by a linear function of the wavenumbers for all phonon momenta as in the interacting Bose gas, we prove that for a large window of momenta below the radiation threshold, the energy-momentum relation at strong coupling is indeed essentially a parabola with semi-latus rectum equal to the Landau–Pekar effective mass, as expected. For the Fröhlich polaron describing electrons in polar crystals where the dispersion relation is of the optical type and the form factor is formally UV–singular due to the nature of the point charge-dipole interaction, we are able to give the corresponding upper bound. In contrast to the regular case, this requires the inclusion of the quantum fluctuations of the phonon field, which makes the problem considerably more difficult. The results are supplemented by studies on the absolute ground-state energy at strong coupling, a proof of the divergence of the effective mass with the coupling constant for a wide class of polaron models, as well as the discussion of the apparent UV singularity of the Fröhlich model and the application of the techniques used for its removal for the energy estimates.

Related papers

New theoretical approaches to Bose polarons

Fabian Grusdt, Eugene Demler

The Fröhlich polaron model describes a ubiquitous class of problems concerned with understanding properties of a single mobile particle interacting with a bosonic reservoir. Originally introduced in the context of electrons interacting with phonons in crystals, this model found applications in such diverse areas as strongly correlated electron systems, quantum information, and high energy physics. In the last few years this model has been applied to describe impurity atoms immersed in Bose-Einstein condensates of ultracold atoms. The tunability of microscopic parameters in ensembles of ultracold atoms and the rich experimental toolbox of atomic physics should allow to test many theoretical predictions and give us new insights into equilibrium and dynamical properties of polarons. In these lecture notes we provide an overview of common theoretical approaches that have been used to study BEC polarons, including Rayleigh-Schrödinger and Green's function perturbation theories, self-consistent Born approximation, mean-field approach, Feynman's variational path integral approach, Monte Carlo simulations, renormalization group calculations, and Gaussian variational ansatz. We focus on the renormalization group approach and provide details of analysis that have not been presented in earlier publications. We show that this method helps to resolve striking discrepancy in polaron energies obtained using mean-field approximation and Monte Carlo simulations. We also discuss applications of this method to the calculation of the effective mass of BEC polarons. As one experimentally relevant example of a non-equililbrium problem we consider Bloch oscillations of Bose polarons and demonstrate that one should find considerable deviations from the commonly accepted phenomenological Esaki-Tsu model. We review which parameter regimes of Bose polarons can be achieved in various atomic mixtures.

arXiv:1510.04934•Oct 16, 2015•60 citations

Bose polarons in ultracold atoms in one dimension: beyond the Fröhlich paradigm

Fabian Grusdt, Gregory E Astrakharchik, Eugene Demler

Mobile impurity atoms immersed in Bose–Einstein condensates provide a new platform for exploring Bose polarons. Recent experimental advances in the field of ultracold atoms make it possible to realize such systems with highly tunable microscopic parameters and to explore equilibrium and dynamical properties of polarons using a rich toolbox of atomic physics. In this paper we present a detailed theoretical analysis of Bose polarons in one-dimensional systems of ultracold atoms. By combining a non-perturbative renormalization group approach with numerically exact diffusion Monte Carlo calculations we obtain not only detailed numerical results over a broad range of parameters but also qualitative understanding of different regimes of the system. We find that an accurate description of Bose polarons requires the inclusion of two-phonon scattering terms which go beyond the commonly used Fröhlich model. Furthermore we show that when the Bose gas is in the strongly interacting regime, one needs to include interactions between the phonon modes. We use several theoretical approaches to calculate the polaron energy and its effective mass. The former can be measured using radio-frequency spectroscopy and the latter can be studied experimentally using impurity oscillations in a harmonic trapping potential. We compare our theoretical results for the effective mass to the experiments by Catani et al (2012 Phys. Rev. A 85 023623). In the weak-to-intermediate coupling regimes we obtain excellent quantitative agreement between theory and experiment, without any free fitting parameter. We supplement our analysis by full dynamical simulations of polaron oscillations in a shallow trapping potential. We also use our renormalization group approach to analyze the full phase diagram and identify regions that support repulsive and attractive polarons, as well as multi-particle bound states.

arXiv:1704.02606•New Journal of Physics•Sep 4, 2017•117 citations

Bose polarons in ultracold atoms in one dimension: beyond the Fröhlich paradigm

Fabian Grusdt, Astrakharchik, Grigori, Eugene Demler

Mobile impurity atoms immersed in Bose–Einstein condensates provide a new platform for exploring Bose polarons. Recent experimental advances in the field of ultracold atoms make it possible to realize such systems with highly tunable microscopic parameters and to explore equilibrium and dynamical properties of polarons using a rich toolbox of atomic physics. In this paper we present a detailed theoretical analysis of Bose polarons in one-dimensional systems of ultracold atoms. By combining a non-perturbative renormalization group approach with numerically exact diffusion Monte Carlo calculations we obtain not only detailed numerical results over a broad range of parameters but also qualitative understanding of different regimes of the system. We find that an accurate description of Bose polarons requires the inclusion of two-phonon scattering terms which go beyond the commonly used Fröhlich model. Furthermore we show that when the Bose gas is in the strongly interacting regime, one needs to include interactions between the phonon modes. We use several theoretical approaches to calculate the polaron energy and its effective mass. The former can be measured using radio-frequency spectroscopy and the latter can be studied experimentally using impurity oscillations in a harmonic trapping potential. We compare our theoretical results for the effective mass to the experiments by Catani et al (2012 Phys. Rev. A 85 023623). In the weak-to-intermediate coupling regimes we obtain excellent quantitative agreement between theory and experiment, without any free fitting parameter. We supplement our analysis by full dynamical simulations of polaron oscillations in a shallow trapping potential. We also use our renormalization group approach to analyze the full phase diagram and identify regions that support repulsive and attractive polarons, as well as multi-particle bound states.

•Jan 1, 2017•106 citations

Energy estimates for the two-dimensional Fermi polaron

Ulrich Linden

This thesis is concerned with the quantum mechanical system of a single particle interacting with an ideal gas of identical fermions by point interaction. In the physics literature this system is often referred to as Fermi polaron. We investigate the two-dimensional Fermi polaron. Unlike the one-dimensional case, point interactions in two or three dimensions cannot be implemented as perturbation of the quadratic form of the Laplacian. Either they are obtained as self-adjoint extensions of the Laplacian restricted to functions that vanish when the coordinates of two particles coincide, or they are constructed by a suitable limiting process. Choosing the second approach, a many-body operator with two-particle point interaction has firstly been rigorously defined by Dell'Antonio, Figari and Teta. We consider the Fermi polaron confined to a box with periodic boundary conditions and we identify a broad class of regularization schemes that approximate the Hamiltonian of the Fermi polaron as limit operator in the strong resolvent sense. The Hamiltonian is not given by a closed form, which could be conveniently used in standard variational principles. We establish a novel variational principle that characterizes all bound states, i.e. all energy eigenstates below the bottom of the spectrum of the kinetic energy. This variational principle turns out to be very useful for the following purposes. The ground state of the Fermi polaron is expected to be well approximated by the polaron and the molecule ansatz in the regime of weak and strong coupling between the impurity and the Fermi gas, respectively. In the physics literature, these two classes of trial states are used for variational computations with the (ultraviolet) regularized Hamiltonian. Although the implicit expressions for the minimal energy of both classes allow for the removal of the ultraviolet cutoff, it remains unclear whether the results are upper bounds to the ground state energy of the Fermi polaron. We show that the minimization of energy over polaron and molecule trial states can be reformulated in a natural way with the help of our variational principle. By doing so, the classes of trial states simplify considerably, and since there is no reference to regularized quantities, we can prove that the expressions for the polaron and the molecule energy in the physics literature are indeed upper bounds to the ground state energy of the Fermi polaron. As a further application of the variational principle, we prove analytically that to first order in a particle-hole expansion the molecule ansatz yields a better approximation to the ground state energy than the polaron ansatz if the coupling between the impurity and the Fermi gas is strong enough. So far, this had only been done numerically. The concluding chapter is devoted to the derivation of a lower bound to the ground state energy of the Fermi polaron in two-dimensional space. We show that the ground state energy can be bounded from below by a quantity that does not depend on the number of fermions in the Fermi gas. This result is correct under the assumption that the ratio of the mass of the impurity and the mass of a fermion exceeds 1.225. We also present a method which might yield a similar result for lower mass ratios. This method gives an estimate for the quadratic form of the regularized Hamiltonian in position space representation. In this connection, we present an inequality that bounds a singular potential of a two-dimensional Fermi gas depending only on the minimal distance between two fermions by the kinetic energy of the Fermi gas uniformly in the number of fermions. This inequality also applies to a potential with singularity 1/r^2, for which the Hardy inequality does not hold in two dimensions. Therefore, the full antisymmetry of the wave function has to be taken into account.

•Jan 1, 2017•5 citations

Polaron physics beyond the Holstein model

D. J. J. Marchand

Many condensed matter problems involve a particle coupled to its environment. The polaron, originally introduced to describe electrons in a polarizable medium, describes a particle coupled to a bosonic field. The Holstein polaron model, although simple, including only optical Einstein phonons and an interaction that couples them to the electron density, captures almost all of the standard polaronic properties. We herein investigate polarons that differ significantly from this behaviour. We study a model with phonon-modulated hopping, and find a radically different behaviour at strong couplings. We report a sharp transition, not a crossover, with a diverging effective mass at the critical coupling. We also look at a model with acoustic phonons, away from the perturbative limit, and again discover unusual polaron properties. Our work relies on the Bold Diagrammatic Monte Carlo (BDMC) method, which samples Feynman diagrammatic expansions efficiently, even those with weak sign problems. Proposed by Prokof'ev and Svistunov, it is extended to lattice polarons for the first time here. We also use the Momentum Average (MA) approximation, an analytical method proposed by Berciu, and find an excellent agreement with the BDMC results. A novel MA approximation able to treat dispersive phonons is also presented, along with a new exact solution for finite systems, inspired by the same formalism.

•Jan 1, 2011•1 citation

Impurities and polarons in bosonic quantum gases: a review on recent progress

Fabian Grusdt, Nader Mostaan, Eugene Demler +1 more

This review describes the field of Bose polarons, arising when mobile impurities are immersed into a bosonic quantum gas. The latter can be realized by a Bose-Einstein condensate of ultracold atoms, or of exciton polaritons in a semiconductor, which has led to a series of experimental observations of Bose polarons near inter-species Feshbach resonances that we survey. Following an introduction to the topic, with references to its historic roots and a presentation of the Bose polaron Hamiltonian, we summarize state-of-the-art experiments. Next we provide a detailed discussion of polaron models, starting from the ubiquitous Fröhlich Hamiltonian that applies at weak couplings. Already this highly simplified model allows insights into ultra-violet divergencies, logarithmic and power-law, that need to be properly regularized. To capture the physics near a Feshbach resonance, two-phonon scattering terms on the impurity as well as phonon-phonon interactions need to be included. We proceed by a survey of concurrent theoretical methods used for solving strongly interacting Bose polaron problems, ranging from Lee-Low-Pines mean-field theory, Chevy-ansatz, Gross-Pitaevskii-equation to diagrammatic Monte Carlo approaches. The subsequent sections are devoted to the large bodies of work investigating strong coupling Bose polarons, including detailed comparisons with radio-frequency spectra obtained in ultracold atom experiments; to investigations of universal few-body and Efimov states associated with a Feshbach resonance in atomic mixtures; to studies of quantum dynamics and polarons out of equilibrium; Bose polarons in low-dimensional 1D and 2D quantum systems; induced interactions among polarons and bipolaron formation; and to Bose polarons at non-zero temperatures. We end our review by detailed discussions of closely related experimental setups and systems, including ionic impurities, systems with strong light-matter interactions, and variations and extensions of the Bose polaron concepts e.g. to baths with topological order or strong interactions relevant for correlated electrons. Finally, an outlook is presented, highlighting possible future research directions and open questions in the field as a whole.

Reports on Progress in Physics•May 15, 2025•37 citations
The polaron model is a basic model of quantum field theory describing a single particle interacting with a bosonic field. It arises in many physical contexts. We are mostly concerned with models applicable in the context of an impurity atom in a Bose-Einstein condensate as well as the problem of electrons moving in polar crystals. The model has a simple structure in which the interaction of the particle with the field is given by a term linear in the field’s creation and annihilation operators. In this work, we investigate the properties of this model by providing rigorous estimates on various energies relevant to the problem. The estimates are obtained, for the most part, by suitable operator techniques which constitute the principal mathematical substance of the thesis. The first application of these techniques is to derive the polaron model rigorously from first principles, i.e., from a full microscopic quantum-mechanical many-body problem involving an impurity in an otherwise homogeneous system. We accomplish this for the N + 1 Bose gas in the mean-field regime by showing that a suitable polaron-type Hamiltonian arises at weak interactions as a low-energy effective theory for this problem. In the second part, we investigate rigorously the ground state of the model at fixed momentum and for large values of the coupling constant. Qualitatively, the system is expected to display a transition from the quasi-particle behavior at small momenta, where the dispersion relation is parabolic and the particle moves through the medium dragging along a cloud of phonons, to the radiative behavior at larger momenta where the polaron decelerates and emits free phonons. At the same time, in the strong coupling regime, the bosonic field is expected to behave purely classically. Accordingly, the effective mass of the polaron at strong coupling is conjectured to be asymptotically equal to the one obtained from the semiclassical counterpart of the problem, first studied by Landau and Pekar in the 1940s. For polaron models with regularized form factors and phonon dispersion relations of superfluid type, i.e., bounded below by a linear function of the wavenumbers for all phonon momenta as in the interacting Bose gas, we prove that for a large window of momenta below the radiation threshold, the energy-momentum relation at strong coupling is indeed essentially a parabola with semi-latus rectum equal to the Landau–Pekar effective mass, as expected. For the Fröhlich polaron describing electrons in polar crystals where the dispersion relation is of the optical type and the form factor is formally UV–singular due to the nature of the point charge-dipole interaction, we are able to give the corresponding upper bound. In contrast to the regular case, this requires the inclusion of the quantum fluctuations of the phonon field, which makes the problem considerably more difficult. The results are supplemented by studies on the absolute ground-state energy at strong coupling, a proof of the divergence of the effective mass with the coupling constant for a wide class of polaron models, as well as the discussion of the apparent UV singularity of the Fröhlich model and the application of the techniques used for its removal for the energy estimates.
The Fröhlich polaron model describes a ubiquitous class of problems concerned with understanding properties of a single mobile particle interacting with a bosonic reservoir. Originally introduced in the context of electrons interacting with phonons in crystals, this model found applications in such diverse areas as strongly correlated electron systems, quantum information, and high energy physics. In the last few years this model has been applied to describe impurity atoms immersed in Bose-Einstein condensates of ultracold atoms. The tunability of microscopic parameters in ensembles of ultracold atoms and the rich experimental toolbox of atomic physics should allow to test many theoretical predictions and give us new insights into equilibrium and dynamical properties of polarons. In these lecture notes we provide an overview of common theoretical approaches that have been used to study BEC polarons, including Rayleigh-Schrödinger and Green's function perturbation theories, self-consistent Born approximation, mean-field approach, Feynman's variational path integral approach, Monte Carlo simulations, renormalization group calculations, and Gaussian variational ansatz. We focus on the renormalization group approach and provide details of analysis that have not been presented in earlier publications. We show that this method helps to resolve striking discrepancy in polaron energies obtained using mean-field approximation and Monte Carlo simulations. We also discuss applications of this method to the calculation of the effective mass of BEC polarons. As one experimentally relevant example of a non-equililbrium problem we consider Bloch oscillations of Bose polarons and demonstrate that one should find considerable deviations from the commonly accepted phenomenological Esaki-Tsu model. We review which parameter regimes of Bose polarons can be achieved in various atomic mixtures.
Mobile impurity atoms immersed in Bose–Einstein condensates provide a new platform for exploring Bose polarons. Recent experimental advances in the field of ultracold atoms make it possible to realize such systems with highly tunable microscopic parameters and to explore equilibrium and dynamical properties of polarons using a rich toolbox of atomic physics. In this paper we present a detailed theoretical analysis of Bose polarons in one-dimensional systems of ultracold atoms. By combining a non-perturbative renormalization group approach with numerically exact diffusion Monte Carlo calculations we obtain not only detailed numerical results over a broad range of parameters but also qualitative understanding of different regimes of the system. We find that an accurate description of Bose polarons requires the inclusion of two-phonon scattering terms which go beyond the commonly used Fröhlich model. Furthermore we show that when the Bose gas is in the strongly interacting regime, one needs to include interactions between the phonon modes. We use several theoretical approaches to calculate the polaron energy and its effective mass. The former can be measured using radio-frequency spectroscopy and the latter can be studied experimentally using impurity oscillations in a harmonic trapping potential. We compare our theoretical results for the effective mass to the experiments by Catani et al (2012 Phys. Rev. A 85 023623). In the weak-to-intermediate coupling regimes we obtain excellent quantitative agreement between theory and experiment, without any free fitting parameter. We supplement our analysis by full dynamical simulations of polaron oscillations in a shallow trapping potential. We also use our renormalization group approach to analyze the full phase diagram and identify regions that support repulsive and attractive polarons, as well as multi-particle bound states.
Mobile impurity atoms immersed in Bose–Einstein condensates provide a new platform for exploring Bose polarons. Recent experimental advances in the field of ultracold atoms make it possible to realize such systems with highly tunable microscopic parameters and to explore equilibrium and dynamical properties of polarons using a rich toolbox of atomic physics. In this paper we present a detailed theoretical analysis of Bose polarons in one-dimensional systems of ultracold atoms. By combining a non-perturbative renormalization group approach with numerically exact diffusion Monte Carlo calculations we obtain not only detailed numerical results over a broad range of parameters but also qualitative understanding of different regimes of the system. We find that an accurate description of Bose polarons requires the inclusion of two-phonon scattering terms which go beyond the commonly used Fröhlich model. Furthermore we show that when the Bose gas is in the strongly interacting regime, one needs to include interactions between the phonon modes. We use several theoretical approaches to calculate the polaron energy and its effective mass. The former can be measured using radio-frequency spectroscopy and the latter can be studied experimentally using impurity oscillations in a harmonic trapping potential. We compare our theoretical results for the effective mass to the experiments by Catani et al (2012 Phys. Rev. A 85 023623). In the weak-to-intermediate coupling regimes we obtain excellent quantitative agreement between theory and experiment, without any free fitting parameter. We supplement our analysis by full dynamical simulations of polaron oscillations in a shallow trapping potential. We also use our renormalization group approach to analyze the full phase diagram and identify regions that support repulsive and attractive polarons, as well as multi-particle bound states.
This thesis is concerned with the quantum mechanical system of a single particle interacting with an ideal gas of identical fermions by point interaction. In the physics literature this system is often referred to as Fermi polaron. We investigate the two-dimensional Fermi polaron. Unlike the one-dimensional case, point interactions in two or three dimensions cannot be implemented as perturbation of the quadratic form of the Laplacian. Either they are obtained as self-adjoint extensions of the Laplacian restricted to functions that vanish when the coordinates of two particles coincide, or they are constructed by a suitable limiting process. Choosing the second approach, a many-body operator with two-particle point interaction has firstly been rigorously defined by Dell'Antonio, Figari and Teta. We consider the Fermi polaron confined to a box with periodic boundary conditions and we identify a broad class of regularization schemes that approximate the Hamiltonian of the Fermi polaron as limit operator in the strong resolvent sense. The Hamiltonian is not given by a closed form, which could be conveniently used in standard variational principles. We establish a novel variational principle that characterizes all bound states, i.e. all energy eigenstates below the bottom of the spectrum of the kinetic energy. This variational principle turns out to be very useful for the following purposes. The ground state of the Fermi polaron is expected to be well approximated by the polaron and the molecule ansatz in the regime of weak and strong coupling between the impurity and the Fermi gas, respectively. In the physics literature, these two classes of trial states are used for variational computations with the (ultraviolet) regularized Hamiltonian. Although the implicit expressions for the minimal energy of both classes allow for the removal of the ultraviolet cutoff, it remains unclear whether the results are upper bounds to the ground state energy of the Fermi polaron. We show that the minimization of energy over polaron and molecule trial states can be reformulated in a natural way with the help of our variational principle. By doing so, the classes of trial states simplify considerably, and since there is no reference to regularized quantities, we can prove that the expressions for the polaron and the molecule energy in the physics literature are indeed upper bounds to the ground state energy of the Fermi polaron. As a further application of the variational principle, we prove analytically that to first order in a particle-hole expansion the molecule ansatz yields a better approximation to the ground state energy than the polaron ansatz if the coupling between the impurity and the Fermi gas is strong enough. So far, this had only been done numerically. The concluding chapter is devoted to the derivation of a lower bound to the ground state energy of the Fermi polaron in two-dimensional space. We show that the ground state energy can be bounded from below by a quantity that does not depend on the number of fermions in the Fermi gas. This result is correct under the assumption that the ratio of the mass of the impurity and the mass of a fermion exceeds 1.225. We also present a method which might yield a similar result for lower mass ratios. This method gives an estimate for the quadratic form of the regularized Hamiltonian in position space representation. In this connection, we present an inequality that bounds a singular potential of a two-dimensional Fermi gas depending only on the minimal distance between two fermions by the kinetic energy of the Fermi gas uniformly in the number of fermions. This inequality also applies to a potential with singularity 1/r^2, for which the Hardy inequality does not hold in two dimensions. Therefore, the full antisymmetry of the wave function has to be taken into account.
Many condensed matter problems involve a particle coupled to its environment. The polaron, originally introduced to describe electrons in a polarizable medium, describes a particle coupled to a bosonic field. The Holstein polaron model, although simple, including only optical Einstein phonons and an interaction that couples them to the electron density, captures almost all of the standard polaronic properties. We herein investigate polarons that differ significantly from this behaviour. We study a model with phonon-modulated hopping, and find a radically different behaviour at strong couplings. We report a sharp transition, not a crossover, with a diverging effective mass at the critical coupling. We also look at a model with acoustic phonons, away from the perturbative limit, and again discover unusual polaron properties. Our work relies on the Bold Diagrammatic Monte Carlo (BDMC) method, which samples Feynman diagrammatic expansions efficiently, even those with weak sign problems. Proposed by Prokof'ev and Svistunov, it is extended to lattice polarons for the first time here. We also use the Momentum Average (MA) approximation, an analytical method proposed by Berciu, and find an excellent agreement with the BDMC results. A novel MA approximation able to treat dispersive phonons is also presented, along with a new exact solution for finite systems, inspired by the same formalism.
This review describes the field of Bose polarons, arising when mobile impurities are immersed into a bosonic quantum gas. The latter can be realized by a Bose-Einstein condensate of ultracold atoms, or of exciton polaritons in a semiconductor, which has led to a series of experimental observations of Bose polarons near inter-species Feshbach resonances that we survey. Following an introduction to the topic, with references to its historic roots and a presentation of the Bose polaron Hamiltonian, we summarize state-of-the-art experiments. Next we provide a detailed discussion of polaron models, starting from the ubiquitous Fröhlich Hamiltonian that applies at weak couplings. Already this highly simplified model allows insights into ultra-violet divergencies, logarithmic and power-law, that need to be properly regularized. To capture the physics near a Feshbach resonance, two-phonon scattering terms on the impurity as well as phonon-phonon interactions need to be included. We proceed by a survey of concurrent theoretical methods used for solving strongly interacting Bose polaron problems, ranging from Lee-Low-Pines mean-field theory, Chevy-ansatz, Gross-Pitaevskii-equation to diagrammatic Monte Carlo approaches. The subsequent sections are devoted to the large bodies of work investigating strong coupling Bose polarons, including detailed comparisons with radio-frequency spectra obtained in ultracold atom experiments; to investigations of universal few-body and Efimov states associated with a Feshbach resonance in atomic mixtures; to studies of quantum dynamics and polarons out of equilibrium; Bose polarons in low-dimensional 1D and 2D quantum systems; induced interactions among polarons and bipolaron formation; and to Bose polarons at non-zero temperatures. We end our review by detailed discussions of closely related experimental setups and systems, including ionic impurities, systems with strong light-matter interactions, and variations and extensions of the Bose polaron concepts e.g. to baths with topological order or strong interactions relevant for correlated electrons. Finally, an outlook is presented, highlighting possible future research directions and open questions in the field as a whole.