{"id":"8f2cb341-5d89-4dc6-b1a6-f5648ead1f12","arxiv_id":"2607.29309","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Across the density-dependent-gauge-field-driven BEC-to-soliton transition, quasiparticle modes harden and a finite mode gap opens, while repulsive interactions push the critical gauge coupling upward.","lead":"This paper studies how the low-energy vibrations of an ultracold-atom cloud in an optical lattice change when a density-dependent gauge field drives the cloud from an extended condensate into a compact soliton state. The authors find that the soliton's excitation modes stiffen, repulsive interactions push the transition to stronger gauge fields, and the soliton resists both Bloch oscillations and trap-quench breathing.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"BdG eigenvalues are plotted as real 'mode energies' despite a non-Hermitian matrix; without a stated extraction rule, mode hardening, the zero-energy mode, and the gap phase diagram are ungrounded.","rationale":"Reader identified the BdG linearization and non-Hermiticity as weakest; I agree. The mode hardening and gap phase diagram are the quantitative core of the paper, and they rest entirely on the eigenvalue computation. The dropped anomalous-density term is a known approximation (Popov-like), but the text does not justify it sufficiently, and the non-Hermitian-eigenvalue handling is undocumented. In addition, Eq. (4) is garbled, so the BdG matrix 'given' in Eq. (8) cannot be audited from the manuscript alone. The γ_c=4.65 perturbative estimate also sits unexplained alongside 3.45, and the abstract's 'stabilizing the lattice soliton' wording conflicts with the body's 'U opposes localization' — but these are secondary to the spectral method unless the eigenvalue handling is clarified. The dynamical claims (no Bloch, no breathing) are less sensitive to this, but they are also supported by the same ground states and are not the primary quantitative claim. Thus a single check on the BdG eigenvalue reality and zero-mode robustness is the decisive test. If it passes, the central claim is likely sound; if it fails, the paper's main quantitative assertions must be revised.","tokens_in":15722,"tokens_out":8867,"duration_ms":81918,"concrete_test":"For a representative soliton state (γ=4.5i, U=0.05, L=110, J=1, Ω=0.001, N=100), obtain c_j via imaginary-time propagation, form the BdG matrix from Eq. (8), and diagonalize with ZGEEV. Count eigenvalues with |Im E|>10^-6; if any, report how Re E and |E| differ and re-plot Fig. 2(b) accordingly. Separately restore the anomalous-density term in the linearized Eq. (4), re-diagonalize, and check whether the minimal |Re E| remains zero within 10^-5. If complex eigenvalues appear or the lowest mode moves off zero, the central spectral claims are not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All spectral claims — mode hardening, the preserved zero-energy mode, the finite mode gap, and the phase diagram — derive from the 2L×2L BdG matrix (Sec. III.A, Eq. 7). The paper itself states that this matrix 'could be non-Hermitian and non-symmetric' and is diagonalized with LAPACK ZGEEV, with no rule given for converting complex eigenvalues to the real E_l displayed in Figs. 2–5 and used for ΔE in Fig. 6. A non-Hermitian BdG matrix can have complex eigenvalues; if any exist in the soliton phase, the plotted 'energies' are ambiguous, and the gap ΔE = E_1 − E_0 is not defined. This concern is compounded by Sec. II's statement that the anomalous density ⟨φφ⟩ was omitted 'in order to obtain the gapless excitation spectrum' — if that omission is what creates the zero eigenvalue, then the 'preserving the zero-energy mode' claim is an artifact. The two issues jointly put the central claim at risk: the hardening curves and phase boundaries would not follow from the stated computations unless an explicit convention (e.g., Re E, |E|, or a symmetrized BdG form) is adopted and the zero mode is shown to be robust to the anomalous-density approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a one-dimensional Bose-Hubbard model with density-dependent tunneling (DDT) in an optical lattice. Using imaginary-time propagation for the ground state and discrete Bogoliubov–de Gennes (BdG) diagonalization, it reports a first-order transition from a plane-wave condensate to a localized soliton as the gauge coupling γ increases. The main claims are: low-lying quasiparticle energies harden and a mode gap opens in the soliton phase while a zero-energy Goldstone mode is preserved; repulsive U shifts the critical γ upward; the dispersion changes qualitatively; and in real-time dynamics the soliton neither Bloch-oscillates nor breathes. A phase diagram in the complex γ plane and a modulational-instability boundary complete the picture.","tokens_in":15941,"tokens_out":10153,"duration_ms":95941,"significance":"If established, the results would provide a useful characterization of a DDT-driven condensate-to-soliton transition in a setting accessible to cold-atom experiments: a finite quasiparticle gap, a U-dependent phase boundary, and unusual dynamical rigidity. The qualitative consistency of the figures and the symmetry-based assignment of Goldstone/dipole/breathing modes lend plausibility. However, the central spectral claims rest on a non-Hermitian BdG diagonalization for which no eigenvalue-extraction rule is stated, and on an explicitly imposed gaplessness; the dispersion relation used for stability also appears inconsistent with the displayed BdG matrix. These issues must be resolved before the conclusions can be relied upon.","major_comments":[{"comment":"The BdG matrix is acknowledged to be non-Hermitian and non-symmetric, and is diagonalized with LAPACK's ZGEEV. However, the paper never states how possibly complex eigenvalues are converted to the real E_l displayed in Figs. 2–5 and used for ΔE in Fig. 6. In the soliton phase a non-Hermitian BdG problem can have complex eigenvalues; if any occur, E_l and ΔE are not defined and the hardening/gap claims are ungrounded. Please either prove the eigenvalues are real in the plotted regime, or give and justify the extraction convention (Re E, |E|, symplectic diagonalization).","section":"Sec. III.A, Eq. (7)"},{"comment":"The anomalous density is omitted \"in order to obtain the gapless excitation spectrum.\" Later the preservation of the zero-energy mode is presented as a finding. If the omission is what forces E_0=0, the zero mode and the gap ΔE=E_1−E_0 are partly artifacts of the approximation. The authors should demonstrate robustness by repeating the BdG calculation with the anomalous term included (or with a controlled Popov-type truncation) and checking that a Goldstone mode still exists.","section":"Sec. II, after Eq. (4)"},{"comment":"Diagonalizing H_BdG = 2J(1−cos ka)I − (4γ_I/L) sin ka σ_y gives branch energies 2J(1−cos ka) ± (4γ_I/L) sin ka, not the square-root expression in Eq. (12). The instability threshold γ_I=JL/2 and the stable-regime diagram in Fig. 6(c) are derived from Eq. (12). This must be reconciled, or the approximation leading to Eq. (12) must be stated and justified.","section":"Eq. (12) vs Eq. (11)"},{"comment":"For U=0.05 the numerical critical couplings are reported as γ_R=3.39 and γ_I=3.45, but the text then asserts a first-order perturbative shift gives γ_c=4.65. No derivation is shown, and the two numbers differ by about 25%. Since the U-induced upward shift of the transition is a central claim, a consistent calculation or an explicit explanation of this discrepancy is required.","section":"Sec. III.A and Figs. 4, 6"}],"minor_comments":[{"comment":"The interaction strength U used for the dispersion curves is not stated. If U=0, γ=3.0i is already in the soliton phase (critical γ_I=1.69), contradicting the label \"BEC phase\"; if U=0.05, this should be stated explicitly in the caption and text.","section":"Fig. 5 and text"},{"comment":"The noninteracting critical imaginary coupling is given as 1.67 in Sec. III.A and 1.69 in Sec. III.B. Please use one consistent value.","section":"Sec. III.A/III.B"},{"comment":"Several typographical issues: \"To ensue such oscillations\" should be \"to ensure\"; \"magnifying evolution\" should be \"magnified evolution\"; and the statement that the transition is determined by both components of γ conflicts with Eq. (12), where only γ_I appears.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The central issues are the missing eigenvalue-extraction rule for a non-Hermitian BdG problem and the apparent circularity of enforcing gaplessness. If the authors can supply the method and demonstrate that the zero mode and hardening are not artifacts, the paper is likely to be of value to the quantum-gas community. The phase diagram and quench dynamics are potentially interesting, but the spectral foundation must be fixed first."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a follow-up to Faugno–Salerno–Ozawa (PRL 2024), which already established the BEC-to-soliton transition driven by density-dependent tunneling. The authors are honest about that lineage, and what they add is real: a mode-resolved view of the transition. The hardening of the dipole and higher modes as γ crosses the critical value, the mode-amplitude analysis that correctly shows Goldstone/dipole/breathing characters, the upward shift of the critical γ under repulsive U, and the dynamical results — suppressed Bloch oscillations and no breathing after a trap quench — are new relative to [20], and the story is internally consistent across figures. The modulational-instability boundary γ_I = JL/2 is cross-checked between the analytic dispersion and the time-domain correlation collapse, which is a nice touch. The physics is plausible: localization reduces effective tunneling, so moving or deforming the soliton costs more energy.\n\nThe soft spots are real but fixable. First, the BdG matrix is non-Hermitian in general, yet the paper never says how possibly complex eigenvalues become the real 'mode energies' plotted everywhere. The stress-test note is right to flag this; one explicit convention sentence would remove the ambiguity. Second, the zero-energy mode is partly by construction — the anomalous density is dropped 'in order to obtain the gapless excitation spectrum.' That is the standard Bogoliubov approximation, but it means 'preserving the zero-energy mode' is not an independent output. The finite-mode hardening survives this, so the central trend stands; the gap-based phase diagram just needs framing accordingly. Third, smaller stuff: the perturbative γ_c = 4.65 at U = 0.05 is quoted without derivation and disagrees with the numerical 3.45; the abstract says repulsion 'stabilizes' the soliton while the body correctly says it opposes localization; Eq. (4) is garbled, so the audit trail runs through Eq. (8); and numerical details like the Gaussian width and Bloch-force strength are underspecified.\n\nThe central qualitative picture holds up. This is for people working on density-dependent tunneling or lattice soliton dynamics — they will get a useful mode map. It deserves a serious referee, and the revision must state the eigenvalue convention and address the enforced gaplessness. I would want to see that revision before citing its numbers.","headline":"Follow-up on the known DDT soliton transition, giving a coherent mode-resolved and dynamical picture; the non-Hermitian BdG diagonalization and the hand-enforced zero mode need explicit treatment before the specific numbers are trusted.","tokens_in":16561,"tokens_out":8783,"would_cite":false,"duration_ms":78304,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Density-dependent tunneling drives a BEC-to-soliton transition in which the low-lying quasiparticle spectrum acquires a gap while the zero-energy Goldstone mode survives.","keywords":["density-dependent tunneling","optical lattice","Bose-Einstein condensate","bright soliton","Bogoliubov theory","quasiparticle spectrum","Bloch oscillations","breathing oscillations"],"falsifier":"A direct numerical check of the BdG matrix at γ=4.5i, U=0.05: if any eigenvalues have non-negligible imaginary parts beyond numerical noise, the plotted 'mode energies' are not well-defined, and the claimed hardening and gap would need re-evaluation. Alternatively, recomputing the spectrum while including the anomalous density term would reveal whether the zero-energy mode survives without the enforced gaplessness.","tokens_in":15456,"feed_emoji":"⚛️","tokens_out":2949,"duration_ms":28018,"temperature":0.7,"pith_summary":"The paper claims that when bosons in a one-dimensional optical lattice hop with an amplitude that depends on the local density difference between sites, the system undergoes a first-order quantum phase transition from a plane-wave condensate to a localized bright soliton. Using discrete Bogoliubov theory, the authors show that as the gauge coupling increases in the soliton phase, the dipole and higher excitation energies rise while a zero-energy mode is preserved; the excitation gap widens with the gauge field. Repulsive on-site interactions increase the critical gauge field needed for the transition and stabilize the soliton. Dynamically, the soliton does not Bloch-oscillate under a constant force and does not breathe after a trap quench, in contrast to the condensate, because the density-dependent gauge field couples center-of-mass motion to internal shape. This matters because it shows that a purely density-dependent hopping term can produce rigid, gap-protected localized states with distinct excitation and dynamical signatures.","feed_headline":"Soliton transition hardens modes and silences Bloch oscillations","feed_subtitle":"Bogoliubov analysis shows density-dependent hopping creates a gapped, rigid soliton that neither Bloch-oscillates nor breathes.","key_machinery":"The central object is the density-difference-dependent hopping term in the Bose-Hubbard Hamiltonian, which couples the local density difference to the phase of hopping and acts as an effective gauge field. The argument is carried by the discrete Bogoliubov–de Gennes equations, linearized around the mean-field ground state, whose eigenvalues give the quasiparticle mode energies; the mode gap and dispersion curves are used to locate the transition and characterize the phases.","core_discovery":"The central discovery is that the condensate-to-soliton transition induced by density-dependent tunneling is accompanied by mode hardening: the energies of the dipole and breathing modes increase with the gauge field in the soliton phase, the spectrum develops a finite gap, and the zero-energy Goldstone mode remains. The repulsive interaction shifts the critical gauge field upward by raising the energy cost of localization. The soliton's effective momentum becomes density-dependent, so it does not undergo Bloch oscillations, and its width is not excited by a trap quench—breathing oscillations are suppressed.","pith_inferences":["The same mechanism—density-dependent hopping—could be engineered via Floquet driving to probe the transition experimentally, which the paper gestures at but does not detail.","The prediction that the effective momentum becomes density-dependent suggests a Hall-like or chiral response if the gauge field has a nonvanishing real part, a testable extension.","The zero-energy mode preservation alongside a finite gap implies spontaneous translational symmetry breaking; measuring local density correlations could confirm this directly.","Because the Bogoliubov treatment enforces gaplessness by dropping the anomalous density, an independent numerical check (e.g., exact diagonalization or quantum Monte Carlo) of the gap in the strongly localized regime would be a worthwhile test of the central claim."],"forward_implications":["If correct, the quasiparticle excitation gap serves as an order parameter for the BEC-to-soliton transition in density-dependent Hubbard models.","The upward shift of the critical gauge field with repulsive interactions is observable as a shift of the mode-hardening point, providing a testable prediction.","The predicted absence of Bloch oscillations for the soliton distinguishes it dynamically from the condensate in a clean, measurable way.","The suppression of breathing after a trap quench offers a second dynamical signature of the soliton phase.","The finite-momentum zero mode suggests the soliton carries an internal phase gradient, analogous to spin-orbit coupled condensates."],"fun_headline_variants":["Density-dependent tunneling creates gapped, rigid soliton","Soliton transition silences Bloch oscillations and breathing","Quasiparticle modes stiffen as condensate becomes soliton","No Bloch or breathing: density-dependent hopping soliton","Density-dependent gauge field hardens soliton modes"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central assumption is that the Bogoliubov treatment—which drops the anomalous density to force a gapless spectrum and diagonalizes a potentially non-Hermitian 2L×2L matrix—remains valid in the strongly localized soliton phase; if the true spectrum is complex or gapped there, the reported mode hardening and phase boundaries do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Density-dependent tunneling creates gapped, rigid soliton","Soliton transition silences Bloch oscillations and breathing","Quasiparticle modes stiffen as condensate becomes soliton","No Bloch or breathing: density-dependent hopping soliton","Density-dependent gauge field hardens soliton modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001345,"raw_usage":{"total_tokens":5314,"prompt_tokens":771,"completion_tokens":4543,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":4463}},"tokens_in":515,"tokens_out":4543,"duration_ms":39057,"temperature":1.0,"reasoning_tokens":4463,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T09:31:02.946787+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical check of the BdG matrix at γ=4.5i, U=0.05: if any eigenvalues have non-negligible imaginary parts beyond numerical noise, the plotted 'mode energies' are not well-defined, and the claimed hardening and gap would need re-evaluation. Alternatively, recomputing the spectrum while including the anomalous density term would reveal whether the zero-energy mode survives without the enforced gaplessness.","supporting_citations":[],"review_version":1}