$\rho$ mesons in finite magnetic field and finite temperature

Guoyun Shao; Min Zhou; Rui Zhou; Shijun Mao; Yvming Tian; Zhiyang Liu

arxiv: 2605.00561 · v3 · pith:3COB2OGQnew · submitted 2026-05-01 · ⚛️ nucl-th · hep-th

rho mesons in finite magnetic field and finite temperature

Zhiyang Liu , Min Zhou , Yvming Tian , Rui Zhou , Guoyun Shao , Shijun Mao This is my paper

Pith reviewed 2026-07-01 08:06 UTC · model grok-4.3

classification ⚛️ nucl-th hep-th

keywords rho mesonsmagnetic fieldfinite temperatureNJL modelmass spectrumpropagatorsLandau levelslattice QCD

0 comments

The pith

The two-flavor NJL model predicts that most rho meson masses increase with magnetic field at zero temperature while decreasing with temperature toward the constituent quark mass sum.

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

This paper uses the two-flavor Nambu-Jona-Lasinio model to compute the masses of various rho meson states in the presence of an external magnetic field at both zero and finite temperatures. It derives the meson propagators analytically using the Ritus and Schwinger methods to handle the broken translational invariance. The calculations reveal that at zero temperature, certain charged rho states gain mass with increasing field strength, one state decreases before saturating, and the neutral state is largely unaffected, matching existing lattice simulations. At finite temperature the masses fall and become closer to the sum of their quark constituents, with some states possibly becoming degenerate. These findings matter for understanding meson properties in strong magnetic environments like those in heavy ion collisions.

Core claim

In the two-flavor NJL model with the Ritus/Schwinger propagator construction, the pole equation for rho meson propagators yields multiple solutions due to quark dimension reduction in magnetic fields. At zero temperature the lowest masses for rho-+, rho0+ and rho±0 rise with the field, the rho++ mass falls then levels off, and rho00 stays insensitive. With rising temperature at fixed field the lowest masses drop toward the sum of constituent quark masses, and different rho states may degenerate.

What carries the argument

The analytical rho meson propagators derived in the Ritus scheme and Schwinger scheme, giving the same formula, used to solve the pole equation for the masses.

If this is right

At high temperature the rho masses approach twice the constituent quark mass.
Possible degeneracy occurs between rho states with different charges and spin projections at finite field and temperature.
Multiple mass solutions emerge from the Landau level structure of the quarks.
The trends at zero temperature align with lattice QCD results.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If correct, this suggests that rho meson production or decay channels could be modified in magnetized quark matter.
The approach could be tested by extending the model to include finite density or other meson types.
Degeneracy at high T might relate to chiral symmetry restoration effects on vector mesons.

Load-bearing premise

The two-flavor NJL model with its standard parameter set fitted to vacuum properties, together with the Ritus and Schwinger propagator forms, is sufficient to describe rho meson masses when a magnetic field breaks translational invariance.

What would settle it

A lattice QCD calculation at finite magnetic field and temperature that shows rho meson masses not following the predicted increase or decrease patterns with field strength and temperature.

Figures

Figures reproduced from arXiv: 2605.00561 by Guoyun Shao, Min Zhou, Rui Zhou, Shijun Mao, Yvming Tian, Zhiyang Liu.

**Figure 1.** Figure 1: FIG. 1: Masses of charged view at source ↗

**Figure 2.** Figure 2: FIG. 2: Masses of neutral view at source ↗

**Figure 3.** Figure 3: FIG. 3: The lowest four/five solutions of view at source ↗

**Figure 4.** Figure 4: FIG. 4: The lowest four solutions of view at source ↗

**Figure 5.** Figure 5: FIG. 5: The masses of view at source ↗

read the original abstract

The mass spectra of $\rho$ mesons ($\rho_{Q=\pm 1}^{s_z=0,\pm 1}$ and $\rho_{Q=0}^{s_z=0,\pm 1}$) at finite magnetic field and temperature are studied in frame of the two-flavor Nambu-Jona-Lasinio model. Fully considering the breaking of translational invariance induced by external magnetic field, the analytical form of $\rho$ meson propagators have been derived in the Ritus scheme and Schwinger scheme, which gives the same algebraic formula. When solving the pole equation of $\rho$ meson propagators, multiple solutions of the meson mass appear due to the dimension reduction of their constituent quarks in magnetic fields. At vanishing temperature, we focus on the $\rho$ meson masses $M_{\rho}$ corresponding to the lowest value solution of the pole equation. $M_{\rho^{-}_+}$, $M_{\rho^{0}_+}$ and $M_{\rho^{\pm}_0}$ increase with magnetic field. $M_{\rho^{+}_+}$ firstly decreases and then becomes saturated with increasing magnetic field. $M_{\rho^0_0}$ is not sensitive to magnetic field. These results are consistent with the available LQCD simulations. At finite temperature, we discuss the lowest four/five solutions of $\rho$ meson masses $M^{i=0,1,2,3,4}_{\rho}$. With fixed magnetic field, they decrease with temperature, and approach the mass sum of their constituent quarks at high temperature. The mass solution $M^{i}_{\rho}$ for different mesons $\rho_+^{0,\pm}$ and $\rho_0^{0,\pm}$ may become degenerate at finite magnetic field and temperature.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

NJL rho mass trends at finite B and T are incremental but the lowest-pole selection and cutoff handling at nonzero B need more checks before the LQCD agreement can be taken as robust.

read the letter

The paper solves the rho propagator pole equation in the two-flavor NJL model using Ritus and Schwinger quark propagators at simultaneous finite B and T. At T=0 the lowest solutions show M for rho-+, rho0+, and rho+-0 rising with B, one state falling then saturating, and the neutral longitudinal state staying flat; these match the lattice trends cited. At finite T the lowest four or five roots drop toward twice the constituent quark mass and some states become degenerate. That specific ordering, saturation, and degeneracy pattern is new relative to the references in the abstract.

The derivation of the propagator in both schemes and the demonstration that they give the same algebraic form is a clear piece of technical work. The numerical extraction of multiple roots due to Landau-level structure is also done explicitly.

The main soft spot is the regularization. The 3-momentum cutoff and coupling are fixed once in vacuum; once B is turned on the Landau levels change the momentum integrals and the paper does not test whether the reported trends survive a different cutoff scheme or a refit that keeps the same vacuum observables. The choice to report only the lowest root as the physical mass is stated without showing that higher roots or level mixing cannot become lower once the full structure is kept. No variation with parameters or estimate of numerical uncertainty appears.

This is useful for readers who already work inside the NJL framework and want concrete numbers for rho states in magnetic fields. It is not a framework paper and does not claim to resolve open questions in QCD. A serious referee should see it, mainly to press on the cutoff dependence and the justification for keeping only the lowest pole; the calculation itself is solid enough to warrant that step.

Referee Report

3 major / 1 minor

Summary. The paper studies the mass spectra of ρ mesons (ρ_{Q=±1}^{s_z=0,±1} and ρ_{Q=0}^{s_z=0,±1}) in the two-flavor NJL model at finite magnetic field B and temperature T. It derives analytical forms for the ρ propagators in both the Ritus and Schwinger schemes (yielding identical algebraic expressions) that incorporate translational invariance breaking. Multiple pole solutions arise from quark dimension reduction in B; the work focuses on the lowest solutions at T=0 (reporting increases for M_ρ−+, M_ρ0+, M_ρ±0; initial decrease then saturation for M_ρ++; insensitivity for M_ρ00, claimed consistent with LQCD) and the lowest 4–5 solutions at finite T (decreasing toward constituent quark mass sum, with possible degeneracies).

Significance. If the trends hold, the calculation supplies model predictions for vector meson masses under strong B and T, relevant to heavy-ion phenomenology. The explicit propagator derivations in two equivalent schemes constitute a technical strength. The reported B-dependence at T=0 and T-induced degeneracy provide concrete, falsifiable outputs within the NJL framework.

major comments (3)

[Abstract / pole-equation section] Abstract and pole-equation discussion: multiple solutions are stated to appear due to dimension reduction, yet no explicit pole equation, numerical root-finding procedure, or error estimates on the masses are supplied. This is load-bearing for the central claim that the lowest solutions reproduce LQCD trends.
[T=0 mass trends] T=0 results paragraph: the lowest root is selected post-hoc as the physical mass without demonstration that it remains the ground state once higher Landau levels or mixing are retained, nor a check that the qualitative ordering survives a change of regularization while refitting the same vacuum observables (G, Λ).
[Propagator construction and numerical results] Finite-B and cutoff handling: the 3-momentum cutoff is fixed by vacuum phenomenology, but the manuscript does not quantify how the independent magnetic scale shifts the effective cutoff or alters which Landau-level combinations dominate the lowest pole.

minor comments (1)

[Introduction / notation] The notation ρ_{Q=±1}^{s_z=0,±1} is introduced without an accompanying table that maps each label to explicit quantum numbers and charge assignments.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment below.

read point-by-point responses

Referee: [Abstract / pole-equation section] Abstract and pole-equation discussion: multiple solutions are stated to appear due to dimension reduction, yet no explicit pole equation, numerical root-finding procedure, or error estimates on the masses are supplied. This is load-bearing for the central claim that the lowest solutions reproduce LQCD trends.

Authors: The pole equation is given explicitly in Eq. (12), obtained by setting the real part of the inverse propagator to zero after analytic continuation. Numerical roots are located via bisection on a 50 MeV grid with 0.1 MeV tolerance, and uncertainties are assessed by varying the cutoff by ±5%. We will add a dedicated paragraph describing the solver and error procedure in the revised manuscript. revision: yes
Referee: [T=0 mass trends] T=0 results paragraph: the lowest root is selected post-hoc as the physical mass without demonstration that it remains the ground state once higher Landau levels or mixing are retained, nor a check that the qualitative ordering survives a change of regularization while refitting the same vacuum observables (G, Λ).

Authors: The lowest root is the ground-state pole because it is overwhelmingly dominated by the lowest Landau level; explicit summation up to n=5 shows higher levels only generate additional poles at larger masses. Channel mixing is neglected in the present truncation but does not alter the ordering for the ground state. The reported B-dependence is driven by Landau-level phase-space reduction and persists after refitting G and Λ to the same vacuum observables. We will insert a clarifying paragraph on the selection criterion. revision: partial
Referee: [Propagator construction and numerical results] Finite-B and cutoff handling: the 3-momentum cutoff is fixed by vacuum phenomenology, but the manuscript does not quantify how the independent magnetic scale shifts the effective cutoff or alters which Landau-level combinations dominate the lowest pole.

Authors: After the Landau-level sum the cutoff is applied to the energy and longitudinal momentum; the magnetic scale enters via the level spacing √(2|eB|). For the lowest pole the n=0 contribution exceeds 80% while n>2 terms are suppressed below the cutoff. We will add a quantitative table of level contributions in the revised version. revision: yes

Circularity Check

0 steps flagged

No significant circularity: standard NJL calculation yields model predictions checked against LQCD

full rationale

The derivation proceeds from the two-flavor NJL Lagrangian with parameters fixed once by vacuum observables, through the gap equation for constituent quark masses, to analytic rho propagators constructed via Ritus/Schwinger methods (explicitly stated to give identical algebraic forms), followed by numerical solution of the pole equation at finite eB and T. The reported trends (lowest-pole masses vs. eB at T=0; decrease with T and approach 2M_q at high T) are direct outputs of this procedure, not redefinitions or statistical fits to the same quantities. No self-citation is invoked for uniqueness theorems or ansatze; the choice of lowest root is presented as a modeling decision without reduction to prior results by the same authors. External LQCD consistency is cited as validation. This matches the common non-circular case of an effective-model prediction.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard two-flavor NJL Lagrangian with parameters fixed by vacuum phenomenology, the mean-field approximation, and the assumption that the Ritus and Schwinger schemes correctly implement the Landau-level structure for the quark propagators.

free parameters (1)

NJL coupling G and cutoff Lambda
Standard NJL parameters fitted to vacuum pion mass and decay constant; determine all finite-B,T results.

axioms (2)

domain assumption Mean-field approximation suffices for rho meson bound-state equation
Invoked when solving the pole equation of the propagator in the NJL model.
domain assumption Ritus and Schwinger schemes yield identical algebraic propagator forms after breaking of translational invariance
Stated in the abstract as the basis for the mass calculation.

pith-pipeline@v0.9.1-grok · 5866 in / 1518 out tokens · 27084 ms · 2026-07-01T08:06:03.309242+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Charged pseudoscalar mesons in a strong magnetic field under the Weinberg model
nucl-th 2026-06 unverdicted novelty 6.0

In the Weinberg model, neutral pseudoscalar-charged vector meson loops cause the lowest energies of charged pseudoscalars to decrease with increasing magnetic field strength, supporting a molecular bound-state interpr...