Excitation of Collective Modes in a Chiral Superfluid by Thermal Quench

J. A. Sauls; Noble Gluscevich

Excitation of Collective Modes in a Chiral Superfluid by Thermal Quench

T0 review · 2 major / 4 minor · reviewed 2026-07-12 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Rapid cooling into chiral superfluid ³He-A excites bosonic collective modes as well as topological defects, with Kibble-Zurek exponents that cross over smoothly from underdamped to overdamped dynamics.

desk verdict Clean computational extension of KZ work that maps the damping-driven z crossover and the PSD of Higgs/clapping modes; the freeze-out threshold is the only real soft spot and it is already flagged. read the letter →

arxiv 2606.28306 v3 pith:BUWA3C6V submitted 2026-06-26 cond-mat.supr-con cond-mat.stat-mech

Noble Gluscevich , J. A. Sauls This is my paper · ORCID

classification cond-mat.supr-concond-mat.stat-mech

keywords chiralsuperfluidKibble-Zurekmechanismtime-dependentGinzburg-LandaucollectivemodesHiggsmodeclappingthermalquenchLangevindamping

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper shows that a fast temperature quench through the second-order transition into the chiral A-phase of superfluid ³He films does more than produce the usual Kibble-Zurek network of vortices and domain walls. Using stochastic time-dependent Ginzburg-Landau simulations, the authors demonstrate that the same quench also launches the bosonic collective modes of the newly formed chiral domains—the Higgs amplitude mode and the clapping modes. In thermal equilibrium the power spectral density of each mode exhibits a sharp mass threshold that is washed out by strong Langevin damping. After a quench the order-parameter amplitude grows, oscillates, and coarsens; the freeze-out time and correlation length obey Kibble-Zurek scaling whose dynamical exponent z interpolates continuously from 1 (weak damping) to 2 (overdamped), while the correlation-length exponent remains ν ≈ ½. The result supplies a concrete laboratory route to study nonequilibrium mode excitation and damping-dependent critical dynamics in a multi-component topological superfluid.

What carries the argument

Stochastic time-dependent Ginzburg-Landau equations for the two-component complex order parameter of a ³He-A film, driven by Gaussian white noise and a linear temperature quench, whose Fourier amplitudes and power spectral densities track both the defect network and the Higgs/clapping modes.

What would settle it

Measure the power spectral density of order-parameter fluctuations immediately after a controlled temperature quench in a ³He film and check whether sharp mass thresholds at 2Δ and √2 Δ appear only in the underdamped regime and whether the freeze-out scaling of correlation length and time interpolates between the predicted underdamped and overdamped exponents as damping is varied.

Watch

Extended reading notes

Core claim

Rapid cooling through the chiral superfluid transition simultaneously generates topological defects by the Kibble-Zurek mechanism and excites the bosonic collective modes of the nascent chiral domains; the dynamical exponent that governs freeze-out scales smoothly from z = 1 to z = 2 with increasing Langevin damping, while ν remains ≈ ½.

Load-bearing premise

The freeze-out time is identified by the moment the spatially averaged order-parameter amplitude first reaches a fixed numerical threshold of 0.01 k_B T_c; any other cut-off or diagnostic can shift the extracted scaling exponents.

Editorial extensions

If this is right

Power spectra of order-parameter fluctuations after a quench should display sharp mass thresholds for the Higgs and clapping modes when damping is weak, and smooth cross-overs when damping is strong.
Kibble-Zurek freeze-out time and correlation length will obey δ_KZM ≈ 1/3, χ_KZM ≈ 1/3 in the underdamped limit and δ_KZM ≈ 1/2, χ_KZM ≈ 1/4 in the overdamped limit, with a continuous crossover in between.
Nonlinear mode coupling after freeze-out generates harmonic sidebands and mass shifts that can be read from the same Fourier spectra.
Large-scale post-quench simulations produce a measurable complex network of domain walls and vortices whose coarsening can be tracked in real and Fourier space.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The same damping-tuned crossover of z should appear in other multi-component superfluids or Bose condensates whose collective-mode spectrum is known, offering a general experimental knob for critical dynamics.
Observation of the predicted mass thresholds in a real ³He film would give a direct spectroscopic signature of the chiral ground state formed by the quench, complementary to NMR or heat-capacity probes.
Because the clapping modes lie well below the pair-breaking continuum, they remain long-lived even when the Higgs mode is broadened, potentially allowing selective excitation and detection of time-reversed chirality fluctuations.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Request a human review

A listed scientist reviews the paper for a fee and the review publishes here regardless of verdict. See the reviewers or get listed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The new piece here is the simultaneous tracking of bosonic collective modes (Higgs and clapping) with topological-defect formation under a temperature quench, plus a clean numerical map of how the dynamical exponent z crosses over from 1 to 2 as Langevin damping is raised while u stays at ½. That is a genuine addition to the authors’ own 2025 defect-statistics paper and to the broader KZ literature on multi-component order parameters.

What they do well is straightforward. The linearized Langevin analysis of the power spectral densities (Eqs. 19–23) is standard and matches the numerical spectra for both weak and strong damping; the sharp mass thresholds and 1/ω tails in the underdamped regime are clear. Large-scale pseudospectral runs (L = 440 ξ_GL) give freeze-out and correlation-length exponents that sit close to the expected mean-field limits, and the continuous exponential fit for δ_KZM(γ) is tidy. The real-space correlation-function analysis is independent of the k-space diagnostics and supports the same picture. Citations are appropriate and the math is inspectable from the text.

The soft spot is exactly the one already noted: freeze-out is defined by the ad-hoc threshold ⟨Δ⟩_t = 0.01 k_B T_c. Different cuts or diagnostics shift the fitted exponents by a few percent, and the reported limiting values already sit slightly below pure mean-field. That is real but not fatal; the paper itself acknowledges it, and the independent PSD and correlation-length results remain consistent. Absence of public code and exhaustive finite-size/noise error bars is a practical limitation, not a conceptual one.

This is for people who work on nonequilibrium phase transitions, quantum turbulence, or collective modes in chiral superfluids. It deserves a serious referee. I would accept it for peer review and would cite the z-crossover and the PSD results if I were writing on KZ scaling or mode excitation in the next year.

Referee Report

2 major / 4 minor

Summary. The manuscript uses stochastic time-dependent Ginzburg-Landau theory for thin-film chiral superfluid 3He-A to study both thermal excitation of bosonic collective modes (Higgs, Nambu-Goldstone, and clapping modes) and their generation by linear temperature quenches through Tc. Analytic power spectral densities for the linearized Langevin dynamics are derived and shown to match large-scale pseudospectral simulations for weak and strong damping. Quench simulations extract Kibble-Zurek freeze-out times and correlation lengths over a range of quench rates and dimensionless damping γ; the dynamical exponent z is found to cross over smoothly from ≈1 (underdamped) to ≈2 (overdamped) while the correlation-length exponent remains ν≈1/2, producing the expected mean-field limits for the KZ exponents δ_KZM and χ_KZM.

Significance. The work cleanly unifies two aspects of nonequilibrium chiral superfluids that are usually treated separately: the spectrum of order-parameter collective modes and Kibble-Zurek defect formation. The analytic PSD formulae (Eqs. 20–23) and their numerical confirmation provide a concrete, falsifiable signature of the mode masses that could be sought in future NMR or acoustic experiments on 3He films. The continuous damping-driven crossover of z, obtained from both freeze-out times and real-space correlation lengths on large (L=440 ξ_GL) grids, is a useful quantitative result for the broader KZ literature. Strengths include the transparent comparison of analytic and numerical spectra, the systematic scan of γ, and the independent real-space check of the correlation-length scaling.

major comments (2)

Sec. IV and Fig. 2: All reported KZ exponents rest on the operational definition that freeze-out occurs when the spatially averaged amplitude first reaches ⟨Δ⟩_t=0.01 k_B T_c. The paper itself notes that alternative diagnostics shift the fitted δ_KZM and χ_KZM by a few percent, and the limiting values already lie slightly below mean-field (0.313 vs 1/3, 0.488 vs 1/2). A short robustness check against at least one other threshold (or against the peak of the structure factor) is needed to confirm that the smooth z(γ) crossover is not an artifact of this particular cut-off.
Sec. V and Fig. 6: The nonlinear homogeneous spectra are shown only for zero damping and two hand-chosen initial amplitudes. Because the central claim is that nonlinearities modify the PSD immediately after freeze-out, a quantitative comparison of the full stochastic PSD (with defects present) just after freeze-out versus the thermal PSD of Sec. III would strengthen the argument that the observed oscillations in Fig. 2 are indeed the Higgs and clapping modes rather than defect-related ringing.

minor comments (4)

Eq. (11) and surrounding text: the dimensionless damping γ is defined with τ_GL and μ, but the numerical values of au_GL and eta_245 used in the simulations are never stated; a single sentence listing the microscopic units would aid reproducibility.
Fig. 1 caption and Sec. III: the analytic curves (dashed) are said to be “Eq. (22)”, yet the slight excess power below the mass thresholds is attributed to nonlinearities without a quantitative residual. A brief statement of the relative L2 deviation would clarify how small the nonlinear corrections actually are.
Appendix A: the strong-coupling eta_sc_i are mentioned but never used; either drop the paragraph or note explicitly that all simulations are performed in the weak-coupling limit.
Typographical: “varia-tons” (App. A), “coicides” (Sec. IV), and inconsistent spacing around “δ_KZM” appear in several places.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: measured KZ exponents and mode PSDs are independent simulation outputs compared to standard mean-field TDGL limits; only minor self-citation of prior defect-counting work.

full rationale

The derivation chain is self-contained. Mode masses (Table I, Eqs. 7–9) follow directly from the quadratic expansion of the standard TDGL Lagrangian (Eq. 3) with the conventional choice μ=β̄1 that places the Higgs mass at 2Δ; this is parameter fixing, not a prediction that reduces to its input. Analytic PSDs (Eqs. 20–23) are obtained from the linearized Langevin oscillator and are then compared to full nonlinear simulations; residual discrepancies are attributed to weak nonlinearities, not forced by construction. Freeze-out time ˆt and correlation length ˆξ are extracted operationally from the stochastic TDGL trajectories (⟨Δ⟩t=0.01 kBTc threshold and HWHM of C(r,t)), after which the scaling exponents δKZM(γ) and χKZM(γ) are fitted; z and ν are then read off via the standard KZ relations. The observed smooth crossover z:1→2 with ν≈1/2 is therefore a numerical result, not an algebraic identity. Self-citations (chiefly Ref. 23 for earlier defect statistics and the authors’ prior TDGL parameter papers) supply background material coefficients and context but are not load-bearing for the new mode-excitation or damping-dependent exponent claims. The only soft spot is the conventional but ad-hoc freeze-out threshold, which can shift fitted numbers by a few percent; that is a methodological choice, not circularity. Score 1 reflects the presence of non-essential self-citation without elevating it to a circular step.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

The central claims rest on the standard TDGL Lagrangian for 2D chiral p-wave pairing, the fluctuation-dissipation relation that fixes the noise correlator to the damping coefficient, and the operational definition of freeze-out. No new particles or forces are introduced; the free parameters are numerical cut-offs and the phenomenological crossover fit.

free parameters (3)

freeze-out amplitude threshold = 0.01 k_B T_c
⟨Δ⟩_t = 0.01 k_B T_c is chosen by hand to mark freeze-out; all scaling exponents are extracted from this cut-off.
crossover fit parameters b, δ_UD, δ_OD = b=7.427, δ_UD=0.313, δ_OD=0.488
Exponential interpolation formula (Eq. 26) fitted to the measured δ_KZM(γ); used only to summarize the crossover, not to define the exponents themselves.
dimensionless damping values γ
Discrete set {1.0, 0.5, 0.1, 0.05, 0.01, 0.005} scanned by hand to map the under- to over-damped regimes.

assumptions (3)

domain assumption Time-dependent Ginzburg-Landau Lagrangian with weak-coupling material coefficients (μ, K_i, β_i) correctly captures the long-wavelength dynamics of the 2D chiral order parameter.
Invoked throughout Sec. II and Appendix A; standard for ³He films but remains an effective theory.
domain assumption Gaussian space-time white noise whose correlator is fixed by the fluctuation-dissipation theorem (Eq. 12) adequately models thermal fluctuations near T_c.
Used to generate all stochastic trajectories; temperature dependence of Γ is left unexplored.
standard math Mean-field values ν=1/2 and the limiting dynamical exponents z=1 (underdamped) / z=2 (overdamped) are the correct asymptotic targets.
Taken from classic Kibble-Zurek scaling theory and used as benchmarks for the numerical exponents.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Excitation of Collective Modes in a Chiral Superfluid by Thermal Quench." pith.science (2026). https://pith.science/paper/BUWA3C6V

@misc{pith2026260628306,
  author       = {Pith},
  title        = {Pith review of: Excitation of Collective Modes in a Chiral Superfluid by Thermal Quench},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BUWA3C6V}},
  note         = {Machine review of arXiv:2606.28306}
}

abstract

Based on time-dependent Ginzburg-Landau field theory we show that rapid cooling through the second-order phase transition into superfluid $^3$He-A excites collective modes of newly formed chiral domains, in addition to topological defects that are formed via the Kibble-Zurek mechanism. Simulations of temperature quenches in the presence of Gaussian space-time white noise generate a highly excited inhomogeneous condensate. Large-scale simulations exhibit a complex network of domain walls and vortices. We report results for the excitation of bosonic collective modes by thermal noise as well as nonequilibrium temperature quenches, followed by coarsening dynamics tracked in terms of the Fourier components of the order parameter amplitudes. For thermal states, the spectrum of bosonic excitations is defined by a power spectral density (PSD) for each mode, which is sensitive to the Langevin damping. For weak damping the PSD onsets sharply at the frequency corresponding to the mass of the bosonic mode, then decays as $1/\omega$. We also track the dynamics of the order parameter following a temperature quench. We report results for the scaling exponents of Kibble-Zurek freeze-out time and correlation length as a function of quench rate for several damping rates. The dynamical exponent $z$ is shown to transition smoothly from $z=1$ to $z=2$ as the damping is increased, while the correlation length exponent, $\nu\approx 1/2$, is independent of damping.