{"id":"f4e39ca1-83de-46d1-b64c-09b2444acac3","arxiv_id":"2607.16820","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"The 2D logarithmic droplet GPE is asymptotically reduced to KP-I and Davey-Stewartson-I models, from which line, lump, ring, and dromion solitary waves are constructed and checked numerically.","lead":"This paper uses a standard approximation technique to turn a model of ultracold quantum droplets into two simpler, exactly solvable wave equations, then builds four kinds of 2D wave shapes — stripes, lumps, rings, and dromions — from their known solutions. Direct simulations of the original droplet equation show the waves persist, with dromions matching theory the least.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dromion branch is not supported by the DS-I reduction as presented: at ρ0=2.25 the residual R[q] is ~56% of the retained nonlinearity, and the R→0 limit coincides with vanishing sound speed (ρ0=e^{-1}), so no simulation is shown in a valid small-perturbation window.","rationale":"Read in good faith, the paper's KP/cKP reductions are standard and plausible; the line, lump, and ring results are probably not affected by my concern. However, the DS/dromion part has an internal consistency problem: the integrable DS-I limit is reached only at lnρ0=-1, where c=0, so the stated 'small-R' window is not inside the modulationally stable regime with finite sound speed. At the only dromion simulation shown, ρ0=2.25, the residual coefficient 1+2A≈0.56 makes the dropped term comparable to the retained nonlinearity. This is a correctness risk within the paper's own asymptotic framework, not a disagreement with external consensus. The paper's Sec. V assertion that deviations are O(ε²) at initial times is therefore unsupported for dromions and likely false. The reader already identified the non-small R at ρ0=2.25 as the weakest assumption; my analysis sharpens this by showing the R→0 limit is degenerate, so the claimed asymptotic window is essentially empty. The appropriate disposition remains CONDITIONAL: the rest of the paper can stand, but the dromion claim needs either a revised theoretical justification or a simulation in a genuinely small-R, non-degenerate regime.","tokens_in":21400,"tokens_out":16239,"duration_ms":138375,"concrete_test":"Compute the t=0 residual of Eq. (2) for the dromion initial condition (37) at fixed ε=0.01, k=1.2 and three backgrounds: ρ0=2.25 (1+2A≈0.56), ρ0=e^{-0.9} (1+2A≈0.15), and ρ0=e^{-0.99} (1+2A≈0.02), using the same spatial/temporal discretization (dx=dy=0.08, dt=10^-4). If the residual at ρ0=2.25 is O(ε) rather than O(ε²), and does not decrease approximately in proportion to 1+2A as ρ0→e^{-1}, the dromion is not a valid asymptotic solution. As a dynamical cross-check, evolve the ρ0=e^{-0.9} dromion and the ρ0=2.25 dromion to t~200 and compare a quantitative error metric (e.g., L2 distance from the corresponding DS-I dromion; also peak-depth/phase-jump errors). If the error does not drop as the residual coefficient shrinks, the DS-I reduction is not controlling the dromion dynamics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing weak point is the DS-I/dromion reduction. In Eq. (31), the residual is R[q]=-(1+2A)|q|²q with A=-(2+lnρ0)/(8+6lnρ0). Hence 1+2A = 2(1+lnρ0)/(4+3lnρ0). Modulation stability requires c²=gρ0(1+lnρ0)>0 (upper Lambert branch), so on the stable branch 1+2A>0. At the simulation point ρ0=2.25, 1+2A≈0.563: the term dropped to reach the integrable DS-I limit is 56% as large as the retained -|q|²q. This is not a small perturbation. The paper's stated small-perturbation window ρ0≳e^{-1} is also ill-posed: at ρ0=e^{-1}, c=0, the linear regime degenerates and the DS coefficients in Appendix A (e.g., c1 and c3) are not well behaved; no stable background with finite sound speed has R small. Therefore the dromion initial condition (37) is not an asymptotic solution of Eq. (2) in the simulated regime, and the observed persistence cannot be attributed to the DS-I reduction. The KP-I/cKP-I line/lump/ring reductions do not share this defect and are not called into question by this concern.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the two-dimensional extended Gross-Pitaevskii equation (eGPE) with a logarithmic nonlinearity, as used for quantum droplets. After analyzing the modulational stability of a homogeneous droplet background, the authors perform multiscale asymptotic reductions to derive effective integrable models: the KP-I equation (Eq. 17), the cylindrical KP-I equation (Eq. 20), and a nearly integrable Davey–Stewartson system (Eqs. 31). From these they construct approximate analytical solutions for line, lump, ring, and dromion solitary waves on the stable background, and test them in direct numerical simulations of the 2D eGPE. The line, lump, and ring solutions are reported to persist with small radiation, while the dromion keeps its envelope but deviates most from the analytical waveform.","tokens_in":21685,"tokens_out":10281,"duration_ms":89842,"significance":"The paper offers a systematic path from a non-integrable 2D eGPE with competing nonlinearities to integrable reductions, and gives explicit waveforms for several 2D solitary-wave families in a droplet environment. The KP/cKP derivations for the line, lump, and ring are carefully presented, the algebra is largely checkable, and the simulations provide qualitative evidence for the persistence of these structures. If fully correct, the work would extend the catalogue of multidimensional solitons in models with Lee–Huang–Yang corrections and could inspire experiments. The main weakness is the DS-I/dromion branch, where the perturbative justification is quantitatively problematic; this issue is analyzed below. The overall contribution would be strengthened by addressing this point and by quantifying the claimed accuracy.","major_comments":[{"comment":"The DS-I reduction for dromions is not justified as a small-perturbation limit at the parameters used in the simulation. The residual in Eq. (31a) is R[q] = -(1+2A)|q|^2 q with A = -(2+ln ρ0)/(8+6 ln ρ0), so 1+2A = 2(1+ln ρ0)/(4+3 ln ρ0). On the modulationally stable branch (c^2 > 0) we have 1+ln ρ0 > 0, hence 1+2A > 0. At the dromion simulation point ρ0 = 2.25, 1+2A ≈ 0.563, so the term dropped to reach the integrable DS-I limit is about 56% as large as the retained -|q|^2 q term. The manuscript's statement that R is small for ρ0 ≳ e^{-1} (Sec. III B) is not useful: as ρ0 → e^{-1}+, the sound speed c → 0, and the coefficients c1, c3, c4 in Appendix A diverge (c3 ~ c^{-3}, c4 ~ (1+ln ρ0)^{-1}), so the long-wavelength DS reduction degenerates. Thus no stable background with finite sound speed has a small R. Consequently, the dromion initial condition (37) is not an asymptotic solution of","section":"§III B, Eq. (31) and §IV C, Fig. 7"},{"comment":"The paper asserts that 'the deviation between our analytical solutions and the time evolved states is of the order of ε^2 at the initial stages' without defining an error norm or providing any numerical measurement. Since this statement is used to characterize the validity of all four approximations, and since the dromion residual is not small at the chosen parameters, a quantitative comparison (e.g., L2 or L∞ difference between the evolved state and the predicted waveform at several early times) should be supplied. Without such data, the claim is unsupported and should be softened or removed.","section":"Sec. V, first paragraph"}],"minor_comments":[{"comment":"The ring-soliton amplitude contains η^2 = κ1^2 (t0/t)^{2/3}, which is singular at t = 0. The text and Fig. 6 do not specify how the t = 0 initial condition is defined (presumably t is set to t0 in the initialization) or what value of t0 is used. Please clarify the initialization procedure and give the parameter value.","section":"Sec. IV B, Eq. (36)"},{"comment":"There is a typographical error: 'the DS-I system of Eqs. (31a)-(31a)' should read 'Eqs. (31a)-(31b)'. This occurs both in Sec. IV C and in the surrounding text in Sec. III B.","section":"Sec. IV C, Eq. (31) references"},{"comment":"The dromion simulation uses k = 1.2 with the long-wavelength approximations for the coefficients c_j (Appendix A). The dispersion ratio k^4/4 vs c^2 k^2 is about 0.09, and higher-order terms in k may not be negligible. It would be useful to state the range of k over which the leading-order coefficients are quantitatively accurate, or to compare with the full coefficient values.","section":"Appendix A and Sec. IV C"}],"recommendation":"major_revision","confidential_remarks":"The KP/cKP parts (line, lump, ring) appear technically sound and are a solid contribution. The dromion/DS-I part is the main problem: the residual is not small in any finite-sound-speed stable regime, and the divergence of the coefficients near ρ0 = e^{-1} makes the claimed small-perturbation window illusory. This is a load-bearing issue for one of the paper's four central results. I would ask the authors to either remove the dromion from the asymptotic analysis or re-frame it as a purely numerical observation, and to quantify the claimed O(ε^2) accuracy. With those changes, the remaining content could be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a real asymptotic-reduction paper: the multiscale derivations of KP-I and cKP-I with explicit coefficients are carefully done, the Lambert-W stability analysis is coherent, and the line, lump, and ring soliton checks look like genuine persistence experiments. Second, the dromion section is the weak part of the manuscript, and the specific simulation shown does not sit in the regime the DS-I reduction claims to describe.\n\nThe genuinely new pieces are the explicit KP-I, cKP-I, and DS-I reductions for this logarithmic eGPE, with coefficients written in terms of rho0, and the first attempt to place dromions on a droplet background. The line and lump reductions are standard in structure but usefully worked out here, and the numerics visually support the claim that those waveforms are approximate coherent structures of the full 2D equation. The citation pattern looks appropriate, and the paper is honest that the dromions deviate the most.\n\nThe soft spots are specific and not all equally serious. The Sec. V claim that deviations are \"of the order of epsilon^2 at the initial stages\" is asserted without showing any error metric or convergence check. That is a fixable omission. More important, the dromion simulation at rho0 = 2.25 gives 1 + 2A approximately 0.56, so the term dropped to reach integrable DS-I is more than half the size of the retained nonlinear term. That is not a small perturbation. The paper's stated small-perturbation window near rho0 = e^{-1} is also ill-posed: at that point the sound speed vanishes and the DS coefficients are not well behaved, so no finite-speed background actually makes R[q] small. This does not destroy the numerical observation that a dromion-like object persists in the eGPE, but it does mean the DS-I reduction as presented does not explain that persistence. The authors should either find a scaling that makes the residual genuinely small, or reframe the dromion as a numerical observation supported only by the simulation. Minor points: no code or data is provided, and the ring-soliton reference time t0 is left unspecified.\n\nThe line, lump, and ring results are not affected by the dromion concern, and the paper is worth engaging with despite that weak spot. I would send it to peer review and ask for revision focused on quantifying the deviations and repairing or reframing the dromion claim.","headline":"Solid asymptotic-reduction paper for the logarithmic eGPE, with credible KP-I/cKP-I line/lump/ring numerics; the dromion branch is oversold because the simulation runs where the discarded DS-I term is not small.","tokens_in":22393,"tokens_out":1638,"would_cite":true,"duration_ms":18687,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35Q53","35Q51","37K40"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that the 2D extended Gross-Pitaevskii equation for quantum droplets reduces, in the long-wavelength limit, to integrable Kadomtsev-Petviashvili and Davey-Stewartson equations, from which line, lump, ring, and dromion solit","keywords":["quantum droplets","Lee-Huang-Yang corrections","extended Gross-Pitaevskii equation","Kadomtsev-Petviashvili equation","Davey-Stewartson equation","dromions","lump solitons","ring solitons"],"falsifier":"Run the dromion initial condition at two background densities: near ρ0 ≈ e^{−1} (where 1 + 2A ≈ 0, so the DS remainder is nearly absent) and at ρ0 = 2.25 (where the remainder is ≈ 0.56). If the dromion disperses or breaks apart far faster at ρ0 = 2.25 than near e^{−1}, then the claim that the DS-I reduction governs the dynamics for the simulated parameters would be refuted; conversely, similar survival times would show the non-integrable term is not the controlling factor. Additionally, computing the Bogoliubov–de Gennes spectrum of the four states would reveal any negative eigenvalues that im","tokens_in":21168,"feed_emoji":"🌀","tokens_out":3638,"duration_ms":35455,"temperature":0.7,"pith_summary":"The paper tries to establish that a 2D quantum droplet model—the extended Gross-Pitaevskii equation with competing mean-field and logarithmic Lee-Huang-Yang nonlinearities—supports four families of previously unseen two-dimensional solitary waves. On the modulationally stable branch of the droplet background, it derives effective integrable models (KP-I, cylindrical KP-I, and DS-I) via multiscale expansions, then reads off approximate analytical solutions for line, lump, ring, and dromion solitons. Direct numerical integration of the full equation shows that all four states propagate over long times, emitting only weak radiation, with line, lump, and ring staying closest to theory and dromions deviating most. A sympathetic reader would care because genuine 2D solitons are usually destroyed by collapse or transverse instabilities in ultracold gases; this work identifies a concrete, experimentally accessible setting where they persist.","feed_headline":"Four 2D soliton families survive in quantum droplet models","feed_subtitle":"Long-wavelength reductions to KP and Davey-Stewartson equations predict line, lump, ring, and dromion states that persist in full simulation","key_machinery":"The multiscale asymptotic reduction is the central mechanism. The density is expanded as ρ = ρ0 + ερ1 + ε²ρ2 + ... and the phase as a similar series with either half-integer or integer powers of ε, with stretched variables X = ε^{1/2}(x − ct), Y = εy, T = ε^{3/2}t. Solvability conditions at successive orders yield the compatibility condition C² = gρ0(1 + lnρ0), then the KP-I equation (Eq. 17), the cylindrical KP-I equation (Eq. 20), and the DS-I system (Eq. 31) with a non-integrable remainder R[q] = −(1 + 2A)|q|²q. These reduced integrable models supply exact soliton solutions that serve as approximate initial conditions for the original problem.","core_discovery":"The central claim is that the 2D eGPE with logarithmic nonlinearity, linearized around a homogeneous droplet background, supports weakly nonlinear 2D solitary waves in the modulationally stable regime. Using density and phase expansions with stretched coordinates, the authors reduce the non-integrable eGPE to the KP-I equation in Cartesian geometry, to Johnson's cylindrical KP-I equation in polar geometry, and to a nearly integrable Davey-Stewartson system for a carrier wave with a mean-flow term. From these reductions they obtain approximate analytical dark line solitons, algebraically decaying lump solitons, expanding ring solitons, and exponentially localized dromions. Simulations of the","pith_inferences":["An editorial extension: the dromion simulation at ρ0 = 2.25 lies far outside the window where the DS remainder R[q] = −(1 + 2A)|q|²q is small (1 + 2A ≈ 0.56 there); testing the dromion near ρ0 ≈ e^{−1}, where 1 + 2A ≈ 0, would cleanly separate the integrable prediction from the non-integrable correction.","The line-to-lump fragmentation observed at moderate perturbation amplitude suggests a striking experimental signature: imprinting a curved or periodically bent dark stripe on a droplet background should spontaneously produce a regular chain of lump solitons, observable in situ.","The paper leaves spectral (Bogoliubov–de Gennes) stability unexamined; if those spectra were computed, one would expect the ring and dromion states to harbor instabilities at longer times than simulated, constraining their true lifetimes.","The same multiscale machinery could be applied to 3D droplet shells or dipolar droplet settings, where analogous reductions to integrable equations might yield spherical or vortex-tangle soliton states."],"forward_implications":["If the reduction is valid, the eGPE with logarithmic nonlinearity is an experimentally relevant platform where 2D dark line, lump, and ring solitons, as well as dromions, can be generated and observed over milliseconds-scale evolution.","The KP-I and DS-I reductions provide quantitative predictions—soliton speeds, depths, phase jumps, and the ring's amplitude decay ∝ (t0/t)^{2/3}—that can be tested against direct imaging of density and phase in ultracold atom experiments.","The simulation showing a perturbed line soliton dissolving into a lump array confirms the KP-I transverse instability scenario in a droplet medium, implying that stripe perturbations are a practical route to generating lump arrays.","All four states are only approximate, so their persistence rests on the smallness of the neglected corrections; the paper bounds initial deviations as O(ε²), with gradual growth due to radiation.","The identified reduction chain (eGPE → KP-I/DS-I → soliton families) suggests the same method can be extended to other nonlinearities that feature competing attractive and repulsive interactions."],"fun_headline_variants":["2D solitons thrive in quantum droplet models","Quantum droplets host four soliton families","Line, lump, ring, dromion: 2D solitons persist","eGPE reductions predict 2D solitons in droplets","Quantum droplet solitons: KP and DS guide the way"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire construction assumes that the neglected O(ε²) density corrections, and the non-integrable Davey-Stewartson remainder R[q], stay small over the simulated timescales—an assumption that is explicitly violated for the dromion parameter choice ρ0 = 2.25, where 1 + 2A ≈ 0.56.","fun_headline_variants_meta":{"raw":{"variants":["2D solitons thrive in quantum droplet models","Quantum droplets host four soliton families","Line, lump, ring, dromion: 2D solitons persist","eGPE reductions predict 2D solitons in droplets","Quantum droplet solitons: KP and DS guide the way"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000802,"raw_usage":{"total_tokens":3355,"prompt_tokens":734,"completion_tokens":2621,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":478,"completion_tokens_details":{"reasoning_tokens":2536}},"tokens_in":478,"tokens_out":2621,"duration_ms":17659,"temperature":1.0,"reasoning_tokens":2536,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T19:52:40.040900+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the dromion initial condition at two background densities: near ρ0 ≈ e^{−1} (where 1 + 2A ≈ 0, so the DS remainder is nearly absent) and at ρ0 = 2.25 (where the remainder is ≈ 0.56). If the dromion disperses or breaks apart far faster at ρ0 = 2.25 than near e^{−1}, then the claim that the DS-I reduction governs the dynamics for the simulated parameters would be refuted; conversely, similar survival times would show the non-integrable term is not the controlling factor. Additionally, computing the Bogoliubov–de Gennes spectrum of the four states would reveal any negative eigenvalues that im","supporting_citations":[],"review_version":1}