Parity Doublet Model in Hadronic Physics

• Parity Doublet Model is an effective hadronic framework characterized by mirror chiral assignments and a nonzero chirally invariant mass term (m₀).
• It explains baryon parity partner splitting through spontaneous chiral symmetry breaking and reproduces key nuclear matter properties.
• The model is applied to finite nuclei, dense matter thermodynamics, and compact stars, with parameters guided by empirical data and lattice QCD.

The parity doublet model (PDM) is an effective hadronic framework in which baryons and their opposite-parity partners are described within linearly realized chiral symmetry. Its defining feature is the presence of a chirally invariant mass term, $m_0$, so that baryons can remain massive even when the chiral condensate vanishes. In the chirally restored limit, the positive- and negative-parity states become degenerate at $m_0$, while the splitting between them is generated by spontaneous chiral symmetry breaking. This structure has been used in descriptions of vacuum hadron properties, nuclear matter saturation, finite nuclei, dense matter thermodynamics, neutron stars, hybrid stars, and baryon spectroscopy (Mukherjee, 8 Dec 2025, Kong et al., 2024).

1. Chiral structure and field content

In the two-flavor formulation, the PDM introduces two Dirac fields, often denoted $\psi_1$ and $\psi_2$, with mirror chiral assignment. A standard form of the baryonic Lagrangian is

$\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$

Under SU(2)$_L\times$SU(2)$_R$, the mirror assignment exchanges the chiral transformation properties of the second field relative to the first, and this assignment permits the chirally-invariant mass term $m_0$; by contrast, the naïve assignment forbids the $m_0$ term (Heinz, 2013, Mukherjee, 8 Dec 2025).

Many nuclear applications supplement the scalar and pion sector with vector mesons. In the extended parity doublet model with hidden local symmetry (HLS), the $\omega$ and $m_0$0 mesons are introduced in a chirally invariant way, and these vector interactions are crucial for nuclear saturation and for isospin physics in nuclei (Mun et al., 2018, Motohiro et al., 2015). In the PDHLS formulation, the mesonic sector is written in terms of the HLS Maurer–Cartan one-forms, and the nucleon sector contains both parity-doublet structure and vector couplings, providing a framework aimed at dense baryonic matter (Paeng et al., 2011).

2. Mass generation, parity doubling, and the chiral-invariant mass

After spontaneous chiral symmetry breaking, the physical positive- and negative-parity baryon masses are obtained by diagonalizing the parity-mixing mass matrix. A standard two-flavor expression is

$m_0$1

with $m_0$2 in the vacuum. In the chirally restored limit, $m_0$3, one has $m_0$4, so parity doubling means degeneracy at nonzero mass rather than vanishing baryon mass (Motohiro et al., 2015, Mukherjee, 8 Dec 2025).

Within nuclear-structure applications, $m_0$5 is interpreted as the part of the nucleon mass that does not arise from spontaneous chiral symmetry breaking. One study notes that its physical origin is possibly related to the QCD trace anomaly, tetraquark or gluon condensates, and summarizes previous estimates ranging from $m_0$6 from $m_0$7 decay up to $m_0$8 from nuclear matter incompressibility (Mun et al., 2018). In the PDHLS model, a tree-level fit to the decay width of the parity doubler $m_0$9 to nucleon-pion and nucleon axial coupling $\psi_1$0 gives $\psi_1$1 (Paeng et al., 2011). A later review of mean-field phenomenology and lattice constraints states that nuclear phenomenology requires $\psi_1$2 and that optimal compressibility and masses are typically found for $\psi_1$3, with lattice QCD also supporting $\psi_1$4 (Mukherjee, 8 Dec 2025).

The size of $\psi_1$5 directly controls medium dependence. In the extended HLS-based model for nuclei, increasing $\psi_1$6 decreases the couplings of the $\psi_1$7 and $\psi_1$8 fields to nucleons and reduces the contribution of spontaneous chiral symmetry breaking to the nucleon mass, leading to weaker attractive scalar and repulsive vector mean fields at saturation (Mun et al., 2018).

3. Mean-field realizations for nuclear matter and finite nuclei

A major line of work treats the meson fields as mean fields and derives thermodynamics or finite-nucleus structure self-consistently. In uniform matter, the model parameters are fixed from free-space hadron properties and from empirical nuclear matter data. One representative fit uses $\psi_1$9 MeV, $\psi_2$0 MeV, $\psi_2$1 MeV, $\psi_2$2 MeV, $\psi_2$3 MeV, $\psi_2$4 MeV, together with the saturation conditions $\psi_2$5 MeV, $\psi_2$6 fm$\psi_2$7, symmetry energy $\psi_2$8 MeV, and incompressibility $\psi_2$9 or $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$0 MeV (Mun et al., 2018).

The finite-nucleus formulation solves meson and Coulomb mean fields self-consistently together with the single-particle Dirac equation,

$\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$1

and computes the total energy $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$2 and binding energy per nucleon $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$3 (Mun et al., 2018).

For stable nuclei from $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$4O to $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$5Pb, the extended parity doublet model with HLS was studied for $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$6 between $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$7 and $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$8 MeV, while calculations do not converge for $\mathcal{L}_B = \overline\psi_{1} i\slashed{\partial} \psi_{1} +\overline\psi_{2} i\slashed{\partial} \psi_{2} -g_{1}\,\overline\psi_{1}(\sigma + i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{1} -g_{2}\,\overline\psi_{2}(\sigma - i\gamma_{5}\vec{\tau}\cdot\vec{\pi})\psi_{2} -m_0 (\overline\psi_1\gamma_5\psi_2 - \overline\psi_2\gamma_5\psi_1)\,.$9 MeV. The results approach the experimental values as $_L\times$0 is increased until $_L\times$1 MeV and start to deviate more from the experiments afterwards with $_L\times$2 larger than $_L\times$3 MeV. With further fine-tuning at $_L\times$4 MeV, the root-mean-square deviation becomes $_L\times$5 MeV for binding energies per nucleon and $_L\times$6 fm for charge radii, and the quality is comparable to relativistic continuum Hartree–Bogoliubov calculations with PC-PK1 (Mun et al., 2018). A 2024 review states in the same direction that finite nuclei studied with RCHB prefer $_L\times$7 MeV and that the symmetry energy is larger for smaller chiral invariant mass (Kong et al., 2024).

The same review also discusses an extension with the iso-vector scalar meson $_L\times$8. In that extension, the inclusion of the $_L\times$9 enlarges the symmetry energy of the infinite nuclear matter; because the $_R$0 contribution is attractive, keeping the empirical symmetry energy at saturation requires a stronger repulsive $_R$1 force, which in turn stiffens the high-density symmetry energy unless additional $_R$2–$_R$3 mixing is introduced (Kong et al., 2024).

4. Thermodynamics and phase structure

In mean-field thermodynamics, the PDM predicts two distinct transitions in nuclear matter: a liquid-gas phase transition at normal nuclear density and a chiral transition at higher density (Sasaki et al., 2010). In the two-flavor parity-doublet nucleon-meson model, the zero-temperature chiral transition in symmetric matter has been analyzed in detail and shown to be driven by a kind of symmetry energy that tends to equilibrate the populations of opposite parity baryons once the phase space for the negative-parity partner is opened (Eser et al., 2023). In that analysis, the parity-doublet model yields a first-order chiral transition at large baryon chemical potential, whereas the related singlet model obtained by disregarding the chiral partner gives a second-order transition in the chiral limit and a smooth crossover for realistic pion mass (Eser et al., 2023).

Isospin asymmetry changes both the liquid-gas and chiral sectors. In the HLS-based extension with a six-point $_R$4 interaction, the first order phase transition for the liquid-gas phase transition disappears in asymmetric matter, and the critical density for the chiral phase transition at non-zero density becomes smaller for larger asymmetry (Motohiro et al., 2015). A later study of asymmetric and neutron matter reformulates the gap equations as ordinary differential equations in baryon density and shows that the threshold for the onset of the population of the chiral partners is exclusively determined by the fermionic parameters, most notably by the chiral-invariant mass of the nucleon. That work also underlines the role of a parity symmetry energy in driving the equilibration of the nucleons and their parity partners across the chiral transition (Eser et al., 2024).

The model also supports inhomogeneous condensates. For the chiral density wave ansatz

$_R$5

the parity doublet model coupled to the linear sigma model including vector mesons admits a homogeneous ground state of nuclear matter at lower baryon chemical potential, but at larger baryon chemical potential the CDW is favored with respect to the homogeneous phase. For physical pion mass, the model exhibits a homogeneous broken phase at low $_R$6, an intermediate homogeneous phase with smaller condensate around $_R$7 MeV, and then a first-order transition to the CDW phase; in the chiral limit, the CDW is always favored at nonzero $_R$8 (Heinz, 2013).

A significant technical refinement concerns baryonic vacuum fluctuations. A multiplicatively renormalizable mean-field treatment that includes the baryonic vacuum contributions in an explicitly renormalization-group invariant form finds that these vacuum contributions smoothen the reduction of the chiral condensate with increasing density or temperature, move the chiral transition to higher $_R$9, and change its order from first-order to crossover for most phenomenologically allowed $m_0$0; only for $m_0$1 MeV does a first-order transition persist in that analysis (Recchi et al., 10 Nov 2025).

5. SU(3) extensions, quarks, and compact stars

The SU(3) parity-doublet framework generalizes the model from nucleons to the baryon octet and their parity partners. A standard SU(3) effective-mass formula is

$m_0$2

so that chiral restoration corresponds to the degeneracy of parity partners while retaining a nonzero invariant mass (Dexheimer et al., 2012, Schramm et al., 2011). In the hadronic SU(3) parity-doublet model, every baryon, including hyperons, has a positive- and negative-parity counterpart; depending on the masses of the chiral partners, the transition to the chirally restored phase can show a first-order line with critical endpoints, in addition to the standard liquid-gas phase transition of self-bound nuclear matter (Steinheimer et al., 2011).

Several works extend the model further to quarks and deconfinement. In the SU(3) parity model with quarks and Polyakov-loop dynamics, chiral restoration and deconfinement are decoupled, and the deconfinement transition is modeled as a crossover with excluded-volume suppression of hadrons at high density and temperature (Steinheimer et al., 2011). In hybrid-star applications, the extended SU(3) parity model containing quark degrees of freedom realizes chiral symmetry restoration inside the star, with chiral partners appearing and their masses becoming degenerate, and deconfinement occurring through the appearance of quark degrees of freedom. One implementation yields a maximum cold neutron-star mass of about $m_0$3 in its stiffer parametrization (Dexheimer et al., 2012).

The quark-hadron chiral parity-doublet model (Q$m_0$4P) pursues a unified equation of state for hadrons and quarks. In the 2017 implementation, the symmetry energy and its slope at saturation are $m_0$5 MeV and $m_0$6 MeV, the mass-radius relation accommodates massive and small neutron stars, and the radius of a $m_0$7 star is about $m_0$8 km. That model also emphasizes early quark appearance, suppression of hyperons, and absence of hadronic direct Urca without fine-tuning (Mukherjee et al., 2017).

A later review of neutron-star matter based on a nucleonic parity-doublet model coupled to $m_0$9, $m_0$0, and $m_0$1 fields, and to strange mesons through the U(1)$m_0$2 anomaly, argues that large $m_0$3 makes the nucleon mass insensitive to the medium and allows nuclear saturation without demanding strong $m_0$4 and $m_0$5 couplings. By confronting hadronic equations of state with nuclear constraints, neutron-star observations, and interpolated quark-matter constraints, that study delineates a range $m_0$6 and finds that the U(1)$m_0$7 anomaly softens equations of state from low to high density (Minamikawa et al., 2023).

6. Spectroscopy, axial structure, and open theoretical issues

The parity-doublet idea has also been extended from thermodynamics to baryon spectroscopy. A central issue is the choice of chiral representations in SU(3)$m_0$8SU(3)$m_0$9. One construction classifies representations by good and bad diquarks, treating $\omega$0 and $\omega$1 as soft and $\omega$2 as hard. In that framework, first order terms in the meson field $\omega$3 do not reproduce the mass hierarchy correctly although the Gell-Mann–Okubo relation is satisfied, while second-order terms reproduce the masses of the positive parity channels well up to the first radial excitations but leave problems in negative-parity mass ordering (Minamikawa et al., 2023). A subsequent octet model including $\omega$4, $\omega$5, and $\omega$6 reports that the ground state baryons are well dominated by the first and third representations, while the excited states require $\omega$7 and its bad-diquark content (Gao et al., 2024). Another SU(3) construction excludes $\omega$8 and states that the $\omega$9 representation containing symmetric bad diquarks is essential for reproducing the correct baryon mass hierarchy, particularly the $m_0$00–$m_0$01 mass ordering (Gao et al., 1 Dec 2025).

The axial sector exposes a separate limitation. In the standard PDM with mass mixing, the nucleon axial charge satisfies

$m_0$02

so finite $m_0$03 implies $m_0$04, whereas phenomenologically it is about $m_0$05. An extended model with kinetic-mixing terms corresponding to meson-baryon derivative couplings introduces two additional parameters, two mixing angles, and axial-charge expressions that allow $m_0$06 and decouple $m_0$07 from $m_0$08 (Kummer et al., 3 Dec 2025).

Beyond mean field, current reviews emphasize several unresolved points: the identification of the negative-parity partner is uncertain; extracting $m_0$09 and vector couplings requires global analyses using lattice, nuclear, and astrophysical data; and functional renormalization group studies indicate that fluctuations smoothen and shift the chiral transition (Mukherjee, 8 Dec 2025). In the PDHLS model, one-loop renormalization-group equations reveal a fixed point corresponding to the dilaton limit, at which vector mesons decouple from nucleons before the vector manifestation fixed point is reached (Paeng et al., 2011). This suggests a distinctive dense-matter trajectory in which parity doubling, vector decoupling, and the survival of a substantial chiral-invariant baryon mass are intertwined.

Within this body of work, a persistent misconception is explicitly excluded: the PDM does not predict that baryon masses vanish when chiral symmetry is restored. Instead, it predicts that the mass splitting generated by the condensate disappears, while a nonzero chirally invariant mass remains. The quantitative value of that invariant mass, and the degree to which it controls nuclear structure, dense-matter thermodynamics, and baryon spectroscopy, remain the central organizing questions of the parity-doublet program (Mukherjee, 8 Dec 2025, Mun et al., 2018).