SciPost Submission Page
The origin of the period$2T/7$ quasibreathing in diskshaped GrossPitaevskii breathers
by J. Torrents, V. Dunjko, M. Gonchenko, G. E. Astrakharchik, M. Olshanii
This is not the latest submitted version.
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  G. E. Astrakharchik · Vanja Dunjko · Maxim Olshanii 
Submission information  

Preprint Link:  scipost_202108_00053v1 (pdf) 
Date submitted:  20210820 22:55 
Submitted by:  Dunjko, Vanja 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approaches:  Theoretical, Computational 
Abstract
We address the origins of the quasiperiodic breathing observed in [Phys. Rev.\ X vol. 9, 021035 (2019)] in diskshaped harmonically trapped twodimensional Bose condensates, where the quasiperiod $T_{\text{quasibreathing}}\sim$~$2T/7$ and $T$ is the period of the harmonic trap. We show that, due to an unexplained coincidence, the first instance of the collapse of the hydrodynamic description, at $t^{*} = \arctan(\sqrt{2})/(2\pi) T \approx T/7$, emerges as a `skillful impostor' of the quasibreathing halfperiod $T_{\text{quasibreathing}}/2$. At the time $t^{*}$, the velocity field almost vanishes, supporting the requisite timereversal invariance. We find that this phenomenon persists for scaleinvariant gases in all spatial dimensions, being exact in one dimension and, likely, approximate in all others. In $\bm{d}$ dimensions, the quasibreathing halfperiod assumes the form $T_{\text{quasibreathing}}/2 \equiv t^{*} = \arctan(\sqrt{d})/(2\pi) T$. Remaining unresolved is the origin of the period$2T$ breathing, reported in the same experiment.
Current status:
Reports on this Submission
Anonymous Report 2 on 2021108 (Invited Report)
 Cite as: Anonymous, Report on arXiv:scipost_202108_00053v1, delivered 20211008, doi: 10.21468/SciPost.Report.3628
Strengths
1 The manuscript provides an elegant theoretical explanation of some aspects of the important experiment of the Jean Dalibard group.
2 The analytical approach employed for the study of the dynamics of an initially discontinuous density profile is highly original.
3 The manuscript provides intriguing results that should inspire further studies.
Weaknesses
1 The manuscript does not discuss differences between the numerical GrossPitaevskii solutions and the hydrodynamic solutions.
Report
The manuscript theoretically addresses an important problem posed by the puzzling experimental work of the Paris group. It (analytically) delivers a quite general, sort of surprising, result for the breathing "halfperiod", which in two spatial dimensions, surprisingly well, reproduces the experimental measurement. It is rather doubtful that such a result could be obtained with less sophisticated analytical techniques. I think that it definitely deserves publication in a highimpact journal such as SciPost after consideration of the remarks from the Requested Changes sections.
Requested changes
1The lefthand side of (2) does not match the expression proportional to $C$. I suppose the definition of $C$ needs to be corrected. The brackets, $\{\cdots\}$, look awkward in (2), etc.
2There is velocity field missing in one place in (5).
3If (8) satisfies (6), after setting $\omega=0$ and $\delta t=t$, only when $d=1$ , then how is (8) found for an arbitrary $d$? Where does it come from? What fixes its form?
4 Below (11): $0+\to0^+$.
5 Should there be "$\text{for} \ r>R_\text{outer}(t)$" in one of the expressions in (12)?
6 What does the statement "the velocity field at $t^*$ is nearly zero" mean (remark right below Statement 2)? In what sense it is nearly zero?
7 I am not sure whether $N$, in the equation for $\sqrt{\langle r^2\rangle(t)}$, is somewhere defined.
8 Could it be explained why singlevaluedness of the velocity field leads to the proof of Statement 5?
9 Should there be $\bar{R}_\text{inner}(t)$ instead of $R_\text{inner}(t)$ in the expression for $1/r$, the one below (23)?
10 It is unclear to me why the empty space, the $\omega t\gtrsim1$ region, is displayed in Fig. 1. Can analytical solutions be put on the plot as, e.g., the dashed line? Change "(34)(28)" in the caption to, e.g., "(28,34)".
11 Should it be said in Statement 7 that (34) holds when $\mu(n)\sim n^\nu$? Similar question applies to Statement 8.
12 Should one explicitely define $V_\text{outer}(t)$ as $dR_\text{outer}(t)/dt$?
13 Figs. 2 and 3: the "dotted curves" seem to me to be rather "dashed" than "dotted".
14 Fig. 3: why the blue, GrossPitaevskii results, are missing? Smaller range of the vertical axis could lead to better presentation of the results.
15 Can GrossPitaevskii equation be written down? Can one comment on how its simulations are carried out (e.g. write explicitly what initial conditions for the numerics are employed)? Are there some specific parameters used for the GrossPitaevskiibased simulations?
Anonymous Report 1 on 2021927 (Invited Report)
 Cite as: Anonymous, Report on arXiv:scipost_202108_00053v1, delivered 20210927, doi: 10.21468/SciPost.Report.3575
Report
The manuscript by Torrents et al. addresses the origin of the 2T/7 quasiperiod of the diskshaped breather that was observed in trapped twodimensional BEC by Dalibard’s group. In this paper, the authors have brilliantly argued that at a certain time t*, the velocity field almost vanishes, which naturally leads to a 2T/7 quasiperiod. The key step of their argument is to focus on the trajectories of the singularities of the density and velocity fields (R_{inner} and R_{outer}) instead of the whole solution to the hydrodynamic equations. By doing this, they are able to extract the exact trajectories of both R_{inner} and R_{outer} and hence prove that the velocity field at these two points vanishes at t*. Since this work solves an experimental puzzle and the method is quite innovative, I think this paper can be published in Scipost.
Requested changes
Minor comments:
1．In eq.(5), a $\mathbf{v}$ is missed on the L.H.S. of the Euler’s equation.
2．Last paragraph on page 12, 2/7 should be around 2.86, instead of 2.96.
3．The converse of statement 1 (If the velocity field vanishes at T/2, then T is a period of the breather.) should also be true. Even though it is not directly related to the following argument in the paper, it will be nice that the authors add it to statement 1 to make things complete.
Author: Vanja Dunjko on 20211116 [id 1946]
(in reply to Report 1 on 20210927)
We thank the referee for this very encouraging review: our replies are below.
Requested changes:
 Corrected.
 Corrected.
 We do not think this is correct in general. At the very least, we do not see a fundamental reason why there should not be more than two timereversalinvariant instants of time per period. Incidentally, for the triangular breathers of [Phys. Rev. X vol. 9, 021035 (2019)], the instants of time $T_{\text{breathing}}/4$ and $3T_{\text{breathing}} /4$ come pretty close, exhibiting an identically vanishing velocity field along six separate rays; the rays are aligned with the corners of the right pentagon formed by the bulk density region at these moments. For this breather, $T_{\text{breathing}} = T /2$.
Author: Vanja Dunjko on 20211116 [id 1945]
(in reply to Report 2 on 20211008)We wholeheartedly thank the referee for the appreciation of our work and for the extremely detailed reading and analysis. Below, we describe the way we are going to incorporate the suggestions that were offered.
Requested changes
Corrected.
Thank you, corrected.
We have added an explanation, after Eq. 12.
Corrected.
Corrected.
We compare the velocity at \(t^{*}\) with the typical velocity inside the shockwave front, in the middle of the time evolution. The relevant velocity scale is the initial speed of sound. We added an explanation, in the second paragraph of the summary.
Corrected.
We substantially lengthened the proof.
Your help is truly invaluable. Corrected.
This is an important point:
We added an explanation to the caption of Fig. 1.
In short, the hydrodynamical solution ceases to exist at \(t^{*}\) . The reason is a conical singularity emerging at the origin, at this instant of time.
If the velocity were vanishing exactly, one might use the timereversal invariants to carry the hydrodynamics through the catastrophe. But since the velocity vanishing is not exact, one needs a theory one level higher to restart hydrodynamics (with new initial conditions) after such a “Big Crunch/Big Bang”. The GrossPitaevslii equation does precisely that.
As the ENS experiment and numerics show, the time evolution after \(t^{*}\) is both regular and erratic: there are regular quasirevivals at integer multiples of \(2T /7\), but also there are quasirandom small modifications of the initial state at the quasirevival points. Nonetheless, at \(t = 7(2T /7) = 2T\) , the revival is exact.
A small remark: formally speaking, the upper curve (outer edge as a function of time) in Fig. 1 could be continued for a bit further, as the singularity at the origin will need some time to propagate to the edges. But we believe that such a continuation would be misleading.
This was a significant omission on our part. Corrected in both places, thank you.
Introduced, between the Eqs. 11 and 12.
Corrected.
We decided against attempting to introduce a quantummechanical analogue of the local velocity. We added a sentence stating this, at the end of the “Summary and numerical results” section.
We added a paragraph at the end of the “Summary and numerical results” section. For this comment, we are especially thankful: it reminded us to mention, in the paper, that our GrossPitaevskii numerics aimed to replicate the ENS experiment [Phys. Rev. X vol. 9, 021035 (2019)] verbatim.