Factorizing quarkonium production matrix elements using effective field theory

This paper utilizes effective field theory and a Hubbard-Stratonovich transformation to factorize quarkonium production matrix elements in NRQCD into state-independent gluon correlators and wavefunctions at the origin, thereby verifying existing relationships for S-wave states, identifying new P-wave contributions, and restoring universality to TMD soft transition functions.

The Gist

The Big Picture: Building a Heavy Metal Ball

Imagine you are trying to understand how a heavy metal ball (called a quarkonium) is formed in a high-speed collision, like two cars crashing. Inside this ball are two very heavy particles (a heavy quark and an antiquark) that are stuck together.

For a long time, physicists used a rulebook called NRQCD to predict how often these balls are made. The rulebook said: "To make the ball, you need to know the probability of the two heavy particles sticking together." These probabilities are called Matrix Elements.

The problem was that the rulebook treated the "glue" holding the particles together as a messy, unseparated blob. It didn't distinguish between the "soft" glue (gentle tugs) and the "ultra-soft" glue (very gentle, long-range whispers). Because of this, the predictions were often vague, and the numbers needed to fit the data didn't always make sense.

The New Tool: The "Hubbard-Stratonovich" Magic Trick

This paper introduces a new way to look at the problem using a mathematical technique called a Hubbard-Stratonovich transformation.

The Analogy:
Imagine you have a group of people (the heavy quarks) trying to hold hands in a crowded room full of wind (gluons).

• Old Way: You tried to track every single person and every single gust of wind simultaneously. It was chaotic and impossible to separate the people from the wind.
• New Way: The authors introduce a "Ghost Team" (composite fields). Instead of tracking the people holding hands directly, they imagine a ghost team representing the finished pair.
• The Trick: They use a mathematical "magic trick" to swap the messy interaction of people + wind for a clean interaction between the Ghost Team and the Wind.

The Big Discovery: Untangling the Knots

The most important finding of this paper is that they proved you can untangle the "soft" wind from the "heavy people" specifically when the ball is being created.

1. The "Zero-Radius" Secret: When the heavy particles are first created in a collision, they are born at the exact same point in space (zero distance apart).
2. The Decoupling: Because they are born at the same point, the "soft wind" (which usually messes things up) cannot stick to the heavy particles in a way that prevents them from forming the final ball. The math shows that the "soft wind" and the "heavy particles" can be separated into two completely independent lists.
3. The Result: The probability of making the ball can now be written as two separate things multiplied together:
   • Part A: How big the ball is (the "wavefunction at the origin").
   • Part B: A universal "glue strength" factor (a vacuum correlator) that is the same for any type of heavy ball, regardless of which specific ball it is.

Why This Matters: The "Universal Glue"

Before this paper, physicists had to measure a different "glue strength" for every single type of heavy ball (J/ψ, ψ(2S), Υ, etc.). It was like needing a different key for every single lock in a house.

This paper proves that the locks are actually the same.

• If you know the "glue strength" for one type of heavy ball, you automatically know it for all the others.
• This reduces the number of unknown variables (free parameters) in the theory from 12 down to 3.
• It makes the theory much more powerful because it connects different experiments. If you measure one type of ball, you can predict the behavior of another type with high confidence.

A New Twist on "P-Waves"

The paper also looked at a specific type of formation called "P-wave" (where the particles have a bit of spin or rotation).

• They found a new type of contribution that was previously overlooked.
• Analogy: Imagine you thought a car engine only had a main piston. They found a smaller, secondary piston that kicks in under specific conditions.
• This new contribution might explain why some current experiments (like those at the LHC) don't quite match the old predictions at low speeds. It suggests the "secondary piston" might be more important than we thought.

The "TMD" Connection: Predicting the Future

Finally, the paper applies this logic to a framework called TMD (Transverse Momentum Dependent), which deals with particles moving sideways.

• In the past, the rules for sideways movement were messy and seemed to depend on the specific experiment (process-dependent).
• By using their new "untangling" method, they showed that even in these sideways scenarios, the "glue strength" is actually universal.
• This means we can now use data from one experiment to predict results in another, completely different experiment, which is a huge step forward for precision physics.

Summary

In short, this paper uses a clever mathematical trick to separate the "messy glue" from the "heavy particles" during the creation of a quarkonium ball. They discovered that the glue is actually universal across different types of balls. This simplifies the rules of the universe, reduces the number of unknowns, and helps physicists make much sharper predictions about how these heavy particles behave.

Technical Summary

Technical Summary: Factorizing Quarkonium Production Matrix Elements Using Effective Field Theory

Problem Statement
The production of heavy quarkonium (bound states of heavy quark-antiquark pairs, $Q\bar{Q}$) is traditionally described using Non-Relativistic QCD (NRQCD). In this framework, the cross-section is factorized into a perturbatively calculable short-distance coefficient and non-perturbative Long-Distance Matrix Elements (LDMEs). These LDMEs describe the probability of a $Q\bar{Q}$ pair in a specific quantum state transitioning into the final quarkonium state.

However, significant challenges remain in the NRQCD formalism:

1. Scale Separation: Traditional NRQCD treats all infrared physics below the hard scale ($2m$) as equivalent, failing to distinguish between "soft" ($mv$) and "ultrasoft" ($mv^2$) scales. This obscures the systematic calculation of corrections and the factorization of these scales.
2. Universality Issues: Recent developments in Transverse Momentum Dependent (TMD) factorization have introduced new production operators (TMD Shape Functions and TMD Soft Transition Functions). These operators encode soft radiation from the entire collision, leading to apparent process dependence that violates the universality of LDMEs.
3. Parameter Proliferation: For S-wave vector quarkonia, there are four independent LDMEs at leading and subleading orders. Global fits to data often yield inconsistent values for these parameters, particularly for color-octet states, suggesting missing theoretical constraints.

While Potential NRQCD (pNRQCD) has successfully derived factorization theorems that relate different LDMEs to universal vacuum correlators of chromo-electric and chromo-magnetic fields, it is unclear how to replicate these results within the Velocity NRQCD (vNRQCD) framework. vNRQCD explicitly retains soft degrees of freedom alongside ultrasoft ones, making it difficult to decouple the soft sector from the heavy quark fields in production matrix elements.

Methodology
The authors employ vNRQCD, an effective field theory that maintains explicit soft and ultrasoft gluon modes and heavy quark/antiquark fields. The core methodology involves:

1. Hubbard-Stratonovich Transformation: The authors apply a Hubbard-Stratonovich transformation to the vNRQCD Lagrangian. This technique introduces auxiliary bosonic fields—composite color-singlet ($S_r$) and color-octet ($O_r$) fields—representing the $Q\bar{Q}$ bound state. This transforms the four-fermion potential interactions into interactions between quarks and these composite fields.
2. Decoupling Scales:
   • Ultrasoft Decoupling: Using a BPS (Bauer-Pirjol-Stewart) field redefinition, the authors decouple the ultrasoft gluons from the heavy quark fields at leading order, separating the ultrasoft Hilbert space from the potential and soft sectors.
   • Soft Decoupling in Production: A critical insight is that quarkonium production operators are matched onto local operators (zero radial separation, $r=0$) due to the multipole expansion of the hard process. The authors demonstrate that for $r=0$, interactions between the composite fields and soft gluons (specifically "seagull" vertices involving two soft gluons) vanish. This allows the soft sector to be factorized from the composite fields in production matrix elements, a step that is generally impossible in bound-state decay calculations where $r \neq 0$.
3. Factorization of Matrix Elements: By integrating out the quark and antiquark fields, the production matrix elements are rewritten in terms of the composite fields evaluated at the origin. The soft gluon interactions are then isolated into vacuum correlators of gauge-invariant chromo-electric ($E$) and chromo-magnetic ($B$) fields.

Key Contributions and Results

• Derivation of Factorization Theorems in vNRQCD: The paper successfully derives factorization theorems for quarkonium production matrix elements within the vNRQCD framework, confirming results previously obtained only in pNRQCD.

  • S-wave Vector States: The color-singlet LDME $\langle O_V(^3S^{[1]}_1) \rangle$ is shown to be proportional to the wavefunction at the origin squared, $|R_V(0)|^2$.
  • Color-Octet S-wave States: The subleading color-octet LDMEs are factorized into the wavefunction at the origin and universal vacuum correlators:
    • $\langle O_V(^1S^{[8]}_0) \rangle \propto |R_V(0)|^2 \langle B \rangle$, where $\langle B \rangle$ is a chromo-magnetic field correlator.
    • $\langle O_V(^3S^{[8]}_1) \rangle \propto |R_V(0)|^2 \langle EE \rangle$, where $\langle EE \rangle$ is a double chromo-electric field correlator.
  • P-wave States: The authors derive new operator contributions for the color-octet P-wave mechanism. They identify a leading transition operator proportional to a single chromo-electric field correlator $\langle E \rangle$ and a subleading "genuine P-wave" operator proportional to the binding energy and a different correlator.

• Universality Relations: The factorization leads to powerful constraints relating the LDMEs of different S-wave vector quarkonium states (e.g., $J/\psi$ and $\Upsilon$). The ratios of LDMEs for different states are determined solely by the ratio of their wavefunctions at the origin and mass-dependent coefficients, reducing the number of independent non-perturbative parameters from 12 to 3 for S-wave vector states.

• TMD Factorization Application: The authors apply these techniques to TMD Soft Transition Functions (TMDSTFs), which are the dominant operators for quarkonium production in Semi-Inclusive Deep Inelastic Scattering (SIDIS) at low transverse momentum.

  • They demonstrate that TMDSTFs can be factorized into a process-dependent vacuum correlator (containing initial-state Wilson lines) and a state-independent wavefunction factor.
  • This establishes a "state-independence" property for TMDSTFs in SIDIS, implying that the TMDSTFs for different quarkonium states are related, thereby restoring a degree of universality and improving the predictive power of the TMD framework.

• Velocity Scaling Analysis: The paper re-evaluates the velocity scaling of the LDMEs. While traditional NRQCD power counting suggests all color-octet LDMEs scale as $v^7$, the authors find that the soft transitions derived here scale differently (e.g., $\langle O_V(^1S^{[8]}_0) \rangle \sim v^6$ and $\langle O_V(^3S^{[8]}_1) \rangle \sim v^7$). They argue that ultrasoft transitions, which would scale as $v^7$, are power-suppressed relative to the soft transitions in the production context due to the time-ordering of scales ($\tau_s \ll \tau_{us}$).

Significance
The paper claims to resolve the apparent tension between the vNRQCD and pNRQCD formalisms regarding the factorization of production matrix elements. By utilizing the Hubbard-Stratonovich transformation and exploiting the local nature of production operators ($r=0$), the authors show that the soft sector can be decoupled in vNRQCD, just as it is in pNRQCD.

The primary significance lies in the restoration of universality for quarkonium production parameters. By expressing LDMEs and TMDSTFs in terms of universal vacuum correlators and the wavefunction at the origin, the number of free parameters in theoretical predictions is drastically reduced. This provides a more robust theoretical foundation for extracting non-perturbative parameters from experimental data and improves the predictive power of NRQCD, particularly in the TMD framework where process dependence had previously been a major obstacle. The identification of new operator contributions for P-wave states also offers a potential explanation for discrepancies between NRQCD predictions and experimental data at low transverse momentum.