{"id":"e3b5825b-c143-4d8a-a310-65260866dc8c","arxiv_id":"2501.02297","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Rotating solitons in self-interacting ultralight dark matter form through a uniform vortex lattice, with a maximum radius about 1.59 times and a maximum rotation rate about 1.34 times the square root of the central density.","lead":"Using analytic theory and 2D simulations, this paper shows that self-interacting ultralight dark matter halos with net angular momentum can relax into a rotating soliton whose spin is carried by a lattice of vortices, not by a single high-angular-momentum quantum state. The result gives concrete limits on soliton radius and spin and predicts observable signatures for this class of dark matter models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability proof for all Ω<Ωmax relies on a finite-disk Bessel spectrum; the soliton's exponential tail invalidates the hard-boundary truncation.","rationale":"The paper's central claim is credible: the analytical TF profile, the vortex-lattice solid-body rotation, and the 2D simulations cohere, and the existence/formation part does not depend on the disputed stability proof. The reader's weakest-assumption identification is the right one. I considered other possible objections—fitted Ω in the profile comparison, no released code/data, and the 3D extrapolation—but these weaken the strength of the numerical check rather than the logic of the central claim. The stability proof is the load-bearing element because the abstract and Sec. IV B claim a stable minimum for all Ω<Ωmax, and the proof's finite-disk Bessel spectrum excludes the very modes that could be unstable. Since the reader's CONDITIONAL verdict already requires strengthening the stability analysis, no additional adjustment is needed.","tokens_in":34993,"tokens_out":18830,"duration_ms":215035,"concrete_test":"Take the numerically relaxed rotating soliton at t≈500 for ϵ=0.01, α=1 (and the ϵ=0.005 and 0.03 runs), add small random perturbations to ψ, and evolve the linearized Gross-Pitaevskii equation in a large box with the true exponential tail; equivalently, diagonalize the Bogoliubov operator about the exact background state and compute complex frequencies without imposing δρ(R)=0. If no mode with Im(ω)<0 is found—especially ℓ=1, ℓ=2, and modes localized beyond RΩ—for Ω<Ωmax, the stability claim survives; if such a mode grows, the bound (70) must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. IV B proves stability by diagonalizing the density Hessian K = 4πΔ^{-1}+λ on a finite disk of radius R>RΩ, imposing δρ(R)=0 for ℓ≥1 and δM(R)=0 for ℓ=0, and assuming the density vanishes identically beyond R. However, App. A (Eq. A11) shows the actual soliton has an exponential WKB tail, so no such hard boundary exists. The eigenvalue equation (69) contains only λ and the domain size, not the background density; therefore choosing a larger effective domain for perturbations that extend into the tail pushes the lowest Bessel eigenvalue κ=(x_1/R)^2 below the Jeans threshold 4π/λ and would reverse the stability sign. The boundary conditions also exclude free-surface modes, including the m=2 bar-like deformation and ℓ=0 surface modes, which are precisely the modes that can destabilize rotating self-gravitating bodies. Thus the claim that every Ω<Ωmax soliton is a stable minimum is not established by the Bessel argument; it depends on an unverified truncation of the perturbation space.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies rotating solitons in two-dimensional ultralight dark matter with repulsive quartic self-interactions, described by a Gross-Pitaevskii equation in the Thomas-Fermi regime. The authors derive a variational continuum-limit configuration at fixed mass and angular momentum, obtaining a solid-body rotation profile, a deformed Bessel density profile, and a uniform vortex lattice; they also derive upper bounds on the rotation rate and radius (Eqs. (64), (65)) and argue, via a second-variation analysis, that all such configurations are dynamically stable. Numerical simulations from stochastic initial conditions with injected angular momentum show the formation of a rotating soliton with a regular lattice of single-charge vortices, and the measured density, velocity, vortex count, and angular momentum profiles are compared with the analytical expressions.","tokens_in":35219,"tokens_out":5776,"duration_ms":64578,"significance":"If the formation and stability claims hold, the paper provides a concrete and physically interesting picture for rotating solitons in self-interacting ultralight dark matter, connecting the vortex-lattice physics of rotating Bose-Einstein condensates with cosmological soliton dynamics. The analytical parts are substantial: the Thomas-Fermi profiles, the matched asymptotic expansions for the soliton boundary layer and vortex core, and the energy decomposition of vortex excess energy are worked out in detail and are internally consistent. The numerical simulations are clean and the diagnostics (density, phase winding, velocity field, vortex trajectories) are well chosen to exhibit the claimed phenomenology. The paper also makes falsifiable predictions, notably the scaling Ωmax ≃ 1.343√ρ0 and the radius bound Rmax ≃ 1.593 R0, which are stated in a parameter-free form. However, the dynamical-stability proof has a technical gap, and several numerical comparisons are partly circular because the rotation rate is measured from the same simulation that is then compared with the formulas.","major_comments":[{"comment":"The stability proof for all Ω < Ωmax relies on a finite-disk truncation: eigenvectors are Bessel functions with hard-wall conditions δρ(R)=0 (ℓ ≥ 1) or δM(R)=0 (ℓ = 0) at a radius R > RΩ, justified by the statement that the density 'identically vanishes' beyond RΩ. This is inconsistent with App. A, Eq. (A11), which shows that the soliton has an exponential WKB tail that never vanishes identically. The eigenvalue equation (69) contains only κ, ℓ, and the domain size, not the background density; therefore a perturbation that extends into the tail can effectively see a larger R, pushing the lowest eigenvalue κ = (x_ℓ1/R)^2 below the Jeans threshold 4π/λ and reversing the sign of ν. The hard-wall boundary conditions also exclude free-surface modes, including ℓ=0 compression modes and ℓ=2 bar-like deformations, which are the modes most likely to destabilize a rotating self-gravitating body. Consequently, the claim that every Ω < Ωmax soliton is a dynamically stable minimum is not established by the present argument; the proof needs a treatment of the exponential tail or a direct variational test of tail and free-surface modes.","section":"Sec. IV B, Eqs. (68)-(70)"},{"comment":"The numerical confirmation is partly circular. In Fig. 3, Ω is obtained by a least-squares fit to the simulated transverse velocity profile, and this same Ω is then inserted into Eqs. (60), (63), and (73) to produce the red curves in Figs. 2 and 5. This demonstrates that the simulated configuration is internally consistent with the ansatz of a solid-body-rotating Thomas-Fermi soliton, but it does not independently confirm the predictive content of the derivation. A stronger test would predict Ω from the initial conditions, for example through the mass-shell angular-momentum estimate used in Sec. V C 3, and then compare the resulting density, vortex number, and angular-momentum profiles with the simulation output. Alternatively, the measured Lz and Nv should be compared with relations that do not share the fitted Ω. The abstract and conclusion currently state that the numerical results agree with the analytical derivations, which overstates the strength of the available evidence.","section":"Sec. V C, Figs. 2, 3, and 5"},{"comment":"The sentence following Eq. (38) states that the uniform-lattice energy is 'about a quarter of the single vortex energy (35), for Nv = |σ|'. For large Nv and ln(R0/ξ) ≫ 1, Eq. (38) is dominated by the term Nv ρ0 ε² ln(R0/ξ), which is smaller than the single-vortex energy Nv² ρ0 ε² ln(R0/ξ) by a factor of order 1/Nv, not by a factor of four. The qualitative conclusion that a high-spin vortex splits into unit vortices is not affected, but the quantitative statement should be corrected or clarified.","section":"Sec. III B, Eq. (38)"}],"minor_comments":[{"comment":"The text says 'at the later times, 84 < t < 84' but the intended interval is evidently 84 < t < 85; this typo appears in the discussion of the late-time vortex trajectories.","section":"Sec. V C 5"},{"comment":"The matched asymptotic vortex profile is derived under the assumption |σ| ≫ 1, and the authors note that for |σ| = 1 the normalization and the transition shape are known only up to a factor of order unity. Since all vortices in the simulations have |σ| = 1, the quantitative vortex-profile predictions of Sec. III A 2 are not directly verified analytically for the simulated case; this limitation should be acknowledged where Eq. (35) and the related energy estimates are used.","section":"App. B"},{"comment":"The anisotropic initial condition in Eq. (99) uses the discontinuous function sign(L) at L = 0, which produces a phase-space distribution that is discontinuous across the L = 0 plane. This is a valid modeling choice for injecting net angular momentum, but its effect on the initial vortex population and on the subsequent relaxation could be commented on, since sign(L) is not smooth.","section":"Sec. V A 3"},{"comment":"The statement that 'most of these properties' extend to 3D with vortex rings is a conjecture rather than a result of this paper; the distinction should be made explicit, since the 2D logarithmic Green function and the stability analysis do not automatically carry over to 3D point-vortex or vortex-ring systems.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for the journal and the core phenomenology—vortex-lattice formation and solid-body rotation in the Thomas-Fermi regime—appears solid and well illustrated. The main technical weakness is the stability proof in Sec. IV B, which is load-bearing for the abstract's claim of dynamical stability; it should be either fixed with a tail-aware argument or substantially softened. The numerical comparisons should also be reframed or strengthened so that they are not described as independent confirmations when Ω is measured from the same simulation. I recommend major revision rather than rejection, because the variational construction, the matched asymptotic expansions, and the numerical demonstration of vortex-lattice formation are valuable and likely correct in their main conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper delivers a genuinely new result—rotating self-gravitating solitons in the Thomas-Fermi regime of self-interacting ULDM are not large-ℓ eigenstates but solid-body rotators built from a uniform lattice of quantized vortices. The profile (60), the maximum radius Rmax ≈ 1.59 R0, and the maximum rotation rate Ωmax ≈ 1.34 sqrt(ρ0) are new and likely to be cited. The asymptotic matching in the appendices (Painlevé II boundary layers for the soliton edge and vortex core) is careful and technically impressive.\n\nThe numerical simulations from stochastic initial conditions are a real plus. The vortex lattice genuinely forms, the density profile matches (60) with a flat center, and the vortex density is uniform as predicted. The authors are honest that the profile comparison uses Ω measured from the same simulation; this makes the agreement a consistency check rather than a prediction. The scaling with ϵ and α gives confidence that the mechanism is right.\n\nThe soft spots are in proportion. The stability proof in Sec. IV B is the weakest link. It diagonalizes K = 4πΔ^{-1}+λ on a finite disk with boundary conditions at R > RΩ, assuming the density vanishes identically beyond R. But App. A shows the density has an exponential tail, so the hard boundary is an artificial truncation. The eigenvalue equation does not contain the background density at all; if you enlarge the domain to admit perturbations that penetrate the tail, the lowest Bessel eigenvalue falls below the Jeans threshold and the quadratic form goes negative. The proof therefore establishes stability only within a restricted perturbation class, not for all Ω < Ωmax as claimed. This is a technical gap, not a demonstrated instability—the simulations run many dynamical times without trouble—but the claim about stability should be softened or the tail modes addressed, e.g., by including quantum pressure in the Hessian for the tail region.\n\nThe 3D generalization is speculative, but the paper says so. No code or data is released, which would help if this becomes a reference result.\n\nWho is this for? People working on self-interacting ULDM solitons, vortex formation in dark matter, and rotation curves in these models. It deserves a serious referee—the core mechanism is plausible, the analysis is mostly clean, and the bounds are useful—but the referee should press on the stability argument and ask for a genuine prediction of Ω.","headline":"A credible 2D demonstration that rotating solitons in self-interacting ULDM are supported by vortex lattices, with useful new bounds, though the stability proof and numerical tests are weaker than the narrative suggests.","tokens_in":35727,"tokens_out":7907,"would_cite":true,"duration_ms":84950,"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":"In self-interacting ultralight dark matter, a halo with nonzero spin collapses into a rotating soliton whose rotation is carried by a uniform lattice of quantum vortices, yielding solid-body motion and a maximum spin set by the central…","keywords":["ultralight dark matter","Gross-Pitaevskii equation","Thomas-Fermi regime","quantum vortices","vortex lattice","rotating soliton","solid-body rotation","angular momentum"],"falsifier":"Solve the second-variation eigenvalue problem for the full solution including the exponential tail, without imposing δρ(R)=0 or δM(R)=0 at a finite radius, and look for a negative eigenvalue; a single negative mode would overturn the claim that all Ω < Ωmax solitons are stable.","tokens_in":34761,"feed_emoji":"🌀","tokens_out":8628,"duration_ms":71273,"temperature":0.7,"pith_summary":"The paper studies halo-sized clouds of ultralight dark matter with repulsive self-interactions, governed by the Gross-Pitaevskii equation, and asks what happens when the cloud has net angular momentum. In 2D simulations starting from stochastic initial conditions, a rotating soliton forms in a few dynamical times. The key claim is that the rotation is not due to a large orbital quantum number of the wave function; it is produced by a regular lattice of quantized vortices that, in the continuum limit, becomes solid-body rotation. These rotating solitons are stable minima of the energy at fixed mass and angular momentum, and for a given central density they cannot rotate faster than about 1.343√ρ0 nor be larger than about 1.593 times the static radius. The analytical profiles, vortex densities, and stability bounds match the simulations.","feed_headline":"Dark-matter solitons rotate via vortex lattices, not wave whirls","feed_subtitle":"A stable rotating core grows from random initial conditions, with spin capped by the central density.","key_machinery":"The key machinery is the multi-vortex wave-function ansatz ψ = $e^{{-iμt/ϵ}}$ √ρ0 ∏_j f(r − r_j) $e^{{iσ_j θ_j}}$, which treats each vortex as a phase singularity carrying circulation 2πϵσ_j. In the limit ϵ→0 the discrete vortices become a smooth vorticity distribution, and extremizing E − μM − ΩL_z yields the solid-body rotation and the balance Φ_N + Φ_I − r²Ω²/2 = μ. The Bessel-function profile and the requirement of finite mass produce the bounds (64)–(65), and the second variation of the energy at fixed angular momentum gives stability.","core_discovery":"The central discovery is that the steady rotating state of a self-gravitating Gross-Pitaevskii condensate in the Thomas-Fermi regime is a uniform lattice of unit-circulation vortices embedded in an axisymmetric soliton. In the continuum limit, where the de Broglie wavelength is much smaller than the system size, the lattice is equivalent to a smooth vorticity field and the velocity becomes solid-body rotation, v = Ω r eθ. The density profile is ρ(r) = (ρ0 − Ω²/2π) J0(z0 r/R0) + Ω²/2π, which is broader than the static Bessel profile; requiring finite mass gives Ωmax ≈ 1.343√ρ0 and Rmax ≈ 1.593 R0. A second-variation calculation shows all such configurations with Ω below Ωmax are dynamically stable. Numerical simulations with three values of the de Broglie parameter and several initial anisotropies confirm the deformed profile, the solid-body rotation, the uniform vortex density, and the vortices moving on circular orbits.","pith_inferences":["Not explored in the paper: the stability proof truncates the density at a finite radius, so a numerical eigenvalue solve on the true exponential tail would test whether any mode living outside that radius is unstable; if one were found, the claim that all Ω < Ωmax solitons are stable would need revision.","An extension the authors leave implicit is that in 3D the point vortices become vortex rings; tracking whether such rings form and stay circular in a 3D simulation would show whether the 2D solid-body picture carries over.","A testable corollary not drawn by the authors: if an observed rotation curve of a dark-matter-dominated dwarf implies a soliton with spin above Ωmax for its inferred central density, this model would be ruled out for that halo."],"forward_implications":["The density profile of a rotating soliton is wider than the static one, so rotation-curve fits that ignore rotation will overestimate the central density or underestimate the radius.","The vortex lattice provides a concrete mechanism for a central soliton to retain a large fraction of the halo's angular momentum instead of expelling it.","There is an upper bound on angular momentum for a given central density; halos with more spin must leave the excess in the outer envelope or form a different configuration.","In 3D, point vortices should become vortex rings, making the prediction testable in full 3D simulations of scalar-field dark matter."],"supporting_citations":[{"why":"gives the critical rotation rate above which creating a vortex lowers the energy in the Thomas-Fermi regime, the baseline this paper extends.","marker":"[32]"},{"why":"derives the Thomas-Fermi soliton hydrostatic equilibrium that the rotating profile generalizes.","marker":"[33]"},{"why":"provides the stochastic initial-condition construction and soliton-formation method used for the 2D simulations here.","marker":"[41]"},{"why":"establishes that vortices from wave interferences have unit winding and dominate the outer halo.","marker":"[53]"},{"why":"supplies the multi-vortex ansatz and the result that each vortex follows the flow of the others.","marker":"[74]"},{"why":"derives the equations of motion for vortices from the action, used to show vortices follow the matter velocity.","marker":"[75]"},{"why":"introduces the Feynman argument that a uniform vortex lattice produces solid-body rotation in a rotating superfluid.","marker":"[78]"}],"fun_headline_variants":["Vortex lattices spin dark matter solitons","Stable spin from vortex lattices in dark matter cores","Max spin for dark matter solitons set by vortex lattices","Rotating dark matter solitons are vortex lattices","Vortex lattices drive rotation of dark matter solitons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The stability proof assumes the soliton density drops exactly to zero at a finite boundary radius, so it does not examine modes living in the exponential tail of the density profile.","fun_headline_variants_meta":{"raw":{"variants":["Vortex lattices spin dark matter solitons","Stable spin from vortex lattices in dark matter cores","Max spin for dark matter solitons set by vortex lattices","Rotating dark matter solitons are vortex lattices","Vortex lattices drive rotation of dark matter solitons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000987,"raw_usage":{"total_tokens":4218,"prompt_tokens":1012,"completion_tokens":3206,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":628,"completion_tokens_details":{"reasoning_tokens":3123}},"tokens_in":628,"tokens_out":3206,"duration_ms":22492,"temperature":1.0,"reasoning_tokens":3123,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:14:58.061824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the second-variation eigenvalue problem for the full solution including the exponential tail, without imposing δρ(R)=0 or δM(R)=0 at a finite radius, and look for a negative eigenvalue; a single negative mode would overturn the claim that all Ω < Ωmax solitons are stable.","supporting_citations":[{"cited_title":"Ultra-light dark matter explanation of NANOGrav observations","cited_arxiv_id":"2311.10148","evidence_quote":"establishes that vortices from wave interferences have unit winding and dominate the outer halo."}],"review_version":1}