SciPost Submission Page
Multicomponent one-dimensional quantum droplets across the mean-field stability regime
by I. A. Englezos, P. Schmelcher, S. I. Mistakidis
This is not the latest submitted version.
Submission summary
| Authors (as registered SciPost users): | Ilias Englezos |
| Submission information | |
|---|---|
| Preprint Link: | https://arxiv.org/abs/2502.08392v3 (pdf) |
| Date submitted: | July 11, 2025, 9:48 a.m. |
| Submitted by: | Ilias Englezos |
| Submitted to: | SciPost Physics |
| Ontological classification | |
|---|---|
| Academic field: | Physics |
| Specialties: |
|
| Approaches: | Theoretical, Computational |
Abstract
The Lee-Huang-Yang (LHY) energy correction at the edge of the mean-field stability regime is known to give rise to beyond mean-field structures in a wide variety of systems. In this work, we analytically derive the LHY energy for two-, three- and four-component one-dimensional bosonic short-range interacting mixtures across the mean-field stability regime. For varying intercomponent attraction in the two-component setting, quantitative deviations from the original LHY treatment emerge being imprinted in the droplet saturation density and width. On the other hand, for repulsive interactions an unseen early onset of phase-separation occurs for both homonuclear and heteronuclear mixtures. Closed LHY expressions for the fully-symmetric three- and four-component mixtures, as well as for mixtures comprised of two identical components coupled to a third independent component are provided and found to host a plethora of mixed droplet states. Our results are expected to inspire future investigations in multicomponent systems for unveiling exotic self-bound states of matter and unravel their nonequilibrium quantum dynamics.
Author indications on fulfilling journal expectations
- Provide a novel and synergetic link between different research areas.
- Open a new pathway in an existing or a new research direction, with clear potential for multi-pronged follow-up work
- Detail a groundbreaking theoretical/experimental/computational discovery
- Present a breakthrough on a previously-identified and long-standing research stumbling block
Author comments upon resubmission
Dear Editor-in-charge of SciPost Physics,
Thank you very much for your communication on the above manuscript submitted to SciPost Physics. We are pleased with the positive receptions of our work by the referee. For instance, the referee remarks: "I see that this work has a clear potential for multi-pronged follow-up work. I recommend publication." Also, the recommendation of the referee is: "Publish (easily meets expectations and criteria for this Journal; among top 50%)".
Yet, the reviewer offers some helpful remarks towards the improvement of the manuscript. In the presently resubmitted revised version of our work, we have taken all of these comments into consideration and performed the corresponding changes. In the following, you can find our detailed reply to all questions and suggestions raised by the referee. A list of changes is appended after the point by point reply to the comments of the referee.
Given the positive outlook of the review, we hope that you will find the revised manuscript suitable for publication in SciPost Physics. We will, thus, be looking forward to your final editorial decision on the manuscript.
On behalf of all the authors, Ilias Englezos
- The referee writes: "The manuscript reports an in-depth analysis of BMF quantum mixtures in 1D. The main merit is the derivation of LHY term valid for a large range of interspecies interactions, not only close to the critical points of MF collapse and immiscibility. The extension to three and four components is valuable, although inevitably limited to a few specific configurations, due to the large number of parameters (masses, interaction strength). In terms of usability for experiments, a limitation could be traced in the restriction to the homogenous case. However, in my opinion this is a necessary and justified choice to keep the number of external parameters at a manageable level. I see that this work has a clear potential for multi-pronged follow-up work. I recommend publication."
Our response is: We thank the referee for their positive evaluation of our work and useful comments which helped us to improve and further clarify our presentation. Below, we address in detail all the issues raised. Corresponding modifications have been performed in the revised manuscript and for convenience a list of changes is also appended at the end of this letter.
- The referee writes: "1- Page 2, in the sentence The emergence of the self-bound liquid-type configurations known as quantum droplets in contact interacting three-dimensional (3D) bosonic mixtures'', the reference [24] should be complemented with [44], as the latter was actually published earlier."
Our response is: We agree with the referee. In the revised version of the manuscript, we have also placed this reference at the suggested location (see also the list of changes).
- The referee writes: "2- Page 4, sec. 2.1: it is stated that 1D configurations with periodic boundary conditions, i.e. a 1D ring,\dots can be experimentally reached with digital micromirror devices''. This is not correct, DMD can make rings, but these are hardly 1D, as the transverse confinement is insufficient. Also, Ref. [69] describing a disk-shaped potential seems quite not appropriate here."
Our response is: We agree with the referee that these DMD setups do not necessarily reach the 1D geometry considered in our work. In the revised manuscript, we have removed this sentence and added a comment about the persistence of our findings with box potentials (see also the list of changes).
- The referee writes: "3- Page 4, sec 2.1: dimensionally, $\hbar^3/(m_Ag_A)$ cannot be a unit of time."
Our response is: We thank the referee for pointing out this typo. We have corrected it to $\hbar^3/(m_A g_A^2)$, which is indeed a unit of time (see also the list the changes).
- The referee writes: " 4- Page 5, below eq. 3: LaTeX, \verb|diag| should be typed as text, I think.
Our response is: We agree with the referee and we have changed ``diag" to text format in the revised manuscript (see the list of changes).
- Finally the referee points out the weakness: "In the three component mixture, the MF stability region is asymmetric: $-g/2\leq G\leq g$. No insight is given about this asymmetry."
Our response is: We thank the referee for their suggestion. We have included a new discussion in Appendix A, to provide an intuitive explanation for this asymmetry. An additional comment has been incorporated in the main text, see also the list of changes.
List of changes
List of changes in the revised manuscript:
1) In the second paragraph of page 2, Section 1, the reference Science 359(6373), 301 (2018) has been placed earlier. 2) In the first paragraph of page 4, Section 2.1, a sentence regarding the possible DMD experimental implementation has been removed and replaced with a comment on the persistence of our results with a box potential. 3) In the second paragraph of page 4, Section 2.1, a typo in the time units used throughout the text has been fixed. 4) In the first paragraph of page 5, Section 2.2, the term "diag" has been changed to text format as suggested. 5) Second paragraph of page 16, Section 4.1, a comment pointing to Appendix A has been added. 6) In the last paragraph of pages 23-24, Appendix A, a short discussion about the asymmetry of the mean-field stability regime in higher-component mixtures has been included. 7) All references have been updated whenever necessary. 8) On page 9, we fixed a typo in the legend of Figure 1. 9) On page 10, a typo has been corrected in the legend of Figure 2. 10) On page 12, a typo has been amended in the legend of Figure 3. 11) On page 18, a typo has been updated in the title of Figure 5.
Current status:
Reports on this Submission
Report #3 by Anonymous (Referee 4) on 2025-8-13 (Invited Report)
- Cite as: Anonymous, Report on arXiv:2502.08392v3, delivered 2025-08-13, doi: 10.21468/SciPost.Report.11741
The referee discloses that the following generative AI tools have been used in the preparation of this report:
I used the AI Assist tool of overleaf to help english readibility.
Strengths
1) The authors provide an analytical derivation of the Lee-Huang-Yang (LHY) energy from two- to four-component one-dimensional bosonic mixtures, valid across the full Mean-Field stability regime.
2) Solving numerically the resulting extended Gross-Pitaevskii equations, they show quantitative deviations from the original MF-boundary LHY approximation.
3) The three- and four-component derivations pave the way for exploring exotic droplet phases.
Weaknesses
1) When considering Hetero-nuclear two-component mixtures, the authors compute numerically the LHY energy density (Eq. 11) to solve the eGPEs. They then compare to the exact equal-mass eGPEs (Eq. 14) with mass-imbalanced inserted by hand. The results perfectly agree (Fig. 3), but the authors do not comment on the reason behind this agreement.
Report
The authors numerically solve the resulting eGPEs to explore the ground-state properties of different interaction regimes and mixtures. They show substantial deviations from the original MF-boundary approximation. As the number of components increases, the parameter space grows drastically and the authors present a few specific scenarios. An interesting result is that the effect of quantum fluctuations is enhanced when the number of components increases.
The manuscript is well written and pedagogical.
Requested changes
1) p.8 "calculations [40–43]or experimental" -> missing space 2) Fig. 1. In the caption it could be specified that g_{AB} is given in unit of g g_{12} should be replaced by g_{AB} 3) p.10 g_{12} is used instead of g_{AB} many times
Recommendation
Publish (easily meets expectations and criteria for this Journal; among top 50%)
Report #2 by Anonymous (Referee 3) on 2025-8-10 (Invited Report)
- Cite as: Anonymous, Report on arXiv:2502.08392v3, delivered 2025-08-10, doi: 10.21468/SciPost.Report.11720
Strengths
1- Thorough and lucid, this paper can be read as a guide to specific calculations or as a reference as is for future work 2- Detailed examples are presented, and the most important scenarios explored
Weaknesses
1- While the content lends itself well as a reference for future work, readers would likely benefit from a succinct summary of results, since it is easy for a reader to get lost in the sheer amount of parameters involved.
2- Only minor points, detailed in Requested Changes section.
Report
Overall, I view this work as a valuable resource for future analytic and numeric studies of mixtures in 1D. Since the LHY correction affects how these systems will behave near instability, it will be a valuable resource in understanding how mixtures behave in 1D, as well as better understanding quantum fluctuations in low dimensions. I recommend publication as it fits within Scipost's expectation of opening a new pathway in an existing research direction, with clear potential for multi-pronged follow-up work.
Requested changes
1- On Pg. 1, the authors write: “In 1D that we focus herein, the LHY term is attractive [53] instead of being repulsive as in 3D. As a result, the 1D droplet parametric regions i) do not feature collapse […] within the respective MF stability regime.” As written, this is confusing, as it seems to imply that an attractive LHY correction leads to stabilization, rather than being stabilized by the 1D geometry (e.g. KE term ~1/L^2) vs. the low powers of |psi| in the corrections. Upon reading the rest of the manuscript, it’s apparent the authors do not intend this meaning, but rather the “early” onset of droplets when compared to the non-LHY case. I would suggest changing these sentences.
2- The ground-state phases including three symmetric droplets (e.g. Fig 2 (c) ;Fig 3 (d), (f); Fig. 5 (e), (f)) give me a slight pause w.r.t. whether they are the true ground states. If the system separates into flat-topped droplets, why would one component split into two separated droplets, heightening its kinetic energy, rather than a single long droplet of the same average density? On P.g 12-13 the authors argue that the heavier component splits (rather than the lighter component) to reduce kinetic energy, but it’s not clear why either should split at all. In the case of even weak confinement [even weaker than mentioned in the footnote 4], then I would agree, the authors remind us we are in free space, I can’t see why we have this symmetric splitting.
Recommendation
Publish (meets expectations and criteria for this Journal)
Report
Recommendation
Publish (meets expectations and criteria for this Journal)
