Magnetic Field Dependent Mass & Coupling Beta Functions

• Magnetic Field Dependent Mass and Coupling Beta Functions describe how effective masses and couplings evolve with an external magnetic field via modified renormalization group flows.
• The framework uses an environmentally friendly RG approach where the magnetic scale replaces traditional UV scales, leading to phenomena like magnetic catalysis and asymptotic freedom-type running.
• Applications in scalar models, QCD, and NJL theories reveal non-monotonic behavior in dynamical masses and effective couplings, with implications for meson properties and quark matter stability.

Magnetic field dependent mass and coupling beta functions describe the behavior of dynamically generated masses and effective interaction couplings as functions of an external magnetic field within quantum field theory. When quantum fields—be they scalars, fermions, or gauge bosons—are exposed to a magnetic background, the standard renormalization group (RG) running with respect to the ultraviolet (UV) scale is supplanted or supplemented by a flow with respect to the magnetic field strength. This environmental dependence modifies the structure of the relevant beta functions, the ordinary differential equations that dictate how couplings and masses evolve as the system's external conditions change. Such magnetic-field-induced running plays a central role in phenomena including magnetic (and inverse magnetic) catalysis, chiral phase transitions, and the stability of exotic matter in astrophysical and heavy ion collision environments.

1. Scalar Models: Environmentally Friendly Renormalization Group Analysis

The prototypical setting for magnetic-field-dependent beta functions is a model with three real scalar fields: a neutral field $\pi_0$ and two charged fields $\pi^+, \pi^-$ with Lagrangian

$$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$

The covariant derivative $D_\mu = \partial_\mu \mp iqA_\mu$ (for $\pi^\pm$), and the external uniform magnetic field $B$ is taken along $z$ with $A_y = Bx$.

The "Environmentally Friendly Renormalization Group" (EFRG) approach employs the magnetic scale $b^2 \equiv |qB|$ as the flow parameter instead of the traditional subtraction scale $\mu$. At one-loop order, Schwinger's proper-time representation is used for charged scalar propagators in the magnetic background. This framework enables computation of the $\pi^+, \pi^-$0-dependent mass and quartic coupling corrections for the neutral scalar via charged-loop vacuum polarization (Ayala et al., 22 Dec 2025).

The one-loop beta functions satisfied by the neutral scalar mass $\pi^+, \pi^-$1 and coupling $\pi^+, \pi^-$2 are:

$\pi^+, \pi^-$3

where $\pi^+, \pi^-$4 and $\pi^+, \pi^-$5 are calculated from proper-time integrals involving gamma and digamma functions; see Table 1 for explicit expressions.

Beta Function	Explicit Form
$\pi^+, \pi^-$6	$\pi^+, \pi^-$7
$\pi^+, \pi^-$8	$\pi^+, \pi^-$9

Numerical integration (e.g., fourth-order Runge-Kutta) shows that $$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$0 increases monotonically with $$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$1—reflecting a bosonic analogue of magnetic catalysis—while $$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$2 decreases, suggestive of asymptotic freedom-type running for the quartic self-interaction. Analytical solutions in the $$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$3 regime yield

$$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$4

and

$$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$5

for large $$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$6 (Ayala et al., 22 Dec 2025).

2. QCD Coupling and Dynamical Mass in Magnetic Backgrounds

In quantum chromodynamics (QCD), the application of a strong magnetic field reorganizes quark and gluon propagators into Landau levels, effectively replacing the standard RG scale $$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$7 by $$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$8 (the product of the elementary charge and magnetic field). The QCD coupling and its beta function become

$$\mathcal{L} = (D_\mu \pi^+)(D^\mu \pi^-) + \frac{1}{2} \partial_\mu \pi_0 \partial^\mu \pi_0 - \frac{1}{2} m^2 (2 \pi^+ \pi^- + \pi_0^2) - \frac{\lambda}{4} (2 \pi^+ \pi^- + \pi_0^2)^2$$9

and

$D_\mu = \partial_\mu \mp iqA_\mu$0

where $D_\mu = \partial_\mu \mp iqA_\mu$1 is the number of colors, $D_\mu = \partial_\mu \mp iqA_\mu$2 is the number of active flavors, and $D_\mu = \partial_\mu \mp iqA_\mu$3 is the QCD scale. This leading-order expression coincides with the vacuum result—except that now the RG flow is with respect to $D_\mu = \partial_\mu \mp iqA_\mu$4 rather than $D_\mu = \partial_\mu \mp iqA_\mu$5 (Li et al., 2016).

Practical model treatments, such as the Nambu–Jona-Lasinio (NJL) framework, implement a $D_\mu = \partial_\mu \mp iqA_\mu$6-dependent running of the four-fermion coupling $D_\mu = \partial_\mu \mp iqA_\mu$7 to match lattice and Dyson–Schwinger results:

$D_\mu = \partial_\mu \mp iqA_\mu$8

for SU(2), or

$D_\mu = \partial_\mu \mp iqA_\mu$9

for SU(3). The dynamical mass $\pi^\pm$0 then follows from a modified gap equation (Li et al., 2016):

$\pi^\pm$1 (GeV$\pi^\pm$2)	$\pi^\pm$3 (MeV)	$\pi^\pm$4 (MeV)	$\pi^\pm$5 (MeV)
0.00	330	330	469
0.05	355	310	460
0.10	370	290	450
0.20	395	260	430

For moderate to strong $\pi^\pm$6, the fall-off in $\pi^\pm$7 dominates over the Landau-level-driven enhancement, driving $\pi^\pm$8 non-monotonically. This produces a crossover from "magnetic catalysis" (mass increases with $\pi^\pm$9 at fixed $B$0) to "inverse magnetic catalysis" (mass decreases for running $B$1 at high fields) (Li et al., 2016).

3. Magnetic Field Beta Functions for Effective Couplings and Masses

One-loop background field calculations in the NJL model rigorously establish the $B$2-dependent corrections to the effective four-fermion coupling $B$3, as well as $B$4-dependent meson masses:

$B$5

For each flavor $B$6, the scalar-channel coupling is

$B$7

with $B$8 given by

$B$9

where $z$0 is determined self-consistently via the modified gap equation (Moreira et al., 2022).

The beta function $z$1 then involves both the explicit $z$2-derivative of the prefactor and the implicit derivative through $z$3, evaluated using the chain rule. For moderate $z$4, $z$5, confirming a monotonic decrease of the effective coupling with increasing field strength.

For neutral pion and kaon masses, the $z$6-dependent mass is obtained from the Bethe–Salpeter equation, and its flow is given by the implicit function theorem. The beta functions for the meson masses are negative for low–to–moderate $z$7: the mass decreases with increasing field, in agreement with lattice QCD for $z$8–$z$9 (Moreira et al., 2022).

4. Physical Interpretation and Phenomenological Consequences

The principal outcome of charged-boson or quark-loop vacuum polarization in a magnetic field is a modification of the effective scalar mass and interaction strength. In scalar models, $A_y = Bx$0 grows monotonically with $A_y = Bx$1 (“magnetic catalysis”), while $A_y = Bx$2 is suppressed, behaving as $A_y = Bx$3 for asymptotically large $A_y = Bx$4 (Ayala et al., 22 Dec 2025). In QCD and NJL settings, increasing $A_y = Bx$5 typically suppresses the four-fermion coupling, and, depending on the interplay with Landau level enhancement, can reverse magnetic catalysis (“inverse magnetic catalysis”). This impacts phase structure, catalysis of symmetry breaking, and stability properties of dense matter.

For oriented applications, such as the stability of strange quark matter and the free energy per baryon in the presence of $A_y = Bx$6:

• A running $A_y = Bx$7 can reduce free energy per baryon by $A_y = Bx$8 MeV at $A_y = Bx$9, potentially stabilizing strange quark matter ($b^2 \equiv |qB|$0 MeV) relevant for compact stars (Li et al., 2016).
• $b^2 \equiv |qB|$1-dependent corrections to meson masses closely track lattice findings at low field, with discrepancies at higher fields presumably arising from omitted higher-loop or non-perturbative physics (Moreira et al., 2022).

5. Extension to Gauge Theories and Lattice Validation

The EFRG and background field methods can be extended to theories with gauge and/or fermionic degrees of freedom:

• In QED, both electron and photon loops contribute; the RG scale $b^2 \equiv |qB|$2 set by $b^2 \equiv |qB|$3 enters Schwinger propagators for electrons and photons, with corresponding $b^2 \equiv |qB|$4 and $b^2 \equiv |qB|$5.
• In QCD, external $b^2 \equiv |qB|$6 drives the separation of color and flavor sectors, and Landau level sums become central for both gluon and quark polarization. The running coupling $b^2 \equiv |qB|$7 satisfies a flow equation $b^2 \equiv |qB|$8. Nonperturbative phenomena such as inverse magnetic catalysis at finite temperature demand going beyond leading order.
• Beta functions obtained from effective field theory or model analysis can be directly compared with lattice QCD computations; for instance, extracted effective $b^2 \equiv |qB|$9 and $\mu$0 agree quantitatively with lattice up to $\mu$1 (Moreira et al., 2022).

Quantity	Trend with $\mu$2	Field/Model
$\mu$3	Increases	Scalar EFRG (Ayala et al., 22 Dec 2025)
$\mu$4	Decreases	Scalar EFRG (Ayala et al., 22 Dec 2025)
$\mu$5	Decreases	NJL background field (Moreira et al., 2022)
$\mu$6	Non-monotonic	NJL gap eq. (Li et al., 2016)
$\mu$7	Decreases (modest)	NJL Bethe‐Salpeter (Moreira et al., 2022)
$\mu$8	Decreases	Quark matter (Li et al., 2016)

6. Limitations and Outlook

The primary limitation of current scalar and NJL-type analyses is the restriction to one-loop order and the neglect of gauge field fluctuations, higher-order corrections, and full gauge invariance. In strongly-coupled QCD, additional effects—such as gluon loop-induced inverse magnetic catalysis and flavor-dependent vacuum polarization—become relevant at large fields. The observed trends in $\mu$9-dependent beta functions, particularly the suppression of effective couplings and non-monotonic mass evolution, motivate further study in the context of gauge theories and at finite temperature or density. Extensions to non-Abelian gauge theories require worldline or proper-time techniques that preserve gauge invariance and systematically incorporate Landau-level dynamics (Ayala et al., 22 Dec 2025, Moreira et al., 2022, Li et al., 2016).

A plausible implication is that environmental RG approaches, using $\pi^+, \pi^-$00 as the flow variable, offer a unifying framework for analyzing quantum field responses to strong external conditions, with direct applicability to lattice QCD, heavy ion collisions, and compact star physics. Quantitative discrepancies at high field strengths and in the pseudoscalar sector highlight the need for higher-order and nonperturbative treatments. The calculable shifts in meson mixing angles, such as the $\pi^+, \pi^-$01–$\pi^+, \pi^-$02 angle, provide potential signatures of magnetic field effects in precision hadronic experiments (Moreira et al., 2022).