{"id":"5eaf5616-d123-4fcc-ac2a-8bc2f334c837","arxiv_id":"2507.10910","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Hopfion quantum droplets with Hopf numbers zero through seven can be stable in a toroidally trapped binary BEC when the Lee-Huang-Yang correction is present.","lead":"Stable ring-shaped knots called hopfions, with topological charges up to seven, are predicted inside trapped quantum droplets made from binary Bose-Einstein condensates. The paper shows that quantum fluctuations, not ordinary attraction, are what keep these twisted structures stable, and it identifies a double-ring pattern as the experimental signature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability is certified only in the Ψ1=Ψ2 subspace; the binary relative (spin) sector is never analyzed, so stable hopfions in binary gases are not yet supported.","rationale":"The reader's weakest assumption is precisely the load-bearing issue I find. The abstract claims stable hopfions in binary atomic gases, but Eq. (5) restricts the problem to the symmetric manifold; Eq. (9) is a scalar GPE with no relative-phase or relative-density degrees of freedom. A standard BdG decomposition around a symmetric stationary point splits into in-phase and out-of-phase sectors, and only the in-phase sector is controlled by the scalar equation. The out-of-phase equation derived above contains terms proportional to the difference between intra- and inter-species scattering lengths; for the chosen 39K parameters this difference is large, so it is not a negligible correction. The paper's direct simulations, VK criterion plots, and stability labels all refer to the scalar equation, so the reported stability cannot be transferred to the binary system without an additional calculation. I do not claim the relative modes are actually unstable; the problem is that the paper supplies no evidence either way. This is exactly the kind of missing check that a conditional acceptance should require. The geometric Hopf-fibration discussion and the scalar numerics are otherwise coherent, and no independent circularity is involved.","tokens_in":16307,"tokens_out":22016,"duration_ms":254829,"concrete_test":"Compute the Bogoliubov-de Gennes spectrum of the full system (2)-(3) around the S=1, M=1 hopfion at the parameters of Fig. 5(a) (a=100a0, a'=-110a0, R0=0.9 µm, Ω=51000 Hz, N=111130), restricting to the antisymmetric branch η=δΨ1-δΨ2. If any eigenfrequency has Im(ω)>0, the scalar-stable solution is unstable in the binary system. As a complementary dynamical check, evolve the same state in the coupled GPEs with an initial antisymmetric perturbation δΨ1=-δΨ2 at 1% amplitude for at least 10 ms and compare the density profiles to Fig. 5(a1). Reporting the antisymmetric BdG spectrum is the minimal check that would settle the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing gap is the reduction made in Eq. (5). All stationary states and all reported stability runs solve the scalar GPE (9), so stability is verified only for perturbations that preserve Ψ1=Ψ2. In the full binary system, the linearization around a symmetric state has an independent antisymmetric sector. With G11=G22=G, G12=G21=G' and background Ψ1=Ψ2=Ψ/√2, the difference η=δΨ1-δΨ2 satisfies iℏ∂Tη = H0η + [(2G+G')|Ψ|2/2 + Γ|Ψ|3]η + (G-G')Ψ2η*/2, where H0=-(ℏ2/2m)∇2+V. The anomalous term (G-G')Ψ2η*/2 is absent from Eq. (9); for the paper's main parameters a=100a0 and a'=-110a0 it has the large coefficient G-G'=840πℏ2a0/m. No Bogoliubov spectrum or coupled two-component simulation addresses this sector; the Vakhitov-Kolokolov criterion quoted in Figs. 2, 6, and 9 constrains only the scalar branch. If any antisymmetric mode has Im ω>0, the scalar-stable hopfions are unstable in the binary gas even though every scalar calculation is correct. The central claim therefore is conditional on a stability check that the manuscript does not provide.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies toroidally trapped binary Bose-Einstein condensates described by coupled Gross-Pitaevskii equations with cubic mean-field terms and the Lee-Huang-Yang quartic correction. Under the symmetric ansatz Ψ1=Ψ2=Ψ/√2 in Eq. (5), the two-component system is reduced to the scalar GPE in Eq. (9), and stationary hopfion solutions with S=1 and M=0,...,7 are computed by a Newton-conjugate-gradient method starting from the twisted toroidal ansatz in Eq. (14). Stability is assessed through the Vakhitov-Kolokolov criterion and finite-time perturbed evolution, leading to the claims that QH=0 states are stable for inner torus radius above a critical value, that QH=1,...,7 states are partly stable, that they are completely unstable without the LHY term, and that a purely LHY-driven ('LHY superfluid') regime supports at least partly stable hopfions. The paper also presents preimage-based geometric representations of the Hopf number and petal structures.","tokens_in":16597,"tokens_out":18253,"duration_ms":218483,"significance":"If the stability results survive a full two-component analysis, the paper would deliver a new class of three-dimensional topological solitons, hopfion quantum droplets with Hopf numbers 1 through 7, in a realistic binary BEC setup. The identification of the LHY term as the stabilizing nonlinearity and the proposed double-ring experimental signature are valuable. The systematic parameter scans in Figs. 2, 6, and 9 and the geometric illustrations in Figs. 4 and 8 are useful and clearly presented. However, the central claim for a binary gas is conditional on the unanalyzed component-antisymmetric sector, and the topological charge is not independently computed for the converged numerical states, so the significance is currently prospective rather than fully established.","major_comments":[{"comment":"All stationary solutions and all reported stability evolutions solve the scalar GPE (9), but the physical model is the binary system (2)-(3). The manuscript never analyzes perturbations that break the condition Ψ1=Ψ2. The relative mode η=δΨ1-δΨ2 obeys a linearized equation that contains an anomalous term proportional to (G-G')Ψ²η*, which is absent from Eq. (9); for the main parameters a=100a0 and a'=-110a0 this coupling is large. Consequently, the Vakhitov-Kolokolov criterion and the perturbed evolutions in Figs. 2, 6, and 9 certify stability only in the symmetric subspace. A coupled two-component Bogoliubov-de Gennes spectrum, or direct binary-GPE simulations seeded with antisymmetric perturbations, is required before the claim of stable hopfions in binary atomic gases is supported.","section":"Sec. II, Eq. (5); Sec. III"},{"comment":"The ansatz (14) explicitly builds in the winding numbers e^{iMθ} and e^{iSϕ}, so the value QH=MS is effectively imposed by construction. The paper does not compute a Hopf invariant of the converged numerical solution, for example through an integral of the Hopf density or through the linking number of preimage curves obtained from an explicitly defined normalized map to S2. The preimage and petal plots therefore illustrate the labels put in by the ansatz rather than independently establish the topological charge of the numerical state. Please define the map Φ:R3→S2, including the boundary condition at the edge of the trap, and report QH computed directly for each converged solution.","section":"Sec. III, Eq. (14); Figs. 4 and 8"},{"comment":"The stability labels are based on the Vakhitov-Kolokolov criterion, which the paper itself notes is only a necessary condition, and on finite-time perturbed evolution lasting roughly 10-40 ms. Weak instabilities with small growth rates can be missed on these timescales, especially for the higher-QH families that are divided into stable and unstable segments. A Bogoliubov-de Gennes spectrum for representative points in each stable and unstable segment would make the labels much more robust. At minimum, the perturbation amplitudes and a convergence check of the finite-time runs should be reported.","section":"Secs. III A and III B; Figs. 2, 6, 9"},{"comment":"The abstract states that true hopfions with QH=1,...,7 form partly stable families 'including the case of the LHY superfluid', and the conclusion states that when only the LHY term is present, hopfions with QH≠0 can stably exist. However, the only LHY-superfluid results shown in Fig. 9 are for QH=1. No data are presented for QH=2,...,7 in the g=0 regime, so the broader statements in the abstract and conclusion are not supported by the reported numerics. Please either provide those families or qualify the claim to QH=1.","section":"Sec. III B, Fig. 9; Abstract and Conclusion"}],"minor_comments":[{"comment":"The last panel labels read '92267 (f1-f3)', but they should refer to the g-column; please correct this typo.","section":"Fig. 5 caption"},{"comment":"The angle ϕ=arctan((r-r0)/z) is singular at r=r0 and z=0; it would be helpful to state how the numerical initial guess regularizes the vortex core in that region.","section":"Eq. (15)"},{"comment":"The symbol R in the caption should be R0 to match the notation introduced in Sec. II.","section":"Fig. 5 caption"},{"comment":"Reference [33] and reference [87] are the same work and should be merged or renumbered.","section":"References"},{"comment":"The manuscript does not report the numerical box size, grid resolution, or convergence tolerance for the Newton-conjugate-gradient method; these details are needed for reproducibility.","section":"Sec. III"},{"comment":"In the sentence 'true hopfions with S=1, M=1~7, which correspond, accordingly, to QH=1~7', the dependence on S=1 is clear, but it would be safer to state explicitly that QH=MS with S=1.","section":"Abstract"},{"comment":"The terms 'vorticity' and 'winding number' are used interchangeably for M; please define the terminology once at first use.","section":"Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and the scalar numerical results are plausible, but the main claim about binary atomic gases is currently conditional on an unanalyzed relative-mode sector. I would request a coupled two-component stability analysis (Bogoliubov-de Gennes or direct binary simulations) before publication, along with an independent computation of the Hopf number and a correction of the LHY-superfluid claim in the abstract and conclusion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here's my read. The paper finds stable hopfion families in a scalar GPE with cubic attraction, LHY repulsion, and a harmonic toroidal trap, for Hopf numbers up to 7, including the pure-LHY superfluid case. That is new and worth knowing: earlier hopfion scenarios used rotating BECs or engineered repulsive media, not quantum droplets. The numerics look careful in what they actually do, and the geometric preimage/petal analysis is a nice illustrative touch.\n\nThe soft spot is structural and load-bearing. The model in Eqs. (2)-(3) is a binary mixture, and the paper reduces it to the scalar Eq. (9) via the symmetric ansatz (5). Every stationary solution and every stability run is in that symmetric subspace. Nowhere do the authors analyze the relative (spin) modes. The stress-test derivation is correct: for η=δΨ1-δΨ2 the linearized equation contains the anomalous term (G11-G12)Ψ²η*/2, which is absent from the scalar equation and is large for the parameters used (a=100a0, a'=-110a0). That term can destabilize the symmetric solution even when the scalar problem is stable. VK and scalar perturbed evolution do not constrain this sector. So the stated claim—stable hopfions in binary gases—is not yet supported. To be fair, the scalar results are real and self-consistent; the conditional language in the abstract ('families ... are stable') is too strong, but the underlying numerics are not sloppy.\n\nMinor points: no Bogoliubov spectra, but direct evolution plus VK is acceptable for a first pass. No code/data, which is typical but would help. The topology section is standard and the QH construction is not circular in a damaging way—the ansatz just seeds the quantum numbers.\n\nWho is this for? Anyone working on topological solitons in quantum droplets or toroidally trapped BECs. It deserves a serious referee, but the referee should push for a two-component stability analysis, at minimum a Bogoliubov spectrum of the relative mode or coupled simulations. If that check comes out stable, this is a solid publication. If not, the scalar-model version still stands as a useful finding, but the paper must be reframed.\n\nMy recommendation: send to peer review with a request for the binary-sector stability check. I would not desk-reject it.","headline":"Promising scalar-model study of LHY-stabilized hopfions; the binary-condensate stability claim outruns the numerics.","tokens_in":17177,"tokens_out":4055,"would_cite":true,"duration_ms":47388,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.75.Lm","67.85.-d"],"model":"deepseek-v4-flash","headline":"Stable hopfion quantum droplets with Hopf numbers 0 through 7 can exist in toroidally trapped binary Bose gases when the Lee-Huang-Yang quantum-fluctuation term is present.","keywords":["hopfions","quantum droplets","Bose-Einstein condensates","Lee-Huang-Yang term","toroidal trap","Gross-Pitaevskii equations","topological solitons","Hopf number"],"falsifier":"Initialize the coupled two-component GPEs (2)–(3) with a stable symmetric hopfion and add a small antisymmetric perturbation, $\\Psi_1=\\Psi/\\sqrt{2}+\\epsilon$ and $\\Psi_2=\\Psi/\\sqrt{2}-\\epsilon$; if the relative mode grows and breaks the torus on the reported timescales, the binary-mixture stability claim fails even though the scalar numerics are correct.","tokens_in":16112,"feed_emoji":"🌀","tokens_out":11816,"duration_ms":121926,"temperature":0.7,"pith_summary":"This paper aims to show that hopfions—three-dimensional knot-like solitons built as twisted vortex tori—can exist as stable quantum droplets in a toroidally trapped binary Bose–Einstein condensate. Working from the coupled Gross–Pitaevskii equations with cubic mean-field self-attraction, quartic Lee–Huang–Yang repulsion, and a harmonic toroidal trap, the authors report stable families with azimuthal winding $M=0$, twist $S=1$ (Hopf number $Q_H=0$) once the trap's inner radius exceeds $R_0 \\simeq 1.1\\,\\mu$m, and partly stable families of true hopfions with $M=1$ through $7$ and $S=1$ ($Q_H=1$ through $7$). The Lee–Huang–Yang term is the load-bearing ingredient: with only the cubic mean-field nonlinearity every $Q_H \\neq 0$ hopfion is unstable, while a purely LHY nonlinearity still supports partly stable hopfions. If correct, the results give concrete, experimentally accessible signatures—double-ring density patterns and petal-like preimage knots—for observing hopfions in ultracold atomic gases.","feed_headline":"Stable hopfion droplets reach Hopf numbers up to 7","feed_subtitle":"The Lee-Huang-Yang quantum-fluctuation term is what keeps these knot-like solitons stable in trapped binary Bose gases.","key_machinery":"The carrying object is the twisted-torus ansatz $\\varphi = (r')^S \\exp[-(r')^2/A + iM\\theta + iS\\phi]$, in which $M$ winds around the vertical axis and $S$ winds in the poloidal $(r,z)$ plane, so the Hopf number is $Q_H = MS$. Stationary solutions are found by Newton–conjugate-gradient iteration of the dimensionless scalar GPE $-\\frac{1}{2}\\nabla^2\\psi + g|\\psi|^2\\psi + \\gamma|\\psi|^3\\psi + \\frac{1}{2}\\omega[(r-r_0)^2+z^2]\\psi = \\mu\\psi$, and stability is checked by the negative-slope criterion $dN/d\\mu < 0$ together with direct simulation of perturbed evolution. The geometric preimage construction—constant values of $(\\mathrm{Re}\\,\\varphi,\\mathrm{Im}\\,\\varphi)$ mapped back to real space through the Hopf map and stereographic projection—turns the abstract Hopf number into visible linking and petal counts.","core_discovery":"The central discovery is that the scalar Gross–Pitaevskii equation obtained from the binary system by the symmetric reduction $\\Psi_1=\\Psi_2=\\Psi/\\sqrt{2}$ possesses stable toroidal soliton solutions with two independent winding numbers: the azimuthal winding $M$ around the vertical axis and the twist $S$ around the torus core, whose product is the Hopf number $Q_H=MS$. Numerically, with parameters for $^{39}$K, the $Q_H=0$ family ($M=0$, $S=1$) is stable for inner torus radii $R_0 \\gtrsim 1.1\\,\\mu$m and ceases to exist below $R_0 \\approx 0.4\\,\\mu$m; the true hopfions with $S=1$, $M=1$–$7$ have shrinking stability regions as $Q_H$ grows. In the absence of the LHY term ($\\gamma=0$) the negative-slope condition on $N(\\mu)$ is violated and these hopfions are unstable, whereas in the LHY-only case ($g=0$) at least the $Q_H=1$ family is mostly stable. The paper also establishes an elementary geometric reading of the Hopf number: curves of constant $(\\mathrm{Re}\\,\\varphi,\\mathrm{Im}\\,\\varphi)$ are non-intersecting concentric circles for $Q_H=0$, and for $Q_H \\ge 1$ they intersect in $Q_H$ points forming a $Q_H$-petal knot in the $Z=0$ plane.","pith_inferences":["Because the stability analysis is confined to the symmetric subspace $\\Psi_1=\\Psi_2$, the binary-mixture claim is an extrapolation from a single-component reduction; the full two-component equations have not been tested for relative-density or relative-phase instabilities.","Since the preimage linking count equals $Q_H$, a phase- and density-resolved reconstruction could in principle read the Hopf number off the number of petals in the horizontal plane.","The sharp contrast between unstable mean-field-only hopfions and partly stable LHY-supported ones suggests the same stabilization mechanism may apply to other three-dimensional topological solitons in quantum droplets, such as vortex knots and skyrmions in binary mixtures.","A systematic scan of the stability boundary in the $(R_0, N, \\delta a)$ parameter space would produce a concrete, experimentally testable phase diagram for hopfion droplets."],"forward_implications":["Hopfion quantum droplets with Hopf numbers up to $Q_H=7$ should be realizable in toroidally trapped binary Bose gases with $^{39}$K parameters, provided the torus inner radius and atom number are above the reported thresholds.","The double-ring density profile in the $Z=0$ plane gives an experimental fingerprint that distinguishes hopfions from ordinary vortex droplets.","The LHY correction is not a minor quantitative shift but the qualitative stabilizer: removing it destabilizes all nonzero-Hopf families, while keeping only LHY restores partial stability.","Stability windows narrow as $Q_H$ grows, so the highest Hopf numbers will be hardest to observe and the lowest ones the best experimental targets.","Without the toroidal trap the hopfions decay on sub-millisecond timescales, so experimental realization must supply the toroidal confinement."],"supporting_citations":[{"why":"Introduces the idea of stable knot-like solitons, the conceptual starting point for hopfions.","marker":"[2]"},{"why":"Establishes static solitons with nonzero Hopf number and the linking-number interpretation used for the preimages.","marker":"[3]"},{"why":"Shows stable Hopf solitons in a rotating BEC, the earlier setting this paper extends and contrasts with the LHY/toroidal trap.","marker":"[59]"},{"why":"Supplies the Lee-Huang-Yang term and the quantum-droplet balance that the paper's model is built on.","marker":"[61]"},{"why":"Provides prior three-dimensional vortex quantum droplets that motivate the toroidal droplet branch studied here.","marker":"[70]"},{"why":"Demonstrates vortex and multipole quantum droplets in a toroidal potential, the trap geometry underlying the present solutions.","marker":"[78]"},{"why":"Defines the Hopf fibration used to map hopfion wavefunctions to geometric preimages.","marker":"[79]"},{"why":"Gives the negative-slope ($dN/d\\mu<0$) necessary stability criterion applied to whole branches.","marker":"[88]"}],"fun_headline_variants":["Quantum hopfions tie knots up to Hopf 7","Hopfion knots stable up to Hopf 7 in trapped gases","LHY term stabilizes hopfion droplets up to Hopf 7","Knot-like quantum droplets stable up to Hopf 7","Hopfions with twists up to 7 remain stable in traps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that setting the two components equal, $\\Psi_1=\\Psi_2=\\Psi/\\sqrt{2}$, faithfully represents the binary system; the paper never simulates the coupled two-component dynamics, so an antisymmetric perturbation between the components could in principle destroy the hopfions even though the scalar calculation is stable.","fun_headline_variants_meta":{"raw":{"variants":["Quantum hopfions tie knots up to Hopf 7","Hopfion knots stable up to Hopf 7 in trapped gases","LHY term stabilizes hopfion droplets up to Hopf 7","Knot-like quantum droplets stable up to Hopf 7","Hopfions with twists up to 7 remain stable in traps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000809,"raw_usage":{"total_tokens":3713,"prompt_tokens":1270,"completion_tokens":2443,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":886,"completion_tokens_details":{"reasoning_tokens":2351}},"tokens_in":886,"tokens_out":2443,"duration_ms":19669,"temperature":1.0,"reasoning_tokens":2351,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:22:18.834418+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Initialize the coupled two-component GPEs (2)–(3) with a stable symmetric hopfion and add a small antisymmetric perturbation, $\\Psi_1=\\Psi/\\sqrt{2}+\\epsilon$ and $\\Psi_2=\\Psi/\\sqrt{2}-\\epsilon$; if the relative mode grows and breaks the torus on the reported timescales, the binary-mixture stability claim fails even though the scalar numerics are correct.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows stable Hopf solitons in a rotating BEC, the earlier setting this paper extends and contrasts with the LHY/toroidal trap."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Lee-Huang-Yang term and the quantum-droplet balance that the paper's model is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides prior three-dimensional vortex quantum droplets that motivate the toroidal droplet branch studied here."},{"cited_title":"Ferrier-Barbut, H","cited_arxiv_id":null,"evidence_quote":"Demonstrates vortex and multipole quantum droplets in a toroidal potential, the trap geometry underlying the present solutions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Hopf fibration used to map hopfion wavefunctions to geometric preimages."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the negative-slope ($dN/d\\mu<0$) necessary stability criterion applied to whole branches."}],"review_version":1}