Multiple dispersive bounds. I) The z-expansion

This paper proposes enhancing the phenomenological application of the Boyd-Grinstein-Lebed $z$-expansion for hadronic form factors by explicitly incorporating a unitarity filter derived from its equivalence to the Dispersion Matrix method and by introducing kernel functions to apply multiple dispersive bounds, thereby enabling more rigorous, model-independent analyses of various physical processes.

The Gist

Imagine you are trying to predict the weather for a specific city next week. You have a bunch of data points: temperature, humidity, wind speed, and pressure readings from various sensors. You also have a set of "laws of physics" that say, for example, "it's impossible for the temperature to be -100°C in July" or "wind speed cannot be infinite."

In the world of particle physics, scientists face a similar challenge. They want to understand Hadronic Form Factors. Think of these as the "weather patterns" of the subatomic world. They describe how particles made of quarks and gluons (like protons or neutrons) interact with each other when hit by a force (like a weak nuclear force or an electromagnetic field).

The problem is that calculating these interactions from scratch is incredibly hard, like trying to predict a hurricane by tracking every single air molecule. So, scientists use a clever mathematical shortcut called the BGL z-expansion. It's like using a standard weather model that fits the data you have, while respecting the laws of physics.

This paper, written by Silvano Simula and Ludovico Vittorio, proposes two major upgrades to this "weather model" to make it more accurate and reliable.

1. The "Reality Check" Filter (The Unitarity Filter)

The Analogy:
Imagine you are filling out a survey about your daily habits. You say you sleep 2 hours, eat 50 meals, and run a marathon every day. Even if your answers are mathematically possible in a vacuum, they violate the basic laws of biology (you can't survive on 2 hours of sleep and 50 meals). A smart analyst would say, "Wait a minute, these numbers don't make sense together. Let's filter out the impossible combinations before we try to predict your future."

The Science:
In particle physics, there is a fundamental rule called Unitarity. It basically says: "Probability must add up to 100%." You can't create energy out of nothing, and you can't have a particle interaction that is more likely than 100%.

The authors point out that when scientists use the BGL model to fit experimental data, they usually check if the final model obeys Unitarity. But they often forget to check if the input data itself makes sense before fitting it.

• The Old Way: Take all the data, fit a curve, and hope the curve doesn't break physics.
• The New Way (The Filter): Before fitting the curve, run the data through a "Unitarity Filter." This filter checks if the data points are consistent with the laws of probability. If a data point (or a group of them) suggests something impossible (like a probability greater than 100%), the filter flags it or removes it.

Why it matters:
Sometimes, experimental data has small errors or "noise" that, when combined, look like they break the laws of physics. If you don't filter this out, your final prediction will be wrong. The authors show that this filter is crucial for getting accurate results, especially when trying to measure fundamental constants of the universe (like the strength of the weak force).

2. The "Multi-Lens" Strategy (Multiple Dispersive Bounds)

The Analogy:
Imagine you are trying to guess the shape of a hidden object in a dark room.

• The Old Way: You shine one giant flashlight from the ceiling. You get a general shadow, but you miss the details. You know the object is roughly "big," but you don't know if it's a ball, a cube, or a pyramid.
• The New Way: You use a set of specialized flashlights with different colored filters and angles. One light highlights the top, another the sides, and another the bottom. By combining these specific views, you get a much sharper, more detailed 3D picture of the object.

The Science:
The BGL model uses a "bound" (a limit) to keep the math from running wild. Traditionally, scientists use one single limit for the entire range of data. It's like saying, "The total energy of all possible outcomes must be less than 100."

The authors propose using Multiple Dispersive Bounds. Instead of one big limit, they break the problem down into smaller pieces using "kernel functions" (think of these as the colored filters on the flashlights).

• They divide the "energy budget" into different categories (e.g., short-distance effects, long-distance effects, or specific mass ranges).
• They apply a limit to each category individually.

Why it matters:
By constraining the model in multiple specific ways, rather than just one general way, the scientists can narrow down the possible answers much more tightly. It's like solving a puzzle where you have more clues. This leads to much more precise predictions about how particles behave, which is essential for testing if our current theories of physics are correct or if there is "New Physics" hiding in the details.

The Big Picture

The authors are essentially saying: "We have a great tool (the BGL expansion) for understanding how particles interact. But to get the best results, we need to be smarter about how we use it."

1. Check your ingredients first: Don't just mix the data; make sure the data itself isn't broken by the laws of physics (The Filter).
2. Look at the details: Don't just look at the whole picture; break it down and apply rules to specific parts to get a clearer image (Multiple Bounds).

By doing this, they can help physicists extract more accurate information from experiments, which is vital for understanding the fundamental building blocks of our universe, from the decay of heavy particles to the behavior of light.

Technical Summary

1. Problem Statement

Hadronic form factors are essential for understanding strong non-perturbative dynamics and are crucial for precision tests of the Standard Model, such as the extraction of the CKM matrix element $|V_{cb}|$ and the calculation of the ratio $R(D^*)$ in semileptonic $B \to D^* \ell \nu$ decays.

The Boyd-Grinstein-Lebed (BGL) $z$-expansion is the standard model-independent parametrization used to describe these form factors, relying on unitarity and analyticity. However, the authors identify two critical limitations in current phenomenological applications:

1. Neglect of Input Data Consistency: Standard BGL fits impose unitarity constraints on the expansion coefficients but often fail to check if the input data itself is consistent with unitarity. Non-unitary data (due to statistical fluctuations, systematic errors, or unphysical extrapolations) can lead to biased results.
2. Single Global Bound: Current methods typically apply a single, global dispersive bound (a total integral constraint) to the form factor. This "averages out" detailed information about the form factor's behavior along the unitarity branch cut, potentially leading to looser constraints than necessary.

2. Methodology

The authors propose a refined framework integrating two main ingredients into the BGL $z$-expansion:

A. The Unitarity Filter (Input Data Consistency)

The authors establish an equivalence between the BGL approach and the Dispersion Matrix (DM) method. They derive an additive representation of the form factor $f(z)$:
$$f(z) = \beta(z) + R(z)$$

• $\beta(z)$: An analytic function constructed solely from the input data $\{f_i\}$ at specific points $z_i$. It has no free parameters.
• $R(z)$: A remainder function representing the uncertainty and higher-order terms.

They derive a specific inequality, $\chi[\beta] \leq \chi_U$, where $\chi[\beta]$ is the susceptibility calculated directly from the input data, and $\chi_U$ is the theoretical upper bound derived from QCD (via polarization functions).

• Function: This acts as a filter. If a dataset yields $\chi[\beta] > \chi_U$, the data is inconsistent with unitarity and must be rejected or re-evaluated before fitting.
• Implementation: The authors propose a procedure where input data is sampled (e.g., via Gaussian distributions), filtered against this inequality, and only the "unitary subset" is used for the final BGL fit.

B. Multiple Dispersive Bounds

Instead of a single global bound $\chi[f] \leq \chi_U$, the authors introduce a set of kernel functions $\tilde{K}_p(z)$ (where $p=1, \dots, P$) that partition the unitarity circle.

• Decomposition: The total susceptibility is split into partial contributions: $\chi[f] = \sum \chi_p[f]$.
• Constraints: Each partial contribution is bounded individually: $\chi_p[f] \leq \chi_U^p$.
• Implementation: These kernels can be defined in the Euclidean time domain (using lattice QCD correlators $V(\tau)$) or as step functions in the energy domain. By applying multiple simultaneous constraints, the allowed region for the form factor coefficients is significantly reduced compared to using a single global bound.

3. Key Contributions

1. Formal Equivalence and Filter Derivation:
   The paper rigorously demonstrates that the BGL expansion with unitarity constraints is mathematically equivalent to the DM method. Crucially, it derives the $\chi[\beta] \leq \chi_U$ condition, proving that input data must satisfy this filter to be physically valid within the unitary framework. This distinguishes between the unitarity of the expansion (coefficients) and the unitarity of the data.

2. Algorithm for Data Filtering:
   The authors provide a concrete algorithm for phenomenological analyses:

   • Generate input data samples.
   • Calculate $\chi[\beta]$ for each sample.
   • Discard samples where $\chi[\beta] > \chi_U$.
   • Quantify the impact of this filtering using metrics $\Delta$ (shift in mean values) and $\epsilon$ (change in uncertainties). They suggest acceptable thresholds ($\Delta \lesssim 1$, $\epsilon \gtrsim 0.5$).

3. Generalization to Multiple Bounds:
   The paper introduces the formalism for multiple dispersive bounds. By introducing kernel functions (e.g., short, intermediate, and long-distance windows in Euclidean time), the authors show how to impose multiple independent constraints on the form factor. This allows for a more granular control over the form factor's behavior, particularly useful for analyzing sub-threshold branch cuts (addressed in detail in the companion paper [6]).

4. Lattice QCD Integration:
   The methodology explicitly connects to Lattice QCD (LQCD) by utilizing Euclidean current-current correlators $V(\tau)$ to compute the susceptibilities $\chi_U$ and the multiple bounds $\chi_U^p$ non-perturbatively.

4. Results and Theoretical Findings

• Filtering Necessity: The authors argue that many existing datasets (e.g., from Lattice QCD collaborations or experimental measurements) may contain non-unitary fluctuations. Applying the filter $\chi[\beta] \leq \chi_U$ is essential to avoid bias. In previous studies (e.g., $B \to D^*$), applying this filter reduced the subset of usable data but significantly tightened the final predictions.
• Impact on Uncertainties: The inclusion of the unitarity filter often leads to a reduction in the uncertainties associated with form factor extrapolations, especially when the number of input data points is large.
• Multiple Bounds Advantage: While the full numerical application is reserved for the companion paper, the formalism shows that multiple bounds are strictly more constraining than a single global bound. They prevent the form factor from satisfying the global bound while violating local physical constraints in specific kinematic regions.
• Robustness: The authors show that the results are robust against the choice of the auxiliary variable $t_0$ in the $z$-mapping, provided the truncation error is under control.

5. Significance

This work represents a significant advancement in the precision analysis of hadronic form factors:

• Rigorous Data Validation: It provides a necessary "sanity check" for input data, ensuring that phenomenological fits do not inadvertently fit unphysical data points that violate fundamental QCD principles.
• Enhanced Precision: By utilizing multiple dispersive bounds, the method extracts more information from the same dataset, leading to tighter constraints on form factors. This is critical for reducing theoretical uncertainties in the determination of $|V_{cb}|$ and $|V_{ub}|$, which are currently sources of tension in flavor physics (the "$V_{cb}$ puzzle").
• Bridge between Lattice and Phenomenology: The framework seamlessly integrates non-perturbative LQCD calculations (via susceptibilities) with experimental data analysis, offering a unified approach to handle both theoretical and experimental inputs consistently.

In summary, Simula and Vittorio propose a more rigorous, two-step unitary analysis: first, filter the input data to ensure physical consistency, and second, apply multiple dispersive bounds to maximize the constraining power of the unitarity principle on the form factor parametrization.