Glauber-Lachs formula-based analysis of three-pion Bose-Einstein correlation data at 7 TeV from the LHCb Collaboration

This paper analyzes 7 TeV LHCb Collaboration data on two- and three-pion Bose-Einstein correlations by combining the Glauber-Lachs formula with a two-component pion production model, utilizing specific pion exchange functions extended into four-dimensional Euclidean configuration space.

The Gist

Imagine you are at a massive, chaotic party (the Large Hadron Collider, or LHC) where billions of particles are smashing into each other. When these particles collide, they don't just disappear; they burst out in a shower of new particles, mostly pions (which are like the "confetti" of the particle world).

This paper is about trying to understand the dance of this confetti. Specifically, the authors are looking at how groups of pions (two at a time, and three at a time) behave when they are identical twins.

Here is the breakdown of their work using simple analogies:

1. The Core Mystery: The "Bose-Einstein" Dance

In the quantum world, identical particles (like two pions) have a weird rule: they love to stick together. If you throw two identical pions out, they tend to land closer to each other than you would expect by pure chance. This is called Bose-Einstein Correlation (BEC).

Think of it like this: If you throw two identical red balloons into a crowd, they might drift toward each other not because they are holding hands, but because the "rules of the universe" make them want to be close. The scientists measure how close they get to figure out the size of the room (the collision area) where they were born.

2. The Old Map vs. The New Map

The LHCb Collaboration (the experimental team) had already looked at this data using a specific mathematical map called the Glauber-Lachs formula.

• The Old Map: They assumed the "room" where pions are born was a single, uniform blob. They also had to add a "correction factor" (like a patch on a map) to make the math fit the data.
• The New Map (This Paper): The authors (Mizoguchi, Matsumoto, and Biyajima) decided to look at the room differently. They used a Two-Component Picture.

The Analogy:
Imagine the collision isn't just one big messy room. Instead, imagine it's a large, fuzzy cloud (the chaotic part) with tiny, sharp laser spots inside it (the coherent part).

• Chaotic Part (The Cloud): This is the messy, random explosion of particles. It's big and spread out.
• Coherent Part (The Laser Spots): This is a more organized, "in-sync" emission of particles. It's small and tight.

The authors realized that the old map tried to describe this whole scene with just one shape. Their new map says, "No, we need two shapes: a big dipole (a smooth, round shape) for the cloud, and a sharper, different shape for the laser spots."

3. The "One-and-a-Half Pole" Secret Sauce

This is the most technical part, but here's the simple version:
In math, shapes are often described by how "steep" their sides are.

• The "Cloud" (chaotic part) has a standard steepness (a "dipole").
• The "Laser Spots" (coherent part) needed a different steepness to fit the data perfectly.

The authors invented a new mathematical shape they call a "one-and-a-half pole."

• Analogy: Imagine a hill. A standard hill has a slope of 1. A very steep cliff has a slope of 2. The authors found that the "laser spots" in the particle collision act like a hill with a slope of 1.5. It's not quite a cliff, but it's steeper than a normal hill. This specific shape allowed their math to fit the experimental data much better than the old methods.

4. What Did They Find?

By using this new "Two-Component" map with the "1.5 slope" shape, they analyzed data from the LHCb experiment at 7 TeV (a very high energy).

• The Size of the Room: They calculated the size of the area where pions are born. Their results matched the LHCb team's previous findings, which is a good sign.
• The Difference: They found that the "cloud" (chaotic part) is quite large (about 1.5 to 1.8 femtometers—imagine a femtometer is a trillionth of a millimeter). The "laser spots" (coherent part) are tiny (about 0.25 femtometers).
• The "Delta R": They noticed that when looking at three pions instead of two, the "cloud" seems to get slightly bigger. This suggests that the way particles are produced changes slightly depending on how many you look at.

5. Why Does This Matter?

• Better Physics: By understanding that the "room" has two different textures (a big fuzzy cloud and tiny sharp spots), we get a clearer picture of how the universe works at the smallest scales.
• Quantum Optics in Particle Physics: The authors borrowed a formula from Quantum Optics (the study of light and lasers) and applied it to Hadron Physics (the study of heavy particles like protons). It's like using a recipe for baking a cake to successfully cook a steak. It turns out the math for light waves works surprisingly well for particle collisions.
• Future Predictions: They even used their new formula to predict what would happen if we looked at four pions at once. They created a graph showing what the data should look like, waiting for the LHCb team to test it in the future.

Summary

The authors took a complex puzzle about how particles stick together. They realized the old way of looking at it was too simple. By splitting the problem into a "big messy part" and a "small organized part," and using a clever new mathematical shape (the 1.5 slope), they created a much clearer picture of the particle collision "dance floor." This helps us understand the fundamental structure of matter better.

Technical Summary

1. Problem Statement

The paper addresses the analysis of Bose–Einstein Correlation (BEC) data for two-pion ($2\pi$) and three-pion ($3\pi$) systems at a center-of-mass energy of 7 TeV, provided by the LHCb Collaboration.

• Context: Previous analyses by LHCb utilized the Glauber–Lachs (GL) formula with a correction factor ($f_c$) to account for the mixture of chaotic and coherent pion production. However, the standard exponential forms used in these analyses often required ad-hoc correction factors and struggled to provide a consistent geometric interpretation of the source size in configuration space.
• Objective: The authors aim to re-analyze this data using a refined theoretical framework that combines the GL formula from quantum optics with a two-component picture (chaotic vs. coherent) and specific geometric forms (dipole and inverse pole) for the pion exchange function. The goal is to improve the fit quality ($\chi^2$), eliminate the need for arbitrary correction factors, and derive consistent source extension parameters ($R$) in 4-dimensional Euclidean space.

2. Methodology

The authors propose a new theoretical formulation for the BEC correlation functions, moving away from simple exponential decays to forms inspired by particle form factors and quantum optics.

A. Theoretical Framework

• Two-Component Picture: The average pion multiplicity $\langle n \rangle$ is decomposed into a chaotic part ($A$) and a coherent part ($|\zeta|^2$), defined by the parameter $p = A/\langle n \rangle$ (chaotic fraction) and $1-p$ (coherent fraction).
• Exchange Functions ($E_{2B}$):
  • Chaotic Term ($p^2$): Modeled using a dipole form: $1/(1 + (R_1 Q)^2)^2$. This is analogous to the proton electric form factor.
  • Mixed Term ($2p(1-p)$): Modeled using a weaker inverse one-and-a-half pole form: $1/(1 + (R_2 Q)^2)^{1.5}$. This accounts for the interference between chaotic and coherent components.
  • Coherent Term: Modeled with a half-pole form.
• Long-Range Correlation (LRC): The authors test two LRC models:
  1. Gaussian: $C/(1 + \alpha e^{-\beta Q^2})$ (decaying function).
  2. Linear: $C(1 + \delta Q)$ (increasing function).
• Configuration Space: The analysis is performed in 4-dimensional Euclidean space ($\xi = \sqrt{|r_1-r_2|^2 + (t_1-t_2)^2}$), utilizing Wick rotation to relate momentum space ($Q$) to configuration space ($\xi$).

B. Data Analysis

• Data Source: LHCb 7 TeV data, categorized by charged particle multiplicity ($N_{ch}$) into low, medium, and high activity bins.
• Fitting Procedure: The authors fit the proposed formulas (Eqs. 5 and 6) to the experimental ratios $N(2\pi)/N_{BG}$ and $N(3\pi)/N_{BG}$.
• Comparative Analysis:
  • Comparison with the original LHCb analysis (Table 1 vs. Table 2).
  • Comparison between the "dipole + 1.5-pole" model (Table 3) and a "dipole + Gaussian" alternative model (Table 7).
  • Comparison between Gaussian and Linear LRC models (Table 3 vs. Table 5).

3. Key Contributions

1. Refined Glauber–Lachs Formulation: The paper introduces a specific geometric ansatz where the chaotic component follows a dipole distribution and the mixed chaotic-coherent component follows an inverse 1.5-pole distribution. This removes the need for the empirical correction factor $f_c$ used in previous LHCb analyses.
2. Two-Component Geometric Picture: The authors successfully interpret the source as having two distinct scales:
   • A large-range component ($R_1$) associated with the chaotic part (Non-Single Diffractive dissociation).
   • A small-range component ($R_2$) associated with the mixed/interference part (Single Diffractive dissociation or localized coherent spots).
3. Configuration Space Extensions: The authors calculate the root-mean-square extensions ($\sqrt{\langle \xi^2 \rangle}$) in 4D space. They find that the large component corresponds to resonance decay lengths (e.g., $\rho$ and $\omega$ mesons), while the small component represents localized coherent emission.
4. Prediction for Four-Pion BEC: Using the derived parameters, the authors provide a theoretical prediction for four-pion Bose–Einstein correlations, extending the applicability of quantum optics to hadron physics.

4. Results

• Fit Quality ($\chi^2$):
  • The new "dipole + 1.5-pole" model (Table 3) yields significantly better $\chi^2$ values compared to the original LHCb analysis (Table 1) and the single-component dipole model (Table 2).
  • For example, in the medium activity $2\pi$ case, $\chi^2/dof$ improved from 623/386 (LHCb) to 441/384 (New Model).
  • The $p$-values for the new model are generally high (e.g., 24%, 2.3%, 25% for $2\pi$), indicating a statistically acceptable fit, whereas the original LHCb fits often had very low $p$-values.
• Source Parameters:
  • $R_1$ (Chaotic): Ranges from $\approx 1.17$ fm (low activity) to $1.55$ fm (high activity).
  • $R_2$ (Mixed): Consistently small, ranging from $\approx 0.24$ fm to $0.29$ fm.
  • Extension ($\sqrt{\langle \xi^2 \rangle}$): The large component extends to $\approx 4.7 - 6.2$ fm, while the small component is $\approx 0.8 - 1.0$ fm.
• LRC Model Comparison:
  • The Gaussian LRC model provides a better description of the data than the Linear LRC model. The linear model results in increasing correlation functions at high $Q$, which is physically less plausible, and yields higher $\chi^2$ values.
• Alternative Models:
  • Replacing the 1.5-pole with a Gaussian distribution for the mixed term (Appendix B) results in slightly worse $\chi^2$ values for most bins, suggesting the 1.5-pole form is the more appropriate description for the mixed term in this context.

5. Significance

• Physical Interpretation: The study provides a robust physical interpretation of BEC parameters, linking the chaotic component to the decay of long-lived resonances (like $\rho$ and $\omega$ mesons) and the coherent/mixed component to localized, short-range emission spots.
• Theoretical Unification: It successfully bridges quantum optics (Glauber–Lachs formula) and high-energy hadron physics, demonstrating that the same mathematical structures used for photon statistics apply to pion production in high-energy collisions.
• Methodological Advancement: By eliminating the need for the correction factor $f_c$ and using a two-component geometric picture, the paper offers a more self-consistent framework for analyzing BEC data.
• Future Directions: The prediction of four-pion BEC correlations sets the stage for future experimental verification by the LHCb collaboration, potentially further validating the quantum optics approach in hadronic systems.

In conclusion, the paper demonstrates that a Glauber–Lachs analysis incorporating a dipole form for the chaotic part and an inverse 1.5-pole for the mixed part provides a superior fit to LHCb 7 TeV data, offering deeper insights into the space-time structure of pion production sources.