Pith. sign in

REVIEW 2 major objections 6 minor 25 references

Non-equilibrium scalar fields at finite temperature and density

T0 review · 2 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

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

Pith's one-line read This paper derives a closed-form expression for the non-equilibrium statistical propagator of a scalar field coupled to a thermal reservoir at finite temperature and chemical potential, valid when the chemical potential is weak, and shows…

desk verdict The paper's advertised non-equilibrium mu propagator is not derived from its own intermediate steps: Eq. (76) does not follow from (72) plus (75), and it contradicts Eq. (42) already at mu=0. read the letter →

arxiv 2505.09104 v1 pith:K7FMGNOG submitted 2025-05-14 hep-th hep-ph

Sebastian Mendizabal This is my paper · ORCID
classification hep-thhep-ph PACS 11.10.Wx
keywords Schwinger-KeldyshformalismKadanoff-Baymequationsfinitetemperaturefieldtheorychemicalpotentialspectralfunctionstatisticalpropagatornon-equilibriumquantumscalarfields
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

The paper establishes a closed-form description of a real scalar field coupled to a thermal reservoir when a chemical potential $\mu$ is present, working in the Schwinger-Keldysh real-time formalism. Starting from the Kadanoff-Baym equations, the author derives the spectral function and the statistical propagator for both equilibrium and non-equilibrium situations. The central result is a compact non-equilibrium statistical propagator, valid for weak chemical potential, in which $\mu$ appears through a phase $e^{i\mu(t_1-t_2)}$, shifted Bose-Einstein factors, and symmetrized hyperbolic cotangents. If the derivation is right, it supplies the two-point building blocks needed to compute particle production, response functions, and transport coefficients in hot and dense scalar systems.

What carries the argument

The load-bearing object is the pair of Kadanoff-Baym integro-differential equations for the spectral function $G^-_{\mu,k}$ and the statistical propagator $G^+_{\mu,k}$. The paper solves the spectral function by Laplace transform and quasiparticle pole identification, defining the renormalized frequency $\Omega_k$ and decay width $\Gamma_k$. The statistical propagator is split into a homogeneous piece fixed by initial conditions and an inhomogeneous memory term built from two spectral-function integrals weighted by the symmetric self-energy. The decisive simplification is the approximation (68): for $\mu \ll \Omega_k$, the two chemical-potential-shifted quasiparticle branches $\Omega^\pm_{\mu,k}$ and $\Gamma^\pm_{\mu,k}$ are replaced by $\Omega_k$ and $\Gamma_k$, collapsing the pole pairs and allowing the frequency integral to be evaluated in closed form. That collapse yields the main formula, Eq. (76).

What would settle it

Evaluate the exact spectral function (56) with the two branches $\Omega^\pm_{\mu,k}$ and $\Gamma^\pm_{\mu,k}$ for a fixed $\mu$ and compare it with the collapsed single-pole form; whenever the two peaks are resolved, the approximation (68) fails. A direct numerical check is to integrate the Kadanoff-Baym equation (59) for a small but finite $\mu$ and compare the full two-time statistical propagator with Eq. (76), watching specifically the $e^{-\Gamma_k(t_1+t_2)/2}$ transient term.

Watch

Extended reading notes

Core claim

The central claim is that for $\mu \ll \Omega_k$ the non-equilibrium statistical propagator takes the closed form of Eq. (76): a cosine envelope $\cos(\Omega_k(t_1-t_2))$ times a phase $e^{i\mu(t_1-t_2)}$, with the transient term weighted by $f_{\rm eq}(\Omega_k+\mu)+f_{\rm eq}(\Omega_k-\mu)$ and decaying as $e^{-\Gamma_k(t_1+t_2)/2}$, plus a memory term weighted by $\frac12[\coth(\beta(\Omega_k+\mu)/2)+\coth(\beta(\Omega_k-\mu)/2)]$ and decaying as $e^{-\Gamma_k|t_1-t_2|/2}$. The spectral function simultaneously reduces to $G^-_{\mu,k}(y)=\frac{e^{i\mu y}}{\Omega_k}\sin(\Omega_k y)e^{-\Gamma_k|y|/2}$, so at this level the chemical potential acts as a phase. As $\Gamma_k(t_1+t_2)\to\infty$ the transient term drops out and the equilibrium statistical propagator computed earlier in the paper is recovered.

Load-bearing premise

The entire closed form rests on the assumption that the chemical potential is small enough that the two shifted quasiparticle energies and decay widths can be replaced by their $\mu=0$ values, and the paper does not quantify how small is small, even though the memory part of the propagator is sensitive to exactly those pole positions.

Editorial extensions

If this is right

  • For weak chemical potential, all two-point functions of a dense scalar field in a thermal bath can be written in closed form, so observables such as particle production rates can be computed without solving integro-differential equations.
  • The equilibrium limit is recovered as the transient term decays, connecting the non-equilibrium result to the standard KMS-constrained propagators.
  • The appearance of $\mu$ in the phase and in the combinations $f_{\rm eq}(\Omega_k\pm\mu)$ gives a concrete recipe for extending existing real-time zero-density computations to finite density.
  • The same Kadanoff-Baym machinery can be applied beyond the weak-$\mu$ regime, using the full two-branch spectral function (56) instead of the collapsed single-pole form.

Reading between the lines

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

  • A natural testable extension is to feed Eq. (76) into the standard expressions for particle number density or energy density and compare the $\mu$-dependent corrections with kinetic-theory or lattice expectations; the paper itself does not perform that comparison.
  • The main limitation shows up at intermediate $\mu$: when the two pole branches of the exact spectral function (56) are resolved, the single-pole approximation (68) is no longer valid, and observable differences should first appear in the transient, not in the $|t_1-t_2|$-dependent, part.
  • If the same phase-with-chemical-potential structure survives in gauge theories, it should modify Debye screening and transport coefficients in dense plasmas, since the statistical propagator enters every loop; this is a direction the paper outlines but does not carry out.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 6 minor

Summary. The manuscript derives non-equilibrium spectral and statistical propagators for a scalar field at finite temperature and chemical potential in the Schwinger-Keldysh real-time formalism. After reviewing the μ=0 Kadanoff-Baym solution from [24], the author extends the construction to a constant chemical potential, obtaining the spectral function G^{-}_{μ,k}(y) in Eq. (56) and, after a 'small chemical potential' pole collapse, the full statistical propagator G^{+}_{μ,k}(t1,t2) in Eq. (76). The advertised central result is a closed-form two-point function with a quasi-particle damping factor, an oscillatory phase e^{iμ(t1-t2)}, and Bose-Einstein/coth coefficients evaluated at Ω_k ± μ. The paper also discusses the analytic structure of the propagators and the role of KMS conditions.

Significance. The paper addresses a relevant problem: real-time non-equilibrium propagators in dense bosonic theories, with potential applications in early-universe and heavy-ion physics. The intended result extends the equilibrium chemical-potential formalism of Weldon [13] and the non-equilibrium μ=0 solution of Anisimov-Buchmüller-Drewes-Mendizabal [24]. Should Eq. (76) be correct, it would provide a compact, falsifiable building block for computing observables. Some internal checks are favorable: at μ=0 Eq. (76) reduces to the earlier result Eq. (42), and the late-time limit regains the equilibrium form. However, these checks do not compensate for the fact that the central formula is not obtained from the stated intermediate equations, and the validity range of the key approximation is unquantified. The paper is therefore not acceptable in its present form.

major comments (2)
  1. [Section IV, Eqs. (72), (75), and (76)] The derivation of the full statistical propagator is internally inconsistent. The text states that Eq. (76) is obtained by adding Eq. (72) and Eq. (75). For μ=0, Eq. (72) reduces to (1/Ω_k) cos(Ω_k y) coth(βΩ_k/2)[exp(-Γ_k(t1+t2)/2) - exp(-Γ_k|t1-t2|/2)] and Eq. (75) reduces to (1/Ω_k) cos(Ω_k y) exp(-Γ_k(t1+t2)/2). The sum therefore has transient coefficient 1+coth(βΩ_k/2)=2+2f_B(Ω_k) and memory coefficient -coth(βΩ_k/2), whereas Eq. (42) -- and Eq. (76) at μ=0 -- have transient coefficient f_B(Ω_k) and memory coefficient +(1/2)coth(βΩ_k/2). No choice of f_B≥0 reconciles these terms; both the sign and magnitude of the |t1-t2| term are wrong. This is a load-bearing algebraic problem, not a typographical nuance: it means Eq. (76) is not derived from the stated intermediate results. The step from Eq. (71) to Eq. (72) is also presented only as a 'lengthy calculation' and must be spelled out and corrected.
  2. [Section IV, Eq. (68)] The closed-form result Eq. (76) relies on the replacement Ω±_{μ,k}→Ω_k and Γ±_{μ,k}→Γ_k in Eq. (68), which collapses the two distinct pole branches of the spectral function (56) into a single pair. No quantitative condition is given for this small-μ limit, and the validity domain is not obvious: Ω_k, Γ_k, μ, and the time variables t1, t2 enter the exponential factors, so the errors induced by Eq. (68) in the transient and memory parts are controlled by different combinations (e.g., μ/Ω_k and μ t). The authors should state the expansion in μ and provide an estimate of the neglected terms in Eq. (76). Without this, the domain of applicability of the headline result is unestablished.
minor comments (6)
  1. [Section II, Eq. (17)] Equation (17) has ω_k without the square; the differential operator should be ∂²_y + ω_k².
  2. [Section II, Eq. (33)] Equation (33) reads coth(βω/3); from the KMS derivation this should be coth(βω/2).
  3. [Section II, Eq. (41)] Equation (41) assigns both ˙Δ_in=0 and ˙Δ_in=Ω_k in the same list; one entry should be ¨Δ_in=Ω_k.
  4. [Section III, Eq. (56)] Equation (56) writes sin(Ω_k) without its argument; it should be sin(Ω_k y).
  5. [Section III, after Eq. (44)] The sentence asserting that the KMS conditions (15)-(16) 'continue to apply' to Δ^{-}_{μ,k} is stated without derivation and appears to conflict with the later appearance of ω_k ± μ in Eqs. (49)-(51); please clarify.
  6. [Throughout] There are numerous LaTeX/rendering errors, including a literal '\cite{weldon}' in Section III, 'IIn' at the start of Section IV, a garbled subscript in Eq. (63), and inconsistent spelling of 'Kadanoff-Baym'; the manuscript needs a careful editorial pass.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the finite-mu propagator is an analytic extension of published Kadanoff-Baym solutions with no fitted inputs; the stated (72)+(75) to (76) algebra is an internal-consistency issue, not circularity.

full rationale

All load-bearing ingredients are either derived in the paper (Laplace-transform solution via Eqs. (19)-(22), spectral function at finite mu via Eqs. (44)-(57)) or taken from prior published work [24,25] as methods, not as conclusions that are then relabeled as predictions. No parameter is fitted to a subset of data, and no 'prediction' is defined in terms of the quantity it is supposed to explain. The small-mu replacement (68) is an approximation whose validity domain is unquantified, and the final expression (76) is presented as the sum of (72) and (75); whether that sum is algebraically correct is an internal-consistency question, not circularity. Self-citations are present but are not used to forbid alternatives or to import a uniqueness theorem, so they do not make the derivation circular.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The derivation relies on the standard Schwinger-Keldysh and Kadanoff-Baym framework and on the quasiparticle approximation; it introduces no new particles, forces, or conserved quantities. The main additional assumptions are a time-translation-invariant equilibrium reservoir, a constant-background chemical potential with unmodified KMS conditions, and weak chemical potential. No numerical fitting to data is performed, and Omega_k and Gamma_k are treated as inputs supplied by the user's self-energy.

free parameters (2)
  • Omega_k (quasiparticle energy)
    Defined in Eq. (21) from Re Pi^R via Omega_k approx omega_k + Re Pi^R(omega_k)/(2 omega_k); no numerical value is needed for the derivation, but the final propagators are functions of it.
  • Gamma_k (quasiparticle damping rate)
    Defined in Eq. (21) from -Im Pi^R(Omega_k)/Omega_k; treated as small and retained to first order in exponents and pole positions.
assumptions (5)
  • standard math Schwinger-Keldysh contour and Kadanoff-Baym equations provide the correct non-equilibrium equations of motion.
    The entire calculation is built on Eqs. (6)-(10), taken from the cited literature [4-7,11].
  • domain assumption The reservoir is time-translation invariant so that Pi_k(t1,t2)=Pi_k(t1-t2) and Delta-_k(t1,t2)=Delta-_k(t1-t2).
    Eqs. (11)-(12) in Section II; this neglects back-reaction and non-stationary reservoirs, but is standard for a thermal bath at leading order.
  • domain assumption Chemical potential enters as a constant background field via partial_0 -> partial_0 + i mu, and the KMS condition remains valid in the same form.
    Section III, after Eq. (44); stated but not derived, and questionable for a real scalar without a conserved charge.
  • domain assumption Weak chemical potential mu much less than Omega_k so that Omega+_mu,k -> Omega_k and Gamma+_mu,k -> Gamma_k.
    Eq. (68) in Section IV; no validity bound is given, and the closed form Eq. (76) depends on this collapse of pole structures.
  • domain assumption Quasiparticle approximation retains only first order in Gamma_k and replaces omega by Omega_k sign(omega) at the poles.
    Eqs. (36)-(38) in Section II; standard Breit-Wigner reduction but neglects off-shell contributions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Non-equilibrium scalar fields at finite temperature and density." pith.science (2026). https://pith.science/paper/K7FMGNOG

@misc{pith2026250509104,
  author       = {Pith},
  title        = {Pith review of: Non-equilibrium scalar fields at finite temperature and density},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/K7FMGNOG}},
  note         = {Machine review of arXiv:2505.09104}
}
read the original abstract

We study propagators in bosonic field theories at finite temperature and chemical potential using the Schwinger-Keldysh real-time formalism. The system is considered in contact with a thermal reservoir, allowing for a consistent treatment of both equilibrium and non-equilibrium situations. The chemical potential, associated with conserved charges, modifies the structure of the propagators and introduces features that require detailed analysis. We focus on how a finite chemical potential affects the analytic structure of the bosonic propagators, including changes in the position of poles and the structure of branch cuts. In our setup, the chemical potential enters the theory as a constant background field, which alters both the dynamics and the boundary conditions. This work provides a basis for understanding the behavior of bosonic fields in thermal and dense environments.

Figures

Figures reproduced from arXiv: 2505.09104 by the authors.

Figure 1
Figure 1. Keldysh contour in the complex plane where Δ >(𝑥1, 𝑥2) ≡ ⟨Φ(𝑥1)Φ(𝑥2)⟩ , (4) Δ <(𝑥1, 𝑥2) ≡ ⟨Φ(𝑥2)Φ(𝑥1)⟩ . (5) One can write the Schwinger-Dyson equation in the contour: (□1 + 𝑚2 )Δ𝐶(𝑥1, 𝑥2) + 𝑖 ∫︁ 𝐶 𝑑 4𝑥 ′Π𝐶(𝑥1, 𝑥′ )Δ𝐶(𝑥 ′ , 𝑥2) = −𝑖𝛿𝐶(𝑥1 − 𝑥2) , (6) where Π𝐶 is defined in the same way as in (3), and □1 = ∂ 2/∂𝑡2 1 − ∇2 . One can also define the symmetric and antisymmetric correlators, fort the Green functions and se… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

25 extracted references · 16 canonical work pages

  1. [24]

    Anisimov, W

    A. Anisimov, W. Buchmuller, M. Drewes, and S. Mendizabal, Annals Phys. 324, 1234 (2009), arXiv:0812.1934 [hep-th]

    arXiv 2009

  2. [13]

    H. A. Weldon, Phys. Rev. D 76, 125029 (2007)

    work page 2007

  3. [1]

    Kolb and M

    E. Kolb and M. S. Turner, The Early Universe (Addison-Wesley Publishing Company, 1994)

    work page 1994

  4. [2]

    Weinberg, Cosmology (2008)

    S. Weinberg, Cosmology (2008)

    2008

  5. [3]

    Le Bellac, Thermal field Theory (Cambridge University Press, 2008)

    M. Le Bellac, Thermal field Theory (Cambridge University Press, 2008)

    work page 2008

  6. [4]

    J. S. Schwinger, J. Math. Phys. 2, 407 (1961)

    1961

  7. [5]

    L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964). 16

    work page 1964

  8. [6]

    P. M. Bakshi and K. T. Mahanthappa, J. Math. Phys. 4, 1 (1963)

    work page 1963

Show all 25 references
  1. [7]

    P. M. Bakshi and K. T. Mahanthappa, J. Math. Phys. 4, 12 (1963)

    1963

  2. [8]

    Matsubara, Progress of Theoretical Physics 14, 351 (1955)

    T. Matsubara, Progress of Theoretical Physics 14, 351 (1955)

    1955

  3. [9]

    R. Kubo, J. Phys. Soc. Jap. 12, 570 (1957)

    1957

  4. [10]

    P. C. Martin and J. S. Schwinger, Phys. Rev. 115, 1342 (1959)

    1959

  5. [11]

    Kadanoff and G

    L. Kadanoff and G. Baym, Quantum Statistical Mechanics (2962)

  6. [12]

    Drewes, S

    M. Drewes, S. Mendizabal, and C. Weniger, Phys. Lett. B 718, 1119 (2013), arXiv:1202.1301 [hep-ph]

    2013 arXiv

  7. [14]

    H. E. Haber and H. A. Weldon, Phys. Rev. Lett. 46, 1497 (1981)

    1981

  8. [15]

    H. E. Haber and H. A. Weldon, Phys. Rev. D 25, 502 (1982)

    1982

  9. [16]

    Bernstein and S

    J. Bernstein and S. Dodelson, Phys. Rev. Lett. 66, 683 (1991)

    1991

  10. [17]

    V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Mod. Phys. Lett. A 19, 1341 (2004), arXiv:hep-ph/0406219

    2004 arXiv

  11. [18]

    K. M. Benson, J. Bernstein, and S. Dodelson, Phys. Rev. D 44, 2480 (1991)

    1991

  12. [19]

    Loewe, S

    M. Loewe, S. Mendizabal, and J. C. Rojas, Phys. Lett. B 617, 87 (2005), arXiv:hep- ph/0412392

    2005

  13. [20]

    Loewe, S

    M. Loewe, S. Mendizabal, and J. C. Rojas, Phys. Lett. B 635, 213 (2006), arXiv:hep- ph/0512042

    2006

  14. [21]

    Mallik and S

    S. Mallik and S. Sarkar, Eur. Phys. J. C 61, 489 (2009), arXiv:0901.1045 [hep-ph]

    2009 arXiv

  15. [22]

    Sasagawa and H

    S. Sasagawa and H. Tanaka, Phys. Rev. C 85, 045201 (2012), arXiv:1106.4379 [hep-ph]

    2012 arXiv

  16. [23]

    Sasagawa and H

    S. Sasagawa and H. Tanaka, (2011), arXiv:1105.1423 [hep-ph]

    2011 arXiv

  17. [25]

    Mendizabal and J

    S. Mendizabal and J. C. Rojas, Int. J. Mod. Phys. A 32, 1750126 (2017), arXiv:1603.01118 [hep-ph]. 17

    2017 arXiv

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.