Pith. Sign in

REVIEW

Local Quantum Cooling for Large Fermi Systems with Pairing

Nuclear Theory · Submitted 14 Dec 2025 · 60 references

Pith's verdict: A solid, checkable extension of local quantum friction to the pairing channel with a useful density-scan protocol; the caveat is that it cools collective modes only, and the spin-imbalanced conclusion leans on that limitation.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 3, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2512.12866 (version 2)
Pith IDQVRHWOWO, Pith's permanent record of this paper and its review
CitationDaniel P\k{e}cak, Gabriel Wlaz{\l}owski, J. E. Alba-Arroyo, Michael McNeil Forbes. "Local Quantum Cooling for Large Fermi Systems with Pairing." arXiv:2512.12866 (2025). https://pith.science/paper/QVRHWOWO
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

A Hermitian cooling potential that acts on both density currents and pairing-field fluctuations can drive large Fermi superfluids toward their ground state while keeping time evolution unitary and conserving particle number.

The reading

This paper extends the local-quantum-friction (LQF) idea to superfluid pairing: by adding a carefully constructed local, Hermitian 'cooling potential' in the pairing channel, together with the existing density-current damper, time-dependent Hartree-Fock-Bogoliubov (TDHFB) evolution remains unitary while steadily removing energy from collective excitations. The central result is an explicit formula for the cooled pairing field, Δ_TD(r)=Δ(r)(1−β sin(θ_B′−θ_κ)), which preserves particle number and damps fluctuations of the pairing field. In tests, pairing-channel cooling often relaxes a system to the true ground state where current-only cooling stalls, and combining both channels is the most efficient. The same construction can be adapted to change particle number in time, enabling single-run 'density scans' of ground-state properties. The authors demonstrate this on spin-imbalanced unitary Fermi gases and on neutron-star crust matter, including nuclear pasta phases, and find large computational savings compared to diagonalization-based static solvers.

What carries the argument

The central object is the Hermitian cooling potential U ∝ −i[H,R], i.e., the negative of the natural unitary flow of the generalized density matrix R. Restricting U to be diagonal in position space gives the density-current damper U_σ^(d) = −(α/mρ_0)∇·j_σ, while a non-diagonal (but still local) block treatment yields the pairing cooling potential Δ_TD(r)=Δ(r)(1−β sin(θ_B′−θ_κ)). This phase-sensitive multiplier suppresses pairing-field fluctuations when the phase of the kinetic pairing flow B′ deviates from the phase θ_κ of the anomalous density, removing energy from the collective channel while preserving particle number and unitarity.

Extended reading notes

Core claim

The paper's central claim is that dissipation can be introduced into large fermionic simulations without breaking unitarity or the convenient local structure of TDHFB. The cooling potential is chosen proportional to the commutator −i[H,R]—the flow the system would naturally follow—then restricted to be diagonal (local) in position space. In the pairing channel this yields the explicit replacement Δ→Δ(1−β sin(θ_B′−θ_κ)), which damps collective oscillations of the superfluid order parameter while conserving particle number because the phase θ_κ of the anomalous density is factored out. The authors show that this pairing-channel cooling is often more effective than damping currents alone, that

Load-bearing premise

The cooling potential removes energy only from collective degrees of freedom—density currents and pairing-field fluctuations—so if the dominant excitations in a target system are intrinsic single-quasiparticle excitations, the evolution will not relax to the true ground state and may stall near a local minimum (as the authors note in Sections III.A and IV.B).

Editorial extensions

If this is right

  • In strongly paired systems, cooling through the pairing channel relaxes states to the true ground state faster than current-only LQF, sometimes reaching the ground state where current-only cooling stalls.
  • The same framework, with a time-dependent imaginary component, can adjust particle number during evolution, enabling a single time-dependent run to map out ground-state properties over a range of densities (a 'density scan').
  • Conservative evolution and locality make the method cheap: for periodic pasta lattices (2×2 and 3×3 grids) the dynamic cooling approach converges in about an order of magnitude less GPU time than diagonalization-based static solvers.
  • Because the method preserves particle number and unitarity, it can be used to prepare high-quality initial states for subsequent time-dependent simulations of e.g. vortex dynamics and turbulence in neutron-star crust matter.
  • The construction is generic enough to be extended with stochastic noise, which the authors identify as a path toward fluctuation–dissipation dynamics and thermalization studies.

Reading between the lines

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

  • If the pairing-channel cooling is as effective as demonstrated, one could use it as a diagnostic: the energy plateau reached after long LQF time signals the energy carried by intrinsic single-quasiparticle excitations, which remain undamped by this method.
  • A testable extension is to apply the cooling potential to a system initialized with a single known quasiparticle excitation (e.g., a local phase twist or a single occupied state above the gap) and check whether the energy stalls, isolating the collective vs. intrinsic contribution.
  • The phase-sensitive form Δ_TD∝(1−β sin(θ_B′−θ_κ)) suggests an interpretation as a superfluid 'velocity-dependent' friction that damps the Goldstone phase mode; this might connect to effective models of superfluid dissipation and to the theory of critical damping in ordered media.
  • The density-scan variant might be adapted to study quantum phase transitions by slowly driving the system across a critical region; the cooling would keep the system near the instantaneous ground state, though traversal of energy barriers (as in pasta transitions) would require an additional mechanism, since the authors note the system can remain trapped near a local minimum.
Original abstract

We present a framework for local quantum cooling that can be efficiently applied to large-scale Fermi systems. The method introduces local Hermitian operators as a cooling potential while strictly preserving the unitarity of time evolution. Our formulation scales favorably with system size and can be seamlessly integrated into time-dependent density-functional theory frameworks. We demonstrate that energy cooling arises from the damping of particle currents and pairing-field fluctuations. Furthermore, we develop a variant of the scheme that allows the particle number to vary in time, enabling controlled density scans. The method is generic and versatile, as illustrated by applications to spin-imbalanced unitary Fermi gases and to nuclear matter in the neutron-star crust. The framework can be naturally extended to include stochastic noise, providing a foundation for studying thermalization in strongly interacting Fermi superfluids.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,448 input tokens.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.08 · human participants 0.01 · controlled comparison 0.02 · randomized assignment 0.01 · states quantitative result 0.04 · compares with prior results 0.06 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

We present a framework for local quantum cooling that can be efficiently applied to large-scale Fermi systems. The method introduces local Hermitian operators as a cooling potential while strictly preserving the unitarity of time evolution. Our formulation scales favorably with system size and can be seamlessly integrated into time-dependent density-functional theory frameworks. We demonstrate that energy cooling arises from the damping of particle currents and pairing-field fluctuations. Furthermore, we develop a variant of the scheme that allows the particle number to vary in time, enabling controlled density scans. The method is generic and versatile, as illustrated by applications to spin-imbalanced unitary Fermi gases and to nuclear matter in the neutron-star crust. The framework can be naturally extended to include stochastic noise, providing a foundation for studying thermalization in strongly interacting Fermi superfluids.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We present a framework for local quantum cooling that can be efficiently applied to large-scale Fermi systems. contribution 1.00 0.92 0.23 asserted 0.98
s1 The method introduces local Hermitian operators as a cooling potential while strictly preserving the unitarity of time evolution. method 0.95 0.89 0.11 asserted 1.00
s2 Our formulation scales favorably with system size and can be seamlessly integrated into time-dependent density-functional theory frameworks. contribution 0.69 0.69 0.10 asserted 0.96
s3 We demonstrate that energy cooling arises from the damping of particle currents and pairing-field fluctuations. result 0.65 0.93 0.29 asserted 1.00
s4 Furthermore, we develop a variant of the scheme that allows the particle number to vary in time, enabling controlled density scans. contribution 0.73 0.95 0.07 asserted 1.00
s5 The method is generic and versatile, as illustrated by applications to spin-imbalanced unitary Fermi gases and to nuclear matter in the neutron-star crust. contribution 0.78 0.79 0.16 asserted 1.00
s6 The framework can be naturally extended to include stochastic noise, providing a foundation for studying thermalization in strongly interacting Fermi superfluids. future work 0.87 0.56 0.09 asserted 0.75

What kind of paper the classifier reports this is: method or tool 1.00 · evidence from simulation or computation 0.80.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1. Particle-conserving time evolution (
Figure 2. FIG. 2. Evolution of the total energy of an initially uniform 2D state in a time-dependent perturbation
Figure 3. FIG. 3. Particle-changing time evolution with
Figure 4. FIG. 4. (a) Decrease in the total energy of a spin-imbalanced unitary
Figure 5. shows an example of such a simulation. We begin with a solution obtained in a cubic box of size 323 with a lattice spacing of d𝑥 = 1.25 fm. We start with 𝑍 = 40 protons and 𝑁0 = 116 neutrons. Since the number of neutrons already exceeds the neutron drip line, the initial configuration corresponds to a Zirconium nucleus immersed in a dilute neutron gas. We then evolve this system in time with lqf (𝛼 = 10, 𝛽 = 5) and include a particle-changing term 𝛾 = 10 to gradually increase the number of neutrons with 𝑁req.(𝑡) = 𝑁0 + 𝑟𝑡 (64) in Eq. (49). The number of protons remains fixed. With a…
Figure 6. FIG. 6. Structural properties of protonic clusters in the neutron-star inner crust obtained from time-dependent
Figure 7. FIG. 7. a) Evolution of the energy difference with the approximate
Figure 8. FIG. 8. Effect of different parameter values
Figure 10. FIG. 10. Comparison of different dissipation strengths for (a) the mean field (𝛽

Showing 9 of the paper's 10 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2512.12866 by the authors.

References
How the literature uses it

What this paper does with the arXiv papers it cites

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

What it cites · 4 sentences without a confident label

Builds on 1 paper
Uses the method or tool of 1 paper
  • For both test cases, and for later applications with spin-imbalanced ultracold Fermi gases, we use the ASLDA [ref] , which defines the intrinsic energy density functional as [equation] where $ρ(r)=ρ_(r)+ρ_(r)$ is the total density, $(r)=(r)^2/2m$ is the local Fermi energy, and $(r)=[3π^2ρ(r)]^1/3$ is the local Fermi momentum.

    The Unitary Fermi Gas: From Monte Carlo to Density Functionals · Tests of dissipative dynamics · 0.98

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.