Magnetized Quark Matter Grows Stiffer When Isospin Imbalance Enters the Equation

Deep inside the most extreme objects in the universe, matter exists in a state that no laboratory on Earth can sustain for more than an instant. In the cores of neutron stars, densities may climb so high that protons and neutrons dissolve into a soup of deconfined quarks, while magnetic fields can reach strengths billions of times greater than anything produced by terrestrial magnets. A new theoretical study published in The European Physical Journal C has now taken a significant step toward understanding how such matter behaves when two additional ingredients are combined: an overwhelming magnetic field and an imbalance between the populations of up and down quarks, known as isospin asymmetry. The work, carried out by Salman Ahamad Khan of Integral University, Sarthak Satapathy of BJB Autonomous College and NISER, and Sumit of the Beijing Institute of Technology, delivers the first one-loop calculation of the thermodynamics of cold quark matter that simultaneously accounts for all three effects within a single resummed perturbative framework.

The quantity at the heart of the calculation is the equation of state, the mathematical relationship that links the pressure of a substance to its density and other thermodynamic variables. For ordinary matter, the equation of state can be measured in a laboratory; for quark matter buried in a neutron star core, it must be computed from first principles using quantum chromodynamics, the theory of the strong force. This is a notoriously difficult task. Lattice simulations, the standard nonperturbative tool of QCD, fail at high densities because of the so-called sign problem, which renders the numerical weight functions oscillatory and unusable. Theorists must therefore rely on controlled approximations, and the new work adopts one of the most systematic of these: hard dense loop perturbation theory, or HDLpt.

HDLpt is the zero-temperature, high-density cousin of hard thermal loop perturbation theory, a resummation technique developed to tame the infrared divergences that plague ordinary perturbation theory in hot or dense gauge theories. The central idea is that at high density, collective excitations of the medium, quarks sloshing collectively near their Fermi surfaces, screen the gluon fields that mediate the strong interaction. By resumming these hard dense loop contributions into the quark and gluon propagators, the theory systematically incorporates screening and damping effects that would otherwise invalidate a naive expansion in the strong coupling constant. In the new paper, the authors derive the resummed propagators for a two-flavor system of up and down quarks, each carrying its own chemical potential, related to the quark chemical potential and the isospin chemical potential through the relation mu_u equals mu_q plus mu_I over two, and mu_d equals mu_q minus mu_I over two.

The magnetic field enters the calculation in a particularly dramatic way. When the field is strong enough, the transverse motion of charged fermions becomes quantized into discrete Landau levels, much like the energy levels of an atom. In the ultra-strong regime considered here, with field strengths between one and two square gigaelectronvolts, corresponding to fields approaching ten billion billion gauss, all quarks are confined to the lowest Landau level. This lowest-Landau-level approximation effectively reduces the dynamics from three spatial dimensions to one, since quarks can move freely only along the direction of the magnetic field. The authors show that this regime is self-consistent when the flavor chemical potential satisfies mu_f squared less than twice the absolute value of q_f B, a condition their chosen parameter ranges, with quark chemical potentials between 0.6 and 0.8 GeV and isospin chemical potentials between 0.5 and 0.7 GeV, carefully respect.

One of the technical triumphs of the paper lies in its careful handling of the quark propagator in the background field. The magnetic field breaks translational invariance in the plane perpendicular to it, introducing the famous Schwinger phase into the coordinate-space propagator and preventing a straightforward Fourier transform to momentum space. The authors demonstrate how a judicious gauge transformation in the symmetric gauge removes this phase, restoring translational invariance of the momentum-space propagator and allowing the self-energy to be decomposed into a set of form factors. These form factors, computed separately for the up and down quarks, encode how the medium modifies the propagation of quarks and gluons, and they carry the entire flavor-asymmetric fingerprint of the isospin imbalance.

With the machinery in place, the team computed the one-loop free energy, split into quark and gluon contributions, and from it derived the thermodynamic observables. The longitudinal pressure, directed along the magnetic field, rises monotonically with both the quark chemical potential and the isospin chemical potential, and it also grows with the magnetic field itself. The dependence on the quark chemical potential reflects the familiar increase of Fermi momentum and degeneracy pressure. The isospin dependence, however, is flavor-asymmetric: raising mu_I boosts the up-quark chemical potential while suppressing the down-quark one, and because the up quark carries the larger electric charge, its larger Landau degeneracy factor dominates in the parameter range studied. The result is a pressure that increases with isospin imbalance, a genuinely new prediction for the cold, magnetized regime.

Perhaps the most striking result concerns the magnetic response of the medium. The computed magnetization turns out to be positive everywhere in the parameter space, signaling that cold quark matter in a strong field is paramagnetic, its spins preferentially aligning with the external field. The authors attribute this to the dominance of spin-polarized quark states in the lowest Landau level over the diamagnetic orbital contribution. Positive magnetization has a profound geometric consequence: because the transverse pressure is given by the longitudinal pressure minus the product of the field strength and the magnetization, it is necessarily suppressed relative to the longitudinal pressure. The pressure of magnetized quark matter is therefore anisotropic, a phenomenon the authors connect to the well-known effect of paramagnetic squeezing, in which the matter is compressed along the field direction under the anisotropic stress.

The study also quantifies how important interactions are compared with an ideal gas of free quarks. The ratio of the HDL-resummed pressure to the ideal pressure exceeds one across the entire range of chemical potentials and field strengths considered, indicating that interactions stiffen the medium. As both chemical potentials grow, however, the ratio approaches unity, a manifestation of asymptotic freedom, the defining feature of QCD whereby the strong coupling weakens at high energy scales and dense quark matter gradually behaves more like a nearly free gas. This crossover behavior is consistent with earlier perturbative calculations of cold quark matter and provides a useful sanity check on the new flavor-asymmetric extension.

The astrophysical implications are considerable. Neutron star surfaces host fields of order ten to the fifteenth gauss, and magnetar cores may conceivably reach far higher values, while recent analyses of NICER X-ray observations of massive pulsars and of the LIGO-Virgo gravitational-wave signal from the binary neutron star merger GW170817 have sharpened the empirical constraints on the equation of state of dense matter. Recent theoretical work has even presented evidence that quark-matter cores may exist inside the most massive neutron stars. An equation of state that incorporates both ultra-strong magnetic fields and isospin asymmetry, the natural condition in stellar interiors where up and down quark densities differ, could feed directly into mass-radius calculations and stability analyses of compact stars. The authors caution that their numerical results are strong-field estimates valid only when higher Landau levels remain energetically suppressed, and that the regime where the down-quark chemical potential changes sign, where pion-condensed phases may appear, lies outside their framework and would require mesonic degrees of freedom.

Looking forward, the authors outline several natural extensions: incorporating nonzero quark masses into the thermodynamic functions, and reusing the quark and gluon structure functions derived here to compute the refractive index of cold quark matter, a quantity that governs how electromagnetic signals propagate through the densest matter in the cosmos. For now, the study stands as the first complete one-loop HDL description of magnetized, isospin-asymmetric cold quark matter, and it delivers a clear qualitative message: strong magnetic fields and flavor imbalance do not merely perturb the equation of state, they reshape it, stiffening the pressure along the field, softening it across, and endowing the medium with an unambiguous paramagnetic character that future neutron star models can no longer afford to ignore.

Subject of Research: One-loop hard-dense-loop thermodynamics of cold, strongly magnetized, isospin-asymmetric quark matter

Article Title: One-loop HDL thermodynamics of a strongly magnetized isospin asymmetric cold quark matter

Article References: Khan, S. A., Satapathy, S., & Sumit (2026). One-loop HDL thermodynamics of a strongly magnetized isospin asymmetric cold quark matter. The European Physical Journal C, 86(9), Article 1050. https://doi.org/10.1140/epjc/s10052-026-16198-7

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16198-7

Keywords: quark matter, quantum chromodynamics, hard dense loop, magnetic field, isospin chemical potential, neutron stars, equation of state, magnetization, pressure anisotropy, Landau quantization, perturbative QCD, compact stars