🔬 physics

Effects of Thomas-Fermi and QED Polarization Approaches on the Equation of State of White Dwarfs

This paper compares three zero-temperature equations of state for white dwarfs—homogeneous relativistic electron gas, nonlinear Poisson–Thomas–Fermi, and QED longitudinal-polarization models—by directly applying them to Tolman–Oppenheimer–Volkoff calculations to quantify how screened models reduce pressure and alter mass-radius sequences relative to the homogeneous gas limit.

Original authors: Saulo Vitor Costa Ramalho, Sergio Barbosa Duarte, Gustavo Santos Vicente

Published 2026-10-13

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Saulo Vitor Costa Ramalho, Sergio Barbosa Duarte, Gustavo Santos Vicente

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Stars are the great furnaces of the universe, burning through their fuel until they can no longer sustain themselves. When a star like our Sun reaches the end of its life, it sheds its outer layers and collapses into a tiny, incredibly dense core known as a white dwarf. These objects are the size of Earth but pack the mass of a star, creating conditions of pressure and density that are impossible to recreate in any laboratory on our planet. Because we cannot build a white dwarf to test our theories, astronomers rely on mathematics to describe how this matter behaves. The key to understanding these stellar remnants lies in the equation of state, a set of rules that tells us how much pressure the material exerts at any given density. This pressure is what holds the star up against its own gravity, preventing it from collapsing further. For decades, scientists have used a standard model based on a uniform, smooth gas of electrons to calculate these rules, assuming the electrons are spread out evenly like a fog.

However, a new study by researchers at the Brazilian Center for Physical Research and the Rio de Janeiro State University suggests that this smooth fog picture might be missing some crucial details. In reality, the electrons in a white dwarf are not perfectly uniform; they swirl around the atomic nuclei, creating a complex, structured environment. The researchers set out to see how this structure changes the pressure inside the star. They compared three different ways of calculating this pressure. The first was the traditional, simple model of a uniform gas. The second was a more detailed approach called the Thomas-Fermi model, which treats the star as a collection of tiny, spherical cells, each containing a nucleus surrounded by electrons that adjust their density to screen out the electric charge. The third approach was the most sophisticated, using quantum field theory to account for how the electrons respond to the nucleus in a way that depends on their momentum, a method that includes effects from the vacuum of space itself.

The team did not try to smooth over the differences between these models or guess what happens at densities they could not calculate. Instead, they took the raw data from each method and fed it directly into the equations that describe how a star holds itself together. They ran these calculations for a theoretical, cold white dwarf made of carbon, a state the star would reach after it has cooled down completely, often called a black dwarf. By keeping the conditions identical for all three models, they could isolate exactly how the different ways of treating the electrons changed the final result. They found that the more complex models, which account for the structure of the atoms, consistently predicted lower pressure than the simple uniform gas model. The electrons, when allowed to arrange themselves around the nuclei, push back less forcefully than the smooth gas model suggests.

This difference in pressure might sound small, but it has a significant impact on the star's size and weight. Because the complex models predict less pressure, the stars they describe are smaller and less massive than those predicted by the simple model. The researchers found that the most advanced quantum field theory approach produced the largest changes. At lower densities, where the star's outer layers are, the difference was most pronounced. For example, at a density of 100 grams per cubic centimeter, the pressure in the quantum field theory model was only about 40 percent of the pressure in the simple uniform gas model. As the density increased toward the center of the star, the models began to agree more closely, but the quantum field theory model still predicted slightly less pressure than the others.

When these pressure differences were applied to the full structure of the star, the results showed a clear pattern. A white dwarf calculated with the quantum field theory approach would be significantly smaller than one calculated with the simple model. For a star with a mass of about 0.4 times that of our Sun, the radius calculated with the advanced model was about 8 percent smaller than the standard prediction. For more massive stars, the difference grew even larger, with the radius shrinking by over 20 percent compared to the simple model. The researchers noted that these differences are not just minor tweaks; they represent a fundamental shift in how we understand the balance of forces inside these dead stars. The study highlights that the way electrons screen the electric charge of the atomic nuclei is a critical factor that cannot be ignored if we want to know the true size and mass of a white dwarf.

The findings are particularly relevant for understanding the final, frozen state of a white dwarf, long after it has cooled into a black dwarf. In this cold state, the internal structure of the matter becomes the dominant factor in determining the star's properties. The study suggests that the simple, uniform gas model, which has been the standard for so long, overestimates the pressure and therefore overestimates the size of these stars. While the researchers did not claim to have solved every mystery of stellar evolution, their work provides a clearer, more accurate map of the relationship between mass and radius for cold, dense stars. It shows that the intricate dance of electrons around atomic nuclei, when described with the full power of quantum physics, leads to a universe of white dwarfs that are slightly more compact and slightly lighter than we previously thought. This level of precision will be essential for future studies that aim to understand the long-term evolution of these stellar remnants and the nuclear reactions that might occur within them as they cool over billions of years.

Technical Summary: Effects of Thomas-Fermi and QED Polarization Approaches on the Equation of State of White Dwarfs

Problem Statement
White dwarf stars are supported against gravitational collapse primarily by electron degeneracy pressure. While the uniform relativistic Fermi gas provides a leading-order description, the presence of a crystalline lattice of fully ionized nuclei necessitates a treatment of non-uniform electron distributions and Coulomb correlations. Standard approaches, such as the Feynman–Metropolis–Teller (FMT) construction, utilize the Thomas–Fermi (TF) approximation, which corresponds to the long-wavelength static limit of the longitudinal response of the degenerate electron gas. However, this local approximation neglects the momentum dependence of the electron-gas response. The authors investigate the structural consequences of replacing the local TF screening with a more rigorous, momentum-dependent static longitudinal polarization derived from Quantum Electrodynamics (QED) and Quantum Field Theory (QFT). The study focuses on the zero-temperature ( $T=0$ ) limit, relevant to the idealized "black dwarf" state of a carbon white dwarf, to isolate the cold semiclassical structural effects without the complications of thermal evolution or atmospheric modeling.

Methodology
The authors compare three distinct equations of state (EoS) for a zero-temperature carbon white dwarf:

Homogeneous Relativistic Electron Gas: A baseline model where the electron density is spatially uniform, ignoring lattice screening.
Nonlinear Poisson–Thomas–Fermi (TF) Model: A self-consistent cell model where the lattice is represented by spherical Wigner–Seitz (WS) cells containing a uniformly charged nucleus. The electron density is determined by solving the nonlinear Poisson equation with a local TF screening potential.
QED Longitudinal-Polarization Model: A construction based on linear-response theory within QFT. This model retains the full momentum dependence of the static longitudinal polarization function, $F(\bar{\mu}_{rel}, k)$ , which includes both vacuum polarization and medium (degenerate electron gas) contributions. The screened potential is derived by inverting the Fourier transform of the interaction kernel, incorporating the nuclear form factor for finite-sized nuclei.

The study employs a direct numerical approach using tabulated EoS data ranging from $\rho = 10^2$ to $10^{10}$ g cm $^{-3}$ . Crucially, no smoothing, lower-density continuation, or high-density extrapolation is applied; the integration stops strictly at the tabulated surface density. The authors solve the Tolman–Oppenheimer–Volkoff (TOV) equations for all three EoS tables using identical numerical settings and a common surface definition ( $\rho_s = 10^2$ g cm $^{-3}$ ). Comparisons are quantified via independent pressure ratios and percentage differences in mass and radius at fixed radius and fixed mass, respectively.

Key Contributions

Direct Comparison of Screening Mechanisms: The paper provides a rigorous, point-by-point comparison between the nonlinear self-consistent TF cell solution and the linear-response QED construction, isolating the effects of momentum-dependent screening.
Quantification of Structural Deviations: By solving the TOV equations with identical boundary conditions, the authors quantify how the different screening treatments propagate from the microphysical EoS to macroscopic stellar structure (mass-radius relations).
Isolation of Cold Limits: The analysis strictly adheres to the $T=0$ limit, avoiding the introduction of atmospheric models or thermal history, thereby focusing purely on the structural implications of the EoS formulation for an idealized black dwarf.

Results

Equation of State Ordering: Throughout the density range ( $10^2$ to $10^{10}$ g cm $^{-3}$ ), the pressure ordering is consistently $P_{QED} < P_{TF} < P_{Hom}$ . The screened models yield lower pressures than the homogeneous gas, with the QED model producing the lowest pressure.
Density Dependence: The relative effect of screening is strongest at lower densities (near the $10^2$ g cm $^{-3}$ cutoff) and progressively decreases as density increases. At $10^{10}$ g cm $^{-3}$ , the TF and QED models reduce the pressure by 1.53% and 2.47%, respectively, relative to the homogeneous model. The QED correction remains systematically larger than the TF correction, though the two descriptions converge at high densities.
Mass-Radius Relations: The EoS ordering propagates directly to the stellar structure. At a fixed radius, the screened models predict smaller masses than the homogeneous model; at a fixed mass, they predict smaller radii. The QED model exhibits the largest deviations.
- At $M \simeq 0.4 M_\odot$ , the TF and QED radii are reduced by 5.34% and 8.25%, respectively, compared to the homogeneous result.
- At $M \simeq 1.3 M_\odot$ , these reductions increase to 12.25% (TF) and 20.53% (QED).
Geometric Sensitivity: The authors note that the large percentage differences in radius near the upper end of the mass range are geometrically amplified by the near-vertical nature of the mass-radius curves in that region, where small displacements in the curve result in large radius changes at fixed mass.

Significance and Claims
The paper claims that the treatment of electronic screening produces measurable changes in the calculated mass-radius relation, particularly when the low-density portion of the EoS significantly influences the stellar radius. The QED approach, by retaining the full momentum dependence of the longitudinal polarization (including both vacuum and medium contributions), yields systematically larger structural corrections than the nonlinear Thomas–Fermi treatment.

The authors position these findings as relevant for applications requiring precise mass-radius relations for cold white dwarfs, specifically in the context of long-term evolution and studies of pycnonuclear processes in the black-dwarf limit. They explicitly state that their model does not constitute a physical atmosphere or a cooling history; rather, the calculated radii are defined at the common density cutoff. The study concludes that while the differences are most pronounced at low densities, the QED-based screening description offers a more complete structural correction within the assumptions of the zero-temperature, static, non-magnetic limit.

Technical Summary: Effects of Thomas-Fermi and QED Polarization Approaches on the Equation of State of White Dwarfs

More like this