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

The new work, published in The European Physical Journal C by Kazem Bitaghsir Fadafan of Shahrood University of Technology, Jesús Cruz Rojas of the National Autonomous University of Mexico, and Jonas Mager of the Vienna University of Technology, applies holography, the correspondence that links strongly coupled quantum field theories to weakly coupled gravitational theories in higher dimensions. The central obstacle in describing quark matter inside neutron stars is that the relevant regime of quantum chromodynamics, or QCD, is strongly coupled. First-principles lattice simulations fail there because of the notorious fermion sign problem at large densities, while traditional phenomenological approaches such as nucleon effective field theories struggle under the extreme conditions found in stellar cores. Holographic models offer a way around this impasse: instead of computing the equation of state of dense quark matter directly, one solves a comparatively tractable problem in a higher-dimensional gravitational dual and reads off the properties of the quantum field theory from the geometry.

The researchers built their model on the well-studied D3/D7-brane configuration of type IIB string theory, a bottom-up construction in which D7-branes embedded in the curved background of D3-branes describe quarks moving in a gauge theory. The key innovation lies in the dilaton, a scalar field that controls the running of the gauge coupling in the dual theory. Rather than keeping the dilaton constant, as in the simplest version of the model, the team adopted a phenomenological profile that interpolates smoothly and monotonically from its ultraviolet value to a finite, regular value in the infrared. Three parameters, denoted A, lambda and kappa, govern the precise shape of this profile, controlling the infrared value of the dilaton, the scale at which conformal symmetry is broken, and the steepness of the transition. This tailored profile captures essential features of QCD such as chiral symmetry breaking and produces a deconfined yet massive quark phase at finite density.

In the holographic setup, the embedding function of the D7-brane encodes the quark mass and condensate, while a gauge field living on the brane describes the quark chemical potential and density. Two distinct phases emerge: a Minkowski, or vacuum, phase with broken chiral symmetry and zero baryon density, and a quark phase in which chiral symmetry is restored and the density is nonzero. By computing the on-shell action of the brane, which the holographic dictionary identifies with the grand canonical potential, the researchers obtained the pressure as a function of chemical potential and hence the equation of state of the quark phase. At asymptotically large densities the model correctly reproduces the perturbative QCD behavior, with the pressure scaling as the fourth power of the chemical potential.

The team also attempted to describe the baryonic phase within the same holographic framework, modeling nucleons as D5-branes wrapped on a five-sphere. The nontrivial dilaton profile stabilizes these wrapped branes against collapse, which is itself a notable achievement. However, in the homogeneous, smeared approximation used here, the resulting nuclear equation of state turned out to be unrealistic, predicting pressures in sharp conflict with phenomenology just above the onset of the baryonic phase. The authors showed analytically that near the transition the chemical potential scales as the one-third power of the density, forcing a fourth-order transition that cannot reproduce the first-order behavior expected of isospin-symmetric QCD. They therefore discarded the holographic baryon phase and instead adopted the phenomenological equations of state of Hebeler and collaborators, derived from chiral effective field theory and constrained by nuclear physics and observation.

With the baryonic phase fixed phenomenologically and the quark phase supplied by holography, the researchers scanned the parameter space of the model, varying the AdS radius and the t’Hooft coupling. They found a region of parameters, with the AdS radius between roughly 0.015 and 0.02 inverse megaelectronvolts and the t’Hooft coupling between 1.9 and 3, that produces rather stiff equations of state with only weakly first-order transitions between baryonic and quark matter. This softness of the transition is crucial: it allows a smooth conversion from nuclear matter to quark matter inside a star rather than a violent discontinuity that would destabilize the configuration. When combined with the stiff Hebeler equation of state, three parameter choices produced equations of state supporting stable quark cores, with phase transitions occurring at chemical potentials around 400 megaelectronvolts.

To connect the microscopic equation of state to observable stellar properties, the team solved the Tolman-Oppenheimer-Volkoff equations of general relativity, which determine the structure of non-rotating compact stars. The resulting mass-radius curves show stars that begin their lives as ordinary nucleonic objects and develop quark cores as the central density increases. Stability requires that the mass increase with central energy density, and the stars with quark cores satisfy this criterion all the way up to the maximum of the mass-radius curve, beyond which a radial mode becomes unstable and the star collapses into a black hole. The maximum masses obtained for the three quark-star-supporting parameter choices are 2.17, 2.13 and 1.90 solar masses, the largest of which is broadly consistent with the heaviest precisely measured neutron star mass of about 2.35 solar masses, albeit somewhat favoring lower values.

The authors also computed the tidal deformability, the quantity that measures how easily a star is deformed by the gravitational field of a companion and that leaves a characteristic imprint on gravitational wave signals from binary neutron star inspirals. Once a star becomes heavy enough to develop a quark core, the tidal deformability drops rapidly, falling by roughly 80 units between the onset of the quark phase and the maximum allowed mass. This rapid decrease is a potentially observable signature of quark matter formation. However, the very stiffness of the baryonic phase required to obtain stable quark cores causes the model to overshoot the tidal deformability constraint from the GW170817 event at 1.4 solar masses, predicting a value near 950 compared with the observational bound of about 190 with large uncertainties. The authors note that a baryonic equation of state intermediate between the stiff and medium phenomenological cases would likely ease this tension while still permitting stable quark cores.

Several other results add nuance to the picture. The speed of sound in the quark phase remains below the speed of light, although the stiff baryonic phase violates a recently derived transport bound. More strikingly, the polytropic index of the holographic quark phase exceeds the value of 1.75 that has been proposed in the literature as a criterion for the onset of quark matter, reaching values as high as 2.5 just after the transition. This suggests that the criterion may be too restrictive, and it contrasts with predictions from the V-QCD holographic model, which places transitions at much larger energy densities with smaller polytropic indices. The authors emphasize that their model cannot definitively establish whether stable quark stars exist in nature, since the answer depends on parameter choices, but it demonstrates that such objects can in principle emerge from holographic models, opening a computational window onto their phenomenology. Future work will explore localized D5-brane configurations, which may improve the low-density baryonic description, along with neutrino transport and rotating stars with quark cores, as multi-messenger observations continue to tighten the constraints on the densest matter in the universe.

Subject of Research: Holographic modeling of stable massive quark stars and the equation of state of dense quark matter in compact star cores

Article Title: Properties of stable massive quark stars in holography

Article References: Bitaghsir Fadafan, K., Cruz Rojas, J., & Mager, J. (2026). Properties of stable massive quark stars in holography. The European Physical Journal C, 86(9), Article 1094. https://doi.org/10.1140/epjc/s10052-026-16338-z

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16338-z

Keywords: quark stars, holographic QCD, neutron stars, D3/D7-branes, equation of state, chiral symmetry breaking, tidal deformability, compact stars, gauge-gravity duality, quark matter, TOV equations, dense nuclear matter

The research, carried out by Khushbu Zala and Sreemoyee Sarkar of SVKM’s NMIMS University in Mumbai, focuses on a phenomenon known as the r-mode instability. R-modes are non-radial oscillation modes in a rotating star that couple to gravity: as the stellar fluid sloshes back and forth, it emits gravitational radiation that carries away angular momentum. In the absence of any counteracting effect, these modes would grow without bound, spinning the star down dramatically in a runaway process. Yet astronomers observe that many millisecond pulsars—stars rotating hundreds of times per second—remain remarkably stable. Something inside them must be damping the oscillations, and the leading candidates are the two forms of viscosity that govern how dense matter dissipates energy: shear viscosity, which resists the sliding of adjacent fluid layers, and bulk viscosity, which dissipates energy when the fluid is periodically compressed and expanded.

What makes the new work distinctive is its treatment of a star whose core contains a mixed phase, a region where hadronic and quark matter coexist in thermodynamic equilibrium. Below a transition density of about 2.23 times nuclear saturation density, the outer layer is modeled as ordinary baryonic matter described by the relativistic mean-field equation of state known as DDME2. Above that threshold, the researchers assume a hybrid interior in which quarks and hadrons mingle. The problem, as the authors emphasize, is inherently ill-posed: the equation of state at supranuclear densities is poorly constrained, quantum chromodynamics cannot yet be solved reliably in the non-perturbative regime relevant to stellar cores, and the transport coefficients—the viscosities that control damping—are even harder to pin down from first principles.

To cut through these uncertainties, the team turned to Bayesian inference, a statistical framework that updates prior physical knowledge with observational data to produce posterior distributions for unknown parameters. The prior distributions for the key dimensionless coefficients—labeled S-tilde, V-tilde, W-tilde and J-tilde, which encode the shear and bulk viscous response of the mixed phase along with the star’s equilibrium structure—were drawn from existing calculations for neutron stars, strange stars and hybrid stars. The likelihood function combined two independent observational constraints: mass-radius measurements from NASA’s Neutron Star Interior Composition Explorer, or NICER, mission, and the spin-frequency and temperature observations of neutron stars in low-mass X-ray binaries. Using the UltraNest nested sampling algorithm, the researchers explored the parameter space efficiently, focusing on regions of high likelihood and estimating the Bayesian evidence for their model.

The formalism at the heart of the analysis describes how the amplitude of an r-mode evolves in time as an exponential whose decay constant is set by the competition between three timescales: gravitational radiation, which drives the instability, and shear and bulk viscosity, which suppress it. At low temperatures, shear viscosity dominates and stabilizes the star; at high temperatures, bulk viscosity takes over. Between these regimes lies a window where damping is least effective and the star is most vulnerable. The minimum of the instability curve—the lowest spin frequency at which the mode can grow—is therefore an exquisitely sensitive probe of the microphysics inside the star, including the equation of state and the weak-interaction processes, such as the direct Urca reaction, that generate bulk viscosity.

Applying this framework to two hybrid star configurations of 1.5 and 1.75 solar masses, the team obtained concrete estimates for the dissipation timescales. The shear viscous damping time came out as approximately 4.99 times ten to the eight, multiplied by the temperature to the five-thirds power, in seconds, while the bulk viscous timescale follows a more complex dependence on both temperature and spin, scaling inversely with the square of the angular velocity. From the inferred coefficients, the researchers calculated the minima of the instability curves: the critical angular velocity reaches its lowest value of about 451.87 hertz at a temperature of 0.259 megaelectronvolts for the 1.5 solar mass star, and 517.47 hertz at 0.234 megaelectronvolts for the 1.75 solar mass star. Normalized to the Kepler frequency, the maximum spin rate a star can sustain before mass shedding, these minima correspond to ratios of roughly 0.069 and 0.071.

Crucially, the resulting instability window does more than produce elegant numbers—it matches what astronomers actually see. When the team compared their inferred instability curves with the observed spin frequencies and temperatures of real millisecond pulsars, they found that the enhanced viscous dissipation from the mixed hadron-quark phase provides sufficient damping to explain the stability of several well-known objects. Among them are the accreting low-mass X-ray binary sources XTE J0929-314 and XTE J1807-294, and the radio millisecond pulsars J0437-4715 and J2124-3358. These stars all spin faster than 100 hertz, placing them squarely in the frequency range where r-mode physics matters, and all of them sit safely outside the region where the instability would grow—precisely as the two-layer hybrid model predicts.

The posterior distributions themselves carried informative structure. The corner plots of the inferred parameters showed moderate correlations among the shear and bulk viscous coefficients, reflecting the coupled role of the two viscosities in setting the instability boundary, while correlations involving the equilibrium parameter remained comparatively weak. Notably, the posterior contours for the more massive 1.74 solar mass configuration were narrower and more tightly localized than those for the lighter star, indicating that frequency-temperature observations constrain the dense-core dissipation properties more stringently in heavier hybrid stars. The team also found that the position of the instability minimum is remarkably robust: it barely shifts when the equilibrium and viscous parameters vary across their full credible intervals, suggesting that the result is not an artifact of statistical noise.

The broader significance of the work lies in its demonstration that r-mode phenomenology, combined with modern statistical inference, can serve as a practical observational tool for probing phase transitions in ultra-dense matter. Because pulsar rotational frequencies and their time derivatives are among the most precisely measured quantities in all of astrophysics, and because NICER continues to deliver mass-radius constraints, the approach links macroscopic observables directly to microscopic transport physics. If the inferred viscous properties of the mixed phase continue to align with observations, it would strengthen the case that some neutron stars genuinely harbor deconfined quark matter in their cores—a question that has remained open since the earliest theoretical speculations about quark stars.

The authors outline several directions for extending the framework. A more realistic three-layer stellar model could better capture the stratification of a hybrid star’s interior, and generalizing the formalism to derive quantitative constraints on the shear and bulk viscosities of each individual layer would sharpen the physical picture. Perhaps most ambitiously, they aim to characterize the nature of the hadron-quark phase transition itself—determining whether it is first order, second order, or a smooth crossover—by performing statistical inference on transport coefficients constrained by the gravitational-wave signatures that r-mode oscillations generate. As gravitational-wave detectors grow more sensitive, the faint hum of a wobbling hybrid star may one day confirm what this Bayesian analysis already hints at: that the universe’s most extreme matter hides its secrets in the way it dissipates motion.

Subject of Research: Bayesian inference of viscous dissipation timescales governing r-mode instability in hybrid stars with hadron-quark mixed phases

Article Title: Modelling dissipative dynamics of r-mode instability in hybrid stars

Article References: Zala, K., & Sarkar, S. (2026). Modelling dissipative dynamics of r-mode instability in hybrid stars. The European Physical Journal C, 86(9), Article 1064. https://doi.org/10.1140/epjc/s10052-026-16199-6

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16199-6

Keywords: hybrid stars, r-mode instability, neutron stars, bulk viscosity, shear viscosity, Bayesian inference, quark matter, millisecond pulsars, NICER, low-mass X-ray binaries, gravitational waves, dense matter