SciPost Submission Page
Nonequilibrium dynamics of dipolar polarons
by Artem G. Volosniev, Giacomo Bighin, Luis Santos, Luis A. Peña Ardila
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Luis Peña Ardila · Artem Volosniev 
Submission information  

Preprint Link:  scipost_202310_00016v1 (pdf) 
Date accepted:  20231122 
Date submitted:  20231017 13:52 
Submitted by:  Peña Ardila, Luis 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approach:  Theoretical 
Abstract
We study the outofequilibrium quantum dynamics of dipolar polarons, i.e., impurities immersed in a dipolar BoseEinstein condensate, after a quench of the impurityboson interaction. We show that the dipolar nature of the condensate and of the impurity results in anisotropic relaxation dynamics, in particular, anisotropic dressing of the polaron. More relevantly for coldatom setups, quench dynamics is strongly affected by the interplay between dipolar anisotropy and trap geometry. Our findings pave the way for simulating impurities in anisotropic media utilizing experiments with dipolar mixtures.
Author comments upon resubmission
Thank you for handling our submission and for the invitation to revise and resubmit the manuscript.
Both Referees found that our work can be suitable for publication provided we address their comments. We revise our manuscript following the Referee's reports and provide a pointbypoint reply.
For convenience, we highlight the main changes in the manuscript. Minor changes (correction of typos and style) are not highlighted.
We hope that our manuscript is now ready for publication in SciPost Physics.
Sincerely,
The Authors
List of changes
\section*{\bf Reply to Reviewer 1}
{\bf Reviewer }: This manuscript presents a perturbation study of the polaron physics of a weakly interacting impurity immersed in a dipolar BoseEinstein condensate. The authors employ the timedependent variational theory on the polaron Hamiltonian in the moving frame of the polaron, obtained through the LeeLowPines transformation. The paper explores various aspects of the polaron, such as steadystate energy, residue, effective mass, and relaxation dynamics. Additionally, the authors investigate the trapped system using a localdensity approximation. Overall, the topic is interesting, and the manuscript is wellwritten and provides clear motivations, explanations, and detailed analysis. However, there are some concerns regarding the formalism that should be addressed before making a final recommendation.
{\bf Our reply:} We thank the Referee for taking the time to review our work, for a positive evaluation of our results and for suggestions that helped us to improve the manuscript considerably. We think we have addressed all of the concerns listed in the report and hope that the improved version of the manuscript is suitable for publication.
{\bf Reviewer }: Derivation of Eqs. (2):
It would be beneficial to include a more detailed derivation of Eqs. (2) from the timedependent variational ansatz that minimizes the action functional. Specifically, it would be helpful to explain why the second term on the righthand side of the second equation in Eqs. (2) is not expressed as $(P  P_B)^2/2m$, given the naive expectation that it should resemble the Hamiltonian in Eq. (S4).
{\bf Our reply:} We thank the Referee for this comment. In the revised manuscript, we have added a detailed derivation of Eq. (2), see new App. C.
{ \bf Reviewer:}
Neglect of $P_B$ in the perturbation analysis:
The authors state that in the LeeLowPines transformation, $P_B$ (which is related to the recoil of the impurity by the phonons) is neglected in their perturbation analysis because it only enters in the nexttoleading order. However, neglecting $P_B$ implies an effective Hamiltonian equivalent to that of an infinitely heavy impurity. This contradicts the assumption of equal masses for the impurity and background bosons, where recoil should play a significant role. It is necessary to provide a more detailed discussion justifying the neglect of $P_B$ and its underlying physical picture.
{\bf Our reply:} We respectfully disagree with the Referee that the limit $\mathbf{P}_{\mathrm{B}}=0$ corresponds to an infinitely heavy impurity. The mass of the impurity enters in our equations even after we set $\mathbf{P}_{\mathrm{B}}=0$, see Eq. (2), which contains the parameter $m$.
We agree with the Referee that the bosonimpurity scattering plays an important role. In particular, it leads to renormalization of the impurity quasiparticle properties, such as energy and mass, which we calculate and benchmark against available results in the literature. It also affects dynamic properties as we discuss in the text.
To further disentangle the limit $\mathbf{P}_{\mathrm{B}}=0$ from the limit of an infinitely heavy impurity, please note that for the ground state with $\mathbf{P}=0$ it is natural to expect (as there is no vector to set a preferred direction) that $\mathbf{P}_{\mathrm{B}}=0$ for any mass of the impurity. However, even in this case, the mass of the impurity enters our calculations as it determines the energymomentum conservation laws for bosonimpurity scattering. In the revised manuscript, we added a clarifying discussion; see the last paragraph in new App.~C and new Eq. (8) in the main text.
{\bf Reviewer :}
Selftrapping effect and earlier works:
Considering that the authors focus on weakly interacting impurities, it would be valuable to mention early works on ultracold polarons that consider the selftrapping effect. The selftrapping effect leads to the spontaneous breaking of spatial translational symmetry, and the momentum P is no longer a good quantum number. Although the meanfield description fails for strongly interacting impurities, it is likely that the selftrapping effect remains relevant for weakly interacting impurities, at least for static state analysis. (I would expect that for dynamical analysis, this is less relevant, as the weaklyinteracting BEC is a slowlyresponded fluid.) This effect becomes particularly interesting in systems with longrange interactions. It would be appropriate to reference earlier works such as [Phys. Rev. Lett. 96, 210401 (2006), Phys. Rev. A 73, 043608 (2006)] and relevant literature in this context.
{\bf Our reply:} We thank the Referee for this comment. Indeed, we have not discussed selflocalization of the impurity in our work. In the revised manuscript, we argue that selflocalization is not relevant for the dilute gas that we consider, see the discussion at the end of App. B. The main reason for that is that the gas parameter, $n a_{11}^3$, is very small for our parameters.
{\bf Reviewer:}
Some Minor Suggestions:
1. Discuss the underlying physics of relaxation:
Since there are no decay channels mentioned, it would be helpful for the authors to comment on the underlying physics of the relaxation process. It could be speculated that the relaxation is caused by some dephasing process. Providing insights into this aspect would enhance the understanding of the dynamics studied.
2. Perturbation argument in droplet or supersolid backgrounds:
While the authors claim that the results can be extrapolated easily to droplet or supersolid backgrounds, it is important to note that the density in such regimes might be significantly higher than in a BoseEinstein condensate. Consequently, the validity of the perturbation argument may no longer hold. It would be fair to mention this limitation to the readers in the Outlook section.
{\bf Our reply:} 1. We thank the Referee for this suggestion. As the Referee correctly pointed out, the underlying mechanism of relaxation is a dephasing process. The systems relaxes even in the absence of bosonboson interactions because the rapid change of the impurityboson interaction populates a continuum of energy states, which then dephase due to the energy mismatch, similar to Refs. [21, 28] of our manuscript.
Importantly, in real experiments, other dephasing relaxation processes are possible. They can
originate due to trap dephasing, decoherence by magneticfield fluctuations, threebody losses and thermal fluctuations. To compare to our results, these decoherence effects should be under control in coldatom experiments (cf.~Ref.~[10] of the revised manuscript).\\
We have added a clarifying discussion, see footnote 2 of the revised manuscript.\\
2. We agree with the Referee that the density of quantum droplets is higher than the one we consider in the paper.
Quantum fluctuations or in our language, the beyondmeanfield (LeeHuangYang) contribution is the main mechanism for the stabilization of quantum droplets. This term may become important for densities one order of magnitude larger than those typical for BECs (namely $\sim10^{21}/\mathrm{m}^{3}$). {\it However}, it should be noted that perturbative approaches such as Gross–Pitaevskiilike equations have been widely used and their accuracy have been confirmed by comparing to experimental data in the droplet and supersolid regimes. See, for example, the reference R.N. Bisset, et.al., Phys. Rev. Lett. 126, 025301 (2021), where the polaron energy was obtained from secondorder perturbation theory and benchmarked against the chemical potential of the impurity immersed in a dipolar quantum droplet.
{\bf Reviewer:} Overall, the manuscript is wellwritten and addresses an interesting topic. By addressing the concerns raised and incorporating the suggestions provided, the authors can further improve the paper and make it suitable for publication.
{\bf Our reply:} We thank these encouraging remarks on our work. We hope that the revised version of our manuscript is ready for publication in SciPost Physics.
\section*{\bf Reply to Reviewer 2}
{\bf Reviewer }:
In this work, the authors study the nonequilibrium dynamics of dipolar Bose polarons after an interaction quench of the impurityboson coupling. Certain dynamical features of this system are calculated using a timedependent variational approach and exploiting the LeeLowPines transformation which allows also for analytical estimates and for distinguishing the relevant timescales. They extract several experimentally relevant observables such as the Ramsey contrast and quasiparticle effective mass. The main result is the manifestation of the polaron anisotropic relaxation dynamics due to the dipolar interactions. Trap effects are also explored within the local density approximation.
{\bf Our reply:} We thank the Referee for taking the time to review our work.
{\bf Reviewer }:
I find these results interesting towards the understanding of polarons in anisotropic media. Moreover, the submission is timely and the setup will be accessible in near future cold atom experiments. The manuscript is also well written. However, I have some comments and questions regarding the presentation and the interpretation of the presented results. If the authors are able to address these concerns which I provide below, the work will be further advanced and be suitable for publication in SciPost Physics. Suggestions for improvement follow
{\bf Our reply:} We thank the Referee for this positive evaluation of our work, and for comments that helped us to improve the manuscript. We hope that the Referee will find that the revised manuscript addresses all of the comments listed in the report, making our work suitable for publication in SciPost Physics.
{\bf Reviewer}:
In the introduction, the statement “collective excitations build the subsequent time evolution” is not clear. What type of collective excitations are meant and in which context? I guess the fewbody correlations should also play a role at long evolution times invalidating the meanfield approximation? Please elaborate at least briefly on this point.
{\bf Our reply:} We thank the Referee for identifying this weak statement. To avoid any confusion, we have rewritten this sentence. In particular, we have decided to refrain from discussing collective excitations that early in the paper. The goal of our introduction is to provide a general motivation for studying nonequilibrium dynamics in general setting and various collective excitations. \\
To explain the importance of few and manybody correlations, we now add a discussion to the beginning of the third paragraph of the Introduction: \textit{Complexity of quench dynamics of impurities stems in particular from the interplay between
few and manybody correlations and the corresponding timescales. Initial dynamics involves
a local excitation of the system, which implies participation of phonons with high momenta that can resolve fewbody physics. After some time, the dynamics is mainly
driven by the propagation of lowenergy manybody excitations}. Of course, the Referee is right that the system always has few and manybody correlations. By discussing few and manybody correlations, we want to stress that the initial time dynamics can be modeled using mainly fewbody physics, while manybody correlations become more important at the later stages of time evolution.
{\bf Reviewer}:
In the introduction the term shorttime polaron physics is stated. This becomes particularly confusing since a sentence later the term initial short time dynamics is used. Please provide an order of magnitude of what is actually meant.
{\bf Our reply:}
We thank the Referee for identifying this possibly confusing issue. In fact, shorttime dynamics and initial shorttime are used in the same context. In order to avoid any confusion, we have rewritten this part of the introduction. In particular, we have decided to avoid introducing timescales in the introduction. They are now introduced later in the main text.\\
{\bf Reviewer }:
In the same paragraph, it would be beneficial for the reader to already explain in a sentence the context of the Ramsey contrast. Otherwise, the comment provided at the end of the Introduction should be moved to this point and explain which interactions are meant in the overlap.
{\bf Our reply:} Following the recommendation of the Referee, we have added a discussion about the Ramseytype scheme to the introduction: {\it In particular, a combination of a RF pulse that coherently populates the impurity state and Ramsey
interferometry, which tracks down the quantum evolution of the impurity coupled coherently to
the BEC, allows one to measure a timedependent Green’s function, which contains information about
nonequilibrium dynamics of the impurity [9, 10]. Its imaginary part gives the contrast, namely the
overlap between the initial state and the polaron state, which is one of the key observables studied
in this work.}
{\bf Reviewer }: Since the authors provide a list of possible techniques that have been used to describe polaron physics it would be fair to also comment on other variational methods that have been extensively employed.
{\bf Our reply:} In the revised manuscript, we have added a number of references in the introduction and throughout the text. In particular, the revised version of the introduction now mentions the Langevin equation, the manybody generalization of WeisskopfWigner theory, and a few papers for an impurity in a Fermi gas. In general, following the recommendation of the Referees, we have extended our bibliography significantly by adding a number of relevant works.
{\bf Reviewer}:
I would also encourage the authors to describe their main specific results in a more clear manner in the Introduction.
{\bf Our reply:} Following this recommendation of the Referee, we have added a corresponding discussion as the last two paragraphs of the introduction.
{\bf Reviewer}:
I wonder whether in the single impurity limit that is considered in Section 2, dipoledipole interactions between the impurity and the bath are indeed experimentally justified. Please comment.
{\bf Our reply:} In experiments, a small number of impurities are created on top of the condensate. In general, the impurity concentration is lower than ten percent. The direct and induced interactions between impurities are (almost) negligible compared to bosonimpurity interactions. Thus, the ultradilute gas of impurities is well described within the framework of a single impurity. For instance, in dipolar mixtures such as $\mathrm{DyDy}$, $\mathrm{ErDy}$ or $\mathrm{ErEr}$, the impurityboson dipolar interaction is proportional to the dipolar length, $d_2$. It is always present, and it prominently enters our analytical expressions.
{\bf Reviewer }:
The fact that the relaxation timescales are less for a headtototail arrangement should be explained. Is it related to the attractive nature of the interaction potential and how? Please comment. The term “relaxation time” should also be explicitly explained since it has different interpretations in the literature.
{\bf Our reply:}
The referee is right that the relaxation timescales are shorter for a headtototail arrangement. In order to explain the previous statement, we added a few sentences in the main text and in the appendix.\\
In the main text, we added:\\
\textit{In the absence of dipolar interactions, polaron formation time, which marks the onset of the final stage of the impurity dynamics, is given by $\xi/c_0$~[25], where $\xi=\hbar/(\sqrt{2}c_0)$ is the coherence length of a nondipolar condensate. If we naively define a directiondependent timescle as $\hbar^2/(\sqrt{2}c(\theta)^2)$ (see App.~D for a derivation of this timescale), then we shall expect slower dynamics in the sidetoside ($\theta=\pi/2$) direction, which features the softest excitations in the system~[50], than in the headtotail ($\theta=0$) one.
This expectation is directly confirmed by our calculations, see Fig.~1 (b).}\\
In App. D, after Eq.~[S32], we added:\\
\textit{For a nondipolar condensate, $t_n$ is the time required for a phonon to travel through the region of the condensate distorted by the impurity. Indeed, $t_n\simeq \xi/c$, where the healing length $\xi$ defines the size of the distortion due the presence of the impurity. For a dipolar condensate, the anisotropy of time evolution is encoded in the prefactor $\bar{\Omega}_{x}$, which effectively leads to a directiondependent speed of sound.}
{\bf Reviewer }:
It would be interesting to comment whether quantum fluctuations (in terms of second order perturbation) play any role in the static (e.g. ground state) polaron energy for the considered interactions. In other words, are there any deviations from the mean field contribution?
{\bf Our reply:} Following the recommendation of the Referee, we have added an estimate of the contribution to the groundstate energy from secondorder perturbation theory. This estimate is about 10\%.
{\bf Reviewer }:
In the contrast defined in Equation (3), I guess the initial wave function refers to the situation where all interactions are switchedoff. Is this correct? Please comment explicitly.
{\bf Our reply:} Indeed, the initial state corresponds to a system where the bosonimpurity interaction is turned off, and the total momentum is vanishing.
We have added a corresponding explanation to the revised manuscript, see Sec. 4.2. It is worth noting also that the most relevant case for current experiments is to consider a nondipolar impurity immersed in a dipolar bath. In this case, only the swave impurityboson scattering length is switched on.
{\bf Reviewer}:
What is the physical context and the formula of the timescale $t_n$ appearing in Equation (4)? How is it compared to $t_\Omega$?
{\bf Our reply:} The parameter $t_n=\frac{m}{8\pi n\hbar a_{11}}$ defines the timescale for the Bose gas at the meanfield level. Following the recommendation of the Referee, we add this explicitly to the manuscript, see the discussion around Eq.~4. For the considered parameters, $t_n\sim 0.1$ms, whereas $t_\Omega\sim 0.2$s, see the caption to Fig.~2. It is worthnoting that the interplay between the parameters $t_\Omega$ and $t_n$ is seen throughout our work, see for example Eq. (5) and Eq. (14) for an experimentally relevant trapped case.
{\bf Reviewer}:
In dipolar systems it is known that loss channels such as threebody recombination play an important role. Can the authors comment on the effect of such processes in their system.
{\bf Our reply:}
We agree with the Referee that
threebody recombination processes play an important role in bosonic systems. However, note that the system under consideration is very weakly interacting, meaning that it survives in a trap for long times that can be of the order of 100ms and even more, (see, e.g., Phys. Rev. Lett. 121(21), 213601 (2018)). As typical time scales for polaron formation are less than 1ms, we do not expect threebody processes to play an important role in our work.
To further address the comment of the Referee, we note that the processes involving recombination of one impurity and two bosons from the condensate yield a distribution of atoms lost from the trap, which is proportional to the number of atoms transferred into the impurity state. These effects can be indeed very important for large interspecies interaction. For example, they are probably relevant for $^{166}$Er$^{162}$Dy mixtures (see, e.g., Phys. Rev. Lett. 121(21), 213601 (2018)). As in the manuscript, we focus on the scattering lengths $\simeq d$, we can neglect these dynamics here.
{\bf Reviewer}:
The decay of the contrast shown in Figure is relatively weak. Would it be possible to be observed in the experiment? Can the authors provide their opinion on this issue?
{\bf Our reply:} We agree with the Referee that the change of the contrast is weak in the uniform case. However, in real experiments, trap dephasing is essential for time dynamics, and therefore, the situation presented in Fig.~3 of our manuscript cannot be observed in current experiments. Note that in the trapped case, the contrast ranges from 0 to 1 and it is accessible via interferometry from very short time scales up to the steady state using already available experimental techniques.
{\bf Reviewer}:
Does the system exhibits any universal characteristics for long evolution times where e.g. the contract is saturated?
{\bf Our reply:} In Eq.~(9), we present an expression for the contrast for longer times. The expression depends on the compressibility of the bath encoded in the gas parameter $na_{11}^{3}$, the dipolar coupling strengths and the swave impurity boson scattering length. Therefore, it seems that there is no universality, as one cannot toss out any of the relevant length scales. Similarly to the nondipolar case, universality is only unveiled for shorttime dynamics, see Ref.~[10] of our work
{\bf Reviewer }:
I expect that depending on the $\epsilon_{dd}$ parameter the phase of the bath is modified. How is this fact be reflected in the contrast?
{\bf Our reply:} Indeed, by changing the parameter $\epsilon_{dd}$ one can modify the phase of the host bath. However, as we focus in our work on the regime with $\epsilon_{dd}<1$, the system is still a BEC with a welldefined global phase. This phase is insensitive to the presence of the impurity. We leave a more detailed investigation on the role of the phase of the bath to future studies.
{\bf Reviewer }:
I would expect that the diffusion process is accompanied by energy redistribution between the impurity and the medium. Can the authors comment on this issue?
{\bf Our reply:} It is true that quench dynamics is accompanied by energy redistribution between the impurity and the medium. To clarify the dynamics of the energy we have modified the text, see the last part of new App. C and new Eq. (8).
{\bf Reviewer }:
Does the fact that the impurity slowdown occurs faster for $\epsilon_{dd} \to1$ is related to the superfluid character of the medium or not? Please clarify.
{\bf Our reply:} The notion of `critical slowdown' requires a superfluid character of the medium  it defines the behavior of the system when the momentum of the impurity is close to the critical momentum, see Ref. [25] of our manuscript.
When $\epsilon_{dd}\rightarrow1$, the speed of sound in one of the directions vanishes, which means that the slowdown occurs for smaller values of the momentum of the impurity.
We have clarified this discussion in the revised manuscript.
{\bf Reviewer}:
I am not sure that the trap frequencies provided before Equation (11) have been used in dipolar condensates. If they have been used please provide relevant references.
{\bf Our reply:} Although, we have chosen these values of frequencies rather arbitrarily, they are not far from what have been seen with cold dipolar atoms (see, e.g., Phys. Rev. Lett. 121(21),
213601 (2018)). In general, we do not foresee any difficulty in achieving experimentally the proposed numbers due to a high tunability of harmonic trapping in coldatom laboratories.
{\bf Reviewer }:
In the outlook section, I would suggest adding at least a brief description of the new results presented in the main text. This will certainly facilitate the understanding of the reader.
{\bf Our reply:} Following this recommendation of the Referee, we have briefly described the main results of the paper in the last section, see the beginning of Sec. 4.6 of the revised manuscript.
We thank the Referee for the comments and hope that the revised version of our manuscript is ready for publication in SciPost Physics.
Published as SciPost Phys. 15, 232 (2023)
Reports on this Submission
Report #2 by Anonymous (Referee 5) on 20231117 (Invited Report)
 Cite as: Anonymous, Report on arXiv:scipost_202310_00016v1, delivered 20231117, doi: 10.21468/SciPost.Report.8135
Report
The authors have addressed the majority of the referees' comments with commendable improvements. However, I would like to clarify my remarks regarding the neglect of P_B in the perturbation analysis. In my comment, I contend that when P_B is zero, the momentum of the polaron aligns with the impurity's momentum. Consequently, in the moving frame, the impurity appears stationary, even though the background atoms recoil from it. This scenario is analogous to having an infinitely heavy impurity in the moving frame, seemingly contradicting the impurity's mass, which still plays a role in the lab frame, as emphasized by the authors.
In response, the authors highlighted that P and P_B can simultaneously be zero. This should be understood as the expectation value of P_B being zero, while P_B, as a general operator, should not. As stated in [PHYSICAL REVIEW X 11, 041015 (2021)], one of the crucial features of LLP transformation is to establish "sufficient entanglement between the impurity and the medium, allowing the HartreeFock approximation to become accurate." The introduction of entanglement between the impurity and the medium by P_B is now recognized as a necessary improvement from the early GPElike treatment of polarons [e.g., PHYSICAL REVIEW A 73, 063604 (2006)]. However, for the scope of this study, the neglect of P_B might still be justifiable, especially if the authors are investigating a regime where the GPElike equation is sufficiently accurate.
Nevertheless, these technical details should not impede the publication of this manuscript, in my opinion. I strongly recommend it for publication.