Generalized Linear Sigma Model

• The generalized linear sigma model is an effective hadronic framework that incorporates quark-antiquark and four‐quark (tetraquark) nonets with distinct U(1)A charges to capture anomaly effects.
• It employs effective Lagrangians with logarithmic anomaly terms to model meson masses, mixings, and decay widths, providing a robust phenomenological tool for low-energy QCD.
• Extensions of the model include glueball and dilaton fields, enabling studies of electromagnetic and trace anomalies, dense-matter effects, and heavy quarkonium–nucleus bound states.

Searching arXiv for the specified paper and closely related generalized linear sigma model work. arxiv_search.query({"search_query":"all:\"generalized linear sigma model\" AND (id:(Fariborz et al., 2017) OR title:\"Electromagnetic trace anomaly in a generalized linear sigma model\" OR title:\"Electromagnetic axial anomaly in a generalized linear sigma model\" OR title:\"Two-body decay widths of lowest lying and next-to-lowest lying scalar and pseudoscalar mesons in generalized linear sigma model\" OR title:\"Generalized linear sigma model with two glueballs\" OR title:\"Spinless glueballs in generalized linear sigma model\" OR title:\"Medium modifications of Heavy Quarkonia masses in a generalized Linear Sigma Model\" OR title:\"Heavy Quarkonium-nuclear bound states within a generalized linear sigma model\")","max_results":10,"sort_by":"relevance"})

Search completed. I’ll use the returned arXiv records corresponding to the specified GLSM papers as citations where needed. The generalized linear sigma model (GLSM) is an effective hadronic framework in which the field content of the linear sigma model is enlarged to accommodate additional QCD degrees of freedom and anomaly structures. In its standard low-energy $SU(3)_L\times SU(3)_R$ realization, the model contains two chiral nonets, one with quark–antiquark structure and one with four-quark content, and is used to analyze masses, mixings, scatterings, and decay widths of light scalar and pseudoscalar mesons. In other realizations, the GLSM is supplemented by scalar and pseudoscalar glueball fields or by a scalar dilaton field $\chi$ that simulates the gluon condensate, allowing the same general framework to address QCD axial and trace anomalies, in-medium heavy-quarkonium mass shifts, and quarkonium–nucleus bound states (Zebarjad et al., 2015, Mondal et al., 2022, Fariborz, 13 Aug 2025).

1. Field content and symmetry assignments

A central $SU(3)$ formulation introduces two $3\times 3$ matrix fields

$M=S+i\Phi,\qquad M'=S'+i\Phi',$

where $S,\Phi$ are the scalar and pseudoscalar nonets built of quark–antiquark degrees of freedom, while $S',\Phi'$ are the corresponding four-quark nonets. Under $SU(3)_L\times SU(3)_R$ they transform identically,

$$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$

but under the anomalous $U(1)_A$ they carry different charges, $\chi$0 and $\chi$1. This distinction is one of the defining structural features of the model, because it permits separate encoding of conventional quarkonium and four-quark components while preserving the same chiral transformation law (Fariborz et al., 2017).

In the low-energy phenomenological applications summarized in the literature, the $\chi$2 version is typically truncated at leading order, corresponding to keeping effective terms that contain no more than eight underlying quarks and antiquarks. This truncation is used in studies of two-body strong decays, two-photon decays, and meson spectroscopy, and it provides the basis for the frequently used “scenario 3I” identifications in the scalar and pseudoscalar isosinglet sectors (Zebarjad et al., 2015, Fariborz et al., 2017).

A distinct $\chi$3 realization is used in dense-matter applications. There the chiral multiplet $\chi$4 is combined with a scalar dilaton field $\chi$5, nucleon fields $\chi$6, and light-vector fields $\chi$7. Heavy quarkonium fields $\chi$8 do not enter at tree level; instead, their medium modifications are induced by the medium dependence of the dilaton field, which simulates the scalar gluon condensate of QCD (Mondal et al., 2024).

2. Effective Lagrangians and vacuum structure

In the two-nonet $\chi$9 model, a representative leading-order effective Lagrangian has the form

$SU(3)$0

with kinetic terms

$SU(3)$1

and covariant derivative

$SU(3)$2

The chiral-invariant potential is written, up to eight quark/antiquark lines, as

$SU(3)$3

while explicit chiral symmetry breaking is introduced through

$SU(3)$4

In equivalent formulations the same structure appears with slightly different normalizations, but always with two chiral nonets, explicit symmetry breaking by a quark-mass spurion, and a logarithmic $SU(3)$5-breaking term (Fariborz et al., 2017, Zebarjad et al., 2015).

The vacuum structure is specified by diagonal condensates,

$SU(3)$6

with the isospin limit $SU(3)$7, $SU(3)$8, $SU(3)$9. The extremum conditions $3\times 3$0 and $3\times 3$1 determine the condensates in terms of the Lagrangian parameters. In the pseudoscalar-anomaly analysis there are initially 12 unknowns, namely six couplings $3\times 3$2, two quark-mass parameters $3\times 3$3, and four condensates $3\times 3$4; the minimum conditions remove four of these, and the remaining eight are fixed by meson masses, $3\times 3$5, $3\times 3$6, and the trace and determinant conditions for the $3\times 3$7 pseudoscalar mass matrix (Fariborz et al., 2017).

In the $3\times 3$8 dense-matter version, the total Lagrangian is written as

$3\times 3$9

where $M=S+i\Phi,\qquad M'=S'+i\Phi',$0 contains the chiral- and scale-invariant kinetic and Yukawa terms, $M=S+i\Phi,\qquad M'=S'+i\Phi',$1 is the scale-anomaly potential, and $M=S+i\Phi,\qquad M'=S'+i\Phi',$2 contains the isovector–vector $M=S+i\Phi,\qquad M'=S'+i\Phi',$3 terms. Here the vacuum expectation values are $M=S+i\Phi,\qquad M'=S'+i\Phi',$4 and $M=S+i\Phi,\qquad M'=S'+i\Phi',$5, with $M=S+i\Phi,\qquad M'=S'+i\Phi',$6 (Mondal et al., 2022).

3. Axial, electromagnetic, and trace anomalies

One of the most distinctive aspects of the GLSM literature is the treatment of anomalies by nonderivative logarithmic terms. For the electromagnetic axial anomaly, the divergence of the $M=S+i\Phi,\qquad M'=S'+i\Phi',$7 current arising from $M=S+i\Phi,\qquad M'=S'+i\Phi',$8 vanishes identically in a linear realization, so an additional term is introduced to reproduce the Adler–Bell–Jackiw anomaly: $M=S+i\Phi,\qquad M'=S'+i\Phi',$9 A convenient form is

$S,\Phi$0

with $S,\Phi$1, and matching gives

$S,\Phi$2

An exactly analogous term $S,\Phi$3 is constructed for the four-quark nonet, and the combined model may be written as $S,\Phi$4 (Fariborz et al., 2017).

For the electromagnetic trace anomaly in the scalar sector, the minimal term is constructed by analogy with the ’t Hooft axial-anomaly term: $S,\Phi$5 with

$S,\Phi$6

This term reproduces $S,\Phi$7 and yields effective $S,\Phi$8 couplings after expansion around the vacuum expectation values (Fariborz et al., 2017).

In the dense-matter and glueball variants, the trace anomaly is instead encoded through a dilaton or scalar-glueball sector. In the $S,\Phi$9 model the scale-anomaly potential is

$S',\Phi'$0

with $S',\Phi'$1 for $S',\Phi'$2, $S',\Phi'$3, and

$S',\Phi'$4

In glueball-extended $S',\Phi'$5 versions, scalar and pseudoscalar glueballs are introduced through their connections with the trace and axial anomalies of QCD, respectively, with explicit anomaly terms $S',\Phi'$6 and $S',\Phi'$7 in the potential (Mondal et al., 2022, Fariborz, 13 Aug 2025).

4. Mixing, mass eigenstates, and spectroscopy

Because $S',\Phi'$8 and $S',\Phi'$9 transform identically under $SU(3)_L\times SU(3)_R$0, physical scalars and pseudoscalars are mixtures of two-quark and four-quark fields. The model therefore employs orthogonal rotations from bare fields to mass eigenstates. In the $SU(3)_L\times SU(3)_R$1 pseudoscalar sector, for example,

$SU(3)_L\times SU(3)_R$2

while in the isosinglet scalar sector one may define

$SU(3)_L\times SU(3)_R$3

and then

$SU(3)_L\times SU(3)_R$4

In the fit of Jora et al. designated “scenario 3I,” one identifies $SU(3)_L\times SU(3)_R$5 and $SU(3)_L\times SU(3)_R$6 (Zebarjad et al., 2015, Fariborz et al., 2017).

The same scenario is used in the pseudoscalar sector. Of six possible assignments for the heavier model-predicted $SU(3)_L\times SU(3)_R$7 states, the “scenario 3I”—identifying $SU(3)_L\times SU(3)_R$8, $SU(3)_L\times SU(3)_R$9 and taking solution I of the resulting quadratic—gives the best overall $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$0 to the mass spectrum and to the decay $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$1. For central inputs $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$2 and $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$3, the model yields $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$4, $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$5, $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$6, $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$7, $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$8, $$M\longrightarrow U_LMU_R^\dagger,\qquad M'\longrightarrow U_LM'U_R^\dagger,$$9, $U(1)_A$0, $U(1)_A$1, $U(1)_A$2, $U(1)_A$3, $U(1)_A$4, and $U(1)_A$5 (Fariborz et al., 2017).

Glueball-extended versions enlarge the isosinglet mass matrices further. In the 2025 spin-0 glueball analysis, the model contains two chiral nonets together with scalar and pseudoscalar glueballs. It is found that in order to satisfy the axial anomaly and at the same time accurately generate all seven eta masses, it is necessary to include at least two pseudoscalar glueballs, a physical one and an unphysical one that gets integrated out and yields an effective instanton-type term. In that analysis, the state that is dominantly made of glue in the pseudoscalar sector is a state with mass above $U(1)_A$6, while in the scalar sector the three isosinglets in the $U(1)_A$7 region can contain substantial glue, with the glue contents determined as functions of the scalar glueball condensate $U(1)_A$8 (Fariborz, 13 Aug 2025).

5. Strong and electromagnetic decays

The GLSM is extensively used to compute tree-level strong two-body widths and, after unitarization, physical pole positions. For a generic scalar decay into two pseudoscalars,

$U(1)_A$9

with

$\chi$00

Loop and final-state-interaction effects are incorporated through K-matrix unitarization of the partial-wave amplitude,

$\chi$01

and the pole positions give the physical masses $\chi$02 and widths $\chi$03. In the decay-width study of lowest lying and next-to-lowest lying scalar and pseudoscalar mesons, the two-body decay widths of lowest lying mesons are well predicted by this model while for the next-to-lowest lying mesons, only some of the decay widths agree with the experimental results (Zebarjad et al., 2015).

The comparison with the ordinary single-nonet linear sigma model is a recurring theme. In the decoupling limit $\chi$04, $\chi$05, or in the ordinary non-renormalizable single-nonet model, the predicted widths are much too large. The explicit comparison quoted for $\chi$06 is

$\chi$07

with corresponding masses $\chi$08, $\chi$09, and $\chi$10. The same analysis states that including two- and four-quark mixing in GLSM brings all lowest-lying widths into the correct ballpark (Zebarjad et al., 2015).

For pseudoscalar two-photon decays, the anomaly term yields

$\chi$11

Using the fitted parameters, the model predicts

$\chi$12

versus $\chi$13 experimentally, and also gives

$\chi$14

Apart from the slight overestimate of $\chi$15, the agreement is at the $\chi$16 level, and the results suggest significant four-quark mixing in the heavier pseudoscalars (Fariborz et al., 2017).

For scalar two-photon decays, the trace-anomaly term gives

$\chi$17

Fitting the well-measured width

$\chi$18

shows that agreement with experiment occurs only if $\chi$19 or $\chi$20, i.e. only when the direct mixing-sensitive term proportional to $\chi$21 is large. Once $\chi$22 is so constrained, the model predicts

$\chi$23

The same analysis states that “turning off” mixing, $\chi$24, renders both decay widths far too small (Fariborz et al., 2017).

6. Glueballs and anomaly-driven extensions

The inclusion of explicit glueball fields extends the GLSM beyond the two-nonet mesonic sector. In the model with two glueballs, one introduces a scalar glueball $\chi$25 and a pseudoscalar glueball $\chi$26 in addition to the nonets $\chi$27 and $\chi$28, and writes

$\chi$29

with

$\chi$30

Here $\chi$31 mocks up the axial $\chi$32 anomaly and $\chi$33 mocks up the QCD trace anomaly. In the decoupling limit, where glueball–quarkonia interactions are switched off, the analysis determines the properties of the pure scalar glueball and yields $\chi$34, stated to be in good agreement with lattice and QCD-sum-rule estimates (Fariborz et al., 2018).

The more detailed 2025 analysis with spinless glueballs shows that decay widths of isosinglet scalars as well as different self consistencies within this framework can be used to probe the glueball condensate. A histogram of successful decays versus $\chi$35 peaks sharply in the high-$\chi$36 band, and the combined mass-spectrum plus width analysis fixes the scalar-glue condensate to $\chi$37. In that treatment, $\chi$38, $\chi$39, and $\chi$40 all carry appreciable glue admixture, but the precise glue contents depend sensitively on $\chi$41 (Fariborz, 13 Aug 2025).

A plausible implication is that, within the GLSM family, the scalar sector above $\chi$42 is best viewed not as a single nonet problem but as a coupled system of quarkonium, four-quark, and glueball configurations whose relative weights depend on anomaly terms and vacuum condensates.

7. Dense matter, heavy quarkonia, and quarkonium–nucleus bound states

In nuclear-matter applications the GLSM is reformulated as an $\chi$43 chiral model with broken scale invariance incorporated through the dilaton field $\chi$44. The model relation

$\chi$45

is inserted into the leading-order heavy-quarkonium mass-shift formula,

$\chi$46

where the momentum-space wave functions are Fourier transforms of harmonic-oscillator eigenfunctions with size parameter $\chi$47 fixed to reproduce the vacuum r.m.s. radius of each quarkonium state. The studies report an appreciable mass drop in the heavy-quarkonium states under consideration and note that the in-medium masses at finite densities should modify the in-medium partial decay widths of heavy quarkonia to open heavy flavor mesons (Mondal et al., 2022).

The 2024 bound-state analysis converts the medium-modified dilaton field into quarkonium–nucleus potentials,

$\chi$48

and then solves the Klein–Gordon equation

$\chi$49

for bound-state eigenvalues $\chi$50. For six nuclei, $\chi$51, $\chi$52, $\chi$53, $\chi$54, $\chi$55, and $\chi$56, the study finds typical $\chi$57 binding energies of approximately $\chi$58, $\chi$59, $\chi$60, $\chi$61, $\chi$62, and $\chi$63 for $\chi$64, while $\chi$65 reaches approximately $\chi$66 in $\chi$67 and approximately $\chi$68 in $\chi$69. For bottomonium, $\chi$70 has no bound states in $\chi$71, $\chi$72, or $\chi$73, but binds in $\chi$74 by approximately $\chi$75 and in $\chi$76 by approximately $\chi$77. The same work states three key features: charmonia bind more deeply than bottomonia for the same nucleus, higher orbital and excited states produce substantially stronger binding than the ground $\chi$78, and binding energies grow with nuclear mass number $\chi$79 and saturate near Pb (Mondal et al., 2024).

These applications lie outside the traditional light-meson domain of the sigma model, but they preserve the same organizing idea: broken chiral symmetry and broken scale invariance are encoded in effective hadronic fields, and anomaly-sensitive condensates determine the observable spectrum. This suggests a broad interpretation of the GLSM as a family of QCD-motivated effective theories rather than a single fixed Lagrangian.