{"id":"d3df0018-7ea7-4762-82fb-ba069d9d6a0b","arxiv_id":"2412.03544","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"TDGL quench simulations of chiral 3He films reproduce Kibble-Zurek vortex scaling and reveal a post-freeze-out population asymmetry between triangular and crescent vortex cores.","lead":"Numerical simulations of a temperature quench in superfluid helium-3 films show vortices and domain walls forming, with the vortex density obeying the Kibble-Zurek scaling law. A smart generalist should read this for testable predictions about chiral superfluids, including an unusual asymmetry between two types of vortex cores.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported δ_KZM=0.436±0.003 rests on a single freeze-out criterion, although the cited Antunes et al. find KZ exponents varying by ~0.14 with criterion; without a sensitivity test the headline scaling exponents may be definition-dependent.","rationale":"The reader’s weakest assumption — that the TDGL/Langevin model remains valid down to 0.5T_c — is legitimate, but it is broad and not easily settled by one targeted check; moreover, for slow quenches the fitted freeze-out times occur near T_c, so the 0.5T_c issue mainly affects saturated fast quenches and late-time asymmetry, not the core scaling fit. A more internal and immediately testable threat is the freeze-out definition. The paper itself cites Antunes et al. reporting that the extracted KZ exponent depends on the freeze-out criterion, with a spread from 0.43 to 0.57. The paper uses one criterion, reports δ_KZM = 0.436 ± 0.003, and counts vortices at the freeze-out time defined by that criterion. Hence the headline exponents and the cumulant slopes in Fig. 2 are not demonstrated to be independent of this choice. If the criterion sensitivity seen in Ref. 29 applies here, the small deviations among β_1, β_2, β_3 — the claimed evidence for deviations from independent-defect full-counting-statistics models — could also change. This is a concrete, falsifiable objection: rerun with alternative freeze-out criteria and compare. If the exponents are stable, the central KZM claim is strengthened; if not, conditional acceptance should require reporting that sensitivity. The generous 20,000-realization sampling and the strong-coupling GL foundation are positive features, but they do not remove the criterion-dependence risk. This concern does not move the verdict away from CONDITIONAL; it sharpens the condition that should be met.","tokens_in":18276,"tokens_out":9314,"duration_ms":99329,"concrete_test":"Re-run the same quench histories (same noise seeds, same parameters, same final temperature) and extract freeze-out times and vortex counts using at least two independent criteria: (i) maximum of ∂_t⟨Δ⟩, as in the paper, and (ii) first time the spatially averaged phase coherence or net winding number exceeds a threshold (or one of the criteria used by Antunes et al., 1999). Refit δ_KZM and β_1, β_2, β_3 over the same τ_Q range for each criterion. If the fitted exponents shift by more than the quoted ±0.003 (or more than the bootstrap uncertainty of the fits), the reported KZ exponents are criterion-dependent and the central claim is not yet verified; if the exponents agree within uncertainty, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s central KZ claim depends on how freeze-out is identified. In Sec. I.2, after Eq. (6), the authors explicitly cite Antunes et al. (Ref. 29) as finding KZ exponents that varied between 0.43 and 0.57 depending on the numerical criterion used to identify freeze-out time. Yet in Sec. III the freeze-out time is defined uniquely as the maximum of ∂_t⟨Δ⟩ (Eq. (24), Fig. 1 center panel), and the fit gives δ_KZM = 0.436 ± 0.003. No alternative freeze-out criterion is tested. The vortex cumulants in Fig. 2 are counted at this same freeze-out time, so the fitted exponents β_1, β_2, β_3 are not independent of the chosen criterion. The four-significant-figure precision therefore measures only the quality of the fit to this definition, not the robustness of the exponent. The concern is compounded by the exclusion of ‘unfilled data points’ that exhibit saturation, without reporting which quench rates were excluded or how the fit range was chosen. A different, equally reasonable criterion — e.g., onset of phase coherence, first appearance of net winding, or one of the criteria used in Ref. 29 — could shift the extracted exponents by an amount comparable to the 0.14 spread cited from the same literature, which would change the interpretation of the mean/variance/skewness scaling and the claimed consistency with KZM. This is an internal sensitivity problem, not an objection to the TDGL model itself.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports numerical simulations of temperature quenches in thin films of superfluid 3He-A described by a time-dependent Ginzburg-Landau (TDGL) theory with strong-coupling coefficients and Langevin noise. The authors identify the freeze-out time from the maximum slope of the spatially averaged order parameter and extract Kibble-Zurek exponents for the freeze-out time and for the first three cumulants of the vortex number distribution. They also analyze late-time vortex core structures and report a triangular/crescent vortex asymmetry. The central claimed results are (i) KZ scaling with delta_KZM = 0.436 +/- 0.003, (ii) cumulant exponents beta_1, beta_2, beta_3 that deviate from the mean-field value 1/2 and from each other by small but measurable amounts, and (iii) a pronounced excess of triangular vortices over crescent vortices.","tokens_in":18557,"tokens_out":6177,"duration_ms":60958,"significance":"If confirmed, these results would provide the first strong-coupling TDGL calculation of KZ dynamics in a chiral superfluid film, including the role of domain walls, and would extend full counting statistics tests to a two-dimensional multicomponent order parameter. The large number of noise realizations (20,000 for most quench rates) and the use of microscopically computed strong-coupling coefficients are clear strengths. However, the headline exponents depend on an untested freeze-out criterion and on fitting choices that are not documented in sufficient detail; the central quantitative claims therefore require additional robustness analysis before the paper can serve as a definitive test of KZ scaling in this system.","major_comments":[{"comment":"The freeze-out time is defined exclusively as the time at which the spatially averaged order parameter has maximum slope, yet the paper itself cites Ref. 29 as finding KZ exponents that varied between 0.43 and 0.57 depending on the numerical criterion for identifying freeze-out. No alternative criterion is tested, and the vortex counts in Fig. 2 are all evaluated at this same freeze-out time, so the exponents beta_1, beta_2, beta_3 are not independent of this definition. The quoted uncertainty delta_KZM = 0.436 +/- 0.003 reflects only the fit to this one definition. Please report the dependence of delta_KZM and the beta_i on at least two alternative freeze-out criteria (e.g., first appearance of net phase winding, or onset of phase coherence across the domain, as in Ref. 29), and document which data points were excluded as 'saturated' and the exact fit range.","section":"III, Eq. (24), Fig. 1"},{"comment":"The claim of 'small, albeit measurable' differences among beta_1, beta_2, and beta_3 is not supported by any displayed uncertainty. The text mentions that delta beta_3/beta_3 ~ 10% for N = 20,000, but the figure has no error bars or confidence intervals, and no statistical test is reported. Without these, the reader cannot judge whether the deviations are significant or within the fit noise. Add error bars (e.g., bootstrap or covariance) and a test of the null hypothesis beta_1 = beta_2 = beta_3, and also specify the fit range for each cumulant, since the scaling range in Fig. 2 appears shorter than that in Fig. 1.","section":"III, Fig. 2"},{"comment":"The triangular/crescent asymmetry is based on only 100 quench realizations and is presented as spline fits without error bars or a significance test. Since one of the suggested explanations is explicitly labeled a conjecture needing further study, the claim of a 'significant asymmetry' must be quantified. Provide standard errors or confidence bands for the vortex populations and, if possible, a test of sensitivity to the vortex-detection thresholds used in the DoG and phase-winding procedure.","section":"III.3, Fig. 6"},{"comment":"The dynamical equations use a constant damping gamma appropriate to the gapless regime near Tc, and strong-coupling beta_p coefficients whose temperature scaling was fitted near the transition. However, the freeze-out and counting occur at or after the quench has reached T = 0.5 Tc (App. A: freeze-out near t = 60 while the final temperature is reached at t = 50). At T = 0.5 Tc the quasiparticle damping and strong-coupling corrections may differ substantially from the near-Tc values used in Eqs. (20)-(21). This is a load-bearing approximation for the final vortex populations and asymmetry. Please either justify the extrapolation to T = 0.5 Tc or test one case with modified damping and strong-coupling coefficients to show that the exponents and asymmetry are insensitive to this approximation.","section":"II.3 and App. A"}],"minor_comments":[{"comment":"The combination K23 is used in Eq. (22) but defined only in the following sentence; define K23 = K2 + K3 at first use.","section":"II.1, Eq. (22)"},{"comment":"The index structure of the mass current expression is unclear; please clarify the summation convention and verify the expression is gauge invariant.","section":"II.2, Eq. (26)"},{"comment":"The video animation is an important part of the presentation, but a clickable link in a manuscript may not be persistent; please state that the animation is available as supplementary material.","section":"App. A"},{"comment":"There are several typographical errors, including 'paramter' and 'transiton'; the manuscript needs a careful proofreading pass.","section":"Abstract and Introduction"},{"comment":"The phrase 'the quench rate (time) is in units of 1/tau_GL (tau_GL)' is confusing; specify clearly the units of the horizontal and vertical axes.","section":"Fig. 2 caption"},{"comment":"The first two columns would be clearer if the headers explicitly said 'A+ amplitude' and 'A- amplitude', with row entries indicating whether that component is dominant or subdominant, to avoid ambiguity with the vortex labels.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of cond-mat.supr-con and would be a useful contribution if the quantitative claims are made robust. The central issue is not that the KZ mechanism is controversial, but that the numerical analysis as presented does not yet establish the claimed exponents and asymmetries. I would encourage the editor to request the sensitivity analyses and uncertainty quantification described in the major comments rather than rejecting the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on the 3He-A film quench paper. The genuinely new thing is the full counting statistics and the triangular/crescent vortex asymmetry in a chiral p-wave film with Z2 ground-state degeneracy. The TDGL machinery with strong-coupling coefficients is appropriate, and the number of noise realizations for the cumulants is generous (20,000). The late-time visualization of domain walls and vortices is also well done.\n\nThe main soft spot is the freeze-out time definition. The authors cite Antunes et al. showing that KZ exponents shift by ~0.14 depending on the numerical criterion, then adopt a single criterion (maximum slope of the averaged order parameter) and report δ=0.436±0.003. That precision is fit error, not robustness. They don't test alternative definitions, and they exclude 'unfilled' data points showing saturation without documenting which quench rates or how the fit range was chosen. This is a legitimate concern for the headline scaling claim. It doesn't kill the paper, but it needs to be addressed.\n\nThe asymmetry result is new and physically plausible, but the statistics are thin: 100 quenches on the large domain, no error bars, and the text admits the mechanism is a conjecture. That's fine as a first report, but it needs more support. Also, the pressure-independence claim is stated without shown data.\n\nA broader caveat: the TDGL coefficients are computed near Tc and used down to 0.5Tc. That's a standard extrapolation in this literature, but it should be acknowledged more explicitly as a model assumption, especially since the final vortex populations are measured at that low temperature.\n\nOverall, the paper is worth taking seriously. The central KZ scaling is plausible, the new asymmetry is a testable prediction, and the simulations are carefully done. The soft spots are fixable: add a freeze-out criterion sensitivity test, provide error bars on the asymmetry, document the excluded data. I'd send it to peer review and ask for those revisions. For a reading group, yes if you're interested in KZ or chiral superfluids.","headline":"First KZ quench simulations for a chiral p-wave film with domain walls; solid work, but the reported exponent rests on a single freeze-out criterion and the asymmetry needs error bars.","tokens_in":19101,"tokens_out":2440,"would_cite":true,"duration_ms":23889,"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 a two-dimensional chiral superfluid 3He film, quench-generated vortices obey Kibble-Zurek scaling, with freeze-out exponent 0.436 ± 0.003.","keywords":["Kibble-Zurek mechanism","chiral superfluid 3He-A","time-dependent Ginzburg-Landau theory","topological defects","vortex counting statistics","domain walls","temperature quench","vortex core structures"],"falsifier":"A concrete calculation that would settle the model-dependence of the central claim is to repeat the same quench protocol with the chiral $Z_2$ degeneracy removed, so the order parameter is a single complex U(1) field in two dimensions. If the three cumulant exponents then coincide, the measured differences are caused by the extra chiral and domain-wall sector; if they remain different, the deviations originate elsewhere. An experimental counterpart would be to measure the vortex-count distribution of a quenched $^3$He-A film as a function of cooling rate and compare the fitted exponents with $0.436\\pm0.003$.","tokens_in":18062,"feed_emoji":"🌀","tokens_out":9942,"duration_ms":96351,"temperature":0.7,"pith_summary":"This paper asks whether the Kibble-Zurek mechanism, the prediction that defect density after a quench is set by the crossing rate of a phase transition, survives when the ordered phase is a chiral p-wave superfluid with an additional discrete degeneracy between two time-reversed states. Using a time-dependent Ginzburg-Landau theory for strong-coupling 3He, with Langevin noise and damping, the authors simulate temperature quenches in a thin film and find that the freeze-out time and the mean number of phase vortices follow the Kibble-Zurek power laws, with freeze-out exponent $\\delta=0.436\\pm0.003$. The variance and skewness of the vortex-number distribution also scale with quench rate, but their exponents differ measurably from one another, which is only qualified support for independent-defect models of full counting statistics. The simulations also show that after freeze-out the two inequivalent vortex core structures, triangular and crescent, are produced in unequal numbers, an asymmetry attributed to vortex interactions with domain walls. If the results hold, quenched $^3$He films become a test bed for Kibble-Zurek statistics in a system whose defects carry topologically nontrivial internal structure.","feed_headline":"Kibble-Zurek scaling survives in chiral 3He films","feed_subtitle":"Simulations find freeze-out exponent 0.436; vortex-count moments deviate slightly from independent-defect models.","key_machinery":"The argument is carried by the time-dependent Ginzburg-Landau (TDGL) field theory for the spin-triplet p-wave order parameter of 3He, reduced to a two-component complex vector in the film plane, with strong-coupling coefficients, Langevin noise, and dissipative damping (Eqs. (20)--(22)). Quenches from $1.5T_c$ to $0.5T_c$ are integrated on a periodic grid; the freeze-out time is defined operationally as the time of maximum growth rate of the spatially averaged order parameter, and vortex numbers are counted from $2\\pi$ phase windings around plaquettes. The defect taxonomy uses the pair of winding numbers $(p,m)$ of the dominant and time-reversed chiral amplitudes, constrained by $p+l=m-l$ with $l$ the Chern number of the chiral ground state, which yields the triangular and crescent core structures. This machinery lets the authors separate Kibble-Zurek scaling of vortex number from the effects of the extra $Z_2$ domain-wall sector and test whether the statistics of defect formation remain those of independent events.","core_discovery":"On the paper's own terms, the central result is that the Kibble-Zurek mechanism remains intact when the symmetry-broken ground state is a doubly degenerate chiral condensate rather than a simple U(1) scalar field. Over about two decades of quench rate, the freeze-out time obeys $\\hat{t}\\sim \\tau_Q^{\\delta}$ with $\\delta=0.436\\pm0.003$, below the mean-field value $\\delta=1/2$. The mean, variance, and skewness of the number of $2\\pi$ phase vortices each scale as powers of $1/\\tau_Q$; the fitted exponents are close to one another but not equal, which the paper reports as small but measurable deviations from the prediction that all cumulants share a single Kibble-Zurek exponent. In the late-stage dynamics, two types of singly quantized vortices form in a given chiral domain, triangular cores with winding numbers $(p,m)=(+1,+3)$ and crescent cores with $(-1,+1)$, and their post-freeze-out populations are asymmetric. The paper argues this asymmetry emerges from the interactions of vortices with domain walls between time-reversed chiral domains and is absent in scalar U(1) theories.","pith_inferences":["The paper does not run a control simulation with the chiral $Z_2$ degeneracy switched off; such a U(1)-only quench in two dimensions would show whether the small differences among $\\beta_1,\\beta_2,\\beta_3$ are caused by the domain-wall sector or by something else.","The triangular-versus-crescent population asymmetry could serve as a dynamical probe of vortex-domain-wall interactions in future experiments, since equilibrium probes such as chiral edge currents do not encode quench history.","If the freeze-out exponent near 0.436 is robust across settings, it hints that the Kibble-Zurek exponent is insensitive to the internal degeneracy structure of the order parameter, while the distributional statistics are not; this extrapolation goes beyond what the paper itself claims."],"forward_implications":["If correct, thin films of superfluid 3He-A provide a quantitative test bed for the Kibble-Zurek mechanism in a chiral condensate, where the predicted freeze-out exponent is 0.436 and higher cumulants scale with quench rate.","The measured power-law scaling of variance and skewness means that quench experiments should see the full distribution of vortex counts narrow and skew in a predictable way as the cooling rate changes.","Because vortices and anti-vortices in a given chiral domain are structurally inequivalent, any Kibble-Zurek prediction for a chiral superfluid must specify which core type is being counted; the triangular-versus-crescent asymmetry is a concrete prediction.","The differences among the cumulant exponents imply that independent-defect binomial models are only approximations for this system, so vortex-domain-wall interactions belong in the theory of full counting statistics for chiral condensates."],"supporting_citations":[{"why":"Introduces the Kibble-Zurek scaling argument for defect formation in a symmetry-breaking phase transition.","marker":"[13]"},{"why":"Supplies the relation between quench rate, freeze-out time, critical exponents, and defect density used in Eq. (6).","marker":"[14]"},{"why":"The earlier 3D U(1) quench simulation whose vortex-string exponent near 0.43--0.57 the present results are compared with.","marker":"[29]"},{"why":"Presents the independent-defect model of topological-defect counting statistics tested against the simulated distributions.","marker":"[30]"},{"why":"Predicts that all cumulants of the defect number share a single Kibble-Zurek exponent, the claim the measured $\\beta_i$ deviations check.","marker":"[31]"},{"why":"Provides the time-dependent Ginzburg-Landau field-theory formulation with Langevin noise and damping used for the dynamics.","marker":"[7]"},{"why":"Supplies the temperature-dependent strong-coupling coefficients in the Ginzburg-Landau functional used in the simulations.","marker":"[8]"},{"why":"Establishes the chiral A-phase ground state in thin 3He films, the physical setting for the quench simulations.","marker":"[33]"},{"why":"Earlier classification of vortex core structures in unconventional superconductors that motivates the triangular and crescent distinctions.","marker":"[36]"},{"why":"Provides the constraint $p+l=m-l$ on vortex winding numbers used to identify the vortex types in chiral domains.","marker":"[57]"}],"fun_headline_variants":["Kibble-Zurek scaling holds in chiral 3He films","Vortex asymmetry in quenched chiral 3He films","Freeze-out exponent 0.436 in chiral superfluid films","Vortex count moments break single-exponent rule in 3He films","Chiral films: Kibble-Zurek freeze-out, asymmetric vortices"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the time-dependent Ginzburg-Landau equations, with coefficients computed near the transition temperature, continue to describe the order parameter faithfully all the way down to half that temperature, where vortex numbers and core structures are counted.","fun_headline_variants_meta":{"raw":{"variants":["Kibble-Zurek scaling holds in chiral 3He films","Vortex asymmetry in quenched chiral 3He films","Freeze-out exponent 0.436 in chiral superfluid films","Vortex count moments break single-exponent rule in 3He films","Chiral films: Kibble-Zurek freeze-out, asymmetric vortices"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000991,"raw_usage":{"total_tokens":4243,"prompt_tokens":1034,"completion_tokens":3209,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":650,"completion_tokens_details":{"reasoning_tokens":3116}},"tokens_in":650,"tokens_out":3209,"duration_ms":24910,"temperature":1.0,"reasoning_tokens":3116,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:17:00.075206+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete calculation that would settle the model-dependence of the central claim is to repeat the same quench protocol with the chiral $Z_2$ degeneracy removed, so the order parameter is a single complex U(1) field in two dimensions. If the three cumulant exponents then coincide, the measured differences are caused by the extra chiral and domain-wall sector; if they remain different, the deviations originate elsewhere. An experimental counterpart would be to measure the vortex-count distribution of a quenched $^3$He-A film as a function of cooling rate and compare the fitted exponents with $0.436\\pm0.003$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the relation between quench rate, freeze-out time, critical exponents, and defect density used in Eq. (6)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier classification of vortex core structures in unconventional superconductors that motivates the triangular and crescent distinctions."}],"review_version":1}