{"id":"0040c0f0-c10c-4252-8da3-28eb8c86982d","arxiv_id":"1909.00992","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A lambda-integration-free Eilenberger free energy functional is derived from Luttinger-Ward theory and generalized to spin-triplet correlations and spin-dependent fields.","lead":"This paper derives a closed-form free energy for superconducting and superfluid systems described by quasiclassical theory, extending Eilenberger's 1968 result to spin-triplet pairing. The result removes a numerical integration step and provides practical expressions for the dirty limit, with applications to superfluid helium-3, Majorana nanowires, and competing superconducting phases.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Real-part subtraction at singularities of [tau_t, g_lambda] is the load-bearing unproven step: Eq. (11) equals the Luttinger-Ward functional only if discarded winding contributions are pure imaginary and Re E selects the correct branch.","rationale":"The paper credibly fills a gap by connecting the Serene-Rainer lambda-integral expression to an Eilenberger-type local functional and extending it to spin-triplet systems. The independent matrix-algebra verification of the variation property (8) in the supplement is a real check, and the reduction to the known diffusive expression (21) is a useful cross-check. I read the central claim as conditional on the treatment of singularities of [tau_t, g]. The authors are explicit about the singularities and about taking the real part, but the proof that the subtracted winding contributions are pure imaginary rests on a local analysis assuming U is nonsingular and that each singularity gives i pi times an integer. This is insufficient to guarantee that Re E[g_1[Sigma]] is independent of the choice of tau_t for all configurations covered by the claim, especially when singularities sweep through the lambda interval or the texture is topologically nontrivial. If the real part is not the correct branch, Eq. (11) and its diffusive limit (21) would not be the Luttinger-Ward free energy. The proposed numerical comparison against the lambda-integral for a case with interior singularities would settle the point. This is the same weakness the reader identified, so I agree with the CONDITIONAL verdict and would not move it without seeing that check or an analytic proof.","tokens_in":11319,"tokens_out":36367,"duration_ms":392671,"concrete_test":"Numerical check in 1D: choose a spin-triplet configuration such as a domain wall in 3He-B or a spin-active interface and a self-energy Sigma(x) for which [tau3, g_lambda] is singular for some lambda in (0,1). Solve the Eilenberger equation for g_lambda on a dense lambda grid and compute Omega from the lambda-integral (3). Then compute Re Omega from Eq. (11) using a tau_t such as tau1 or U0^{-1} tau1 U0 that is nonsingular at lambda = 1, and repeat with a second inequivalent tau_t that also avoids endpoint singularities. If Re Omega_11 disagrees with Omega_3 by more than the lambda-grid error, or if the two tau_t choices disagree in their real parts, the real-part subtraction does not reproduce the Luttinger-Ward functional and the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on replacing the lambda-integral in Eq. (3) with the endpoint functional (11) using E[g] from Eq. (12). This substitution is valid only if E has the variation (8) throughout the lambda-sweep and if the Stokes boundary terms are under control. The paper itself notes that [tau_t, g_lambda] need not be invertible on M = [0,1] x [-inf, inf] and that singular points generate winding contributions discarded by taking the real part. This is the load-bearing point: no proof is given that these contributions are always purely imaginary multiples of i pi, that U remains nonsingular at the poles of [tau,g]^{-1} as assumed around Eq. (S10), or that taking Re E does not select a branch depending on the chosen tau_t and on how singularities cross the lambda path. Because E is fixed by Eq. (8) only up to topological terms, the claim that Re Omega_11 equals the Luttinger-Ward functional at the saddle point depends exactly on this unproved subtraction. This matters most for the advertised spin-triplet/textured applications, where zeros of [tau_t, g_lambda] may not be avoidable by a single global tau_t. The derivation is plausible and the local loop-integral calculation is suggestive, but the global real-part statement is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a free-energy functional for quasiclassical superconductors and superfluids starting from the Luttinger–Ward formulation. The main result is Eq. (11), with the gradient functional E[g] given by Eq. (12), which avoids the λ-integration of the earlier Serene–Rainer expression (3). The functional is shown to reduce to Eilenberger's free energy in the spin-diagonal case and is generalized to spin-triplet correlations and spin-dependent fields by allowing an arbitrary matrix τt in Eq. (12). The paper also derives a diffusive-limit (Usadel-type) free energy in terms of momentum-averaged propagators, Eqs. (20)–(21), and provides a supplement containing the variation calculation, a Riccati parametrization, and a discussion of singularities of [τ, g].","tokens_in":11559,"tokens_out":12706,"duration_ms":113566,"significance":"If the central claim is correct, the paper fills a genuine gap by providing a closed-form, λ-integration-free free energy for nonuniform superconducting and superfluid states with arbitrary spin structure, including spin-triplet pairing. The construction of E[g] is explicit, the supplement contains a self-contained derivation of its variation, and the reduction to the known Eilenberger and Usadel results in appropriate limits is demonstrated. These are concrete, checkable contributions that would be useful for analyzing competing phases, vortices, and topological systems. The main weakness is that the equality between the real part of the proposed functional and the Luttinger–Ward functional at the saddle point relies on an unproven treatment of singularities of [τ, g], which is especially relevant for the advertised spin-triplet applications.","major_comments":[{"comment":"The central claim that Re Ω_{11} equals the Luttinger–Ward functional depends on the assertion that singularities of [τ, g_λ] produce only pure-imaginary winding contributions that are removed by taking the real part. The supporting calculation, Eq. (S10), assumes that U in the representation g = U τ3 U^{-1} is nonsingular at the zeros of [τ, g]. The manuscript does not prove this assumption, nor does it state a condition on g or τ that would guarantee it. For the advertised spin-triplet and textured applications, where the order parameter has zeros, it is not evident that a single global τ exists for which U is nonsingular at all singular points of [τ, g_λ] along the λ-sweep. Without such a proof, the equality Re Ω_{11} = Eq. (3) is not established; at best the paper shows that Eq. (12) has the correct variation away from singularities. Please provide a proof or a precise, physically motivated condition under which the real-part subtraction is exact, or weaken the claim accordingly.","section":"Supplement, Eq. (S10); main text below Eq. (12)"},{"comment":"The paper states that one should choose τ to avoid singularities, for instance τ1 near the normal state or τ = U0^{-1}τ1 U0 for a reference g0 ≈ g. For spin-triplet pairing, the anomalous component f can vanish on the Fermi surface, and for inhomogeneous textures with nontrivial topology a global smooth U0(x) may not exist. The manuscript does not show that these recommended choices make [τ, g_λ] invertible for all λ ∈ [0,1] and all x, nor does it quantify the error if they do not. Since the spin-triplet generalization is the paper's main advertised advance, this gap is load-bearing rather than cosmetic.","section":"Eq. (12) and the paragraph following Eq. (17)"},{"comment":"The manuscript correctly notes that E[g] is defined by Eq. (8) only up to topological terms and that the presence of singularities depends on τ. It then asserts that the free energy is real-valued and takes the real part. However, the real part of E[g] itself could depend on the choice of τ if the discarded winding contributions are accompanied by real branch-cut terms, or if different τ choices shift singularities across the λ-path. The paper should either prove that Re E[g] is independent of τ up to total derivatives, or state the conditions under which it is; otherwise the proposed free energy is not manifestly unique.","section":"Discussion of non-uniqueness of E[g] after Eq. (12)"}],"minor_comments":[{"comment":"The reference to 'Supplementary information at ???' is a placeholder and must be filled with the actual supplementary material identifier.","section":"Reference [20]"},{"comment":"There is a typo: 'The the general form' should read 'The general form'.","section":"Text after Eq. (16)"},{"comment":"The title contains a spacing artifact: 'supercon ducting' should be 'superconducting'.","section":"Title"},{"comment":"For spin-triplet states, the anomalous component f may fail to be invertible; the text should note explicitly that Eq. (17) is only valid when the chosen τ (here τ3) avoids such zeros, and that Eq. (12) with a different τ should be used otherwise.","section":"Eq. (17)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the variational construction is plausible, but the singularity handling is the crux of the proof that the real part of Eq. (11) equals the Luttinger–Ward functional. If the authors can supply a rigorous treatment of the singularities, or carefully restrict the claim to configurations where the chosen τ avoids them, the paper would be suitable for publication. As presented, the advertised generality for spin-triplet and topological applications is not fully established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper delivers what it promises: a λ-integration-free free energy functional for quasiclassical superconductors and superfluids with arbitrary spin structure, derived systematically from the Luttinger–Ward functional. The central result, Eqs. (11)–(12), is new and fills a gap that has been open since Eilenberger's 1968 paper. The spin-singlet limit reduces to the known expression, and the triplet generalization is explicit and usable. The diffusive-limit expression, Eq. (21), will be handy for the Usadel community. The derivation is coherent; the supplement actually shows the main variational calculation, and the Riccati parametrization provides a concrete check. This is not a hand-waving paper.\n\nThe stress-test note flags the real-part subtraction at singularities of [τ_t, g_λ]. I read the supplement carefully, and the authors do address this: they compute the singularity contributions and find they are winding numbers of the form i π m, hence purely imaginary. That is more than most papers would do. The remaining gap is narrower than the stress-test suggests: they assume U stays nonsingular at the poles of [τ,g]^{-1}, and they do not prove that taking Re E always selects the correct branch for arbitrary textures. That is a legitimate technical caveat, but it is not a load-bearing flaw. For most practical applications, including the advertised triplet and textured cases, one can choose τ to avoid the worst singularities, and the real-part prescription is physically sensible. I would want the authors to state the conditions more carefully in a final version, but I would not block publication on it.\n\nThe diffusive-limit derivation uses a leading-order spherical-harmonic expansion with the standard normalization argument; that is fine. The citation pattern is appropriate, including the authors' own prior work, which is directly relevant. No hidden parameters, no fitting, nothing circular.\n\nWho is this for? Anyone working on superfluid 3He phases, FFLO states, vortex physics, or superconductor/ferromagnet hybrids. This will become a standard citation for free-energy calculations in those areas. I would bring it to our reading group and would cite it. It deserves serious peer review and, in my view, acceptance after minor revisions that clarify the singularity conditions and the branch choice.","headline":"A genuinely useful closed-form free energy for quasiclassical superconductors with general spin structure; the singularity handling is a real caveat but not a fatal one.","tokens_in":12112,"tokens_out":1805,"would_cite":true,"duration_ms":21876,"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":"The paper derives a closed-form quasiclassical free-energy functional covering spin-singlet and spin-triplet superconducting and superfluid states, avoiding coupling-constant integration.","keywords":["quasiclassical theory","Eilenberger equation","Luttinger-Ward functional","spin-triplet pairing","superfluid helium-3","diffusive limit","Usadel equation","free energy functional"],"falsifier":"Compute a nonuniform spin-triplet or spin-singlet configuration where $[\\hat\\tau,\\hat g_\\lambda]$ has a zero for some $\\lambda\\in[0,1]$, evaluate the real part of (19), and compare it with the exact $\\lambda$-integrated Luttinger-Ward expression (3) at the same saddle point: if the real parts differ by a nonzero winding contribution, the real-part prescription fails.","tokens_in":2106,"feed_emoji":"⚛️","tokens_out":2076,"duration_ms":74112,"temperature":0.7,"pith_summary":"Superconducting and superfluid states are described by quasiclassical Green's functions, but computing their free energy has required a coupling-constant ($\\lambda$) integration over auxiliary propagators. The paper removes that integration by constructing a closed functional of the full quasiclassical propagator whose saddle point gives the Eilenberger equations and whose value at that saddle point is the free energy. The construction holds for arbitrary spin structure, including spin-triplet pairing generated by exchange fields or spin-orbit coupling, generalizing the 1968 Eilenberger expression. In the dirty limit the functional reduces to a simple expression in momentum-averaged propagators that yields the Usadel equation. This closes a gap between the fundamental Luttinger-Ward free energy and practical quasiclassical functionals.","feed_headline":"Superconducting free energy unified for singlet and triplet pairing","feed_subtitle":"A closed-form Eilenberger functional now covers arbitrary spin structure, including superfluid 3He.","key_machinery":"The load-bearing object is the gradient functional $E[\\hat g] = \\frac12\\mathrm{Tr}\\bigl(\\hat g[\\hat\\tau_t,\\hat g]\\, v_F\\cdot\\check\\nabla[\\hat\\tau_t,\\hat g]^{-1}\\bigr)$, defined on the quasiclassical propagator $\\hat g$ (a $4\\times4$ matrix in spin and Nambu space) and an auxiliary matrix $\\hat\\tau_t$ with $\\hat\\tau_t^2=1$. It turns a $\\lambda$-integration over auxiliary propagators into a surface term, so the free energy becomes a direct functional of $\\hat g$; choosing $\\hat\\tau$ appropriately also avoids most singularities in numerical work. The same object, written in Riccati amplitudes, provides the explicit variational functional whose saddle point reproduces the Riccati form of the Eilenberger equations.","core_discovery":"Starting from the Luttinger-Ward formulation, the paper evaluates the $\\lambda$-integral analytically and obtains the free-energy functional $\\Omega[\\hat g,\\hat\\Sigma] = \\frac12 E[\\hat g_1[\\hat\\Sigma]] + \\Phi[\\hat g] + \\frac12\\mathrm{Tr}[\\hat\\Lambda(\\hat g_n - \\hat g_1[\\hat\\Sigma])]$ before minimization, with the saddle-point value $\\Omega = \\frac12 E[\\hat g_*] + \\Phi[\\hat g_*] + \\frac12\\mathrm{Tr}[\\hat\\Lambda(\\hat g_n - \\hat g_*)]$. The central ingredient is the gradient functional $E[\\hat g] = \\frac12\\mathrm{Tr}\\bigl(\\hat g[\\hat\\tau_t,\\hat g]\\, v_F\\cdot\\check\\nabla[\\hat\\tau_t,\\hat g]^{-1}\\bigr)$, where $\\hat\\tau_t$ is any matrix field with $\\hat\\tau_t^2=1$. This functional is constructed so that its variation reproduces the gradient term of the Eilenberger equation, making the saddle point of the full functional equivalent to the quasiclassical equations plus the self-consistency relation. In the spin-singlet, spin-diagonal case it reduces to Eilenberger's original expression, while the same formula covers spin-triplet correlations; the paper also derives a Riccati-amplitude form and the diffusive-limit free energy in terms of momentum-averaged propagators. Where $[\\hat\\tau_t,\\hat g]$ is not invertible, the functional produces imaginary winding contributions, and the paper identifies the physical free energy with the real part of the expression.","pith_inferences":["The freedom in choosing $\\hat\\tau_t$ suggests a gauge-like tool: one can pick $\\hat\\tau_t$ locally so that $[\\hat\\tau_t,\\hat g]$ stays invertible, and a natural choice is to rotate $\\hat\\tau_1$ by the local transformation diagonalizing $\\hat g$, which the paper only sketches in the normal-state limit.","The Berry/Wess-Zumino reading of the gradient term hints that free-energy differences between quasiclassical configurations might be computable as geometric phases, which could simplify topological-state comparisons beyond the saddle-point regime.","A direct stress test would evaluate the real part of expression (19) for a nonuniform spin-triplet texture and compare it with the exact $\\lambda$-integrated expression (3) at the same saddle point, checking the real-part prescription exactly where $[\\hat\\tau,\\hat g]$ nearly vanishes.","The dirty-limit functional could be applied to Majorana nanowire models to see whether vortex-formation energetics change when using this microscopically derived free energy instead of the commonly used approximate forms."],"forward_implications":["For spin-singlet, spin-diagonal systems the new functional reproduces Eilenberger's original free energy, settling the earlier expression's origin as the saddle-point reduction of the Luttinger-Ward functional.","For spin-triplet states, including superfluid $^3$He and systems with spatially inhomogeneous exchange fields or spin-orbit coupling, the same closed form supplies a workable thermodynamic potential at the same level of rigor as singlet calculations.","In the diffusive limit, the derived expression (21) has the Usadel equation and the gap self-consistency relation as its saddle point, giving a microscopically derived free energy for dirty superconductors.","Previous variational functionals from nonlinear sigma-model approaches and the Serene-Rainer construction coincide with the Eilenberger-type expression derived here, unifying these formulations.","Because no $\\lambda$-integration remains, numerical minimization no longer requires solving the auxiliary propagator for many coupling constants, which simplifies phase-competition studies such as $0$-$\\pi$ junctions, FFLO states, and vortex structures."],"supporting_citations":[{"why":"Provides the original Eilenberger free-energy expression for spin-singlet superconductors that the present result generalizes.","marker":"[1]"},{"why":"Supplies the Serene-Rainer quasiclassical free-energy formulation with the $\\lambda$-integral that serves as the starting point.","marker":"[5]"},{"why":"Gives the Luttinger-Ward functional from which the derivation begins.","marker":"[16]"},{"why":"Provides the $\\lambda$-integrated form (3) and the auxiliary propagator equation (4) whose analytic evaluation is the paper's main step.","marker":"[11]"},{"why":"Introduces the Riccati parametrization used to find a candidate gradient functional and to derive the closed form (12).","marker":"[21]"},{"why":"Identifies the gradient term as a Berry/Wess-Zumino term, connecting the functional to a geometric phase.","marker":"[25]"},{"why":"Supplies the Usadel dirty-limit equation that the diffusive free-energy expression (21) has as its saddle point.","marker":"[3]"},{"why":"Represents the nonlinear sigma-model action whose relation to the variational functionals is clarified by the derivation.","marker":"[14]"}],"fun_headline_variants":["Free energy for all superconducting pairings, now unified","Quasiclassical free energy extended to spin-triplet pairing","From Luttinger-Ward to a universal superconductor free energy","Superfluid 3He included in new quasiclassical free energy","One free energy functional for singlet and triplet pairing"],"cache_read_input_tokens":14208,"weakest_assumption_plain":"The derivation assumes that the physical free energy equals the real part of the constructed functional, discarding imaginary winding contributions that arise when $[\\hat\\tau_t,\\hat g_\\lambda]$ fails to be invertible inside the $\\lambda$-integration domain.","fun_headline_variants_meta":{"raw":{"variants":["Free energy for all superconducting pairings, now unified","Quasiclassical free energy extended to spin-triplet pairing","From Luttinger-Ward to a universal superconductor free energy","Superfluid 3He included in new quasiclassical free energy","One free energy functional for singlet and triplet pairing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000335,"raw_usage":{"total_tokens":1898,"prompt_tokens":1026,"completion_tokens":872,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":642,"completion_tokens_details":{"reasoning_tokens":786}},"tokens_in":642,"tokens_out":872,"duration_ms":7835,"temperature":1.0,"reasoning_tokens":786,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:29:41.549028+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute a nonuniform spin-triplet or spin-singlet configuration where $[\\hat\\tau,\\hat g_\\lambda]$ has a zero for some $\\lambda\\in[0,1]$, evaluate the real part of (19), and compare it with the exact $\\lambda$-integrated Luttinger-Ward expression (3) at the same saddle point: if the real parts differ by a nonzero winding contribution, the real-part prescription fails.","supporting_citations":[{"cited_title":"The generalized trace operator in Eq","cited_arxiv_id":null,"evidence_quote":"Provides the original Eilenberger free-energy expression for spin-singlet superconductors that the present result generalizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Serene-Rainer quasiclassical free-energy formulation with the $\\lambda$-integral that serves as the starting point."},{"cited_title":"The func- tional generalizes the well-known Eilenberger free energy for the systems with arbitrary type of pairing and inter- acting with spin-dependent fields","cited_arxiv_id":null,"evidence_quote":"Gives the Luttinger-Ward functional from which the derivation begins."},{"cited_title":"( 12) with ˆτt = ˆτ3","cited_arxiv_id":null,"evidence_quote":"Provides the $\\lambda$-integrated form (3) and the auxiliary propagator equation (4) whose analytic evaluation is the paper's main step."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the gradient term as a Berry/Wess-Zumino term, connecting the functional to a geometric phase."},{"cited_title":"However, the remaining λ-integration necessitates solving Eq","cited_arxiv_id":null,"evidence_quote":"Supplies the Usadel dirty-limit equation that the diffusive free-energy expression (21) has as its saddle point."},{"cited_title":"Eilenberger, Z","cited_arxiv_id":null,"evidence_quote":"Represents the nonlinear sigma-model action whose relation to the variational functionals is clarified by the derivation."}],"review_version":1}