{"id":"aecaf009-247b-4bc6-949b-81561f5bc20b","arxiv_id":"2505.04079","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The vorticity structure factor in phase-separating superfluid mixtures follows a dynamic scaling law with scale exponent 1, producing a vortex energy decay K ~ t^{-2/3} ln t.","lead":"This paper uses computer simulations to show that the swirling flow, or vorticity, in phase-separating superfluid mixtures becomes statistically similar at different times when rescaled by the growing domain size. The result yields a new decay law for the energy of that flow, which can be tested in cold-atom experiments on lithium gases.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The t^{-2/3} ln t decay law rests on the unmeasured exponent b=-1; the paper itself says b is set by physical interpretation and only partially explains the data, and any b<-1 removes the log correction.","rationale":"The reader's weakest assumption is the same one I would flag. The alpha=1 scaling of S_v has independent support from the collapse in Fig. 2(d) and the near-constancy of E0, so it is not the softest point. The b=-1 assumption is explicitly acknowledged as interpretation, controls the logarithm in Eq. (12), and has a direct numerical test. Because the paper is otherwise a plausible numerical discovery and the concern is checkable rather than demonstrably wrong, I would keep the CONDITIONAL verdict rather than reject. The condition should be: show sensitivity of K(t) to b and report the measured b with uncertainty.","tokens_in":9176,"tokens_out":8851,"duration_ms":92148,"concrete_test":"Recompute the vortex energy from Eq. (10) with the same fitted E0, a, and measured l(t), but replace the imposed b=-1 with the actual late-time slope measured in 4<tilde k<20, including its uncertainty (e.g., b=-1.0±0.1 and b=-1.3). If K(t) is better described by a pure power law t^{-2/3} and the log factor appears only for the single imposed value b=-1, then the decay-law claim is not robust. Report the measured b and its error at l/xi >= 30 from the same data used in Fig. 3(a).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central analytic result, Eq. (12), K ~ t^{-2/3} ln t, is obtained by setting b=-1 in the ansatz (9). The log factor is a direct consequence of the high-k integral being effectively ∫ dk/k for b=-1. If the true late-time microscopic slope is b=-1.1 or closer to -2, the integral saturates and K(t) ~ t^{-2/3} without the logarithmic factor (for b<-1), so the claimed distinction from freely decaying turbulence disappears. The paper does not measure b: after Eq. (12) it states 'the exponent b is determined based on a physical interpretation,' and the measured slope 'only partially explains the time evolution of slope in Fig. 3(a).' Appendix D simply imposes b=-1. This is load-bearing because the abstract advertises a decay law, and Fig. 3(b) validates it with this imposed value; a sensitivity test is therefore required before the log-corrected law is accepted.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the vorticity structure factor in phase-separating binary superfluid mixtures using two-dimensional Gross-Pitaevskii simulations. The authors observe that the peak of the energy spectrum E(k,t) is nearly time-independent in the inertial hydrodynamic stage, which implies a scale exponent α=1 in the proposed dynamic scaling form S_v(k,t) l(t)^α = \\tilde S_v(kl). They demonstrate a data collapse for the vorticity structure factor and the energy spectrum in the late stage. They further introduce an interpolating ansatz for the energy spectrum with exponents a and b, set b=-1 from the assumption of an uncorrelated point-vortex distribution, and derive a decay law K(t) ~ t^{-2/3} ln t for the vortex energy, which they compare with numerical data.","tokens_in":9472,"tokens_out":12331,"duration_ms":112218,"significance":"If the scaling law holds, this is an interesting extension of dynamic scaling to vorticity in phase-separating superfluids, connecting phase-ordering kinetics with turbulence theory and directly relevant to recent experiments on 7Li BECs. The numerical collision of the structure factor and the connection to known spectral exponents are the paper's main strengths. The proposed decay law is falsifiable but currently rests on the unmeasured exponent b=-1; the logarithmic correction is sensitive to that value. The manuscript also contains a normalization inconsistency in Eq. (9) that affects the quantitative comparison. The core scaling-law result is plausible and worth publishing after the quantitative issues are resolved.","major_comments":[{"comment":"The exponent b is imposed as -1 rather than measured from the late-time spectra. The logarithmic correction in Eq. (12) follows only for b=-1; for any b<-1 the high-wavenumber integral saturates and K(t) ~ t^{-2/3} without the log factor. Because the theoretical curve in Fig. 3(b) uses b=-1, the comparison does not independently test the log-corrected decay law. Please measure b from the late-stage slope in the microscopic regime (e.g., the local exponent in 4 < \\tilde k < 20 at the largest times) and provide a sensitivity test for b near -1. If b is not accurately determined, the abstract and summary should state the decay law as conditional on this exponent.","section":"Energy Decay Law / Eq. (12) / Appendix D"},{"comment":"The ansatz (9) is described as taking its maximum value \\tilde E_0, but for the quoted fit a=2.35 and b=-1, the maximum of 2 \\tilde k^a/(1+\\tilde k^{a-b}) is about 1.09 \\tilde E_0 at \\tilde k ≈ 1.29, not \\tilde E_0. Since Appendix D sets \\tilde E_0 to the observed max(\\tilde E), the ansatz normalization is inconsistent with the data. This affects the prefactor in Eq. (10) and the quantitative comparison in Fig. 3(b). Please correct the normalization of Eq. (9) or redefine \\tilde E_0 accordingly, and re-evaluate the agreement.","section":"Dynamic Scaling Plots / Eq. (9) / Appendix D"}],"minor_comments":[{"comment":"The Fig. 1 caption states a system size of 512ξ on each side, while Appendix A gives L/Δx = 4096 and Δx/ξ = 0.5, implying L = 2048ξ. Please reconcile these values.","section":"Fig. 1 caption / Appendix A"},{"comment":"The quantity K(t) is called the vortex energy, but the system conserves total energy; K(t) is the spectral energy in the wavenumber band [1/L, 1/ξ] and decays because the spectral peak moves to lower wavenumbers. Please clarify this distinction in the main text.","section":"Energy Decay Law"},{"comment":"The collapse in Fig. 2(d) is asserted visually. A quantitative measure of the spread of the rescaled curves (e.g., a residual or a collapse metric) over the claimed scaling window would strengthen the evidence for the α=1 law.","section":"Dynamic Scaling Plots"},{"comment":"The sentence 'This prediction partially explains the time evolution of slope in Fig. 3(a)' is vague. Please specify which features of the measured slope evolution are captured by the picture of a crossover from vortex sheets to point vortices, and which are not.","section":"Energy Decay Law / Fig. 3(a)"},{"comment":"The late-time discrepancy in Fig. 3(b) is attributed to noise accumulation at the resolution scale. Please provide supporting evidence, such as a plot of the spectrum at high k becoming flat or white-noise-like at the final times.","section":"Fig. 3(b) and Summary"}],"recommendation":"major_revision","confidential_remarks":"The normalization error in Eq. (9) is easily verified and should be fixed. The larger scientific risk is that the decay law (12) depends critically on the imposed b=-1, which is not measured; a sensitivity analysis or a direct measurement of b is needed before the log-corrected decay law can be accepted. The dynamic scaling law itself appears plausible and is the most valuable contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the target quantity, not the machinery. Prior work on dynamic scaling in these mixtures focused on the magnetization structure factor; this paper studies the vorticity structure factor and its associated energy spectrum. The numerical observation that the rescaled energy spectrum collapses with alpha = 1 in the late stage is plausible and well supported by Fig. 2(b,d). The energy-decay law K ~ t^{-2/3} ln t is a concrete, falsifiable prediction that distinguishes this system from ordinary free decay, and it is clearly presented. The paper is honest about its limitations: it notes early-stage data do not collapse, it notes the microscopic exponent b is assigned by physical interpretation rather than measured, and it attributes the late-time discrepancy to noise accumulation. Those self-critical statements are worth taking at face value.\n\nThe soft spot is the log correction. Equation (12) depends on setting b = -1, which the paper justifies by invoking an uncorrelated distribution of point vortices. But the measured slope approaches -1 only gradually, and the text admits the prediction 'partially explains' the data. If b is actually less than -1, the integral saturates and the log factor disappears, leaving an ordinary t^{-2/3} law that would not distinguish this from freely decaying turbulence. The authors do not provide a sensitivity analysis. This is a load-bearing weakness, but not a fatal one, because the scaling collapse itself and the qualitative decay trend do not hinge on the precise value of b. What is missing is a statement of how robust Eq. (12) is to b in [-1.2, -0.8] or an attempt to measure b more precisely in the late stage.\n\nThe citation pattern looks fair. The paper builds on the authors' own prior work, but the relevant prior results are there and the experiment [17] is cited. No obvious stretching. The absence of code or data archive is a minor complaint for a Letter, though the community would benefit from it.\n\nWho should read this: anyone working on phase-ordering kinetics in superfluids, especially the 7Li experiments. It deserves a serious referee — not a desk reject — because the new quantity, the scaling observation, and the analytic decay law are worth checking, and the sensitivity issue can be fixed in revision. My own recommendation is conditional acceptance with the request to add a sensitivity analysis for b and a clearer statement of the uncertainty in the log factor.","headline":"A credible first look at vorticity scaling in phase-separating superfluids, with the advertised t^{-2/3} ln t decay resting on an imposed exponent rather than a measured one.","tokens_in":9895,"tokens_out":617,"would_cite":true,"duration_ms":7897,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that the vorticity structure factor in phase-separating superfluid mixtures obeys a dynamic scaling law with exponent $\\alpha=1$, and that the vortex energy decays as $K(t)\\sim t^{-2/3}\\ln t$, linking phase-ordering…","keywords":["dynamic scaling","vorticity structure factor","phase separation","superfluid mixtures","energy spectrum","quantum turbulence","Gross-Pitaevskii","Bose-Einstein condensate"],"falsifier":"Measure the high-wavenumber slope $b$ of the energy spectrum at $l/\\xi$ well beyond 30 in a higher-resolution simulation, or extract the vortex energy decay from a $^7$Li phase-separation experiment over a full decade of time: if $|b|$ saturates near 2 instead of 1, or the energy decay is better fit by a pure power law $t^{-p}$ with no logarithmic factor, the central scaling and decay predictions fail.","tokens_in":9007,"feed_emoji":"🌀","tokens_out":7126,"duration_ms":66439,"temperature":0.7,"pith_summary":"This paper aims to show that vorticity, not just the magnetization domains, in a phase-separating two-component superfluid follows a dynamic scaling law in the late, inertia-dominated stage. The authors find that the vorticity structure factor $S_v(k,t)$ collapses onto a single time-independent curve when scaled by the characteristic domain size $l(t)$ with exponent $\\alpha=1$. They further derive that the total vortex energy decays as $K(t)\\sim t^{-2/3}\\ln t$, a law that is distinguishable from ordinary freely decaying turbulence. Because the system they simulate is the same regime realized in recent $^7$Li Bose-Einstein condensate experiments, the prediction is directly testable. The result connects phase-ordering kinetics—usually studied through order-parameter structure factors—with fluid-turbulence statistics.","feed_headline":"Superfluid vorticity obeys one universal scaling curve","feed_subtitle":"A single exponent rescales the vorticity spectrum, linking phase-ordering kinetics to turbulence and predicting a concrete energy decay.","key_machinery":"The central object is the vorticity structure factor $S_v(k,t)=\\langle|\\hat{\\omega}(k,t)|^2\\rangle$, the angular-averaged power spectrum of the two-dimensional vorticity field, together with its companion energy spectrum $E(k,t)=L^2 S_v/(4\\pi k)$ borrowed from fluid turbulence. The argument runs through the dynamic hierarchy hypothesis: on scales larger than the domain size $l$ the spectrum has one power-law exponent $a$, and on scales between $l$ and the healing length $\\xi$ it has another exponent $b$. The load-bearing ansatz is $\\widetilde{E}(\\tilde{k})=2\\widetilde{E}_0\\tilde{k}^{\\,a}/(1+\\tilde{k}^{\\,a-b})$ for the universal function, which interpolates between the two regimes, keeps the peak pinned at $\\widetilde{E}_0$, and lets the decay law be computed analytically.","core_discovery":"Using numerical simulations of the Gross-Pitaevskii model, the paper establishes that in the inertial hydrodynamic stage ($l/\\xi\\gtrsim 30$) the vorticity structure factor satisfies the dynamic scaling law $S_v(k,t) l(t)^{\\alpha}$ equals a universal function of $kl$ with $\\alpha=1$. The central observation is that the peak $E_0$ of the energy spectrum $E(k,t)=L^2 S_v(k,t)/(4\\pi k)$ stays constant in time, so the rescaled spectrum collapses onto one curve. The universal spectrum has two power-law regimes: macroscopic scales behave as $\\tilde{k}^{\\,a}$ with $a=2.35(4)$, while microscopic scales approach $\\tilde{k}^{-1}$, the signature of an uncorrelated distribution of point vortices. Combining these regimes through the ansatz $\\widetilde{E}(\\tilde{k})=2\\widetilde{E}_0\\tilde{k}^{\\,a}/(1+\\tilde{k}^{\\,a-b})$ with $b=-1$ yields the vortex energy decay $K(t)\\sim t^{-2/3}\\ln t$, which the numerics support in the late stage.","pith_inferences":["If the $\\alpha=1$ collapse is generic, the same rescaling should hold in other immiscible two-component BECs and in binary classical fluids with vortex-sheet-like interfaces; looking for a log-corrected energy decay in those systems would test the universality class.","The $b=-1$ assumption could be probed by measuring the microscopic-regime slope at larger $l/\\xi$ than simulated; a gradual drift toward $-2$ at higher resolution would indicate that vortex sheets, not point vortices, dominate the smallest scales and would erase the logarithmic factor.","A practical consequence for cold-atom experiments: the constancy of the spectral peak $E_0$ could be used as a calibration-free way to identify the inertial hydrodynamic stage, even when the domain-wall length is difficult to extract from absorption images.","Connecting to turbulence theory, the logarithmic correction may be the observable remnant of the finite healing-length cutoff; experiments with different interaction strengths (hence different $\\xi$) could vary the crossover window in which the $t^{-2/3}\\ln t$ law is visible."],"forward_implications":["The vorticity statistics of a phase-separating superfluid mixture become statistically self-similar in the late stage, so the full two-point vorticity distribution is determined by the domain-size length $l(t)$ up to one scale exponent $\\alpha=1$.","The peak of the energy spectrum staying constant gives an experimental marker for the inertial hydrodynamic stage that does not require measuring the domain-wall length.","The vortex energy decays as $K(t)\\sim t^{-2/3}\\ln t$, which is slower than the power laws $K(t)\\sim t^{-1}$ to $t^{-1.2}$ of freely decaying classical turbulence, providing a clean distinguishing prediction.","Late-stage small-scale vorticity behaves as an uncorrelated set of point vortices, giving $E(k)\\sim k^{-1}$; if vortex sheets dominated instead, the slope would be $-2$.","The theory applies directly to $^7$Li ferromagnetic condensates under strong quadratic Zeeman shift, where the predicted scaling can be checked against the observed $1/z=2/3$ phase-separation dynamics."],"supporting_citations":[{"why":"Reports the $^7$Li experiment whose phase-separation dynamics define the inertial hydrodynamic stage ($1/z=2/3$) that this theory targets.","marker":"[17]"},{"why":"Introduces the dynamic hierarchy hypothesis for the microscopic regime of phase-separating superfluids that the paper extends to vorticity spectra.","marker":"[13]"},{"why":"Supplies the freely decaying turbulence framework and the $K(t)\\sim t^{-p}$ decay law that the derived $t^{-2/3}\\ln t$ law is compared against.","marker":"[36]"},{"why":"Gives the constraint on the macroscopic exponent $a$ and the classical turbulent energy-spectrum analysis used to formulate the ansatz.","marker":"[38]"},{"why":"Establishes quantum Kelvin-Helmholtz instability at domain walls as the mechanism generating vortices, justifying the late-stage point-vortex picture.","marker":"[31]"},{"why":"Provides the result that an uncorrelated distribution of point vortices gives a $k^{-1}$ energy-spectrum slope, the basis for setting $b=-1$.","marker":"[55]"},{"why":"Determines $g_{+-}/g=3$ as the interaction ratio realized in $^7$Li gases, fixing the parameter regime of the simulations.","marker":"[47]"}],"fun_headline_variants":["Universal vorticity scaling unites phase-ordering and turbulence","Superfluid vorticity collapses onto one universal curve","Single exponent rescales vorticity spectrum in superfluids","Dynamic scaling law predicts vortex energy decay","Vorticity structure factor shows one scaling curve"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that in the late stage the high-wavenumber vorticity is an uncorrelated distribution of point vortices, so the microscopic spectral exponent is exactly $b=-1$; the numerics only show the slope gradually approaching $-1$, and if it were closer to $-2$ the derived logarithmic correction in the decay law would disappear.","fun_headline_variants_meta":{"raw":{"variants":["Universal vorticity scaling unites phase-ordering and turbulence","Superfluid vorticity collapses onto one universal curve","Single exponent rescales vorticity spectrum in superfluids","Dynamic scaling law predicts vortex energy decay","Vorticity structure factor shows one scaling curve"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000203,"raw_usage":{"total_tokens":1378,"prompt_tokens":931,"completion_tokens":447,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":375}},"tokens_in":547,"tokens_out":447,"duration_ms":5259,"temperature":1.0,"reasoning_tokens":375,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:38:21.051888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the high-wavenumber slope $b$ of the energy spectrum at $l/\\xi$ well beyond 30 in a higher-resolution simulation, or extract the vortex energy decay from a $^7$Li phase-separation experiment over a full decade of time: if $|b|$ saturates near 2 instead of 1, or the energy decay is better fit by a pure power law $t^{-p}$ with no logarithmic factor, the central scaling and decay predictions fail.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the $^7$Li experiment whose phase-separation dynamics define the inertial hydrodynamic stage ($1/z=2/3$) that this theory targets."},{"cited_title":"Takeuchi, Phys","cited_arxiv_id":null,"evidence_quote":"Introduces the dynamic hierarchy hypothesis for the microscopic regime of phase-separating superfluids that the paper extends to vorticity spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the freely decaying turbulence framework and the $K(t)\\sim t^{-p}$ decay law that the derived $t^{-2/3}\\ln t$ law is compared against."},{"cited_title":"Kida, Phys","cited_arxiv_id":null,"evidence_quote":"Gives the constraint on the macroscopic exponent $a$ and the classical turbulent energy-spectrum analysis used to formulate the ansatz."},{"cited_title":"Takeuchi, Phys","cited_arxiv_id":null,"evidence_quote":"Establishes quantum Kelvin-Helmholtz instability at domain walls as the mechanism generating vortices, justifying the late-stage point-vortex picture."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the result that an uncorrelated distribution of point vortices gives a $k^{-1}$ energy-spectrum slope, the basis for setting $b=-1$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Determines $g_{+-}/g=3$ as the interaction ratio realized in $^7$Li gases, fixing the parameter regime of the simulations."}],"review_version":1}