Energy-Momentum Tensor and D-term of Baryons in Top-down Holographic QCD

This paper numerically computes the energy-momentum tensor of baryons in top-down holographic QCD, revealing a D-term value of approximately -2.05 that is significantly larger in magnitude than previous findings.

The Gist

Imagine the proton (a building block of atoms) not as a tiny, solid marble, but as a bustling, invisible city. Inside this city, there are tiny particles zooming around, pushing and pulling against each other. For a long time, physicists have been able to map out the "population" (where the mass is) and the "traffic" (how the particles spin). But there was one missing piece of the map: the pressure and the internal forces holding the city together.

This paper is like a team of cartographers using a very advanced, futuristic telescope to finally map out those invisible forces inside the proton.

Here is the breakdown of what they did, using some everyday analogies:

1. The "Holographic" Telescope

The scientists used a theory called Holographic QCD. Think of our 3D world (where the proton lives) as a flat piece of paper. Now, imagine that the complex physics happening on that paper is actually a shadow cast by a 3D object standing in a higher dimension.

Instead of trying to solve the incredibly messy math of the "flat paper" (our real world), the authors stepped into the "3D room" (a 5-dimensional universe in their math model). In this higher dimension, the messy quantum particles turn into smooth, stable shapes called solitons (think of them like a permanent, swirling knot in a rope). It's much easier to study the shape of a knot than to track every single fiber of the rope individually.

2. The "D-Term": The Proton's Internal Stress

The main thing they wanted to measure is called the D-term.

• The Analogy: Imagine a water balloon. The water inside pushes out (pressure), and the rubber skin pulls in (tension). The D-term is a single number that tells you how strong that internal tug-of-war is.
• Why it matters: If the D-term is too weak, the balloon pops (the proton falls apart). If it's too strong, it implodes. It's the "glue" that keeps the proton stable.

3. The Mistake in the Old Map

In a previous study (referenced as [1] in the paper), scientists tried to draw this map. However, they took a shortcut.

• The Analogy: Imagine trying to draw a mountain range. The previous team looked at the peak (the center of the proton) and the base (the edge), and they just drew a smooth, straight line connecting them. They assumed the middle was just a gentle slope.
• The Problem: They missed the jagged cliffs and deep valleys in the middle. Because of this, their calculation of the internal stress (the D-term) was way off. They thought the stress was very low (about -0.14).

4. The New, Accurate Map

The authors of this paper decided to do the hard work. Instead of guessing the middle, they used a supercomputer to solve the equations of motion step-by-step, filling in every jagged cliff and valley of the "mountain."

• The Result: They found that the middle of the proton is much more turbulent than they thought. The internal forces are much stronger.
• The New Number: Their new calculation for the D-term is -2.05.

Why is this a big deal?
The absolute value of their number is 15 times larger than the previous guess. It's like realizing that the water balloon isn't just gently taut; it's being squeezed by a hydraulic press. This changes our understanding of how the proton holds itself together.

5. What Else Did They Find?

Besides the stress (D-term), they mapped out other features of this "proton city":

• Energy Density: Where the "heavy stuff" is packed.
• Pressure: Where the particles are pushing out (positive pressure) vs. pulling in (negative pressure). They found that the center pushes out, but the edges pull in, creating a stable balance.
• Size: They calculated the "radius" of the proton based on these forces, finding it to be about 0.94 femtometers (a femtometer is a quadrillionth of a meter).

The Bottom Line

This paper is a "quality control" update for our understanding of the proton. By using a more rigorous mathematical method (solving the equations properly instead of approximating them), the authors discovered that the internal forces holding a proton together are significantly stronger and more complex than we previously believed.

It's a reminder that in physics, the devil is in the details, and sometimes, the most important forces are hiding in the middle of the "mountain," not just at the peak or the base.

Technical Summary

Here is a detailed technical summary of the paper "Energy-Momentum Tensor and D-term of Baryons in Top-down Holographic QCD" by Shigeki Sugimoto and Taichi Tsukamoto.

1. Problem Statement

The paper addresses the calculation of the gravitational form factors (GFFs) of baryons, specifically focusing on the D-term (also known as the "druck term").

• Context: GFFs describe the internal distribution of mass, spin, pressure, and shear forces within hadrons. The D-term is a fundamental quantity related to the mechanical stability and internal force distribution of the particle, often called the "last unknown global property" of the nucleon.
• Previous Limitations: In a previous study by the same group (Fujita et al., 2022), the D-term for the nucleon in top-down holographic QCD (Sakai-Sugimoto model) was estimated to be approximately $D \approx -0.14$. However, this calculation relied on an approximation where the soliton configuration was constructed by smoothly connecting solutions near the center and the boundary, rather than solving the full equations of motion (EOMs).
• Goal: The authors aim to improve the precision of this calculation by numerically solving the full EOMs for the soliton solution. They hypothesize that the D-term is highly sensitive to the gauge field configuration in the intermediate region, which the previous approximation failed to capture accurately.

2. Methodology

The study utilizes the Sakai-Sugimoto model, a top-down holographic QCD model derived from Type IIA string theory with $N_c$ D4-branes and $N_f$ D8-$\overline{\text{D8}}$ branes.

• Theoretical Framework:

  • Baryons are modeled as topological solitons (instantons) in a 5-dimensional $U(N_f)$ Yang-Mills-Chern-Simons theory.
  • The low-energy effective action includes a Yang-Mills term and a Chern-Simons term on a curved 5D background.
  • The Energy-Momentum Tensor (EMT) is derived via the variation of the action with respect to the metric (AdS/CFT prescription).

• Ansatz and Reduction:

  • The authors employ Witten's ansatz to reduce the 5D non-Abelian gauge field to a 2D Abelian-Higgs system on the $(r, z)$ plane (where $r$ is the radial coordinate in 3D space and $z$ is the holographic coordinate).
  • The fields are decomposed into an $SU(2)$ part (instanton-like) and a $U(1)$ part (related to the baryon number current).
  • The system is reduced to solving coupled non-linear partial differential equations for the scalar field $\Phi$, the gauge fields $A_r, A_z$, and the scalar potential $a_0$.

• Numerical Implementation:

  • Coordinate Transformation: The holographic coordinate $z$ is transformed to $w = \arctan(z)$ to handle the infinite domain, mapping $z \in (-\infty, \infty)$ to $w \in (-\pi/2, \pi/2)$.
  • Discretization: The $(r, w)$ plane is discretized on a $30 \times 300$ lattice.
  • Solver: The equations of motion are solved numerically using the Gauss-Seidel method.
  • Boundary Conditions: Strict boundary conditions are imposed at $r=0$, $r \to \infty$, and $w = \pm \pi/2$ to ensure regularity, finite energy, and the correct topological charge (baryon number $B=1$).
  • Parameters: The model parameters are fixed to reproduce experimental values: $M_{KK} = 949$ MeV and 't Hooft coupling $\lambda = 16.65$.

• Extraction of Quantities:

  • Once the soliton solution is obtained, the authors compute the energy density $\epsilon(r)$, pressure $p(r)$, and shear force $s(r)$ by integrating the 5D EMT components over the holographic coordinate $z$.
  • The D-term is calculated via two independent integrals: one using the shear force $s(r)$ (Eq. 2.14) and one using the pressure $p(r)$ (Eq. 2.17).

3. Key Contributions

1. First Full Numerical Solution: The paper provides the first rigorous numerical solution of the full equations of motion for the baryon soliton in the Sakai-Sugimoto model, replacing the previous "smooth connection" approximation.
2. Resolution of Intermediate Region Sensitivity: The authors demonstrate that the D-term is critically sensitive to the gauge field configuration in the intermediate region ($w \neq 0$), which was inadequately treated in previous work.
3. Decomposition of Contributions: The study explicitly separates and analyzes the contributions of the $SU(2)$ (non-Abelian) and $U(1)$ (Abelian) sectors to the total mass, pressure, and D-term.
4. Validation of Conservation Laws: The numerical solution satisfies the mechanical stability condition (von Laue condition) and the conservation of the EMT to high precision ($\sim 1\%$ error relative to central pressure).

4. Results

The numerical analysis yields the following key results (summarized in Table 1 of the paper):

• D-term Value:

  • The calculated D-term is $D \approx -2.05$ (specifically $D_s = -2.05$ and $D_p = -2.06$).
  • This value is significantly larger in magnitude (about 15 times larger) than the previous estimate of $-0.14$.
  • The result falls within the theoretically expected range for nucleons ($-5 < D < -1$) derived from experimental data and other models (Bag model, Lattice QCD, Skyrme model).

• Component Analysis:

  • $SU(2)$ Contribution: The dominant contribution comes from the $SU(2)$ sector, which is negative ($D_{SU(2)} \approx -2.54$). The magnitude is roughly 4 times larger than in the previous approximation. This is attributed to the deviation of the gauge field from the self-dual instanton solution in the intermediate $w$ region.
  • $U(1)$ Contribution: The $U(1)$ sector contributes positively ($D_{U(1)} \approx 0.49$), partially canceling the negative $SU(2)$ term. This contribution remained relatively stable compared to the previous work.

• Mass and Radii:

  • Total Mass: $M_{sol} \approx 1.18$ GeV ($M_{SU(2)} \approx 950$ MeV, $M_{U(1)} \approx 232$ MeV).
  • Mean Square Radii:
    • Energy density radius: $\langle r^2 \rangle_\epsilon \approx (0.662 \text{ fm})^2$.
    • Mechanical radius: $\langle r^2 \rangle_{mech} \approx (0.938 \text{ fm})^2$.
  • These radii are consistent with recent results from other models and experimental extractions.

• Pressure and Shear Force Profiles:

  • The pressure $p(r)$ is positive in the core ($r \lesssim 0.64$ fm) and negative in the outer region, satisfying the von Laue condition ($\int r^2 p(r) dr = 0$).
  • The shear force $s(r)$ is positive, consistent with the stability requirements of the system.

5. Significance

• Validation of Holographic QCD: The result brings the top-down holographic QCD prediction for the D-term into much closer alignment with experimental data and other theoretical approaches (like Lattice QCD and the Skyrme model), validating the model's capability to describe the mechanical structure of hadrons.
• Importance of Numerical Precision: The study highlights that approximations in the intermediate region of the holographic coordinate can lead to order-of-magnitude errors in derived quantities like the D-term. Accurate numerical solutions of the full EOMs are essential for reliable predictions in holographic models.
• Mechanical Structure: The paper provides a detailed map of the internal pressure and shear forces of the nucleon within the holographic framework, offering insights into the "mechanical radius" and the balance of forces that stabilize the baryon.
• Future Directions: The authors identify that future work should include the quantization of collective coordinates (to account for spin and the $J(t)$ form factor), the inclusion of quark masses (crucial for hyperons and IR behavior), and higher-order corrections ($\alpha'$ and $1/N_c$) to refine the EMT expression at higher momentum transfers.

In conclusion, this paper significantly advances the understanding of baryon structure in holographic QCD by correcting a major numerical approximation, resulting in a D-term value that is physically realistic and consistent with the broader landscape of hadron physics.