BHLS$_2$, a New Breaking of the HLS Model and its Phenomenology

TL;DR

This paper introduces BHLS$_2$, an improved Hidden Local Symmetry model that better describes light meson processes and reduces uncertainties in the theoretical calculation of the muon anomalous magnetic moment.

Contribution

BHLS$_2$ extends the BHLS framework with new symmetry breaking schemes, enhancing phenomenological accuracy and predictive power for light meson processes and muon g-2 calculations.

Findings

BHLS$_2$ accurately predicts the lattice QCD pion form factor.

It significantly reduces systematic uncertainties in muon g-2 theoretical estimates.

The leading order HVP contribution is precisely evaluated as 686.65 ± 3.01 (stat) + syst.

Abstract

Previous studies have shown that the Hidden Local Symmetry (HLS) Model, supplied with appropriate symmetry breaking, provides an Effective Lagrangian (BHLS) able to encompass a large number of processes within a unified framework. A global fit provides a simultaneous description of the $e^+ e^-$ annihilation into $\pi^+\pi^-$, $\pi^0\gamma$, $\eta \gamma$, $\pi^+\pi^-\pi^0$, $K^+K^-$, $K_L K_S$, the dipion spectrum in the $\tau$ decay and other light meson partial widths. Additional breaking schemes extend its scope to spacelike processes within a new framework (BHLS$_2$). The phenomenology previously explored with BHLS is revisited in the BHLS$_2$ context with emphasis on the $\Phi$ mass region. It is shown that BHLS$_2$ addresses the close spacelike region covered by NA7 and Fermilab data; also,the Lattice QCD pion form factor is accurately predicted by BHLS$_2$ using fits to annihilation data only. These channels contribution to $a_\mu$ over the range of validity of BHLS$_2$ is shown to strongly reduce the former BHLS systematics. The uncertainty on $a_\mu^{\rm th}(\sqrt{s}< 1.05$ GeV) is much improved compared to standard approaches relying on direct integration methods of measured spectra. Using the BHLS$_2$ results, the leading order HVP contribution to the muon anomalous moment is $a_\mu^{HVP-LO}= 686.65 \pm 3.01 +(+1.16,-0.75)_{syst.}$ in units of $10^{-10}$. Using a conservative estimate for the light--by--light contribution, our evaluation is $a_\mu^{\rm th}=\left [11\,659\,175.96 \pm 4.17 +(+1.16,-0.75)_{syst.}\right]\times 10^{-10}$. The $\rho^0\gamma$ mixing in dispersive and LQCD approaches is discussed and may amount to a shift $\delta a_\mu[\pi\pi]_{\rho\gamma} = +(3.10\pm0.31)\times 10^{-10}$ at LO+NLO, presently treated as additional systematics. Taking this shift into account, $a_\mu^{\rm th}-a_\mu^{\rm BNL}$ exhibits a significance not smaller than $3.8 \sigma$.

Tables

Table 1. Table 1: BHLS 2 fit properties in the kaon sectors. The kaon data from NSK, BaBar and CMD-3 are fitted in isolation. Running BaBar data ”alone” is complemented by the NSK data for the K L ​ K S subscript 𝐾 𝐿 subscript 𝐾 𝑆 K_{L}K_{S} channel. The data and channels other than kaon are identical in all cases; the full χ 2 superscript 𝜒 2 \chi^{2} value and the total number of data points ( N pts subscript 𝑁 pts N_{\rm pts} ) refer to the global fit. The last line displays the shift values δ ​ E CMD3 𝛿 subscript 𝐸 CMD3 \delta E_{\rm CMD3} or δ ​ E BaBar 𝛿 subscript 𝐸 BaBar \delta E_{\rm BaBar} relative to the NSK data energy calibration.

${\chi }^2$ ($K\overline{K}$ Sample Set)	NSK	BaBar	CMD3
NSK $K_L{K}_S$ (92)	$107$	$100$	$-$
NSK $K^{+}{K}^{-}$(49)	$51$	$-$	$-$
BaBar $K^{+}{K}^{-}$ (27)	$-$	$46$	$-$
CMD3 $K_L{K}_S$ (17)	$-$	$-$	$120$
CMD3 $K^{+}{K}^{-}$(16)	$-$	$-$	$88$
${\chi }^2/N_{\mathrm{pts}}$	$1087/1210$	$1078/1188$	1155/1102
Probability	97.5%	95.2 %	3.8 %
$\delta E$ (keV)	$-$	$-125.8\pm 19.0$	$-51.4\pm 13.7$

Table 2. Table 2: BHLS 2 fit properties in the kaon sectors. The kaon data from NSK, BaBar and CMD-3 are fitted pairwise combined. The χ 2 superscript 𝜒 2 \chi^{2} values are displayed for individual channels or for the total number of data points N pts subscript 𝑁 pts N_{\rm pts} . In the two rightmost data columns, the top entry is δ ​ E CMD3 𝛿 subscript 𝐸 CMD3 \delta E_{\rm CMD3} and the bottom one δ ​ E BaBar 𝛿 subscript 𝐸 BaBar \delta E_{\rm BaBar} . The rightmost column displays the fit results when treating the CMD-3 systematics as uncorrelated.

$K\overline{K}$ Sample Set	NSK+BaBar	NSK+CMD3	BaBar+CMD3	NSK +BaBar
(# data points)		(CMD3 corr.)	(CMD3 corr.)	+ CMD3 uncorr.
NSK $K_L{K}_S$ (92)	$106$	$105$	$\left(108\right)$	$113$
NSK $K^{+}{K}^{-}$(49)	$50$	$52$	$\left(50\right)$	$50$
BaBar $K^{+}{K}^{-}$ (27)	$46$	$-$	$45$	$56$
CMD3 $K_L{K}_S$ (17)	$-$	$121$	$122$	$18$
CMD3 $K^{+}{K}^{-}$(16)	$-$	$98$	$97$	$23$
${\chi }^2/N_{\mathrm{pts}}$	1133/1237	$1331/1243$	$1212/1129$	$1193/1270$
Probability	93.3%	1.1 %	0.8 %	81.5 %
$\delta E_{\mathrm{CMD3}}$ (keV)	$-$	$-46.5\pm 6.7$	$-45.9\pm 3.9$	$-38.2\pm 11.7$
$\delta E_{\mathrm{BaBar}}$ (keV)	$-105.8\pm 17.2$	$-$	$-105.5\pm 16.1$	$-103.2\pm 16.7$

Table 3. Table 3: Global fit properties of BHLS and of the BS and RS variants of BHLS 2  : χ 2 superscript 𝜒 2 \chi^{2} values yielded by each group of data samples for the indicated channel. Numbers within parentheses in the first column indicate the numbers of data points submitted to fits; for BHLS, the actual number of fitted data points in the 3 π 𝜋 \pi channel (model B) is shown in the entry as χ 2 / N pts superscript 𝜒 2 subscript 𝑁 pts \chi^{2}/N_{\rm pts} .

Channels	BHLS [25]	BHLS2 (BS)	BHLS2 (BS)	BHLS2 (RS)
		space, NO $\tau$	space, $\tau$	space, $\tau$
NSK ${\pi }^{+}{\pi }^{-}$ (127)	$137$	$140$	$139$	$137$
KLOE ${\pi }^{+}{\pi }^{-}$ (135)	$140$	$138$	$154$	$137$
$\tau$ (ABC) (85)	$85$	$-$	$98$	$92$
${\pi }^0\gamma$ (112)	$91$	$88$	$89$	$86$
$\eta \gamma$( 182)	$120$	$121$	$122$	$120$
${\pi }^{+}{\pi }^{-}{\pi }^0$ (158)	$79/96$	$148$	$153$	$141$
NSK $K_L{K}_S$ (92)	$97$	$107$	$106$	$106$
NSK $K^{+}{K}^{-}$(49)	$51$	$49$	$48$	$50$
BaBar $K^{+}{K}^{-}$ (27)	$49$	$45$	$44$	$46$
space $K^{+}{K}^{-}$ (25)	$-$	$19$	$19$	$18$
space ${\pi }^{+}{\pi }^{-}$(59)	$-$	$58$	$63$	$55$
${\chi }^2/N_{\mathrm{pts}}$	$949/1056$	$1054/1152$	1185/1237	1133/1237
Probability	96.7%	93.9 %	67.6%	93.2%

Table 4. Table 4: BS variant global fits excluding τ 𝜏 \tau data: The leftmost data column shows the partial χ 2 superscript 𝜒 2 \chi^{2} ’s in the fit including both kinds of spacelike data, the rightmost one shows the corresponding results in the fit excluding all spacelike data. Global fit information is displayed in the last lines. The numbers of fitted data points in each channel are given within parentheses.

BHLS2 : BS fit	incl. spacelike data	excl. spacelike data **
space ${\pi }^{+}{\pi }^{-}$ (59)	$57.6$	$58.6$
space $K^{+}{K}^{-}$ (25)	$18.8$	$19.0$
${\chi }^2/N_{\mathrm{dof}}$	$1054/1152$	$977/1068$
Probability	93.9%	92.8 %

Table 5. Table 5: Pion and kaon charged radii in units of fm 2 . See also [ 118 ] for additional LQCD evaluations of < r π ± 2 > expectation superscript subscript 𝑟 superscript 𝜋 plus-or-minus 2 <r_{\pi^{\pm}}^{2}> .

BHLS2 fit	$0.430\pm {0.002}_{\mathrm{mod}}\pm {0.001}_{\mathrm{fit}}$	$0.268\pm {0.004}_{\mathrm{mod}}\pm {0.001}_{\mathrm{fit}}$
NA7 [47, 48]	$0.439\pm 0.008$	$0.40\pm 0.11$
Fermilab [93, 94]	$0.439\pm 0.030$	$0.28\pm 0.05$
CHS [114]	$0.429\pm {0.001}_{\mathrm{stat}}\pm {0.004}_{\mathrm{syst}}$	–
ACD [115]	$0.432\pm 0.004$	–
ChPT 2 flavors [116]	$0.437\pm 0.016$	–
ChPT 3 flavors [117]	$0.452\pm 0.013$	$0.363\pm 0.072$
HPQCD [111]	$0.403\pm {0.018}_{\mathrm{stat}}\pm {0.006}_{\mathrm{syst}}$	–
ETM [113]	$0.443\pm {0.021}_{\mathrm{stat}}\pm {0.020}_{\mathrm{syst}}$	–

Table 6. Table 6: Model parameter values in BHLS and in the BS and RS variants of BHLS 2  ; second line in the Table title indicates the running conditions for what concerns the data samples submitted to fit. The last lines display the χ 2 superscript 𝜒 2 \chi^{2} , N pts subscript 𝑁 pts N_{\rm pts} and the probability of the fits.

Channels	BHLS [25]	BHLS2 (BS)	BHLS2 (RS)	BHLS2 (RS)
	(model B)	(no $\tau$)	($z_V$ free)	($z_V\equiv 1$)
$a_{\mathrm{HLS}}$	$2.620\pm 0.001$	$1.970\pm 0.002$	$1.590\pm 0.001$	$2.658\pm 0.002$
$g$	$5.570\pm 0.001$	$6.855\pm 0.003$	$6.707\pm 0.002$	$6.702\pm 0.003$
$\left(c_3+c_4\right)/2$	$0.886\pm 0.003$	$0.762\pm 0.003$	$0.756\pm 0.003$	$0.755\pm 0.003$
$c_1-c_2$	$1.130\pm 0.022$	$0.771\pm 0.020$	$0.833\pm 0.021$	$0.789\pm 0.018$
${10}^2\times z_3$	$\times$	$-0.373\pm 0.004$	$-0.424\pm 0.003$	$-0.428\pm 0.003$
$z_A$	$1.57\pm 0.01$	$1.55\pm 0.01$	$1.55\pm .01$	$1.48\pm .01$
$z_V$	$1.204\pm 0.001$	$1.403\pm 0.002$	$1.674\pm 0.002$	$1.00$
${10}^2\times {\Delta }_A$	$1.955\pm 0.517$	$0.11\pm 0.51$	$0.83\pm 0.51$	$0.82\pm 0.50$
${10}^2\times {\Sigma }_V$	$-15.63\pm 0.06$	$-7.69\pm 14.8$	$10.42\pm 0.15$	$-34.53\pm 0.09$
${10}^2\times {\xi }_0$	$\times$	$-6.40\pm 0.04$	$-3.60\pm 0.03$	$-1.85\pm 0.04$
${10}^2\times {\xi }_3$	$\times$	$2.34\pm 0.16$	$2.92\pm 0.15$	$2.02\pm 0.080$
${10}^2\times {\psi }_0$	$\times$	$\times$	$-3.91\pm 0.27$	$-4.24\pm 0.24$
${10}^2\times {\psi }_{\varphi }$	$\times$	$\times$	$-0.80\pm 0.09$	$-0.63\pm 0.09$
${10}^2\times {\psi }_{\omega }$	$\times$	$\times$	$-3.20\pm 0.11$	$-3.39\pm 0.12$
${\chi }^2/N_{\mathrm{pts}}$	$949.1/1056$	$1054/1152$	1133/1237	1137/1237
Probability	96.7%	93.9 %	93.3%	92.4%

Table 7. Table 7: HLS contributions to 10 10 × a HVP − LO superscript 10 10 superscript 𝑎 HVP LO 10^{10}\times a^{\rm HVP-LO} integrated up to 1.05 GeV, including FSR. The first data column displays the results using the former BHLS [ 25 , 27 ] and, the second one, those derived from the Basic Solution for BHLS 2 , the τ 𝜏 \tau decay data being discarded. The next two data columns refer to the results obtained using the Reference Solution for BHLS 2 using the largest set of data samples, keeping or discarding the τ 𝜏 \tau data. The last data column refers to the numerical integration for each channel of the same set of data which are used in the BHLS/BHLS 2 fits.

Channel	BHLS	BHLS2 (BS)	BHLS2 (RS)	BHLS2 (RS)	Data Direct
		(excl. $\tau$)	(incl. $\tau$)	(excl. $\tau$)	Integration
${\pi }^{+}{\pi }^{-}$	$493.73\pm 0.70$	$494.04\pm 0.76$	$494.08\pm 0.90$	$494.10\pm 1.0$	$497.82\pm 2.80$
${\pi }^0\gamma$	$4.42\pm 0.03$	$4.46\pm 0.03$	$4.45\pm 0.03$	$4.45\pm 0.03$	$3.47\pm 0.11$
$\eta \gamma$	$0.63\pm 0.01$	$0.63\pm 0.01$	$0.64\pm 0.01$	$0.64\pm 0.01$	$0.55\pm 0.02$
${\pi }^{+}{\pi }^{-}{\pi }^0$	$42.56\pm 0.54$	$42.97\pm 0.53$	$43.02\pm 0.54$	$43.08\pm 0.48$	$41.38\pm 1.28$
$K^{+}{K}^{-}$	$18.10\pm 0.14$	$18.04\pm 0.14$	$18.01\pm 0.14$	$17.99\pm 0.14$	$17.37\pm 0.55$
$K_L{K}_S$	$11.53\pm 0.08$	$11.70\pm 0.08$	$11.71\pm 0.08$	$11.72\pm 0.09$	$11.98\pm 0.36$
HLS Sum	$570.97\pm 0.92$	$571.84\pm 1.06$	$571.90\pm 1.10$	$571.98\pm 1.20$	$572.57\pm 3.15$
${\chi }^2/N_{\mathrm{pts}}$	$949/1056$	$1054/1152$	1133/1237	1052/1152	$\times$
Probability	96.7%	93.9 %	93.3 %	93.6 %	$\times$

Table 8. Table 8: LO- and NLO-HVP results in the standard dispersive approach based on e + ​ e − superscript 𝑒 superscript 𝑒 e^{+}e^{-} data. The NLO diagrams are shown in Fig.  12 .

		$a_{\mu }\times {10}^{-10}$	stat	syst	tot
LO		690.93	0.71	3.83	[3.90]
NLO	diagram a	-20.71	0.03	0.13	[0.13]
NLO	b	10.39	0.01	0.06	[0.06]
NLO	c	0.34	0.00	0.01	[0.01]
NLO	sum	-9.98	0.03	0.07	[0.08]
LO+NLO		680.95	-	-	[3.82]

Table 9. Table 9: LO- and NLO-HVP results in the dispersive approach based on e + ​ e − superscript 𝑒 superscript 𝑒 e^{+}e^{-} data, after having removed the ρ 0 − γ superscript 𝜌 0 𝛾 \rho^{0}-\gamma mixing from the hadronic blob.

		$a_{\mu }\times {10}^{-10}$	stat	syst	tot	shift
LO		694.11	0.71	3.86	[3.93]	+3.18
NLO	diagram a	-20.85	0.03	0.13	[0.13]	-0.14
NLO	b	10.45	0.01	0.06	[0.06]	+0.05
NLO	c	0.34	0.00	0.01	[0.01]	+0.00
NLO	sum	-10.06	0.03	0.07	[0.08]	-0.09
LO+NLO		684.05	-	-	[3.85]	3.10

Table 10. Table 10: The left-side Table displays the contributions to a μ HVP − LO superscript subscript 𝑎 𝜇 HVP LO a_{\mu}^{\rm HVP-LO} from the various energy regions, including the contribution of the non-HLS channels from the s < 1.05 𝑠 1.05 \sqrt{s}<1.05 GeV region; only total errors are shown. The right-side Table provides the various contributions to a μ subscript 𝑎 𝜇 a_{\mu} and information for Δ ​ a μ = a μ ​ ( exp ) − a μ ​ ( th ) Δ subscript 𝑎 𝜇 subscript 𝑎 𝜇 exp subscript 𝑎 𝜇 th \Delta a_{\mu}=a_{\mu}({\rm exp})-a_{\mu}({\rm th}) for the case of the RS variant of BHLS 2 . The statistical significance of Δ ​ a μ Δ subscript 𝑎 𝜇 \Delta a_{\mu} is given at central value (see text).

Contribution from	Energy Range	${10}^{10}\times a_{\mu }^{\mathrm{HVP}-\mathrm{LO}}$	Contribution from	${10}^{10}\times a_{\mu }$
missing channels	$\sqrt{s}\le 1.05$	$1.28\pm 0.17$	LO-HVP	$686.65\pm 3.01+{[{}_{-0.75}{}^{+1.16}]}_{\mathrm{syst}}+{[{}_{-0.0}{}^{+3.10}]}_{\rho \gamma }$
$J/\psi$		$8.94\pm 0.59$	NLO HVP[43]	$-9.927\pm 0.072$
${\rm Y}$		$0.11\pm 0.01$	NNLO HVP [131]	$1.24\pm 0.01$
hadronic	(1.05, 2.00)	$62.21\pm 2.53$	LBL [132]	$10.34\pm 2.88$
hadronic	(2.00, 3.20)	$21.63\pm 0.93$	NLO-LBL [133]	$0.3\pm 0.2$
hadronic	(3.20, 3.60)	$3.81\pm 0.07$	QED [134, 135]	$11658471.8851\pm 0.0036$
hadronic	(3.60, 5.20)	$7.59\pm 0.07$	EW [136]	$15.36\pm 0.11$
pQCD	(5.20, 9.46)	$6.27\pm 0.01$	Total Theor.	$11659175.96\pm 4.17+{[{}_{-0.75}{}^{+1.16}]}_{\mathrm{syst}}+{[{}_{-0.0}{}^{+3.10}]}_{\rho \gamma }$
hadronic	(9.46, 11.50)	$0.87\pm 0.05$	Exper. Aver.	$11659209.1\pm 6.3$
pQCD	(11.50,$\infty$)	$1.96\pm 0.00$	${10}^{10}\times \Delta a_{\mu }$	$33.14\pm 7.56-{[{}_{-0.75}{}^{+1.16}]}_{\mathrm{syst}}-{[{}_{-0.0}{}^{+3.10}]}_{\rho \gamma }$
Total	1.05 $\to \infty$	$114.67\pm 2.76$	Significance ($n\sigma$)	$4.38\sigma$
	+ missing chann.

Equations

$$a_{\mu }^{\mathrm{had}}={\left[\frac{\alpha m_{\mu }}{3\pi }\right]}^2{\int }_{s_0}^{\infty }\frac{ds}{s^2}\hat{K}(s)R(s),$$

$${\mathcal{L}}_{\mathrm{chiral}}=\frac{f_{\pi }^2}{4}\mathrm{Tr}\left[{\partial }_{\mu }U{\partial }^{\mu }{U}^{†}\right]=-\frac{f_{\pi }^2}{4}\mathrm{Tr}{\left[{\partial }_{\mu }{\xi }_L{\xi }_L^{†}-{\partial }_{\mu }{\xi }_R{\xi }_R^{†}\right]}^2,$$

$$\begin{array}[]{lll}\displaystyle U(x)=\exp{\left[2iP(x)/f_{\pi}\right]}~{}~{},&\xi_{R/L}(x)=\exp{\left[\pm iP(x)/f_{\pi}\right]}~{}~{}\Longrightarrow U(x)=\xi_{L}^{\dagger}(x)\xi_{R}(x)\,,\end{array}$$

$$D_{\mu }{\xi }_{R/L}={\partial }_{\mu }{\xi }_{R/L}-ig{V}_{\mu }{\xi }_{R/L}+ie{\xi }_{R/L}{A}_{\mu }Q,$$

$$\displaystyle\begin{array}[]{lll}\displaystyle{\cal L}_{A}=-\frac{f_{\pi}^{2}}{4}{\rm Tr}\left[D_{\mu}\xi_{L}~{}\xi_{L}^{\dagger}-D_{\mu}\xi_{R}~{}\xi_{R}^{\dagger}\right]^{2}~{}~{},&\displaystyle\displaystyle{\cal L}_{V}=-\frac{f_{\pi}^{2}}{4}{\rm Tr}\left[D_{\mu}\xi_{L}~{}\xi_{L}^{\dagger}+D_{\mu}\xi_{R}~{}\xi_{R}^{\dagger}\right]^{2}\;.\end{array}$$

$${\mathcal{L}}_{\mathrm{HLS}}={\mathcal{L}}_A+a{\mathcal{L}}_V,$$

$${\mathcal{L}}_{{}^{\prime }\mathrm{tHooft}}=\frac{{\mu }^2}{2}{\eta }_0^2+\frac{\lambda }{2}{\partial }_{\mu }{\eta }_0{\partial }^{\mu }{\eta }_0.$$

$$\displaystyle\begin{array}[]{lll}\displaystyle{\cal L}_{A}=-\frac{f_{\pi}^{2}}{4}{\rm Tr}\left[(L-R)X_{A}\right]^{2}~{}~{},&\displaystyle{\cal L}_{V}=-\frac{f_{\pi}^{2}}{4}{\rm Tr}\left[(L+R)X_{V}\right]^{2}~{}~{}~{},\end{array}$$

$${\mathcal{L}}_V^{\mathrm{mass}}=\frac{m^2}{2}\left[(1+{\Sigma }_V)({\rho }_I^2+{\omega }_I^2)+2{\Delta }_V{\rho }_I\cdot {\omega }_I+z_V{\varphi }_I^2+2(1+{\Sigma }_V){\rho }^{+}{\rho }^{-}\right],$$

$$\left(\begin{array}[]{l}\displaystyle\rho_{I}\\[14.22636pt] \displaystyle\omega_{I}\end{array}\right)=\left(\begin{array}[]{l}\displaystyle\rho_{R}\\[14.22636pt] \displaystyle\omega_{R}\end{array}\right)-\Delta_{V}\left(\begin{array}[]{c}\displaystyle h_{V}\omega_{R}\\[14.22636pt] \displaystyle(1-h_{V})\rho_{R}\end{array}\right)\,,$$

$$D_{\mu }{\xi }_{R/L}={\partial }_{\mu }{\xi }_{R/L}-ig{V}_{\mu }^I{\xi }_{R/L}+ie{\xi }_{R/L}{A}_{\mu }Q.$$

$$D_{\mu }{\xi }_{R/L}={\partial }_{\mu }{\xi }_{R/L}-ig\left[V_{\mu }^I+\delta V_{\mu }\right]{\xi }_{R/L}+ie{\xi }_{R/L}{A}_{\mu }Q,$$

$$\left\{\displaystyle\begin{array}[]{l}\displaystyle\delta V_{\mu}^{0}=\frac{\xi_{0}}{\sqrt{2}}\left[\frac{\sqrt{2}\omega_{\mu}^{I}+\phi_{\mu}^{I}}{3}\right]{\rm Diag}[1,1,1]\,,\\[14.22636pt] \displaystyle\delta V_{\mu}^{8}=\frac{\xi_{8}}{\sqrt{2}}\left[\frac{\omega_{\mu}^{I}-\sqrt{2}\phi_{\mu}^{I}}{3\sqrt{2}}\right]{\rm Diag}[1,1,-2]\,,\\[14.22636pt] \displaystyle\delta V_{\mu}^{3}=\frac{\xi_{3}}{\sqrt{2}}\left[\frac{\rho_{I}^{0}}{\sqrt{2}}\right]{\rm Diag}[1,-1,0]\,,\end{array}\right.$$

$$\left\{\begin{array}[]{ll}\displaystyle\frac{\rho_{I}}{\sqrt{2}}\Rightarrow&\displaystyle\frac{\rho_{I}}{\sqrt{2}}[1+\xi_{3}]\,,\\[14.22636pt] \displaystyle\frac{\omega_{I}}{\sqrt{2}}\Rightarrow&\displaystyle\frac{\omega_{I}}{\sqrt{2}}\left[1+\frac{2\xi_{0}+\xi_{8}}{3}\right]+\frac{[\xi_{0}-\xi_{8}]}{3}~{}\phi_{I}\,,\\[14.22636pt] \displaystyle\phi_{I}\Rightarrow&\displaystyle\phi_{I}\left[1+\frac{\xi_{0}+2\xi_{8}}{3}\right]+\frac{\sqrt{2}[\xi_{0}-\xi_{8}]}{3}~{}\omega_{I}\;.\end{array}\right.$$

$$\left\{\begin{array}[]{ll}\displaystyle D_{\mu}\xi_{L}=\partial_{\mu}\xi_{L}-ig\left[V_{\mu}^{I}+\delta V_{\mu}\right]\xi_{L}+i\xi_{L}{\cal L}_{\mu}\,,\\[14.22636pt] D_{\mu}\xi_{R}=\partial_{\mu}\xi_{R}-ig\left[V_{\mu}^{I}+\delta V_{\mu}\right]\xi_{R}+i\xi_{R}{\cal R}_{\mu}\,,\end{array}\right.$$

$$\begin{array}[]{ll}\displaystyle{\cal L}_{\mu}=\displaystyle eQA_{\mu}+\frac{g_{2}}{\sqrt{2}}(W^{+}_{\mu}T_{+}+W^{-}_{\mu}T_{-})\hskip 28.45274pt{\rm and}&\hskip 28.45274pt{\cal R}_{\mu}=\displaystyle eQA_{\mu}~{}~{}~{},\end{array}$$

$$\left\{\begin{array}[]{ll}\displaystyle{\cal L}_{\mu,\nu}=&\displaystyle\partial_{\mu}{\cal L}_{\nu}-\partial_{\nu}{\cal L}_{\mu}-i\left[{\cal L}_{\mu},{\cal L}_{\nu}\right]\,,\\[14.22636pt] \displaystyle{\cal R}_{\mu,\nu}=&\displaystyle\partial_{\mu}{\cal R}_{\nu}-\partial_{\nu}{\cal R}_{\mu}-i\left[{\cal R}_{\mu},{\cal R}_{\nu}\right]\,,\end{array}\right.$$

$$\begin{array}[]{ll}\displaystyle\widehat{{\cal L}}_{\mu,\nu}=\xi_{L}{\cal L}_{\mu,\nu}\xi_{L}^{\dagger}~{}~{},{}{}{}&\widehat{{\cal R}}_{\mu,\nu}=\xi_{R}{\cal R}_{\mu,\nu}\xi_{R}^{\dagger}~{}~{},\end{array}$$

$$\begin{array}[]{ll}\displaystyle\widehat{{\cal V}}_{\mu,\nu}=\frac{1}{2}\left[\widehat{{\cal R}}_{\mu,\nu}+\widehat{{\cal L}}_{\mu,\nu}\right]~{}~{}~{},&\displaystyle\widehat{{\cal A}}_{\mu,\nu}=\frac{1}{2}\left[\widehat{{\cal R}}_{\mu,\nu}-\widehat{{\cal L}}_{\mu,\nu}\right]~{}~{}~{},\end{array}$$

$${\mathcal{L}}_z=z_1\mathrm{Tr}\left[{\hat{\mathcal{V}}}_{\mu ,\nu }{\hat{\mathcal{V}}}^{\mu ,\nu }\right]+z_2\mathrm{Tr}\left[{\hat{\mathcal{A}}}_{\mu ,\nu }{\hat{\mathcal{A}}}^{\mu ,\nu }\right]+z_3\mathrm{Tr}\left[{\hat{\mathcal{V}}}_{\mu ,\nu }{V}^{\mu ,\nu }\right],$$

$$X_{\mu ,\nu }={\partial }_{\mu }{X}_{\nu }-{\partial }_{\nu }{X}_{\mu }-i[X_{\mu },X_{\nu }],$$

$${\mathcal{L}}_{z_3}=\frac{ge}{2}{z}_3{A}_{\mu ,\nu }\left[{\rho }_{\mu ,\nu }^0+\frac{1}{3}{\omega }_{\mu ,\nu }-\frac{\sqrt{2}}{3}{\varphi }_{\mu ,\nu }\right]+\frac{g{g}_2}{4}{z}_3\left[{\overline{V}}_{ud}{W}_{\mu ,\nu }^{-}{\rho }_{\mu ,\nu }^{+}+V_{ud}{W}_{\mu ,\nu }^{+}{\rho }_{\mu ,\nu }^{-}\right]$$

$${\mathcal{L}}_{z_3}=eg{z}_{3}s{A}_{\mu }\left[{\rho }_{\mu }^0+\frac{1}{3}{\omega }_{\mu }-\frac{\sqrt{2}}{3}{\varphi }_{\mu }\right]+\frac{g_{2}g{z}_{3}s}{2}\left[{\overline{V}}_{ud}{\rho }_{\mu }^{+}{W}_{\mu }^{-}+V_{ud}{\rho }_{\mu }^{-}{W}_{\mu }^{+}\right].$$

$${\mathcal{L}}_{\mathrm{mass}}=a{f}_{\pi }^2\mathrm{Tr}{\left[X_V(V_{\mu }^I+\delta V_{\mu })\right]}^2.$$

$$\begin{array}[]{ll}{\cal L}_{\rm mass}=&\displaystyle\frac{m^{2}}{2}\left[(1+\Sigma_{V}+2\xi_{3})~{}[\rho^{0}_{I}]^{2}+(1+\Sigma_{V}+\frac{4}{3}\xi_{0}+\frac{2}{3}\xi_{8})~{}[\omega_{I}]^{2}+z_{V}(1+\frac{2}{3}\xi_{0}+\frac{4}{3}\xi_{8})~{}[\phi_{I}]^{2}\right.\\[14.22636pt] {}{}&\left.\displaystyle+2\Delta_{V}~{}\rho^{0}_{I}\cdot\omega_{I}+\frac{2\sqrt{2}}{3}(1+z_{V})(\xi_{0}-\xi_{8})~{}\omega_{I}\cdot\phi_{I}\right]+m^{2}(1+\Sigma_{V})\,\rho^{+}\cdot\rho^{-}\;.\end{array}$$

$$\begin{array}[]{ll}\displaystyle\Delta_{V}=0~{}~{}~{},&\displaystyle\xi_{0}=\xi_{8}~{}~{}.\end{array}$$

$$\left\{\begin{array}[]{ll}\displaystyle m_{\rho^{\pm}}^{2}=\displaystyle m^{2}~{}(1+\Sigma_{V})\,,\\[14.22636pt] m_{\rho^{0}}^{2}=\displaystyle m^{2}\left[1+\Sigma_{V}+2~{}\xi_{3}\right]\,,\\[14.22636pt] \displaystyle m_{\omega}^{2}=\displaystyle m^{2}\left[1+\Sigma_{V}+\frac{4}{3}~{}\xi_{0}+\frac{2}{3}~{}\xi_{8}\right]=m^{2}\left[1+\Sigma_{V}+2~{}\xi_{0}\right]\,,\\[14.22636pt] \displaystyle m_{\phi}^{2}=\displaystyle m^{2}~{}z_{V}\left[1+\frac{2}{3}~{}\xi_{0}+\frac{4}{3}~{}\xi_{8}\right]=m^{2}~{}z_{V}\left[1+2~{}\xi_{0}\right]\;.\\[14.22636pt] \end{array}\right.$$

$$\left(\begin{array}[]{ll}\rho_{I}\\[14.22636pt] \omega_{I}\\[14.22636pt] \phi_{I}\\[14.22636pt] \end{array}\right)=\left(\begin{array}[]{ccc}\displaystyle 1&\displaystyle-\psi_{\omega}&\psi_{\phi}\\[14.22636pt] \displaystyle\psi_{\omega}&\displaystyle 1&\displaystyle\psi_{0}\\[14.22636pt] \displaystyle-\psi_{\phi}&\displaystyle-\psi_{0}&\displaystyle 1\end{array}\right)\left(\begin{array}[]{ll}\rho_{R}\\[14.22636pt] \omega_{R}\\[14.22636pt] \phi_{R}\\[14.22636pt] \end{array}\right)~{}~{}~{},~{}~{}{\rm also~{}denoted}~{}~{}~{}{\cal V}_{I}=R(U_{3})\cdot{\cal V}_{R}~{}.$$

$$\left\{\begin{array}[]{ll}\displaystyle\rho_{R}-\omega_{R}=0\Rightarrow&\displaystyle\Delta_{V}+\left[\frac{4}{3}\left\{\xi_{0}+\frac{1}{2}\xi_{8}-\frac{3}{2}\xi_{3}\right\}\psi_{\omega}\right]=0\,,\\[14.22636pt] \rho_{R}-\phi_{R}=0\Rightarrow&\displaystyle[1-z_{V}]~{}\psi_{\phi}+\frac{\sqrt{2}}{3}\left\{1+z_{V}\right\}(\xi_{0}-\xi_{8})+\left[~{}\psi_{\omega}~{}\psi_{0}\right]=0\,,\\[14.22636pt] \omega_{R}-\phi_{R}=0\Rightarrow&\displaystyle 3[1-z_{V}]~{}\psi_{0}+[1+z_{V}]\sqrt{2}~{}~{}(\xi_{0}-\xi_{8})-3\left[\psi_{\omega}~{}\psi_{\phi}\right]=0\;.\end{array}\right.$$

$$D_{\rho }(s)=s-m_{\rho }^2-\Pi (s),$$

Full text

BHLS2, a New Breaking of the HLS Model and its Phenomenology

M. Benayouna, L. DelBuonoa, F. Jegerlehnerb,c

a LPNHE des Universités Paris VI et Paris VII, IN2P3/CNRS, F–75252 Paris, France

b Humboldt–Universität zu Berlin, Institut für Physik, Newtonstrasse 15, D–12489 Berlin, Germany

c Deutsches Elektronen–Synchrotron (DESY), Platanenallee 6, D–15738 Zeuthen, Germany

Previous studies have shown that the Hidden Local Symmetry (HLS) Model, supplied with appropriate symmetry breaking mechanisms, provides an Effective Lagrangian (BHLS) able to encompass a large number of processes within a unified framework. This allowed to design a global fit procedure which provides a fair simultaneous description of the $e^{+}e^{-}$ annihilation into 6 final states ($\pi^{+}\pi^{-}$, $\pi^{0}\gamma$, $\eta\gamma$, $\pi^{+}\pi^{-}\pi^{0}$, $K^{+}K^{-}$, $K_{L}K_{S}$), the dipion spectrum in the $\tau$ decay and some more light meson decay partial widths. In this paper, additional breaking schemes are defined which improve the BHLS working and extend its scope so as to absorb spacelike processes within a new framework (BHLS2). The phenomenology previously explored with BHLS is fully revisited in the BHLS2 context with special emphasis on the $\phi$ mass region using all available data samples. It is shown that BHLS2 addresses perfectly the close spacelike region covered by NA7 and Fermilab data; it is also shown that the recent Lattice QCD (LQCD) information on the pion form factor are accurately predicted by the BHLS2 fit functions derived from fits to only annihilation data. The contribution to the muon anomalous magnetic moment $a_{\mu}^{\rm th}$ of these annihilation channels over the range of validity of BHLS2 (up to $\simeq$ 1.05 GeV) is updated within the new BHLS2 framework and shown to strongly reduce the former BHLS systematics. The uncertainty on $a_{\mu}^{\rm th}(\sqrt{s}<1.05$ GeV) is much improved compared to standard approaches relying on direct integration methods of measured spectra. Using the BHLS2 results, the leading order HVP contribution to the muon anomalous moment is $a_{\mu}^{\rm HVP-LO}=686.65\pm 3.01+(+1.16,-0.75)_{\rm syst}$ in units of $10^{-10}$. Using a conservative estimate for the light–by–light contribution, our evaluation for the muon anomalous magnetic moment is $a_{\mu}^{\rm th}=\left[11\,659\,175.96\pm 4.17+(+1.16,-0.75)_{\rm syst}\right]\times 10^{-10}$. The relationship between the dispersive and LQCD approaches to the $\rho^{0}-\gamma$ mixing is also discussed which may amount to a shift of $\delta a_{\mu}[\pi\pi]_{\rho\gamma}=+(3.10\pm 0.31)\times 10^{-10}$ at LO+NLO, presently treated as additional systematics. Taking also this shift into account, the difference $a_{\mu}^{\rm th}-a_{\mu}^{\rm BNL}$ exhibits a significance not smaller than $3.8\sigma$.

1 Introduction

The Standard Model is widely recognized as the (gauge) theory which unifies the whole realm of weak, electromagnetic and strong interactions among quarks, leptons and the various gauge bosons (gluons, photons, $W^{\pm}$, $Z^{0}$). For the physics processes – and quantities – involving strong interactions, QCD is at work under two different regimes tightly connected with the energy involved and the onset of the perturbative regime is a priori expected to occur at high energies. However, as clearly shown by the data recently collected by KEDR [1, 2] on the ratio $R(s)=\sigma(e^{+}e^{-}\rightarrow{\rm hadrons})/\sigma(e^{+}e^{-}\rightarrow\mu^{+}\mu^{-})$, the $\{u,d,s\}$ sector of this ratio reaches the perturbative regime at energies as low as $\simeq 2.0$ GeV (see [3, 4], for instance); above this energy, the observed departures from perturbative QCD predictions – including, of course, the $c\overline{c}$ and $b\overline{b}$ threshold effects therein – are only spikes and narrow bumps associated with the charmonium and bottomonium states which require an additional specific treatment in the mass range where they are located.

However, in the low energy region where the non-perturbative regime of QCD is involved, getting theoretical predictions able to compete with the accuracy of some important experimental measurements may be challenging. This is, in particular, the issue met with the photon hadronic vacuum polarization (HVP), which contributes importantly to the muon anomalous magnetic moment $a_{\mu}$, one of the best measured particle properties. The hadronic part of $a_{\mu}$ is related to the so-called $R(s)$ ratio and is given (at leading order) by :

[TABLE]

where $s_{0}=m_{\pi^{0}}^{2}$ is the lowest hadronic threshold and $\hat{K}(s)$ is a known smoothly varying positive function [5]. The $1/s^{2}$ factor strongly enhances the low energy contribution of the spectrum, i.e. just in the validity domain of non-perturbative QCD. Fortunately, as obvious from the KEDR data [1, 2], one can undoubtedly consider that the non-perturbative regime does not extend to energies larger than $\simeq 2.0$ GeV.

So, the real issue with predictions for objects such as $a_{\mu}$ is to get precise estimates of the effects covered by the non-perturbative regime of QCD. Chiral Perturbation Theory (ChPT)[6, 7], the low energy limit of QCD, is of limited help for the present purpose as its realm does not extend much beyond the 400 $\div$ 500 MeV region and misses the quite important meson resonance mass region.

The most promising approach to the non-perturbative regime of QCD is certainly Lattice QCD (LQCD) which already provides valuable information at low energies [8]. The recently derived LQCD evaluations of $a_{\mu}^{\rm had}$ [9, 10, 11, 12] have been found in accord with the experimental measurement performed at BNL [13, 14]. However, the magnitude of the reported uncertainties is by far too large to fruitfully compare with the already existing BNL datum and, a fortiori, with the measurements expected from the Fermilab experiment [15, 16], already running, or from the experiment planned to start later on at J-PARC [17], as both are expected to improve the uncertainties by a factor of 4. So much progress remains to be done before getting satisfactory uncertainties from LQCD.

This leaves room for low energy Effective Resonance Lagrangian Approaches, which can contribute to improve the knowledge of $a_{\mu}^{\rm had}$ by providing a good description of the rich amount of experimental data collected in the timelike region. Among the richest possible Lagrangians, the elegant Hidden Local Symmetry (HLS) Model [18] is worth to be considered; it has been proven to be equivalent to R$\chi$PT [19, 20], provided consistency with the QCD asymptotic behavior is incorporated. It thus follows that the HLS model is a motivated and constraining QCD rooted framework, moreover, easy enough to handle in phenomenological applications.

The original HLS Model deals with the lowest vector meson nonet and provides a framework for hadron production in $e^{+}e^{-}$ annihilation, naturally bounded by the $\phi$ mass region – i.e. up to $\simeq 1.05$ GeV. It thus represents a tool giving a handle on a mass region contributing for $\simeq 83\%$ of the total muon hadronic VP. The region extending from just above the $\phi$ meson mass to 2 GeV only contributes for 7%, slightly less than the $[2{\rm~{}GeV},\infty]$ region (10%).

As such, the non-anomalous HLS Lagrangian [21] sets up a unified framework which encompasses the $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}$, $e^{+}e^{-}\rightarrow K^{+}K^{-}$ and $e^{+}e^{-}\rightarrow K^{0}\overline{K}^{0}$ annihilation channels; the $\tau^{\pm}\rightarrow\pi^{\pm}\pi^{0}\nu$ spectrum also belongs to the same framework, allowing naturally a simultaneous treatment of all the mentioned annihilation and decay processes. On the other hand, the HLS Model possesses an anomalous sector [22, 18] which also brings the $e^{+}e^{-}\rightarrow\pi^{0}\gamma$, $e^{+}e^{-}\rightarrow\eta\gamma$ and $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ annihilation channels inside the same framework, together with some radiative decay modes and the dipion spectra in the $\eta/{\eta^{\prime}}\rightarrow\pi^{+}\pi^{-}\gamma$ decays as shown in former studies [23, 24]. The annihilation channels just listed exhaust almost completely the processes contributing to $a_{\mu}^{\rm had}$ below 1.05 GeV; indeed the missing channels111The $e^{+}e^{-}\rightarrow{\eta^{\prime}}\gamma$ channel has been considered in [25, 26] and found to contributes to $a_{\mu}^{\rm had}$ at the $10^{-12}$ level only and can thus be safely discarded. ($4\pi$, $2\pi\eta$, …) contributes only $\simeq 2\permil$ of the full HVP. Such a broken HLS (BHLS) Model has already been built up [25, 26] and shown to provide a pretty good simultaneous description of almost all data samples covering the six channels listed above up to $\simeq 1.05$ GeV. This work has also proved that there was no contradiction between the $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}$ and $\tau^{\pm}\rightarrow\pi^{\pm}\pi^{0}\nu$ spectra.

The non-anomalous [21] and anomalous [22] sectors of the HLS Model thus open a unified framework able to encompass a large corpus of data and physics processes. However, as such, the HLS framework – with only the universal vector coupling $g$ as free parameter – cannot pretend to provide a satisfactory simultaneous description of the wide ensemble of high statistics data samples collected by several sophisticated experiments in several annihilation channels.

In order to achieve such a program, the HLS model must be supplied with appropriate symmetry breaking mechanisms not present in its original formulation [18]. A first successful attempt has been done in [25, 26, 27] which has set up a model (BHLS) based on a breaking mechanism originally proposed in [28], hereafter named BKY.

The present work reports on a new breaking scheme for the HLS model which aims at improving the behavior of the form factors in the dipion threshold region and also in the $\phi$ region where the original BHLS [25] meets some difficulty leading to introduce additional systematic uncertainties in the evaluation of the muon $g-2$. The model (BHLS2) presented here will be shown to widen the scope of BHLS, in particular to the spacelike domain and improve the BHLS prediction for the muon $g-2$. Indeed, using variants of BHLS2, together with the original BHLS should allow evaluating the model dependence effects in the estimates of derived physics quantities, noticeably the photon HVP.

The layout of the paper is as follows. Sections 2 gives a brief reminder of the original HLS model while Section 3 briefly reminds the BKY mechanism, still involved in BHLS2. Sections 4 and 5 introduce a new breaking mechanism at the level of the covariant derivative which is the basic object leading to the original HLS model. Section 6 describes the ${\bf O}(p^{4})$ terms [18] of the HLS framework which are included in the BHLS2 Lagrangian.

At this step, one has at hand the first variant of the BHLS2 model, named Basic Solution (BS), and, relying on its vector meson mass term, Section 7 defines its parameter properties. Remarking that there is no fundamental reason why the neutral vector fields involved in physics processes should be their ideal combinations ($\rho_{0}^{I}$, $\omega^{I}$, $\phi^{I}$), one allows for combinations of these via a mechanism called Primordial Mixing described in Section 8. This gives rise to a second variant of the BHLS2 model, named Reference Solution (RS). In Section 9, one first reminds the parametrization of the propagators for the $\omega$ and $\phi$ fields [25] in connection with form factor behaviors at $s=0$; this part of the work also shows that the breaking of nonet symmetry in the vector meson sector must be accompanied by breaking of SU(3) of the same intensity.

The dynamical breaking of vector meson [29, 25] is revisited in Section 10 and the role of anomalous loop corrections is emphasized, especially in order to connect smoothly the timelike and spacelike branches of the pion form factor. The derivation of the pion and kaon form factors is the subject of Section 11. The anomalous sector is fully analyzed in Section 12 and the cross sections for the $e^{+}e^{-}\rightarrow(\pi^{0}/\eta)\gamma$ and $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}\pi^{0}$ annihilations are derived.

Altogether, the above-listed Sections and the Appendix fully report on the BHLS2 model properties and tools. Section 13 provides a comprehensive discussion of the available data samples falling into the BHLS2 scope and their peculiarities. Relying on several preliminary studies, one comments here on the 3 discarded data samples; one should note that the number of data samples found to accommodate with each other within the BHLS2 framework exceeds now 50 and the total number of data points reaches 1237.

Sections 15 and 16 give a full description of all aspects of the fits to the 6 annihilation channels performed within the global BHLS2 framework. Moreover, Section 16 also provides a comprehensive analysis within our global framework of the existing $K\overline{K}$ data samples. The aim is to substantiate the issue between the BaBar, CMD-3 and CMD-2 spectra revealed by fits. A summary of the various fit properties is presented in Section 17.

The description of the spacelike region for the pion and kaon form factors as coming from BHLS2 is the purpose of Section 18. The existing model-independent data for the $\pi^{\pm}$ and $K^{\pm}$ form factors are found to naturally accommodate our global framework, giving support to the low energy behavior predicted by BHLS2 for all meson form factors. This fair agreement extends to the $\pi^{\pm}$ form factor data provided by several Lattice QCD (LQCD) groups. One may infer from these data and from other model-dependent data that BHLS2 should almost certainly apply down to $\simeq-1$ GeV2. The prediction for the $K^{0}$ form factor is shown to cope with expectations at $s\simeq 0$. Information on the charge radii of the $\pi^{\pm}$ and $K^{\pm}$ mesons is provided. Section 19 examines the $P$-wave $\pi\pi$ phase shift prediction and shows its agreement with data and other predictions, including our former BHLS, over an energy range extending up to the $\phi$ mass region.

The values found for some parameters, noticeably the $a$ HLS parameter and some BKY breaking parameters, deserve a specific account given in Section 20.

The main motivation of the present work is the study of the BHLS2 predictions and their consequences on the muon $g-2$. This topic is addressed in Section 21. The part ($\simeq 83$%) of the muon HVP contribution to $a_{\mu}$ covered by BHLS2 is estimated in several contexts to derive $a_{\mu}^{\rm HVP-LO}(\sqrt{s}<1.05~{}{\rm GeV})$ and estimate additional systematics possibly due to modeling effects and to observed tension among some data samples. Complementing this piece by the non-HLS part of the HVP derived by more usual means, one gives our best evaluation of the full HVP and our estimate of the muon anomalous magnetic moment. We also discuss an issue in the relationship between the LO-HVP results derived by dispersive approaches – in particular BHLS/BHLS2 – and the LQCD results, where the $\rho^{0}-\gamma$ mixing appears as a NLO effect. Comparison is done with other published evaluations for $a_{\mu}^{\rm HVP-LO}$. Finally, Section 22 summarizes our results and conclusions.

2 Outline of the Hidden Local Symmetry (HLS) Model

Let us briefly outline the derivation of the HLS model222See, for instance, [18] for a full derivation.. The usual ChPT Lagrangian [6, 7] can be written in two different manners [18, 30] :

[TABLE]

where $f_{\pi}$ (=92.42 MeV) is the pion decay constant and :

[TABLE]

when working in the so-called unitary gauge which removes a scalar field term in the definition of $\xi_{R/L}(x)$; $P(x)$ is the pseudoscalar (PS) field matrix. The hidden local symmetry which gives its name to the HLS Lagrangian has transformation properties which affect $\xi_{R/L}(x)$ letting $U$ unchanged. Ignoring for the moment the weak sector to ease the discussion, the HLS Lagrangian is derived by replacing in Equation (2) the usual derivative by the covariant derivative :

[TABLE]

where $A_{\mu}$ is the photon field, $Q={\rm Diag}[2/3,-1/3,-1/3]$ the quark charge matrix and $V_{\mu}$ is the vector field matrix; the expressions for $P$ and $V$ are the usual ones – fulfilling the $U(3)$ flavor symmetry – and can be found in [18, 30, 25]. In the expressions for $D_{\mu}\xi_{R/L}$, the universal vector coupling $g$ occurs beside the unit electric charge $e$. The $\omega$ and $\phi$ fields occurring in the diagonal of the $V$ matrix, also denoted $V^{I}$ below, are the so-called ideal combinations generally denoted $\omega_{I}$ and $\phi_{I}$.

Substituting the covariant derivative Equation (4) to the usual derivative in Equation (2), ${\cal L}_{\rm chiral}$ becomes the first HLS Lagrangian piece denoted ${\cal L}_{A}$ while another piece ${\cal L}_{V}$ shows up which vanishes in the reversed substitution $D_{\mu}\Rightarrow\partial_{\mu}$ :

[TABLE]

The full HLS Lagrangian is then defined by :

[TABLE]

where $a$ is a parameter specific of the HLS approach [18]. The usual VMD framework is obtained by setting $a=2$ [30] while the presently reported phenomenology prefers larger values [24, 29]. The explicit expression for this unbroken ${\cal L}_{\rm HLS}$ can be found fully developed in [30]. Let us note that ${\cal L}_{\rm HLS}$ contains a photon mass term; however this Reference also showed that the loop dressing (in the HLS context) cancels out the physical photon mass and thus, this mass term can be safely ignored in phenomenological studies.

So, the unperturbed (i.e. unbroken) HLS Lagrangian depends on only two parameters ($g$ and $a$) which can be adjusted using data; it is obviously unrealistic to expect this basic ${\cal L}_{\rm HLS}$ to describe precisely the large amount of data covering its scope with so few parameters. To have any chance to account for real data, additional input should enter the HLS Lagrangian to make it more flexible; this is what motivates the introduction of breaking procedures within the HLS framework.

As already noted, both HLS Lagrangian pieces fulfill a $U(N_{f})\times U(N_{f})$ symmetry and not $SU(N_{f})\times SU(N_{f})$. To reduce this symmetry which introduces a ninth PS meson, one includes [31, 25] the ’t Hooft determinant terms [32] which break the axial $U(1)$ symmetry; this turns out to complement the ${\cal L}_{A}$ Lagrangian piece with :

[TABLE]

The singlet mass term manifestly breaks nonet symmetry in the PS sector. Nothing analogous affects the vector sector.

3 Breaking the HLS Lagrangian I : The BKY Mechanism

Up to now, a single breaking mechanism for the HLS Lagrangian has been proposed [28]; it is named here BKY after its proponent names. However, undesirable properties of the original BKY proposal have led to modifications of this breaking scheme [33, 34] originally proposed to break (only) the flavor symmetry of the Lagrangian. The ”new scheme” variant from [34] has been analyzed with data and its results shown to compare fairly well to ChPT expectations [31]. It has been extended to include isospin breaking effects in [25] following the lines of [35]. This (modified) BKY breaking scheme, intended to cover isospin and $SU(3)$ breaking effects within the same framework, has been proved fairly successful, allowing for a high quality global fit of (almost) all available data covering the validity range of the HLS model [25, 26, 27] : $e^{+}e^{-}$ annihilations, $\tau$ decay spectra and radiative decays of light flavor mesons; in particular, the dipion spectra in the $\eta/{\eta^{\prime}}\rightarrow\pi^{+}\pi^{-}\gamma$ decays are fairly well predicted [23, 24], especially the detailed form of the drop-off in the $\rho-\omega$ interference region as recently measured at BESIII [36, 37]. Nevertheless, as already stated in [25, 26, 27], the description of the threshold region in the dipion spectrum and of the $\phi$ mass region in the three pion annihilation channel deserves improvements; this issue is addressed by the present paper.

As BKY is one of the breaking mechanisms used in the present work, let us briefly remind how it is implemented. To lighten writing, we use the notations $L=D_{\mu}\xi_{L}~{}\xi_{L}^{\dagger}$ and $R=D_{\mu}\xi_{R}~{}\xi_{R}^{\dagger}$. The (modified and extended) BKY breaking which also underlies the broken HLS (BHLS) model developed in [25] turns out to modify Equations (5) as follows :

[TABLE]

where $X_{A/V}=$Diag$(q_{A/V},y_{A/V},z_{A/V})$ are constant real matrices. In practice, one prefers setting $q_{A/V}=1+(\Sigma_{A/V}+\Delta_{A/V})/2$ and $y_{A/V}=1+(\Sigma_{A/V}-\Delta_{A/V})/2$. As $z_{A/V}$ are affecting the $s\overline{s}$ entries, their departure from 1 can be large compared to $q_{A/V}$ and $y_{A/V}$ – which refer to resp. the $u\overline{u}$ and $d\overline{d}$ entries; so, the $\Sigma$’s and $\Delta$’s are expected small. Within the previous broken HLS framework – hereafter named BHLS – [25, 26], one got $z_{A}\simeq[f_{K}/f_{\pi}]\simeq 1.5$, $z_{V}\simeq 1.2$ while the $\Sigma$’s and $\Delta$’s were found at the few percent level.

Once these breakings are applied, the PS kinetic energy term contained in ${\cal L}_{A}$ is no longer diagonal and a PS field redefinition has to be performed in order to restore the kinetic energy term to canonical form. The procedure is fully described in Section 4 of [25] and is valid unchanged in the present work; it is not repeated here.

In the BHLS model [25], after the $X_{V}$ breaking, the vector meson mass term is no longer diagonal and writes ($m^{2}=ag^{2}f_{\pi}^{2}$) :

[TABLE]

discarding the unessential $K^{*}$ sector and neglecting terms of degree greater than 1 in the breaking parameters. This mass term can be diagonalized by a 45∘ rotation, a quite unacceptable solution as it does not give a smooth limit in the symmetry limit $\Delta_{V}\rightarrow 0$. A smooth solution is obtained by the following transform to renormalized (R) fields :

[TABLE]

where $h_{V}$ is a parameter submitted to fit, together with $z_{V}$, $\Sigma_{V}$ and $\Delta_{V}$. This solution, adopted in [25, 26, 27], provides a quite satisfactory description of the data. This transformation introduces a $\Delta_{V}\partial\rho_{R}\partial\omega_{R}$ term in the kinetic energy term of the vector mesons; this issue is known to imply the occurrence of wave-function renormalization factors absorbed in the effective couplings [38]; in the case of BHLS, they can be considered absorbed in the breaking parameters.

The full account of the BHLS model also requires the dynamical breaking [29, 25] generated by PS loop effects which also calls for an additional $s$-dependent renormalization step. It is rediscussed below within the context of the present work.

4 Breaking the HLS Lagrangian II : The Covariant Derivative (CD) Handle

As clear in Equations (8), BKY breaks the symmetry of the HLS Model in the very definition of its Lagrangian. This is quite legitimate and strongly validated by its remarkable description of a considerable amount of data [25, 26, 27]. Nevertheless, it is of interest to explore other possibilities to improve the description of the data in some specific mass regions and, then, reduce or cancel out the additional systematic uncertainties reported in [27].

Let us focus on Equation (4) which defines the original HLS covariant derivative :

[TABLE]

The superscript $I$ assigned to the vector field matrix $V$ aims at reminding that the Isospin zero mesons are the so-called ideal combinations $\omega^{I}$ and $\phi^{I}$. These naturally occur in the $U(3)$ symmetric expression of $V$.

The vector meson term in $D_{\mu}\xi_{R/L}$ can be written $gV^{I}=g\sum_{a=0,8}V_{a}T_{a}$ where the $V_{a}$’s are the vector meson fields and $T_{a},~{}(a=1,\cdots 8)$ are the Gell-Mann matrices normalized such that ${\rm Tr}[T_{a}T_{b}]=\delta_{ab}/2$; $T_{0}$ is the unit matrix appropriately normalized : $T_{0}=I/\sqrt{6}$. Thus, the (nonet) $U(3)_{V}$ symmetry of HLS is obtained by :

• •

Plugging the appropriately constructed $V^{I}$ into the covariant derivative $D_{\mu}\xi_{R/L}$,

• •

Assuming the universality of the vector coupling $g$ to the external world.

Therefore, we propose a direct breaking of the covariant derivative, a new tool independent of the BKY mechanism. Such a breaking mechanism is expressed by the modified covariant derivative :

[TABLE]

where $\delta V_{\mu}$ can be chosen to break the $U(3)_{V}$ symmetry in a controlled way. For instance, breaking solely the nonet symmetry of $V^{I}$ turns out to allow the coupling to the singlet $V^{0}_{\mu}$ to differ from those of the octet fields, preserving in this way the $SU(3)_{V}$ symmetry. It will become clearer below why a more systematic approach must be preferred.

Identifying the field combinations associated with each of the canonical $T_{a}$ matrices, one is led to define the following components which can participate to $\delta V_{\mu}$ separately or together :

[TABLE]

in terms of the usual ideal field combinations.

The (free) breaking parameters $\xi_{0}$, $\xi_{8}$ and $\xi_{3}$ are only requested to be real in order that $\delta V_{\mu}$ is hermitian as $V_{\mu}^{I}$ itself. Clearly, $\delta V_{\mu}^{0}$ defines a breaking of the nonet symmetry down to $SU(3)_{V}\times U(1)_{V}$, $\delta V_{\mu}^{8}$ rather expresses the breaking of the $SU(3)_{V}$ symmetry, while $\delta V_{\mu}^{3}$ is related to a direct breaking of Isospin symmetry in the vector sector.

One could also introduce some breaking affecting the $\rho^{\pm}$ entries of the $V^{I}$ matrix333For the sake of simplicity, the case for the $K^{*}$ sectors is left aside; substantially, this would introduce two additional real parameters $\xi_{4,5}$ and $\xi_{6,7}$ which may be fixed by the two radiative $K^{*}$ decays; additional constraints may come from the $\tau$ decay to $K\pi\nu$. . The $\rho^{\pm}$ entries are given by the sum of the Gell-Mann matrix terms $T_{1}$ and $T_{2}$; forcing a breaking for these entries requires two real parameters which should be equal ($\xi_{1,2}$) in order to preserve hermiticity.

So, $\xi_{3}$ and this $\xi_{1,2}$ could summarize the whole isospin breaking effects in the vector meson side; however, one can choose $\xi_{1,2}=0$ as no theoretical value for the (unbroken) universal coupling $g$ is presently available. Indeed, all vector couplings in the HLS Lagrangian could then be reexpressed in terms of $g^{\prime}=g(1+\xi_{1,2})$; in this case, all physics quantities will depend on $g^{\prime}$ and $\xi_{i}^{\prime}=\xi_{i}/(1+\xi_{1,2})$ ($i=0,8,3$) without any other occurrence of $\xi_{1,2}$ dependency.

Therefore, phenomenologically, the $\xi_{i}$’s and $g$ itself are defined up to a normalization factor presently out of reach. This scaling property has obviously no consequence on physics observables like cross-sections or form factors.

So, from now on, one assumes the maximal breaking of $U(3)_{V}$ experimentally accessible : $\delta V_{\mu}=\delta V_{\mu}^{0}+\delta V_{\mu}^{8}+\delta V_{\mu}^{3}$. Compared with the original $V^{I}$ entries, this turns out to modify only the diagonal entries of $V^{I}$ by the following substitutions444 Remind that a $1/\sqrt{2}$ term is factored out in the definition of $V^{I}$. :

[TABLE]

Then, the $U(3)_{V}$ breaking of the covariant derivative generates a breaking of the vector coupling universality. For this purpose, one should note that a vector mixing is generated, except if $\xi_{0}=\xi_{8}$. For the sake of conciseness, this mechanism is referred to below as CD breaking.

5 BHLS2 : A New Broken HLS Lagrangian

Here, we define a new version of the broken HLS Lagrangian, merging the two breaking schemes just presented – BKY and CD – as these two ways of breaking are conceptually independent. Their contributions are thus complementing each other. In order to fully take into account the electroweak sector, one should modify correspondingly the covariant derivative to :

[TABLE]

where, as usual[18] :

[TABLE]

discarding the $Z^{0}$ boson terms of no concern for the phenomenology we address. The $T_{\pm}$ matrices are constructed out of the matrix elements $V_{ud}$ and $V_{us}$ of the Cabibbo-Kobayashi-Maskawa matrix and can be found in [18]. Finally, the weak coupling $g_{2}$ is related to the Fermi constant by $g_{2}=2m_{W}\sqrt{G_{F}\sqrt{2}}$.

The expression for the non-anomalous HLS Lagrangian pieces given in Section 2 remains formally valid, being understood that the covariant derivatives are modified according to Equations (14–15). As a result, the ${\cal L}_{V}$ piece substantially differs from its partner in [25] while the ${\cal L}_{A}$ piece in the present scheme is strictly identical to those given in [25]. The non-anomalous BHLS2 Lagrangian is given expanded in Appendix A.

6 The Order ${\cal O}(p^{4})$ Terms of the HLS Lagrangian

Beside the ${\cal L}_{A}$ and ${\cal L}_{V}$ pieces which are ${\cal O}(p^{2})$, the HLS approach also possesses terms of order ${\cal O}(p^{4})$(see Section 4.3 in [18]), which modify the $V-\gamma/W$ couplings in a specific way. As the role of such terms has never been really examined in phenomenology555See, however, the discussion in [23] which is revisited here., it looks worthwhile examining their relevance when dealing with data of high accuracy.

Using ${\cal L}_{\mu}$ and ${\cal R}_{\mu}$ just given, one first defines :

[TABLE]

and also :

[TABLE]

where $\xi_{R}=\xi_{L}^{\dagger}=\exp{[iP/f_{\pi}]}$. Furthermore defining :

[TABLE]

the ${\cal O}(p^{4})$ Lagrangian writes [18] :

[TABLE]

where one has generically defined :

[TABLE]

and where the $z_{i}$ are constants not theoretically constrained. The $z_{3}$ term, the most involved in the phenomenology one addresses, writes :

[TABLE]

at lowest order in the expansion of the $\xi_{L/R}$ fields; here one has kept the unbroken $V_{\mu}$ matrix for clarity. By integrating by part and fixing the gauge condition to $\partial_{\mu}X^{\mu}=0$ for all vector fields, this piece becomes :

[TABLE]

Thus, the ${\cal L}_{z_{3}}$ piece exhibits $s$-dependent parts for the $V-\gamma$ and $V-W$ transition amplitudes which complement their constant parts provided by the usual ${\cal O}(p^{2})$ HLS Lagrangian. The ’broken’ version of Equation (22) is given in Appendix A.3.

7 The Basic Solution (BS) and the Vector Mass Term

In the context where both the BKY and CD breaking mechanisms are at work, the vector mass term in ${\cal L}_{V}$ becomes :

[TABLE]

Ignoring the $K^{*}$ sector and setting $m^{2}=ag^{2}f_{\pi}^{2}$, it writes :

[TABLE]

Therefore the CD breaking, via its $\delta V_{3}$ component, allows for a $\rho^{0}-\rho^{\pm}$ (HK) mass difference, provided fits favor non-vanishing $\xi_{3}$ values. Moreover, all the components of $\delta V$ defined above contribute to generate different HK masses for the $\rho^{0}$ and $\omega$ mesons as $\Delta m^{2}_{\rho^{0}-\omega}=m^{2}[\xi_{3}-(2\xi_{0}+\xi_{8})/3]$. This contrasts with the previous BHLS model [25] where $m^{2}_{\rho^{0}}=m^{2}_{\rho^{\pm}}=m^{2}_{\omega}$ is fulfilled. Because of the non-vanishing HK mass difference between the $\rho^{0}$ and $\omega$, one can already expect an improved treatment of the dipion threshold and spacelike regions in the BHLS2 framework compared to BHLS.

Equation (24) is manifestly diagonalized by setting :

[TABLE]

This solution – which lets $z_{V}$ unconstrained – defines our Basic Solution (BS). Within the framework of this solution, the breaking parameters to be derived from fits are $\Sigma_{V}$, $z_{V}$, $\xi_{3}$, $\xi_{0}~{}(=\xi_{8})$. Here one feels the issue with assuming solely nonet symmetry breaking; indeed, as this turns out to state $\xi_{8}=0$, the CD breaking implies $\xi_{0}=0$ and thus breaking HLS via the covariant derivative intrinsically implies that nonet symmetry cannot be solely broken.

At leading order in breaking parameters, the vector meson mass term in ${\cal L}_{V}$ becomes diagonal and one has :

[TABLE]

So, BHLS2 a priori yields different HK masses for all the vector mesons.

8 The Primordial Mixing (PM) and the Reference Solution (RS)

Another unused mechanism can be invoked; indeed, independently of the BKY and CD mechanisms just defined, one can always consider that the neutral vector fields $\rho^{0},~{}\omega,~{}\phi$ involved in physical processes are not directly the ideal ones but combinations of these. For this purpose, let us define an infinitesimal rotation matrix $R(U_{3})~{}={\bf 1}+{\cal O}(\epsilon)$, which associates to the ideal field vector ${\cal V}_{I}=(\rho_{I},\omega_{I},\phi_{I})$ a (first step) renormalized field vector ${\cal V}_{R}=(\rho_{R},\omega_{R},\phi_{R})$ via :

[TABLE]

This is also a quite legitimate tool to extend the flexibility of BHLS2.

The matrix $R(U_{3})$ fulfills $R(U_{3})\widetilde{R}(U_{3})=\widetilde{R}(U_{3})R(U_{3})=\bf{1}$, up to terms of order ${\cal O}(\epsilon^{2})$ which are discarded in our ${\cal O}(\epsilon)$ approach. As it is a rotation, this transformation preserves the canonical structure of the vector field kinetic term provided by the Yang-Mills Lagrangian (up to ${\cal O}(\epsilon^{2})$ terms).

However, the (real) $\psi$ (Euler) angles from transformation in Equation (27) should be chosen in such a way that ${\cal L}_{\rm mass}({\cal V}_{I})$ remains canonical in the change of fields ${\cal V}_{I}\rightarrow{\cal V}_{R}$, i.e. the crossed terms in ${\cal L}_{\rm mass}({\cal V}_{R})$ should be canceled out. Restarting from the mass term given by Equation (24), three conditions should be fulfilled :

[TABLE]

The last bracketed terms in these expressions can be discarded as clearly of second order in the breaking parameters $\{\xi_{i},~{}\psi_{j}\}$; this already implies that $\Delta_{V}=0$ and that $\psi_{\omega}$ becomes unconstrained.

As for the parameter $z_{V}$ – generated by the BKY breaking in its $s\overline{s}$ entry – the situation deserves further comments :

• •

Either $z_{V}$ is found such that666 One may think that a reference magnitude for any generic $\epsilon$ is ${\cal O}(\epsilon)\simeq e=\sqrt{4\pi\alpha_{\rm em}}\simeq 0.3$. ${1-z_{V}}\simeq{\cal O}(\epsilon)$, then none of the $\psi$ angles is constrained at order ${\cal O}(\epsilon)$ and, moreover, the last two Equations (28) then imposes $\xi_{0}=\xi_{8}$.

• •

Or $z_{V}$ is returned by fits such that $1-z_{V}$ is not ${\cal O}(\epsilon)$. Then, at first order in breakings, the last two equations imply $\psi_{\phi}=\psi_{0}$ and that they are proportional to $\xi_{0}-\xi_{8}$.

Preliminary fits, performed using this mechanism – named from now on Primordial Mixing (PM) – indicate that $1-z_{V}$ is in the range $(1\div 2)\epsilon$ and thus, in this case also $\xi_{0}=\xi_{8}$ can be imposed. Therefore, besides the Basic Solution (BS), one gets an additional one, we name it Reference Solution (RS), which also includes the Primordial Mixing (PM). Also considering the former BHLS model [25, 26, 27], one thus has at disposal three different models. This allows for a better insight into the systematics and model dependence effects.

9 Vector Meson Propagators and Form Factors

Let us anticipate on form factor calculations using the Lagrangian given in Appendix A and have a look at the form factor values at $s=0$ – before introducing mixing effects due to loop corrections. For this purpose, let us start with a preliminary digression on vector meson propagators, especially those for the $\omega$ and $\phi$ mesons, as these play an important role at the chiral point and in the physics of the (close) spacelike region which is also addressed in the present work.

For the $\rho$ meson, following the pioneering work [39], the inverse propagator at one loop can be written (also see [23, 24, 25]) :

[TABLE]

where the self-energy $\Pi(s)$ is the sum of the loops allowed by the non-anomalous Effective Lagrangian given in Appendix A; in our context, these are essentially pion and kaon loops777 Loops generated by the anomalous HLS Lagrangian pieces [18, 22] also contribute but can be discarded [25] as reminded in the next Section; however, some anomalous loops, suppressed by a factor of $e^{2}$, can have to be kept as they may play some role in specific expressions as will be emphasized below. . $\Pi(s)$ is a real analytic function ($\Pi(s^{*})=[\Pi(s)]^{*}$) which vanishes at $s=0$ because of current conservation, and is real for negative $s$ – in fact, this property holds already below $m_{\pi^{0}}^{2}$, the lowest energy hadronic threshold. Hence, this implies that $D_{\rho}(0)=-m_{\rho}^{2}$, where $m_{\rho}^{2}$ is displayed in Equations (26).

In principle, this also applies to the $\omega$ and $\phi$ propagators; however, taking into account their narrow character, it looks unmotivated to enter into such complications when fitting objects like $e^{+}e^{-}$ annihilation cross sections. In this case, phenomenology has lengthily experienced a broad success using Breit-Wigner (BW) lineshapes. However, as discussed in [25], some physics quantities like the contribution of the $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}$ cross-section to the estimate for the muon HVP, may require some care about the behavior of approximations close to the chiral point, as this region gives an enhanced contribution to the muon HVP. For this purpose, [25] proposed the following modified BW lineshape for the $\omega$ and $\phi$ mesons :

[TABLE]

where $(\widetilde{m}_{V},\widetilde{\Gamma}_{V})$ are parameters to be fitted and the $m_{V}^{2}$’s are the relevant HK square masses as displayed in Equations (26). Numerically, this expression for the $D_{V}(s)$’s is close to usual BW lineshapes888Setting $\widetilde{m}_{V}=m_{V}$ in this approximation gives $D_{V}(s)=s-m_{V}^{2}+is\widetilde{\Gamma}_{V}/m_{V}$; more common choices for the imaginary part here are $m_{V}\widetilde{\Gamma}_{V}$ or $\sqrt{s}\widetilde{\Gamma}_{V}$, quite analogous to our own choice.. This BW-like parametrization imposes the $\omega$ and $\phi$ meson inverse propagators to fulfill $D_{V}(0)=-m^{2}_{V}$ as already does the $\rho$ meson inverse propagator; this property of the propagators at $s=0$ is essential to recover the expected values for the pion and kaon form factors at the origin within the HLS framework.

On the other hand, if one wishes to examine the analytic continuation of the $\omega$ and $\phi$ propagators somewhat inside the spacelike region, it might be desirable to make it real there; this can be achieved by simply replacing $\widetilde{\Gamma}_{V}$ in Equation (29) by $\theta(s)\widetilde{\Gamma}_{V}$ – or rather, $\theta(s-m_{\pi^{0}}^{2})\widetilde{\Gamma}_{V}$.

Using the Basic Solution (BS) or the Reference Solution (RS) to redefine the physical vector meson fields, one can show, using the Lagrangian given in the Appendix A, that :

[TABLE]

up to terms of order ${\cal O}(\epsilon^{2})$. One can check that the conditions $\Delta_{V}=0$ and $\xi_{0}=\xi_{8}$ are essential in the derivation of these constraints.

So, the breaking parameters to be derived from fits are $\Sigma_{V}$, $\xi_{3}$, $\xi_{0}~{}(=\xi_{8})$ and $z_{V}$ when working within the BS framework; one should also include the $\psi$ rotation angles when extending BHLS2 to the Reference Solution.

It is worthwhile reminding that the effects of $\delta V_{\mu}^{0}$ and $\delta V_{\mu}^{8}$ are intimately intricated within our new Lagrangian frameworks. Here, indeed, a breaking of solely nonet symmetry (i.e. $\xi_{0}\neq 0$) cannot occur if not accompanied by a breaking of the $SU(3)_{V}$ symmetry of the same intensity.

In conclusion, the present model (BHLS2) is not a trivial variant of BHLS [25]. Some important properties of BHLS2 versus BHLS mentioned below, will further substantiate this statement.

10 Dynamical Mixing of Vector Mesons

10.1 The Squared Mass Matrix at One Loop …

As previously noted [29, 25] and reminded above, all variants of the HLS Model, including BHLS2 (see Appendix A.1), exhibit $\rho^{0}/\omega/\phi\rightarrow K\overline{K}$ couplings. This implies that, at one loop, the squared mass matrix of the $\rho^{0}/\omega/\phi$ system receives non-diagonal entries, i.e. the vector fields occurring in the Lagrangian Eq (114) are no longer mass eigenstates. Mass eigenstates are constructed using perturbative methods as performed in [29, 25]. As the loops are (real analytic) functions of $s$, the mass eigenstates become also $s$-dependent. Within the BHLS and BHLS2 frameworks, the kaon loops produce a $s$-dependent difference between the $\rho^{\pm}-W^{\pm}$ and $\rho^{0}-\gamma$ transition amplitudes which has provided the first solution [29, 25] to the long standing $e^{+}e^{-}$ versus $\tau$ puzzle [40, 41, 42]. Nevertheless, in another Effective Lagrangian context, this $s$-dependent difference can be successfully generated by other means [43].

Basically, the Dynamical (or Loop) Mixing of Vector Mesons has been first defined999See also [25]. and motivated in [29]. In order to ease the reading, we remind it and emphasize the new features provided by the BHLS2 context.

The (squared) mass matrix of the $\rho^{0}/\omega/\phi$ sector can be written :

[TABLE]

where :

[TABLE]

The Higgs-Kibble masses $m_{V}$ are displayed in Equations (26), and $\Pi_{\pi\pi}(s)$ is the pion loop weighted by the square of the $\rho^{0}\pi^{+}\pi^{-}$ coupling constant. Because all the loop functions are real analytic function of $s$, $M_{0}^{2}(s)$, $\delta M^{2}(s)$, and hence $M^{2}(s)$, are hermitian analytic matrices (e.g. fulfilling $[X(s)]^{\dagger}=X(s^{*})$).

In order to construct explicitly $\delta M^{2}(s)$, it is worth reexpressing some $VPP$ coupling constants in a suitable manner. One can write :

[TABLE]

having defined $\widetilde{G}=ag/(4z_{A})$. Using Equations (25), common to both the BS and RS variants, one finds :

[TABLE]

and :

[TABLE]

When dealing specifically with the BS variant, these two sets of equations are used by simply dropping out the $\psi_{\alpha}$ parameters.

Denoting by resp. $\Pi_{\pm}(s)$ and $\Pi_{0}(s)$ the amputated charged and neutral kaon loops, the ${\cal V}_{R}^{i}\rightarrow{\cal V}_{R}^{j}$ transition amplitudes ($i,j=\rho^{0},~{}\omega,~{}\phi$) are given by :

[TABLE]

using the notations just defined. Then, the elements of the $\delta M^{2}(s)$ matrix are :

[TABLE]

As in the former BHLS [25], $\delta M^{2}(s)$ is non-diagonal for non-zero values of $s$. Therefore, at one-loop order, the field combinations defined in both the BS and RS variants are not mass eigenstates, as in the previous BHLS release. This calls for a mass-dependent diagonalization which is reminded in the next Subsection.

Let us make a few more remarks about the non-diagonal entries of $\delta M^{2}$. In the no breaking limit, where the neutral and charged kaon loops coincide, the $\rho^{0}-\omega$ and $\rho^{0}-\phi$ entries identically vanish; however, the $\omega-\phi$ entry, proportional to the sum of the neutral and charged kaon loops, does not vanish, indicating that the $\omega-\phi$ mixing is always at work at one loop order within the HLS framework.

In the former BHLS, however, if no breaking is applied in the ${\cal L}_{A}$ sector and if there is no mass breaking among the kaons, there was no $\rho^{0}-\omega$ and $\rho^{0}-\phi$ transitions. This remains true within the BS framework defined above. However, within the RS framework, the $\psi_{\alpha}$’s come in such a way that none of the entries of $\delta M^{2}(s)$ vanishes.

10.2 …And its Diagonalization

The diagonalization of ${\cal L}_{\rm mass}$ at one loop order is performed by means of another ($s$-dependent) rotation matrix to a second step of vector field renormalization $R\Rightarrow R^{\prime}$ :

[TABLE]

so that the full renormalization procedure is obtained by :

[TABLE]

The $R^{\prime}$ fields are the physical vector meson fields. When working within the Basic Solution variant, ${\cal V}_{R}\equiv{\cal V}_{I}$ and $R(U_{3})$ can be dropped out from Equation (39); otherwise, within the RS framework, $R(U_{3})$ is given in Equation (27).

The complex ”angles” in Equation (38) can be derived from the $\delta M^{2}(s)$ matrix elements :

[TABLE]

where the numerators manifestly depend on the kaon loops only and the $\lambda$’s are the eigenvalues of $M^{2}(s)$ matrix. At leading order in breaking parameters, these write [25] :

[TABLE]

the vector meson masses being those displayed in Equations (26). As the three $m_{V}^{2}$’s occurring here are different, the three angles defined above vanish at $s=0$; in contrast, within the former BHLS context [25], because the HK masses for $\rho^{0}$ and $\omega$ were equal, $\alpha(0)$ has a non-zero limit at $s=0$, compromising a proper analytic continuation of the form factors downwards into the spacelike region101010See, in particular, the comments in Section 6.3 of [25]..

In summary the full transform from ideal to fully renormalized $\rho_{0}$, $\omega$, $\phi$ fields is given either by (BS solution) :

[TABLE]

or by (RS solution) :

[TABLE]

where breaking terms of order greater than 1 have been discarded. We will refer from now on to Equation (43) for both the BS and RS variants, being understood that in the former case, the $\psi_{i}$’s are zero. When relevant, the entries in the matrix just above will be denoted $\widetilde{\alpha}(s)$, $\widetilde{\beta}(s)$ and $\widetilde{\gamma}((s)$, using obvious notations. Moreover, in the following, the fully renormalized fields can be either indexed by $R^{\prime}$ – if useful – or deprived of indexation to lighten expressions. For clarity, in all displayed Lagrangian pieces, one always writes the couplings for the ideal $\rho_{0}$, $\omega$, $\phi$ fields. To go to physical fields, one has to apply the appropriate transformation $R(I\Rightarrow R^{\prime})$ and collects the contributions in order to yield the physical $\rho_{0}$, $\omega$, $\phi$ couplings.

10.3 The Effects of Anomalous Loops

When defining the loops contributing to $\delta M^{2}(s)$, we have only considered those generated by the non-anomalous HLS Lagrangian. However, the anomalous (FKTUY) HLS Lagrangian pieces [22, 18] allow for (anomalous) couplings also generating loop contributions111111The Yang-Mills part of the full Lagrangian may also contribute with other loops like $K^{*}\overline{K}^{*}$. The $VPPP$ Lagrangian generate 2-loop contributions with coefficients of ${\cal O}(\epsilon)$ or ${\cal O}(\epsilon^{2})$ to the transitions among the vector mesons; they are presently ignored.

to $\delta M^{2}(s)$ as, for instance, $K^{*}\overline{K}$ or 3-pion loops. As all the two particle channels have thresholds above our fitting region, their loops – of order ${\cal O}(g^{2})$ – are real in our fitting range and, thus, can be discarded, assuming their effects numerically absorbed by the subtraction polynomials of the other (pion and kaon) loops involved [25].

However, one should also note that the anomalous (FKTUY) sector of the HLS Lagrangian generates couplings of the neutral vector mesons to the $\pi^{0}\gamma$, $\eta\gamma$ and ${\eta^{\prime}}\gamma$ final states, even when setting $c_{3}=c_{4}$ [25, 27]. Then the FKTUY sector generates the corresponding loops, which have couplings of order ${\cal O}(g^{2}e^{2})$ and develop tiny imaginary parts far inside our fitting region for the first two and close to its upper border for the third one. These contributions, and their imaginary parts, are higher order and can generally be discarded; nevertheless, their tiny imaginary parts were accounted for in our previous studies [25, 26, 27] by simply adding some fixed $i\varepsilon$ to the eigenvalue differences $\lambda_{i}(s)-\lambda_{j}(s)$ which occur in the denominators of the mixing ”angles” shown in Equations (40).

So, these (complex) mixing ”angles” exhibit a dependence upon differences of the $M^{2}(s)$ eigenvalues. However, these eigenvalues – see Equations (41) – being $s$-dependent, their differences may exhibit zeros for real values of $s$ which, accordingly, generate real poles for the mixing ”angles”. Poles occurring at real negative $s$ are not a real issue as they can be dropped out by means of customary dispersive techniques. For the others, the $i\varepsilon$ trick emphasized just above, permits to avoid unwanted singularities on the physical region $s\geq m_{\pi^{0}}$.

In the previous BHLS version[25], the main purpose of this ad hoc $i\varepsilon$ was to force121212As clear from the first Equations (40), the numerator of $\alpha(s)$ behaves as $s$ close to the origin and, as $m_{\rho^{0}}=m_{\omega}$ in the previous BHLS [25], $\lambda_{\rho}(s)-\lambda_{\omega}(s)$ exhibits a similar $s$-behavior. $\alpha(s=0)=0$, which is no longer necessary in the present version of BHLS. A more natural way to proceed – giving the same results – is, however, to take into account the $\pi^{0}\gamma$ loop which introduces a small imaginary part to the eigenvalue differences, for $s_{0}\geq m_{\pi^{0}}^{2}$ upwards. For this purpose, we use the loop expression derived in [44] (see its Appendix C) which can be written ($s_{0}=m_{\pi^{0}}^{2}$) :

[TABLE]

Along the real axis in the complex $s$-plane, one has :

[TABLE]

The derivation of the $g_{V\pi^{0}\gamma}^{R}$ couplings is given in Subsection 12.2. The constant term inside the bracket in Equation (44) ensures that $\epsilon_{\pi^{0}}(0)=0$. Full consistency would impose to add (using obvious notations) :

[TABLE]

to $\delta M^{2}(s)$ as expressed above; however, because of the $e^{2}$ factor, they are higher order and can be neglected, except in the diagonal where they avoid possible real poles in the physical region. Additionally, the $\epsilon_{\eta}(s)$ and $\epsilon_{{\eta^{\prime}}}(s)$ contributions displayed for completeness will also be discarded.

In order to substantiate the effect of the $\pi^{0}\gamma$ loops, let us anticipate on the global fit information. The main effect of these loops is to prevent naturally the occurrence of poles on the physical region. This effect is prominent for the ”angle” $\alpha(s)$ which produces a $\rho^{0}-\omega$ mixing; Figure 1 displays the behavior of Re$[\alpha(s)]$ for the BHLS2 BS and RS variants and also, for comparison, its behavior within the former BHLS, where a $i\epsilon$ trick prevents a real pole close to the 2-pion threshold [25]. Even if the overall lineshapes are similar, one observes significant differences between the various solutions; for instance, $\alpha(s)$ vanishes at $s=0$ within the two variants of BHLS2, as the denominator is non-zero at the chiral point, allowing for a smooth connection between the spacelike and timelike regions.

Over the range $[-3.5\div 1.2]$ GeV2, the behavior observed for the two other mixing angles $\beta(s)$ and $\gamma(s)$ is much smoother (see also [23]); nevertheless, specific parametrizations may lead to a real pole131313One should note that a pole at such a location is exhibited by the usual Gounaris-Sakurai formula [45]. for $\beta(s)$ at negative $s$ below $\simeq-1.5$ GeV2. As noted above, such poles can be eliminated by means of usual dispersive methods.

11 The Model Pion and Kaon Form Factors

Using the non-anomalous sectors of the HLS model, broken as pointed out in the previous Sections (see Appendix A), one can already derive some of the form factors expressions needed for our global fit framework.

11.1 The Pion Form Factor in the $\tau$ Decay

The $\tau$ sector of BHLS2 is almost identical to that of BHLS [25]; indeed the main ${\cal L}_{\tau}$ piece reminded in Appendix A.2 is unchanged, but the newly introduced ${\cal L}_{z_{3}}$ contributes. This yields :

[TABLE]

with :

[TABLE]

where $z_{3}$ is a new HLS parameter [18] introduced by the ${\cal O}(p^{4})$ terms of the HLS Lagrangian (see Appendix A.3). $\Pi_{\rho\rho}^{\tau}(s)$ is the loop correction to the $\rho^{\pm}$ propagator and $\Pi_{W}(s)$ the loop correction to the $W^{\pm}-\rho^{\pm}$ transition amplitude; both are defined just below.

Using the following short-hand notations :

[TABLE]

one first defines the pion and kaon loop contributions to the $\rho^{\pm}$ self-energy $\Pi_{\rho\rho}^{\tau}(s)$ :

[TABLE]

where $\ell_{\pi}^{\tau}(s)$ and $\ell_{K}^{\tau}(s)$ denote the amputated loop functions for resp. $\pi^{\pm}\pi^{0}$ and $K^{\pm}K^{0}$, each having absorbed only a factor of resp. $G_{\pi}^{2}$ and $G_{K}^{2}$; in this way, only the breaking parameters affecting the $V$ sector of BHLS2 are factored out and displayed. $P^{\tau}_{\pi}(s)$ and $P^{\tau}_{K}(s)$ are subtraction polynomials chosen to vanish at $s=0$ because of current conservation. Then one has141414Here and in the following, one can be led to keep breaking parameter dependencies higher than first order to simplify the expression writing. :

[TABLE]

As the loop functions vanish at $s=0$, one clearly has $F_{\pi}^{\tau}(0)=1$. One should note that the same subtraction polynomials occur in the $W-\rho$ transition amplitude $F_{\rho}^{\tau}(s)$ and in the inverse propagator $D_{\rho}(s)$ but differently weighted.

The expression for the loop functions $\ell_{\pi}^{\tau}(s)$ and $\ell_{K}^{\tau}(s)$ along the real $s$-axis can be found in Subsection F.2 of the Appendix to [29]; the subtraction polynomial needed to ensure that they vanish at $s=0$ has coefficients given in Equations (122) of this reference. The coefficients for the constant term and for the second degree one given therein are correct, but the coefficient for the term linear in $s$ should be corrected to :

[TABLE]

using the notations151515$m_{0}^{2}$ and $m_{c}^{2}$ are the $s$ values of the direct and crossed thresholds energies. of [29].

11.2 Parametrization of Loops in $e^{+}e^{-}$ Annihilations and $\tau$ Decays

The $\rho^{0}$ self-energy and the $\gamma\Rightarrow(\rho^{0}/\omega/\phi)$ transition amplitudes involve important loop corrections; these already play a central role in the dynamical mixing of vector mesons as shown in Section 10. Let us define a parametrization common to the computation of the ”angles” $\alpha(s)$, $\beta(s)$ and $\gamma(s)$, to the self-energies and to the $A\Rightarrow V$ transition amplitudes as these are involved in all form factors addressed in our global fitting code.

One first defines the following loops  :

[TABLE]

where $\ell_{\pi}^{e}(s)$, $\ell_{K^{\pm}}^{e}(s)$ and $\ell_{K^{0}}^{e}(s)$ are resp. the $\pi^{+}\pi^{-}$, $K^{+}K^{-}$ and $K^{0}\overline{K}^{0}$ loops. These are defined as the amputated loops, including only a factor of $G_{\pi}^{2}$ for the former loop and of $G_{K}^{2}$ for the two others (see Equations (49)). The expression for these loops along the real $s$-axis can be found in Section F.1 of the Appendix to [29].

The loop parts being defined this way, $P_{\pi}^{e}(s)$, $P^{e}_{K^{\pm}}(s)$ and $P^{e}_{K^{0}}(s)$ are the subtraction polynomials, chosen of the second degree and vanishing at $s=0$.

However, as already observed in Section 10, besides the pion loop, it is rather the sum and difference of the kaon loops which are relevant :

[TABLE]

where $P^{e}_{S}(s)$ and $P^{e}_{D}(s)$ are denoting their respective subtraction polynomials. As already noted, one can consider $\Pi_{K_{\rm diff}}^{e}(s)$ and $P^{e}_{D}(s)$ to be first order in breaking (${\cal O}(\epsilon)$). So, it is certainly more appropriate to parametrize $P^{e}_{S}(s)$ and $P^{e}_{D}(s)$ for fits and propagate them into the subtraction polynomials of $\Pi_{K^{\pm}}^{e}(s)$ and $\Pi_{K^{0}}^{e}(s)$ by stating :

[TABLE]

keeping in mind that $P^{e}_{D}(s)$ is ${\cal O}(\epsilon)$, while $P^{e}_{S}(s)$ is ${\cal O}(1)$.

On the other hand, as one certainly has :

[TABLE]

and it is appropriate to impose :

[TABLE]

and submit to fit the coefficients of the $\delta P^{\tau}(s)$ polynomial rather than those for $P_{\pi}^{\tau}(s)$ directly. This turns out to attribute the full breaking in $e^{+}e^{-}$ versus $\tau$ to the pion loop subtraction term; as our fitting region is mostly below the $K\overline{K}$ thresholds, this looks a safe assumption. The second relation, additionally, is an assumption which allows to reduce the number of free parameters in the minimization procedure by two units without any degradation of the fit quality.

As a general statement, all our subtraction polynomials have been chosen of second degree and vanishing at $s=0$. An exception is made, however, by choosing the third degree for $\delta P^{\tau}(s)$ which is found to provide a significantly better description of the $\tau$ spectra, especially for the Belle dipion spectrum [46].

11.3 The Pion Form Factor in the $e^{+}e^{-}$ Annihilations

Because of the vector meson mixing, deriving the pion form factor in the $e^{+}e^{-}$ annihilation is slightly more involved than in the $\tau$ decay. To start, one needs to propagate the transformation in Equation (43) into the Lagrangian in Appendix A.1 and collect all contributions to the fully renormalized $\rho$, $\omega$ and $\phi$ field couplings. The derived pion form factor in $e^{+}e^{-}$ annihilations $F_{\pi}^{e}(s)$ writes :

[TABLE]

Using the $\rho^{I}\pi^{+}\pi^{-}$ coupling in the Lagrangian displayed in Appendix A.1 and the $I\Rightarrow R$ transform given by Equation (43), one readily gets :

[TABLE]

keeping only the leading ${\cal O}(\epsilon)$ terms in breaking. So, BHLS2 generates $\omega$ and $\phi$ direct couplings to $\pi\pi$. The inverse $\rho^{0}$ propagator writes :

[TABLE]

where $m_{\rho_{0}}^{2}$ is defined in Equations (26) and $\Pi_{\rho\rho}^{e}(s)$ is the self-energy of the physical (i.e. fully renormalized) $\rho^{0}$. This can be defined in terms of the loop functions constructed in the Subsection just above. $\Pi_{\rho\rho}^{e}(s)$ is made up of two pieces :

[TABLE]

where the couplings $\widetilde{g}^{\rho}_{K\pm}$ and $\widetilde{g}^{\rho}_{K0}$ are given in Appendix A.4. A term$\left[\left[\widetilde{g}^{\rho}_{K\pm}\right]^{2}-\left[\widetilde{g}^{\rho}_{K0}\right]^{2}\right]P^{e}_{D}(s)/2$ has been omitted in the expression for $\Pi_{K}^{e}(s)$ as being ${\cal O}(\epsilon^{2})$.

As for the inverse propagators $D_{\omega}(s)$ and $D_{\phi}(s)$, one uses the expression given by Equation (29) with the appropriate values for HK masses $m_{V}^{2}$, and having ([$\widetilde{m}_{V}$, $\widetilde{\Gamma}_{V}$], $V=\omega,~{}\phi$) as parameters to be determined by fit.

The loop corrected $V-\gamma$ transitions amplitudes $F_{V\gamma}^{e}(s)$ are given by :

[TABLE]

where the new-comer $z_{3}s$ terms occur as in $F_{\rho}^{\tau}(s)$ above. Using together the ideal couplings Equations (A.3), the renormalization matrix Equation (43) and the RS conditions (see Section 8), one derives :

[TABLE]

which become $s$-dependent thanks to the dynamical vector mixing. Similarly, using the couplings which can be read off ${\cal L}_{z_{3}}$ in Equation (A.3) and the transformation matrix given by Equation (43), one also derives :

[TABLE]

neglecting terms of breaking order greater than 1. Finally, the relevant loop transition terms $\Pi_{V\gamma}(s)$ are collected in Appendix A.5.

11.4 The Charged and Neutral Kaon Form Factors

The charged and neutral kaon electromagnetic form-factors can also be derived from the BHLS2 Lagrangian given in Appendix A.1 :

[TABLE]

where the $\gamma-V$ transition amplitudes $F_{V\gamma}$ and the propagators have been already defined and the $VK\overline{K}$ couplings are displayed in the Appendix A.4.

The experimental spectra measured in the spacelike region, in particular by the NA7 Collaboration [47, 48], being squared form factors, the modulus squared of the form factors given in Equation (57) and in Equations (11.4) apply directly. In the timelike region, experimentalists rather publish cross-sections :

[TABLE]

where $q_{P}=\sqrt{s-4m_{P}^{2}}/2$ is the meson momentum in the center-of-mass system ($P=K^{+},K^{0}$). In the case of the $K^{+}K^{-}$ final state, the cross-section should take into account the significant interaction between the charged kaons emerging from the $e^{+}e^{-}$ annihilation. Traditionally, this is done by multiplying the $e^{+}e^{-}\rightarrow K^{+}K^{-}$ cross-section by the leading term describing the Coulomb interaction as formulated by [49, 50] :

[TABLE]

As the data have become more and more precise, it looks more convenient to use the exponentiated expression derived in [51] which includes also the Final State Radiation (FSR) correction factor; indeed, the $K^{+}K^{-}$ spectra provided by the various experiments are not amputated from FSR effects.

11.5 Form Factors in BHLS versus BHLS2

It looks worth pointing out the difference between the form factor expressions in BHLS and in BHLS2. As already noted, the ${\cal L}_{A}$ Lagrangian is common to both frameworks. In contrast, the ${\cal L}_{V}$ Lagrangian is significantly different. Obvious differences for ${\cal L}_{V}$ have already been noted. For instance, within BHLS2, $\Delta_{V}=0$ and its sharing parameter $h_{V}$ drop out whereas they are a key ingredient in the diagonalization procedure of BHLS.

Three parameters161616Additional constraints arising in the diagonalization procedure of the vector mass term, have reduced this number to two as $\xi_{0}=\xi_{8}$ is requested. arise from the Covariant Derivative breaking ($\xi_{3}$, $\xi_{0}$, $\xi_{8}$). Moreover, a key role is played by the $\psi$ rotation angles within the RS solution of BHLS2. None of these parameters occur in the BHLS framework of [25].

The effect of a significant $\xi_{3}$ value is obviously equivalent to shifting apart the Higgs-Kibble masses of both $\rho$’s and their couplings to $\pi\pi$ in a correlated way. As this kind of exercise performed within BHLS concluded to a small effect [25, 26], one expects a small value for $\xi_{3}$. Also, BHLS2 deals with different subtraction polynomials – the $\delta P^{\tau}(s)$ function already referred to above – for the $\pi^{\pm}\pi^{0}$ and $\pi^{+}\pi^{-}$ loops, a new (loop) mechanism able alone to contribute to physical mass and width differences between the $\rho^{0}$ and $\rho^{\pm}$ mesons.

One should also note another difference : In the present framework, as already noted in Subsection 11.1, the pion and kaon loops entering the amplitudes $F_{\rho}^{\tau}(s)$ and $F_{\rho}^{e}(s)$ carry fitted subtraction polynomials already involved in the charged and neutral $\rho$ inverse propagators $D_{\rho}(s)$. In BHLS, one allowed the subtraction polynomial in $F_{\rho}^{e/\tau}(s)$ to carry a piece independent of both $D_{\rho}(s)$’s. Because of the ${\cal O}(p^{4})$ term still introduced, it means that a second-degree term in $s$ has been removed from $F_{\rho}^{e/\tau}(s)$ in the BHLS2 framework. We should see below the effect of this on the $\pi\pi$ $P$-wave phase-shift above the $\phi$ mass.

12 The Anomalous Sector of BHLS2

As in [25], the anomalous FKTUY Lagrangian [22, 18] is used to address the $e^{+}e^{-}\rightarrow(\pi^{0}/\eta)\gamma$ and $e^{+}e^{-}\rightarrow\pi^{0}\pi^{+}\pi^{-}$ processes, and also the radiative decays with couplings of the form $VP\gamma$ and $P\gamma\gamma$ which are parts of our global fitting framework. One can display its various pieces :

[TABLE]

where, temporarily, the weight of each piece has been extracted out to exhibit its dependence upon the unconstrained $c_{i}$ constants of the anomalous HLS Lagrangian [22]. If one imposes to yield the Wess-Zumino-Witten terms [52, 53] at the chiral point – the so-called triangle and box anomalies – the condition $c_{3}=c_{4}$ is known to be mandatory [21, 18]. Reference [25] having also shown that this condition is indeed favored by fits to the annihilation data, one endorses the $c_{3}=c_{4}$ constraint and, thus, cancels out the last term in Equation (66). Hence, the $VP\gamma$ couplings become entirely generated by the combination of $VVP$ couplings with the $V-\gamma$ transitions provided by the non-anomalous HLS Lagrangian.

The surviving Lagrangian pieces displayed in Equation (66) are given by :

[TABLE]

where $N_{c}$ is the number of colors and the $c_{i}$ factors have been reabsorbed in the normalization of the various Lagrangian pieces. $P$ and $V$ denote the bare pseudoscalar and vector field matrices. In the present context, one has, however :

[TABLE]

where $V_{\mu}^{I}$ is the (usual) $U(3)$ symmetric vector field matrix and $\delta V_{\mu}=\delta V_{\mu}^{3}+\delta V_{\mu}^{0}+\delta V_{\mu}^{8}$ which has been expressed in term of the ideal field combinations in Equations (12).

For the reader’s convenience, Appendix B reminds briefly the renormalization procedure which leads from the bare PS fields to their renormalized partners [25].

12.1 Renormalization of the $VVP$ Couplings

As displayed in Equation (39), referring to :

[TABLE]

the relation between the fully renormalized vector fields (denoted $R$ here and from now on) and their ideal combination is given by :

[TABLE]

where the product $R(U_{3})\cdot R(Loop)$, given by Equation (43), will be denoted from now on :

[TABLE]

Let us point out that ${\cal R}(s)$ is a real analytic $3\times 3$ matrix function and fulfills at order ${\cal O}(\epsilon)$ :

[TABLE]

along the right-hand cut on the physical sheet [29]. The $VVP$ Lagrangian displayed in Appendix C can be symbolically written :

[TABLE]

in terms of ideal field combinations, $P_{0}$ being any of the $\pi^{0}$, $\eta$ or ${\eta^{\prime}}$ fields; summation over greek and latin indices is understood. The 3-vector $G^{I}(\pi^{\pm})$ and the 3$\times$3 matrix $G^{I}(P_{0})$ can easily be constructed using the relevant pieces of information given in Appendix C. One thus has :

[TABLE]

with the obvious $(\rho,~{}\omega,~{}\phi)$ component indexing. On the other hand, one has defined :

[TABLE]

where, for each $P_{0}$, $\widetilde{g}_{ij}$ is the coupling of the form $P_{0}{\cal V}_{i}{\cal V}_{j}$ which can be read off the relevant expressions given in Appendix C.

Performing the replacement ${\cal V}^{I}={\cal R}(s){\cal V}^{R}$ in Equation (71), one can derive the fully renormalized $VVP$ Lagrangian :

[TABLE]

where, the coupling vector $G^{R}(\pi^{\pm})$ and the coupling matrix $G^{R}(P_{0})$ inherit from ${\cal R}$ a $s$-dependence; they are related to their ideal partners by :

[TABLE]

which is a concise way to express the renormalized couplings.

12.2 The Renormalized AVP Effective Lagrangian

As already noted, setting $c_{3}=c_{4}$ turns out to cancel the direct FKTUY $AVP$ Lagrangian piece; hence, all $AVP$ couplings inside BHLS2 are generated by $VVP$ couplings followed by one $V\rightarrow\gamma$ transition. Let us consider the $V-\gamma$ transition term of the non-anomalous BHLS2 Lagrangian (see Equation (114)). Having defined the transposed vector $\widetilde{f_{V\gamma}^{I}}=(f_{\rho\gamma}^{I},~{}f_{\omega\gamma}^{I},~{}f_{\phi\gamma}^{I})$, one can rewrite it :

[TABLE]

which defines the 3-vector $f_{V\gamma}^{R}$ as :

[TABLE]

The effective ${\cal L}_{AVP}$ Lagrangian can be derived from Equation (74) by replacing one neutral vector meson, say ${\cal V}_{i}^{R}$, using the rule :

[TABLE]

where $m_{i}^{2}$ is the renormalized squared mass of the ${\cal V}_{i}$ vector meson as given in Equations (26); let us note that each $H_{i~{}\gamma}$ carries a hidden factor of $1/g$. Starting with the charged pion part of Equation (74), one gets :

[TABLE]

where $H_{\gamma V}$ is the 3-vector constructed from the $H_{i~{}\gamma}$ defined just above :

[TABLE]

For the $P_{0}$ part of the Lagrangian Equation (74), the algebra is slightly more involved and gives :

[TABLE]

with :

[TABLE]

This effective Lagrangian is the tool to parametrize the $e^{+}e^{-}\rightarrow(\pi^{0}/\eta)\gamma$ cross-sections and the $VP\gamma$ radiative decays.

12.3 The $e^{+}e^{-}\rightarrow(\pi^{0}/\eta)\gamma$ Cross-Sections

Using the $AAP$ and $AVP$ Lagrangians derived above in terms of renormalized vector and pseudoscalar fields, one can construct the amplitudes for $\gamma^{*}\rightarrow(\pi^{0}/\eta)\gamma$ transitions. These can be written :

[TABLE]

where the (co)vector $\widetilde{H}_{\gamma V}$ is defined in Equation (80), the $G^{R}(P_{0})$ matrix in Equation (75) (and Equation (73)). One has also defined :

[TABLE]

the functions $F_{\gamma V}^{R}(s)$ and $D_{V}(s)$ having been defined in Subsection 11.3. Finally, one has :

[TABLE]

with $g_{P_{0}\gamma\gamma}$, derived in [25], and reminded in Appendix D.

In the case of the $e^{+}e^{-}\rightarrow(\pi^{0}/\eta)\gamma$ annihilations, the available data are always cross-sections which are related to $F_{P_{0}\gamma}(s)$ just defined by :

[TABLE]

12.4 The $e^{+}e^{-}\rightarrow\pi^{0}\pi^{+}\pi^{-}$ Cross-section

The amplitude for the $\gamma^{*}\rightarrow\pi^{0}\pi^{+}\pi^{-}$ transition involves most of the FKTUY Lagrangian pieces; it can be written :

[TABLE]

labeling each term by the particular piece of the FKTUY Lagrangian from which it originates. As already noted, because $c_{3}=c_{4}$ is assumed, there is no $T_{AVP}$ piece.

The $T_{APPP}$ contribution to the full $T(\gamma^{*}\rightarrow\pi^{0}\pi^{+}\pi^{-})$ amplitude can be derived from the information given in Appendix D; it writes :

[TABLE]

where one has defined :

[TABLE]

Three pieces are coming from the $VPPP$ Lagrangian also displayed in Appendix E; they are collected in :

[TABLE]

where the renormalized vector couplings $g_{V\pi}^{R}(s)$ to 3 pions have to be derived from using Equations (145) and Equation (148), ${\cal R}(s)$ being given by Equation (43). The $V-\gamma$ amplitudes $F^{e}_{V\gamma}(s)$ have been constructed in Subsection 11.3 and the inverse propagators in Subsection 11.3 for the $\rho^{0}$ meson. The inverse propagators for the $\omega$ and $\phi$ mesons have been discussed and defined in Section 7. One has also defined :

[TABLE]

The $VVP$ Lagrangian piece in Equation (137) has been rewritten in terms of renormalized fields in Equation (74) with (renormalized) couplings derivable by means of Equation (75). The simplest way to write down $T(\gamma^{*}\rightarrow\pi^{0}\pi^{+}\pi^{-})$ in a way easy to code within our global fit procedure is displayed just below.

One first defines the $H_{i}(s)$ functions :

[TABLE]

where $s$ is the incoming squared energy and the $s_{ij}$’s indicate the invariant mass squared of the corresponding outgoing $(i,j)$ pairs; the primed mixing angles have been defined in Subsection 10.2. $T_{VVP}$ depends on the 3 functions $(H_{i}(s),~{}i=1~{}\cdots 3$ with $s$-dependent coefficients $F_{i}(s)$ given below.

Collecting all terms, the full amplitude writes :

[TABLE]

with :

[TABLE]

In this way, to write down the full amplitude, the various $F_{i}(s)$ functions only depend on the incoming off-shell photon energy squared $s$; the dependence upon the various sub-energies $s_{ij}$ is, instead, only carried by the $H_{i}(s)$ functions as clear from Equations (92). One has :

[TABLE]

where $\widetilde{\alpha}(s)=\alpha(s)+\psi_{\omega}$ and $\widetilde{\gamma}(s)=\gamma(s)+\psi_{0}$.

One should note, as in BHLS [25] and earlier in [23], that all parameters already fitted with the (five) $e^{+}e^{-}\rightarrow(\pi\pi/K\overline{K}/\pi^{0}\gamma/\eta\gamma)$ processes fully determine all the $(F_{i}(s),~{}i=1~{}\cdots 3)$. The only new parameter coming with $e^{+}e^{-}\rightarrow\pi^{0}\pi^{+}\pi^{-}$ is $c_{1}-c_{2}$, only affecting $F_{0}(s)$; for practical purposes, it is handled in a specific manner, as in the previous versions of the (broken) HLS model just quoted.

Indeed, a global fit to all cross-sections but $e^{+}e^{-}\rightarrow\pi^{0}\pi^{+}\pi^{-}$ allows to yield the relevant parameters with a good approximation; thus, having at hand all ingredients defining the $F_{i}(s)$’s, a first minimization run including the $e^{+}e^{-}\rightarrow\pi^{0}\pi^{+}\pi^{-}$ cross-section can be performed to also derive a first estimate for $c_{1}-c_{2}$. The output of this minuit minimization run is then used as input for a next minimization step. This initiates an iteration procedure involving all cross-sections which is carried on up to convergence – when some criterion is met.

This method has been proven to converge in a very few minimization steps [25]. What makes such a minimization procedure unavoidable is that the $e^{+}e^{-}\rightarrow\pi^{0}\pi^{+}\pi^{-}$ cross-section expression implies to integrate over the $s_{ij}$’s at each step. This is obviously prohibitively computer time consuming for a negligible gain. Hence, at each step, one starts by tabulating coefficient functions, exhibited in the next expression between curly brackets :

[TABLE]

and these tables are used all along this step. Equation (12.4) uses the Kuraev-Silagadze parametrization [54] and its kernel $G(x,y)$ function; these are reminded in Appendix H of [25].

13 The Data Samples Submitted to Global Fits

Basically, 48 experimental data samples are presently available which enter the scope of the HLS Model. Relying on the global fits performed in [25, 26], one has been led to discard a few of them from our fitting procedure; this situation happens again within the present framework. Before going on, let us remind the most important part of the data samples and list the newcomers. Newcomers encompass either newly collected samples or existing form factor spectra covering spacelike momenta not treated within the former BHLS framework [25].

Substantially, the set of available data samples includes :

• •

i/ For the $\pi^{+}\pi^{-}$ annihilation channel :

The earliest data samples date back to the beginning of the eighties [55, 56]. These have been followed by quite precise data samples collected in scan mode by CMD-2 [57, 58, 59] and SND [60] on the VEPP-2M collider at Novosibirsk. The published spectra cover the energy region from 370 MeV to 970 MeV; presently, there is still no published spectrum covering the $\phi(1020)$ energy region.

These measurements, referred to globally as NSK, have been followed by large statistics data collected using the so-called initial state radiation (ISR) mechanism (see [61], for instance). In this way, three data samples [62, 63, 64] have been collected with the KLOE detector running on the DAPHNE collider up to $\simeq 1$ GeV, just below the $\phi$ meson mass. In the same period of time, another $\pi^{+}\pi^{-}$ sample has been collected by BaBar [65, 66] covering the energy range up to 3 GeV. By 2015, BESIII has also published a $\pi^{+}\pi^{-}$ sample [67] and, by end of 2017, a group from the CLEO-c Collaboration has published the latest to date $\pi^{+}\pi^{-}$ sample [68] collected with the CESR detector.

As a whole, the most precise data samples represent now 8 spectra. One has been led to discard from our global framework KLOE08 [62] and BaBar [65, 66] for different reasons171717As for the BaBar spectrum, a recent comprehensive study [69] analyzing the $\omega\rightarrow\pi^{+}\pi^{-}$ branching ratio reaches conclusions similar to ours. One also has to keep in mind that the BaBar $\pi\pi$ spectrum is in conflict with the $\tau$ data corrected for the isospin breakings. This is what we learn from Figure 1 of Ref. [27]. discussed in [26, 27]. A reanalysis unifying the three KLOE data samples has been recently performed [70, 71] and an 85 data-point spectrum has been derived which will be commented at the appropriate place. Moreover, the full error covariance matrix including the correlations between KLOE08, KLOE10 and KLOE12 is now available; this is used here to account for the (weak) correlations between KLOE10 and KLOE12.

• •

ii/ For the $K^{+}K^{-}$ annihilation channel :

Up to very recently, the available data samples covering this channel were the three scans collected by CMD-2 [72, 73] and the two scans of SND [74]. They were included within the global fits performed in the (previous) BHLS framework (see [25, 27]). Quite recently, the CMD-3 Collaboration has found that the two CMD-2 scans [73] should undergo a rescaling by [75] $1.094\pm 0.04$ which represents an important correction of the spectrum absolute normalization and of its (correlated) systematics ($2.2\%\rightarrow 4.45\%$). The influence of this correction on the physics information previously derived within BHLS in [25, 26, 27] imposes the quite significant update performed in the present study.

BaBar has published in 2013 the first measurement [76] of this cross-section performed in the ISR mode. This high statistics spectrum covers the full energy range from threshold to 3 GeV. Moreover, very recently, CMD-3 has also produced a new high statistics measurement of this spectrum in scan mode [77]. However, as emphasized in the CMD-3 publication, this measurement is inconsistent with the BaBar data sample [76]. This issue is addressed below.

• •

iii/ For the $K_{L}K_{S}$ annihilation channel :

To our knowledge, the first reported measurement of the $e^{+}e^{-}\rightarrow K_{L}K_{S}$ cross section has been provided by CMD-2 in 1995 [72] and has been followed in 2004 by the 4 scans reported in [78]. In the meantime SND also produced 4 data samples, 2 in the charged decay mode of the $K_{S}$ meson, 2 in its neutral mode [74]. These were the data encompassed in our former studies. The present work also includes the CMD-3 data sample [79] recently published.

• •

iv/ For the $(\pi^{0}/\eta)\gamma$ annihilation channels :

As for the $e^{+}e^{-}\rightarrow\pi^{0}\gamma$ cross section, in our energy range, the available data samples are scarce : One sample has been provided by CMD-2 [80] and another one by SND [81] – which supersedes the former [82] used in our previous studies[25, 26, 27]; the older SND $\pi^{0}\gamma$ spectrum given in [83] is also considered.

Regarding the $e^{+}e^{-}\rightarrow\eta\gamma$ cross section, the situation is more favorable as six data samples corresponding to different $\eta$ decay modes have been collected by CMD-2 [84, 85, 80] and SND [83, 86]. Nevertheless, no newly collected data sample for this channel has been reported.

• •

v/ For the $\pi^{+}\pi^{-}\pi^{0}$ annihilation channel :

As no newly collected sample covering the energy region up to the $\phi$ mass has been reported, we are left with only the samples already examined in our previous analyses [25, 26, 27]. These have been collected by CMD-2 [72, 78, 87, 88] and SND [89, 90] and cover either the $\omega$ peak region or the $\phi$ peak region. We also consider the 3-pion spectrum given in [91] as it mostly deals with the region in between the $\omega$ and $\phi$ peaks, allowing for a valuable constraint on the $\rho^{0}$ meson physics background.

However, preliminary fits showing that the sample in [87] exhibits an average $\chi^{2}$ per data point much above 2 – as also yielded by [92] in a different theoretical context – we have been led to discard it from our fits. Overall, one is left with nine data samples, including an old sample from the former CMD Collaboration.

• •

vi/ The pion and kaon form factors in the spacelike region :

One of the motivations to develop a new breaking procedure for the HLS Model is to design it in such a way that the resulting broken model (BHLS2) could encompass accurately the close spacelike region of the pion and kaon form factors. Indeed, a fair account of $F_{\pi}(s)$ in a region which stretches on both sides of $s=0$ is the best way to ascertain the low energy behavior of the $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}$ cross section which provides an enhanced contribution to $(g-2)_{\mu}$.

The NA7 Collaboration has measured the pion form factor [47] in the spacelike range $Q^{2}\in[-0.253,-0.015]$ GeV2 with good statistics and a small reported normalization error (0.9%); two years later, NA7 has also published a measurement of the charged kaon form factor [48] over the range $Q^{2}\in[-0.0950,-0.0175]$ GeV2 with also a small reported normalization uncertainty of 1%. On the other hand, an experiment at Fermilab had previously published less precise spectra on the pion and kaon form factors [93, 94] over resp. $Q^{2}\in[-0.092,-0.039]$ GeV2 and $Q^{2}\in[-0.1145,-0.0409]$ GeV2. Unlike former experiments181818 Reference [95] has collected the results of several of them. at higher virtualities and, more modern experiments191919Former references can be found in the quoted papers. [96, 97, 98, 99, 100, 101], Refs. [47, 93] – and [48, 94] as well – report on results derived in resp. $\pi$ and $K$ scattering on atomic electrons and are, thus, model independent.

As a matter of fact, one should also note that there is no $K^{0}$ form factor spectrum measurement in the spacelike region, which would be a real challenge; nevertheless, measurements of $<r_{K^{0}}^{2}>$ exist [102, 103] which report negative values.

A common feature of the data samples we have to deal with is that their systematic errors are quite generally reported as correlated; this implies that the corresponding error covariance matrices $W$ can be of large dimension (especially for the BaBar and KLOE samples). As the input to the $\chi^{2}$ expressions which feed our fitting code are just their inverses $W^{-1}$, the issue has been to ascertain the quality of these. For this purpose, we have chosen using the subroutines of the NAG library. Numerical checks have nevertheless been performed on the eigenvalues, the diagonal and non–diagonal elements of the products $W\cdot W^{-1}$ to verify that numerical departures from unit matrix properties were indeed negligible.

14 Comments on the Model Parameters

Before going on, it is worthwhile commenting on the fit parameters of the broken HLS Models. Basically, the number of parameters of the original (unbroken) HLS Model [21, 18, 22] to be fixed from data is small. As for its non-anomalous sector, besides $f_{\pi}$, these are the universal vector coupling $g$ and the specific HLS parameter $a$, generally found to slightly depart from 2, the value expected from standard VMD approaches. Its anomalous sector [22] involves a few more parameters already displayed in the Lagrangian pieces presented in Section 12 : $c_{1}-c_{2}$, $c_{3}$ and $c_{4}$. There are strong motivations [21, 18] to impose $c_{3}=c_{4}$. So, when unbroken, the HLS model depends on only four unconstrained parameters.

There are two parameter categories introduced in our approach: the coefficients of the subtraction polynomials which supplement the corresponding loops and, on the other hand, the breaking parameters inherently affecting the present Lagrangian – and those defined in [25].

Traditionally, as illustrated by the Gounaris-Sakurai formula [45], the coefficients of the subtraction polynomials are reexpressed in terms of the resonance mass and width, which are determined by a fit to the data. In a global approach like ours, the subtraction polynomials affect loops which come simultaneously in the various amplitudes associated with the various annihilation processes simultaneously under examination. In this case, as there is no reason to favor a specific channel to fix the meaning and the value of parameters which come in all channels, the strategy adopted is to let the data, as a whole, determine the subtraction parameter values via a global fit. Even if the parameter content may look less intuitive, this approach – also adopted in BHLS [25] – looks the most motivated. This does not prevent to extract usual physics quantities, relying on the fit results, as done since [24, 25], for instance.

In the second category of parameters, besides the parameters of the unbroken HLS model reminded just above, one finds all those introduced by the various aspects of the breaking procedure at work. As for the parameters generated by the BKY mechanism and reminded in Section 3, six in total, the condition $\Sigma_{A}=0$ is already stated since the original BHLS Model [25]. It has been shown in Section 7 that the diagonalization procedure of BHLS2 imposes $\Delta_{V}=0$. Therefore, four BKY parameters only ($\Sigma_{V}$, $\Delta_{A}$, $z_{A}$ and $z_{V}$) remain free.

In Sections 7 and 9, it was also shown that $\xi_{0}=\xi_{8}$ comes as a natural constraint on the parameters of the CD breaking defined in Section 4. Therefore the CD breaking involves two more parameters202020Rigorously speaking, as already noted, two more breaking parameters can be defined addressing the charged and neutral $K^{*}$ mesons. $\xi_{0}(=\xi_{8})$ and $\xi_{3}$. Finally, switching on the Primordial Mixing introduces three more parameters ($\psi_{0}$, $\psi_{\omega}$ and $\psi_{\phi}$) into the fit procedure.

15 The BHLS2 Global Fits : General Features

Various kinds of global fits have been performed using the data samples and channels listed in Section 13, with or without the spacelike data, with or without the $\tau$ data, fixing some of the model parameters (the $\psi_{\alpha}$’s angles and $z_{V}$) or leaving all of them free.

When mixing as many data samples carrying important – sometimes dominant – overall normalization uncertainties, the effect of the so-called Peelle’s Pertinent Puzzle (PPP) [104] cannot be ignored to avoid getting biased estimates for physics quantities. The use of $\chi^{2}$ minimization methods, first questioned, has been finally justified [105, 106]; however, when dealing with distributions or form factors, the pertinent method turns out to invoke the underlying (theoretical) model function which is just what fits are supposed to provide. Following [107], [27] has defined an iterative procedure proved to cancel out biases (see the Appendix in [27]).

It is worthwhile briefly sketching the conclusions of212121See, in particular, its Section 4. [27] concerning the issue just reminded : 1/ The normalization uncertainty (denoted $\sigma$, possibly $s$-dependent) of any spectrum $T$ is absorbed within the total covariance matrix contributing to the general $\chi^{2}$, 2/ The correction which affects the normalization of $T$ is derived from the data, the theoretical function provided by the fit, $\sigma$ and the statistical covariance matrix – which may absorb the uncorrelated systematics.

Let us make a few general statements about fit results and data before focusing on specific topics, mainly the $K\overline{K}$ and $\tau$ sectors and the spacelike region behavior of the pion and kaon form factors. The detailed properties will be emphasized in the next Sections.

• •

As for the behavior of the data samples within our global fits, one can state that the issues met with the dipion spectra from BaBar [65, 66] and KLOE08 [62] within the BHLS framework [25] are confirmed within the BS and RS variants of the BHLS2 framework studied here; they are also discarded in the present study.

Whatever the specific conditions of the various fits performed, the $\pi^{+}\pi^{-}$ data samples from NSK (i.e. CMD-2 and SND), KLOE10/12, BESIII [67] and CLEO-c [68] are as well described within the BHLS2 framework as they were already within BHLS [25]; more precisely, NSK and KLOE10/12 contribute an average $\chi^{2}$ per data point of $\simeq 1$, while the BESIII and CLEO-c samples yield $\chi^{2}/N_{\rm pts}\simeq 0.6$. Figure 2 displays the spectra and fits for the pion form factor in the $e^{+}e^{-}$ annihilation and in the $\tau$ decay.

• •

In order to account for the difficulty to fully address the $\pi^{+}\pi^{-}\pi^{0}$ channel in the $\phi$ mass region within BHLS, a so-called B model was defined which simply turns out to discard from fit this mass region for (only) this channel [25]. In contrast, BHLS2 does not meet any issue222222As for the discarded data sample [87], see the above Section 13 and also [92]. and achieves a fair description of the $\pi^{+}\pi^{-}\pi^{0}$ data up to, and including, the $\phi$ mass region. The spectra are shown in Figure 3.

• •

Finally, the $e^{+}e^{-}\rightarrow(\pi^{0}/\eta)\gamma$ spectra are also nicely fitted within BHLS2 – as well as in the former BHLS – as illustrated by Figure 4. Quite generally, one obtains $\chi^{2}/N_{\rm pts}\simeq 90/112$ and $\chi^{2}/N_{\rm pts}\simeq 120/182$ for resp. the $\pi^{0}\gamma$ and $\eta\gamma$ channels.

In order to perform global fits, one obviously should make an assumption on the energy calibration of each of the data samples considered. For the dipion spectra, one has gathered into the fit data coming from the VEPP-2M machine at Novosibirsk (CMD-2 and SND), data from DAPHNE (KLOE10 and KLOE12), from BESIII and from CESR. A significant mismatch of the relative energy calibrations of these would have produced a failure of their common fit because of the $\rho-\omega$ drop off region; this is clearly not seen neither in the $\chi^{2}$ contributions of the various samples, nor in Figure 2 (see the inset in its leftmost panel). Moreover, comparing both panels in this Figure excludes a significant energy calibration mismatch between $e^{+}e^{-}$ and $\tau$ data.

Moreover, the data collected by CMD-2 and SND around the $\omega$ and $\phi$ peaks in the 3-pion channel as well as in the $(\pi^{0}/\eta)\gamma$ final states indicate a good consistency between the energy calibrations of the various data samples collected in different runs distributed over several years. This is visible in the fit properties as well as in the Figures where the fit function and the various data are displayed.

16 The BHLS2 Global Fits : The $e^{+}e^{-}\rightarrow K\overline{K}$ Channels

The available data covering the $K^{+}K^{-}$ channel232323In the following, as already noted in Section 13, the two CMD-2 scans [73] are corrected as stated in [75]. As for the BHLS predictions in [25, 26, 27], they are updated correspondingly. come from SND, CMD-2 and BaBar. Recently CMD-3 has published a new $K^{+}K^{-}$ measurement [77]. As for the neutral channel, most data come from CMD-2 and SND and have been listed in Section 13. A new data sample covering our fitting range has also been produced by CMD-3 [79]; unfortunately, the published $K_{L}K_{S}$ spectrum from BaBar starts at 1.06 GeV, just above our fitting range [108].

16.1 Preliminary Remarks

As for the $K^{+}K^{-}$ channel : In its [77], the CMD-3 Collaboration points toward disagreements with both the BaBar [76] and the original CMD-2 data samples [73]. Once the correction factor ($1.094\pm 0.040$) for the trigger efficiency of CMD-2 is applied, the discrepancy between the CMD-3 and CMD-2 $K^{+}K^{-}$ samples looks much reduced [75] – possibly washed out. Anyway, a full examination of the various $e^{+}e^{-}\rightarrow K\overline{K}$ data samples242424As for the contribution to the muon HVP of the energy region from $s\in[4m_{K}^{2},1.06$ GeV2] to $a_{\mu}(K^{+}K^{-})$, BaBar obtains $18.64\pm 0.16_{\rm stat}\pm 0.13_{\rm syst}\pm 0.03_{VP}$ while CMD-3 gets $19.33\pm 0.040$ in units of 10*-10*, a difference of $\simeq 0.7\times 10^{-10}$ between the central values, a $3.25\sigma$ effect. is worth addressing.

On the other hand, preliminary BHLS2 analyses based on the CMD-3 and BaBar data, made us aware of a visible degradation of both $K\overline{K}$ channel data description above $1.025\div 1.030$ GeV; therefore, in the present work, one limits the fit region for the $K\overline{K}$ data to the $\sqrt{s}\in[2m_{K},1.025$ GeV] range. The explanation for this effect is unclear. It could be an intrinsic weakness of the BHLS/BHLS2 models252525The standard HLS framework [18], which underlies the BHLS/BHLS2 models, does not account for vector particles not belonging to the fundamental vector nonet. which, as it stands, does address a possible onset of the high mass vector meson contributions; it could also reflect an experimental issue. For the sake of consistency, we have chosen to impose the same fitting range also to the former CMD-2 and SND data samples in both the charged and neutral modes,

16.2 A Specific Issue of the $K\overline{K}$ Data

As already pointed out, one does not detect any indication for a mismatch in the energy calibration of the various data samples covering the 3-pion and $(\pi^{0}/\eta)\gamma$ channels, all collected with the CMD-2 and SND detectors on the same VEPP-2M collider. One also does not observe any mismatch between the energy calibration of the 3-pion and $(\pi^{0}/\eta)\gamma$ data and those collected by the same collaborations in the $K\overline{K}$ channels. Nevertheless262626E. Solodov, Budker Institute, Novosibirsk, private communication., CMD-2 and SND carry a (correlated) absolute energy uncertainty in the range of $30\div 40$ keV, while for CMD-3, the absolute energy uncertainty is $\simeq 60$ keV, statistically independent of that of CMD-2 & SND. Moreover, BaBar also reports a 50 keV correlated energy scale uncertainty[76]. When dealing with an object as narrow as the $\phi$ meson within a global framework, these uncertainties should be accounted for otherwise, the fits cannot succeed.

As the bulk of data samples for both $K\overline{K}$ channels has been provided by CMD-2 and SND, it looks quite natural to choose their common energy calibration as a reference; this will be denoted $E_{\rm NSK}$ in the following. For this purpose, one defines $E_{\rm CMD3}=E_{\rm NSK}+\delta E_{\rm CMD3}$ and $E_{\rm BaBar}=E_{\rm NSK}+\delta E_{\rm BaBar}$ to recalibrate the CMD-3 and BaBar data point energies in a fully correlated way and derive their $E_{\rm NSK}$ equivalent energy.

Therefore, in order to deal with the $K\overline{K}$ data, we are faced with a case where two kinds of correlated uncertainties have to be accounted for: Energy calibration and absolute normalization of the cross-sections. Hence, one merges in the minimization procedure the method developed in [27] and the two calibration parameters $\delta E_{\rm CMD3}$ and $\delta E_{\rm BaBar}$ to be also determined within the same fits.

For this purpose, one assumes that there is no functional relation between the normalization uncertainty $\lambda_{i}$ – of an experiment labeled $i$ – and the energy scale uncertainty $\delta E_{i}$; then, as $\partial\lambda_{i}/\partial[\delta E_{i}]=0$ for each experiment, substituting, as emphasized in [27], the value for $\lambda_{i}$ derived from $\partial\chi^{2}/\partial\lambda_{i}=0$ into the expression for $\chi^{2}$ is unchanged compared to [27].

16.3 Fitting the $K\overline{K}$ Channels : NSK, BaBar and CMD-3 Separately

In the following, the non-kaon sector data and channels used in our fits include, unless otherwise stated, the $\pi^{+}\pi^{-}$, $\pi^{+}\pi^{-}\pi^{0}$, $\pi^{0}\gamma$, $\eta\gamma$ annihilations, the $\tau$ decay to $\pi^{\pm}\pi^{0}$ and a few decay partial widths (see [25]). One starts by discussing the fits where the CMD-2 & SND (NSK), BaBar and CMD-3 kaon data are not influencing each other within the minimization procedure. Table 1 gathers most of our results in this case; these are commented right now.

• •

With regard to NSK : $<\chi^{2}>$, the average (partial) $\chi^{2}$ per NSK data points, is of order 1 in both the neutral and charged kaon channels; thanks to the fair description of the other channels involved in the fit, the global fit probability exceeds the 90% level.

As reminded above, the procedure defined in [27] allows us to derive from the fit procedure the normalization (absolute scale) correction to be applied to each spectrum272727One should stress that such a ”renormalization” is the natural outcome to the scale uncertainty problem and that the main advantage of any kind of global fit is to allow determining such a correction in consistency with a (large) set of data samples coming from different sources.. For the charged channel, one gets +5.0% for SND (expect. value 7.1%), a $0.7\sigma$ effect, and –1.6% for CMD-2 (updated expectation 4.4%), a $0.3\sigma$ effect; in the neutral channels, one correspondingly obtains –1.8% (expect. value 4.2%) for SND and +0.6% (expect. 1.7%) for CMD-2, resp. $0.4\sigma$ and $0.3\sigma$ effects.

Therefore the CMD-2 and SND data samples successfully fit the BHLS2 framework and, moreover, the derived scale corrections meet the reported scale uncertainty expectations.

• •

With regard to BaBar : One has chosen to supply the NSK neutral data to the fits in order to populate each dikaon channel as the BaBar $K_{L}K_{S}$ data [108] do not cover our fitting energy range. As the BaBar charged spectrum [76] has been collected in the ISR mode, it looks appropriate to compare, within the fit, each spectrum datum in a bin to an average value of the model spectrum over the bin width. The second data column in Table 1 shows that the BaBar sample yields $<\chi^{2}>=1.7$, somewhat large, and a global fit probability also above the 90% level.

The minimization procedure returns as energy shift vs NSK $\delta E_{\rm BaBar}=-125.8\pm 19.0$ keV; this looks fair when compared to the independent $30\div 40$ keV (NSK) and 50 keV (BaBar) energy calibration uncertainties. Moreover, the spectrum scale correction is at the $10^{-4}$ level – so that the absolute scale of the BaBar spectrum has not to be corrected. On the other hand, the scale corrections for the neutral NSK data included in the fit are –2.1% (SND) and +0.7% (CMD-2), almost identical to what has been obtained in the fit of the NSK data in isolation mode.

Interestingly one can also derive the scale correction for the charged data of SND (+5.2%) and CMD-2 (–3.7 %) even if they are not submitted to the fit; these predictions are clearly in the expected ballpark. This indicates, beforehand, that NSK and BaBar data in the kaon sectors are consistent with each other.

• •

As for CMD-3 : This Collaboration has provided data samples in both dikaon final states, When fitting these data, the correspondence is made – as with NSK – between each datum and the theoretical model computed at the nominal energy. As for NSK and BaBar, the systematic uncertainties are treated as fully correlated and serve to build up the systematic covariance matrices to be added up to the statistical ones; these are used in the fitting of the neutral and charged CMD-3 data within the global BHLS2 framework. The fit results are displayed in Table 1.

The $<\chi^{2}>$ for both CMD-3 data samples are clearly very large and the global fit probability poor. The scale corrections found are $+1.2$% (charged channel) and $-1.2$% (neutral channel) to be compared to the reported r.m.s., resp. 2.0% and 1.8%. The possibility of treating the bin information as in ISR experiments282828Comparing the datum information to the fit function averaged over the bin width, as just sketched for BaBar. has been examined – specifically, only in the neutral channel – and in view of the outcome, we gave up.

For the sake of completeness, we have nevertheless also performed the fit of the CMD-3 data in isolation by discarding the non-diagonal part of the total covariance matrix292929This turns out to consider that the systematic uncertainties are fully uncorrelated; some configurations intermediate between full correlation and no correlation can lead to similar results., as seemingly done in [76] for the BaBar data; the fit results change dramatically as one reaches $<\chi^{2}>=14/16$ and $<\chi^{2}>=17/17$ for resp. the charged and neutral modes of the CMD-3 data and with a fit probability exceeding 90%. Moreover, in this case (CMD-3 data in isolation) the corrections to the original normalization of the cross-sections become negligible ($\simeq 10^{-4}$) for both the charged and neutral channels303030The $<\chi^{2}>$’s  of the fit prediction to NSK data remain in reasonable correspondence with the NSK fit in isolation (55/49 and 105/92) while the distance of this prediction to BaBar data degrades more severely (76/27).. As for the energy shift vs NSK it moves from $-51.4\pm 13.7$ keV (see Table 1) to $-44.6\pm 14.8$ keV, still in the expected ballpark.

16.4 Fitting the $K\overline{K}$ Channels : NSK, BaBar and CMD-3 Combined

This Subsection is devoted to analyzing the consistency of the three available groups of data. Quite generally, the covariance matrices used within the fit procedure assume full correlation for the systematics, as in the previous paragraph. Nevertheless, one also briefly reports on additional information about a fit performed assuming the systematics uncorrelated for CMD-3 data. It is worthwhile proceeding with pairwise combinations.

• •

As clear from the first data column in Table 2 – and as could be inferred from the results already reported – one observes a fair consistency between NSK and BaBar samples within the global BHLS2 framework. Comparing the present NSK and BaBar contributions to the (total) $\chi^{2}$ with their analogs in their fits in isolation clearly shows that there is no significant tension between them. The only noticeable difference is the central value for $\delta E_{\rm BaBar}$ found closer to 100 keV by $1\sigma$.

Other pieces of information concerning the absolute spectrum normalization313131 For SND and CMD-2, as their corrections are computed run by run and sample by sample, here we only give average values, for simplicity. are worth to be mentioned. The fit outcome indicates that the BaBar data normalization should be downscaled by $\simeq 0.4\%$, the (updated) CMD-2 charged spectrum by 1% while the SND charged data normalization should be increased by 8.1%. These corrections are always consistent with the respective expectations reminded above.

• •

The second data column in Table 2 displays results from fits combining NSK and CMD-3. Here also one observes a striking resemblance between the individual partial $\chi^{2}$’s in the combined fit and in their fits in isolation, i.e. the NSK data samples go on yielding $<\chi^{2}>_{{\rm NSK}}\simeq 1$ while one gets $<\chi^{2}>_{{\rm CMD3}}>6$; the global fit probability is also poor.

Besides the energy shift $\delta E_{\rm CMD3}=-46.5\pm 6.7$ keV, the spectrum rescalings derived from the fit are interesting pieces of information. In this way, the fit returns an increase of the charged CMD-3 data normalization by only 0.2% while the neutral CMD-3 ones should be downscaled by 2.6% (a $1.5\sigma$ effect). One also finds that the SND charged data normalization should undergo an increase by 9.2%, while for CMD-2 the increase is only 1.5%. As for the neutral kaon normalizations, SND should be downscaled by 2.1%, while CMD-2 is increased by 0.5%. So, one does not observe major changes in the NSK normalizations.

• •

The third data column in Table 2 displays results from fits combining BaBar and CMD-3 kaon data, excluding their NSK analogs. The global fit probability clearly indicates that such a combination is not really favored; one should remark, nevertheless, that all properties exhibited by CMD-2 and SND are in fair correspondence with their analogs yielded in their fits in isolation; additionally, even the $\chi^{2}$ distances of the CMD-3+BaBar fit function to the NSK data (not under fit) – given within parentheses – are analog to their values derived in NSK alone or NSK+BaBar fits.

Quite interestingly, the tension observed between CMD-3 on the one hand, and BaBar and NSK on the other hand, is not manifest in plots, as soon as the energy calibration of CMD-3 vs NSK is applied. For instance Figure (5), shows the $K\overline{K}$ cross sections derived from the fit to NSK + BaBar; the corresponding data are rescaled and appropriately shifted as emphasized in the previous item; one has superimposed the (unfitted) data from CMD-3 shifted by $-45.9$ keV and rescaled as predicted by the NSK + BaBar fit : $-4.3$% and -2.4% for resp. the charged and neutral CMD-3 data samples.

On the other hand, one has found worth to also perform the fit sketched in the rightmost data column in Table 2, i.e. treating the NSK and BaBar uncertainties as fully correlated and both CMD-3 data samples as carrying uncorrelated uncertainties. Using the corrections to the absolute normalization of each sample as derived from this fit323232 This amounts to the following rescalings w.r.t. the original normalizations : –1.3% (CMD-3 $K^{+}K^{-}$), –0.5% (CMD-3 $K_{L}K_{S}$), +8.1% and –0.3% ($K^{+}K^{-}$ from resp. SND & CMD-2), +0.4% and –0.5% ($K_{L}K_{S}$ from resp. SND & CMD-2) and $\simeq-0.3\%$ for BaBar $K^{+}K^{-}$., one gets a figure which superimposes exactly to Figure (5), in the smallest detail.

This clearly indicates that the observed tension is essentially due to the non-diagonal parts of the CMD-3 error covariance matrices constructed from experimental information, as also done for the other data samples, CMD-2 in particular.

17 Fit Properties of the Various Broken HLS Models

All fits referred to in this Section involve the collection of data samples and channels already described and listed in Section 13; as for the $K\overline{K}$ channels, one limits oneself to using only the samples from CMD-2 (updated), SND and BaBar which do not show any obvious tension. The motivation for this choice is to avoid disentangling model dependence effects, on the one hand, and identified tensions among data samples333333The issue between BaBar and CMD-3 charged kaon data has already been noted in [77]. on the other hand. There is indeed no gain in increasing the statistics if this turns out to increase the systematics in an uncontrollable way.

The BHLS2 fits presented up to now have been derived using the so-called Reference Solution (RS); this solution combines the so-called Basic Solution (BS) defined in Section 7 and the Primordial Mixing defined in Section 8.

Nevertheless, in order to explore model dependence effects, analyzing also the fit results derived using BS is certainly relevant. On the other hand, for the same concern, we have rerun the BHLS model [27] using the SND together with the updated CMD-2 data; it then looked appropriate to also include the BaBar dikaon sample.

Therefore, with the BS and RS variants of BHLS2 and with BHLS, one actually has at disposal three approaches to the same data; as for BHLS, one concentrates on the so-called model B which circumvents the issue met at the $\phi$ mass in the three pion annihilation channel which is one of the motivations having led to BHLS2. The other one was to ascertain the description of the dipion threshold region and to improve the accountability of the spacelike region. The main fit properties are sketched in Table 3 and deserve some comments.

The leftmost data column – dealing with BHLS – shows the various $\chi^{2}$ for the displayed channels and the sum of the contributions of all data submitted to the fit. In all channels the average $\chi^{2}$ is of the order 1 or better343434One should note the smaller number of data points considered in the fit procedure for the $\pi^{+}\pi^{-}\pi^{0}$ annihilation channel, because of excluding the data points from the $\phi$ region. and the (global) fit probability is fairly good.

The rightmost data column displays the corresponding information derived with the RS variant of BHLS2. Here also, the fit quality is good and one yields a fair probability. One should note the quite favorable $<\chi^{2}>$’s for the pion and kaon form factors in the spacelike energy region in conjunction with good $\chi^{2}$’s in their timelike regions; the spacelike data samples – coming from NA7 and Fermilab – are discussed in detail in the next Section. As far as we know, the fits presented here are the first ones involving the pion and kaon form factors simultaneously in their timelike and spacelike regions with such quality.

The third data column displays the fit properties derived from the BS variant of BHLS2 using the same data as those just discussed for the RS variant. The fit probability is lower than usual in the HLS context. Scrutinizing the various items displayed in Table 3, one clearly observes a significant tension in the BS fit absent from the RS fit results when the $\tau$ data are involved. Indeed, when going from RS to BS $\chi^{2}_{\tau}$ increases by 6 units, but $\chi^{2}_{3\pi}$ increases by 12 and $\chi^{2}_{{\rm space}\pi}$ by 8; moreover353535The different behavior of the NSK and KLOE dipion data under BS global fits reflects the different magnitude of the uncertainties which makes the NSK dipion data – and likewise for those from BESIII and CESR – more permissive than their KLOE analogs., $\chi^{2}_{KLOE~{}2\pi}$ increases by 17 units whereas $\chi^{2}_{NSK~{}2\pi}$ is marginally affected (a 2 unit increase).

Instead, the BS fit performed when discarding the $\tau$ data leads to a picture in good correspondence with the RS fit results already commented; indeed, the partial $\chi^{2}$’s reported in the second data column of Table 3 are similar to those obtained by fitting with the RS variant of the BHLS2 model. Nevertheless, some tension survives in the 3-pion sector while the dipion spectra and the spacelike data almost recover their optimum fit quality; so, the BS variant of BHLS2 exhibits a quite acceptable behavior once the $\tau$ sector is discarded.

In conclusion, despite their different structures – and taking into account the peculiarities of the BS variant of BHLS2 – these models/variants open the possibility to examine model dependence effects in the evaluation of some reconstructed physical quantities; this may allow to compare these with the systematics generated by dealing with contradictory samples.

18 The Pion and Kaon Form Factors in the Spacelike Region

By allowing different HK masses for the $\rho^{0}$ and $\omega$ mesons, BHLS2 naturally permits a smooth connection between the spacelike and timelike sectors of the pion and kaon form factors. Moreover, applying in this context the vector meson dynamical mixing (see Section 10), does not break this smooth connection because $\alpha(s\!=\!0)\!=\!0$ becomes automatic and supplements the (already) vanishing properties of the other angles : $\beta(s\!=\!0)\!=\!\gamma(s\!=\!0)\!=\!0$.

18.1 The Pion and Kaon Form Factors in the Close Spacelike Region

Figure 6 displays the pion and kaon form factors in the spacelike region coming out from the fit. The spacelike data supplied to the global fit procedure are only the (model-independent) data samples collected by the NA7 [47, 48] and Fermilab experiments [93, 94]. Altogether, these experiments cover photon virtualities down to $\simeq-0.25$ GeV2 and $\simeq-0.12$ GeV2 for resp. the pion and kaon form factors.

The pion and kaon form factors derived from the global fit including the spacelike data are the red curves in Figure 6. Global fits have also been performed using the BS and RS variants of BHLS2, excluding the spacelike data. The relevant fit results obtained using, for instance, the BS variant are displayed in Table 4. They correspond to the green curve in Figure 6; as for the kaon form factor the predicted (green) curve superimposes exactly to the fitted (red) one.

Table 4 and Figure 6 obviously prove that the spacelike data are not a real constraint in the BHLS2 global fit approach; stated otherwise, the BHLS2 model supplied with only timelike data fairly well predicts the spacelike behaviors far inside the spacelike region for both the pion and kaon form factors. It is worthwhile to note that these analytic continuations are performed across energy squared distances which can be as large as 1 GeV2. For instance, the kaon form factor displayed here is the analytic continuation of the charged kaon form factor fitted in the $\phi$ mass region, i.e. the curve shown on the rightmost panel of Figure 5. Such a property is rare enough to be underlined.

Figure 6 gives the opportunity to comment on the representation of spectra affected by a normalization (scale) uncertainty; the NA7 spacelike data are a quite simple illustration of handling the rescaling effects. When a scale uncertainty affects a spectrum, there are obviously two possible ways to display the result : Either rescaling the data, keeping the fit/prediction untouched or performing the converse.

So, in Figure 6 one has chosen to display the NA7 data rescaled, keeping unmodified the fit (and prediction) function(s), the plotted uncertainties being the statistical ones. As for the Fermilab data samples, taking into account their low accuracy, no rescaling has been performed and the plotted uncertainties are just the reported total errors.

One can, for once, go into a few plotting details in simple wording. Following the method given in [27], the scale correction is given by :

[TABLE]

where (NA7) $\sigma=0.9\%$ or 1.0%, $f_{i}=|F_{P}(s_{i})|^{2}$ ($P=\pi^{\pm},~{}K^{\pm}$), $V$ is the statistical error covariance matrix and $m_{i}$ is the measured datum at $s=s_{i}$. The plotted NA7 data $m^{\prime}$ are related with the original ones $m$ by :

[TABLE]

Looking at Figure 6, the precise data (NA7 pion form factor) are indeed observed to tightly follow the fit function, once the normalization correction is applied.

It is relevant to provide the precise values for the $\lambda$’s derived from the fitted parameters. One thus finds :

[TABLE]

[TABLE]

which exhibit noticeable accuracies. Subsection 18.3 below – and Figure 8 – illustrates that the predicted scale correction for the pion form factor is supported by the LQCD data which should be intrinsically free of scale uncertainties.

18.2 The Pion Form Factor : The Model-Dependent Data

As one cannot measure the pion form factor in a model-independent way for photon virtualities $Q^{2}=-s>0.3$ GeV2, this physics is performed using the scattering process $pe\rightarrow e\pi^{+}n$, referred to in the literature as 1H$(e,e^{\prime}\pi^{+})n$. In this approach, two questions should be addressed :

i/ Separating out the longitudinal cross section $\sigma_{L}$,

ii/ Extracting $F_{\pi}$ from $\sigma_{L}$.

One can refer to [97, 98, 99] for comprehensive accounts of the issues. As for item i/, the method has evolved and improved in the course of time from the early experiments (also reported in [95]) to the more recent ones [100, 101, 96]. Item ii/ addresses more deeply the nuclear physics content of the process in order to account at best for the $\pi^{+}pn$ vertex and for the off-shell character of the intermediate $\pi^{+}$. Presently, the favored approach is the Regge Model developed by Vanderhaeghen, Guidal and Laget [109, 110]. The data extracted for $F_{\pi}(s<0)$ have thus improved in time but remain unavoidably model-dependent. These data samples have not been included in our fitting approach, but simply compared with the extrapolation of our fit function at $s<0$, a pure prediction below $s\simeq-0.25$ GeV2.

In Figure 6, besides the fitted data from NA7 and Fermilab, one has displayed the (unfitted) data from the DESY experiments [100, 101] and the lowest photon virtuality data collected at Jefferson Lab by the $F_{\pi}$ Collaboration. The data points flagged by JLAB 2008 and DESY have been extracted from the latest – to our knowledge – publication of the $F_{\pi}$ Collaboration [99]. One can observe that the updated [99] DESY data points fall exactly on the $F_{\pi}(s)$ function derived from a fit using also the NA7 [47] and Fermilab [93] data. The three JLab plotted data points are found systematically below the fit function; however, most of these JLab data points can also be found in a previous JLab publication [96]. The difference between the values in [99] and [96] seems mostly coming from using different extraction methods of $\sigma_{L}(Q^{2})$ [97]. The data points from [96] seem in closer agreement with the continuation of our $F_{\pi}(s)$ than the estimates given in [99]. They are also in better agreement with what can be expected from the DESY data points as updated in [99]; this updating seems to deal with the extraction of $F_{\pi}$ from the longitudinal cross section $\sigma_{L}$. The behavior of (the extrapolated) $F_{\pi}(s)$ beyond $-1$ GeV2 is shown in Figure 7 together with the other reported data.

The content of Figure 6 gives a hint that the extrapolation of the BHLS2 $F_{\pi}(s)$ could well be valid down to $s\simeq-1$ GeV2. This also reflects a $\rho$ dominance well inside the non-perturbative spacelike region. This dominance is illustrated by fit properties : When fitting the pion spacelike data, beside the $\rho$ contribution, the fit function contains also those from the $\omega$ and $\phi$ mesons; the fit returns $<\chi^{2}>/N=55/59$; however, computing the $\chi^{2}$ distance of the fit function amputated from the $\omega$ and $\phi$ contributions to the pion data from NA7 and Fermilab, one obtains $<\chi^{2}>/N=56/59$. This indicates that the (tiny) contributions of the $\omega$ and $\phi$ mesons introduce a (tiny) noise and that nothing else than the $\rho$ contribution is actually required.

This $\rho$ dominance phenomenon does not happen for the kaon form factor as, truncating the prediction for $F_{K}(s)$ from its $\omega$ and $\phi$ contributions returns a disastrous $\chi^{2}$ distance; for instance, the full $\chi^{2}$ distance of the NA7 kaon data to $F_{K}(s)$ which is $\simeq 14$, becomes $\simeq 2300$ when $F_{K}(s)$ is amputated from its $\omega$ and $\phi$ contributions. Such property is quite unexpected as the parametrization for the $\omega$ and $\phi$ [25], reminded in Section 9, is nothing but a Breit-Wigner lineshape modified in such a way that the inverse propagators fulfill $D_{V}(0)=-m_{V}^{2}$, these HK masses being those in Equations (26). That the analytic continuation, over more than $\simeq 1$ GeV2, of such a simple parametrization, with parameter values fixed at the $\phi$ mass, can provide a good description of the existing kaon spacelike data may look quite amazing.

18.3 The Pion Form Factor : The Lattice QCD Data

Some Collaborations have published information on the pion form factor derived from Lattice QCD (LQCD) data at pion masses close to the physical pion mass value [3]. We compare these LQCD form factors to our fit function $F_{\pi}(s<0)$ which, already, fairly well describes the accurate NA7 experimental data as obvious from Figure 6 and Tables 3 or 4.

The HPQCD Collaboration has thus produced a parametrization363636We thank L. Lellouch for having drawn our attention to this paper. of the pion form factor [111] $F_{\pi}(Q^{2})$ for virtualities in the range $Q^{2}\in[-0.1,0.0]$ GeV2. In order to perform a numerical comparison with our fit $F_{\pi}(s<0)$, we have plotted their pole parametrization, setting the lattice spacing $a$ at zero and using their pion radius squared $<r_{\pi}^{2}>=0.403(18)(6)$ fm2. In Figure 8, our fit function for $F_{\pi}(s<0)$ is shown by the red curve; the blue curve is the HPQCD parametrization at the central value for $<r_{\pi}^{2}>$. The red and blue curves are clearly very close to each other, if not overlapping. The green curves display the HPQCD parametrization with $<r_{\pi}^{2}>$ at $\pm 1\sigma$ from its central value. This indicates that the LQCD parametrization [111] is in fair agreement with the BHLS2 fit function $F_{\pi}(s<0)$ much beyond the expected HPQCD range of validity and validated by the DESY data points [100, 101]; Figure 8 clearly shows that the agreement between HPQCD [111] and BHLS2 extends to the prediction, shown in this Figure by the green curve, derived when fitting with discarding the spacelike data.

Six spectra for $F_{\pi}(Q^{2})$ have been published by the ETM Collaboration [112]; two of them correspond to a pion mass close to its physical value, namely the gauge ensembles $cA2.09.48$ and $cA2.09.64$. These spectra are given by the form factor values $F_{\pi}(Q^{2})$ and uncertainties as functions of $Q^{2}/M^{2}_{\pi}$. In order to compare with phenomenology and experimental data, one has to restore the values for $Q^{2}$ from the provided values for $Q^{2}/M^{2}_{\pi}$. Ref. [113], which sketches the generation of these ensembles, indicates that, numerically, $M_{\pi}\simeq 134.98$ MeV, almost identical to the RPP $\pi^{0}$ mass.

Figure 8 shows the $|F_{\pi}(Q^{2})|^{2}$ spectra derived from the $cA2.09.48$ and $cA2.09.64$ ensembles in [112]. They are displayed with their reported statistical errors373737It is worth noting that the original (unscaled) NA7 stay below the ETM data in the spectrum range where they are the most precise; this should be compared with Figure 6 in [112], especially its right panel, which exhibits the same trend.. The $\chi^{2}$ distances383838Reference [112] does not give special information about systematics and performs its own fits using uncorrelated $\chi^{2}$’s. So, following this line, we define the $\chi^{2}$ distance of the ETM data to our prediction by the sum of squared terms $d_{i}^{2}$ where : $d_{i}=(F_{ETM}(Q^{2}_{i})-F_{BHLS_{2}}(Q^{2}_{i}))/\sigma_{ETM}(Q^{2}_{i})$, using obvious notations. The quantity $\sum d_{i}^{2}$ is not submitted to the minimization procedure. of $cA2.09.48$ and $cA2.09.64$ to the BHLS2 form factor squared $|F_{\pi}(s)|^{2}$ have been computed, assuming the systematics negligible compared to the statistical errors [112]; they are displayed inside Figure 8. One can remark that the average $\chi^{2}$ distance per datum is 0.67 for $cA2.09.48$, quite a good value for data outside the fit procedure; it is only 0.06 for the $cA2.09.64$ form factor. This can be considered as an independent confirmation concerning the smallness of the systematics versus the statistical errors for the considered ETM spectra.

Stated otherwise, the data points of the $cA2.09.64$ form factor are almost exactly on the fit function. Out of curiosity, we have redone the fit in the RS variant of BHLS2 discarding the spacelike experimental data in order to get, online, the predicted form factor for any $s<0$ using only timelike data; the average $\chi^{2}$ distances become 0.59 and 0.10 for resp. $cA2.09.48$ and $cA2.09.64$. So, the $\chi^{2}$ distance of the ETM samples to the BHLS2 prediction indicates that the form factor lineshape is quite consistent with LQCD expectations; the inset in Figure 8, by magnifying the energy region close to the chiral point, illustrates the tight agreement of the prediction and of the ETM data. In other words, the ETM spectra provide strong support to the BHLS2 fits; this BHLS2 model spectrum is also intimately related to the dealing with the reported correlated systematics as promoted in [27].

For illustration, one has also displayed in Figure 8 the closest among the two data points from DESY and a part of the NA7 spectrum rescaled and not rescaled. As expected from their $\chi^{2}$ distances to the fit function, the ETM spectra are observed closer to the rescaled NA7 spectrum than to the original one.

18.4 The Neutral Kaon Form Factor at $s<0$

As seen in the previous Subsections, BHLS2 predicts the charged pion and kaon form factors in the close spacelike region with fair accuracy. Likewise, BHLS2 also predicts the neutral kaon form factor $F_{K0}(s)$; it is displayed in Figure 9. This is the analytic continuation of the form factor given in Subsection 11.4 and numerically determined within the fit procedure by the $e^{+}e^{-}\rightarrow K_{L}K_{S}$ annihilation data. One clearly observes a behavior quite different of those for $F_{K^{\pm}}(s)$, with a negative slope down to almost $s=0$, when coming from the spacelike region.

Actually, as already stated in Section 9 one has $F_{K0}(s\!=0)\!=\!0$ up to contributions of second order in breaking parameters. What is displayed in Figure 9 is $F_{K0}(s)$ computed at the fitted central value of the fit parameters393939This is equivalent to performing only one sampling on the multidimensional error covariance matrix of the fitted parameters. . To take into account their uncertainties, a Gaussian Monte-Carlo has been run which returns the averaged central value $F_{K0}(0)\simeq 10^{-6}$ or less, depending on the fit conditions. This indicates that the negative excursion close to zero observed in Figure 9 could well be a numerical artifact.

As for the charge radius squared of the neutral kaon, both BS and RS variants of BHLS2 return an average value :

[TABLE]

consistent with all reported measurements, in particular, those of the NA48 Collaboration [102] ($<r^{2}>_{K^{0}}=-0.090\pm 0.021$ fm2), or those of KTEV [103] ($<r^{2}>_{K^{0}}=-0.077\pm 0.007_{\rm stat}\pm 0.011_{\rm syst}$ fm2). It also agrees with the ChPT result obtained by Bijnens and Talavera [117] ($<r^{2}>_{K^{0}}=-0.042\pm 0.012$ fm2) and with the former measurements collected in the Review of Particle Properties [3]. We don’t know about predictions for $F_{K0}(s<0)$ at larger virtualities to which one could compare.

18.5 Pion and Kaon Charge Radii

The fit parameter central values and the parameter error covariance matrix can serve to generate several samplings ($\alpha=1,\cdots n$) and, for each, the corresponding pion and kaon form factors $F_{P}^{\alpha}(s)$. These can be used to evaluate derived physics quantities by computing the average and the r.m.s. values of the various estimates obtained from the set $\{F_{P}^{\alpha}(s),\alpha=1,\cdots n\}$.

Among the relevant physical properties of the pion and kaon form factors, their behavior close to $s=0$ deserves special interest. As clear from Section 11, the charged pion and kaon vector form factors are analytic functions of $s$ and, so, can be expanded in Taylor series around the origin. Neglecting terms of degree higher than 2 in $s$, one can write :

[TABLE]

where the $a_{P}$’s are expected equal to 1 (up to second order terms in breakings) and the $b_{P}$’s are related with the so-called charged radii by $b_{P}=<r_{P}^{2}>/6$ for which several values have been reported.

In order to have a reasonable lever arm, one has chosen a $s$ interval extending on both sides of $s=0$ and bounded404040 Considering the real part of $F_{P}(s)$ permits to avoid (minor) issues with the $\pi^{0}\gamma$ threshold effects. by $\pm 0.05$ GeV2. Using the values of each $F_{P}^{\alpha}(s)$ at three $s$ values ($0,\pm 0.05)$ GeV2, one constructs the sets $\{[a_{P}^{\alpha},b_{P}^{\alpha},c_{P}^{\alpha}],\alpha=1,\cdots n\}$, and derives their averages and r.m.s.’s. Working within the RS variant of BHLS2 and using the largest set of experimental data samples, one derives :

[TABLE]

and :

[TABLE]

As it is more customary to express the slope term in units of fm2, one displays our values together with external information in Table 5. Our results are given with two uncertainties; one is derived from the sampling on the error covariance matrix and is indexed by ”fit”. One has also varied the content of the data sample set submitted to the global fit and compared the results to those derived with the maximal sample set. It is observed that the central values for $<r_{\pi^{\pm}}^{2}>$ and $<r_{K^{\pm}}^{2}>$ get shifted while the uncertainties are almost unchanged. One has also switched from the RS to the BS variant of BHLS2; in this case, the central values get marginally shifted. Merged together, these shifts contribute an uncertainty denoted ”mod” for ”model” which, however, reflects as much the (weak) tensions between the various data samples of the largest sample set – defined in Section 17 – than the model properties by themselves.

As for $<r_{\pi^{\pm}}^{2}>$, Table 5 shows an agreement of our result with the experimental data from NA7 and Fermilab at the $1\sigma$ level or better. The result of G. Colangelo et al. [114] (CHS) is derived from a global fit of the timelike and spacelike pion form factor data based on a parametrization through the Omnès representation of $F_{\pi}(s)$ taking into account isospin breaking effects ($\omega$) and prescriptions for the asymptotics of the $\pi\pi$ $P$-wave phase shift; this estimate is in perfect agreement with our own fit result. Using an elaborate method mixing the low energy $\pi\pi$ phase information and some form factor modulus information, B. Ananthanarayan et al. [115] recently derived the value flagged by ACD in Table 5 which is very close to ours. The two flavor result of J. Bijnens et al. [116] flagged as ”ChPT 2 flavors” as well as the three flavor evaluation in [116] and the most recent Lattice QCD derivations of the HPQCD [111] and ETM [113] Collaborations are also in good correspondence with ours. Former LQCD evaluations for $<r_{\pi^{\pm}}^{2}>$ can also be found in Table 22 of [118].

As for the curvature term, varying the fit conditions as just indicated for $<r_{\pi^{\pm}}^{2}>$, one can assess the effects observed when switching from the RS to the BS variant of BHLS2 and when modifying the set of data samples submitted to fit. Altogether, this leads to our final result :

[TABLE]

The 2 flavor ChPT estimate based on all data samples available at that time (1998), [116] reports $c_{\pi}=3.85\pm 0.60$ GeV2 and the 3 flavor evaluation [117] ($c_{\pi}=4.49\pm 0.28$ GeV2) look in fair agreement with the BHLS2 evaluation given in Equation (101).

As the pion and kaon form factors, in the close spacelike region, are frequently parametrized by a pole expression :

[TABLE]

we performed alike without, however, introducing additional normalization factors as commonly done as indeed we believe in the relevance of the absolute scale of the $F_{P}(s)$’s returned by the BHLS2 fits. Such a parametrization turns out to impose the constraint $c_{P}=b_{P}^{2}$ to the curvature term. This leads to : $<r_{\pi^{\pm}}^{2}>=0.424\pm 0.002_{\rm mod}\pm 0.001_{\rm fit}$ and $<r_{K^{\pm}}^{2}>=0.268\pm 0.003_{\rm mod}\pm 0.001_{\rm fit}$.

For what concerns $<r_{K^{\pm}}^{2}>$, the numbers reported in Table 5, show that our evaluation is at $1\sigma$ from the NA7 measurement [48] and coincides with those from [94]. As for the ChPT prediction from Bijnens and Talavera [117], it differs from the value returned by BHLS2 at the $1\sigma$ level. There is no reported evaluation of the curvature term of the kaon form factor to which one could compare our own evaluation in Equations (100).

Therefore, concerning the low energy behavior of the pion and kaon form factors, the overall picture is a fair agreement of BHLS2 with the corresponding experimental measurements as well as with the existing LQCD and ChPT evaluations.

19 The Isospin 1 $P$-wave $\pi\pi$ Phase Shift

Another relevant piece of information is the $I=1$ dipion $P$-wave phase shift $\delta_{11}$. This physics has not been included within our HLS frameworks414141A relative energy scale calibration issue between phase shift and annihilation data may have been identified [114], probably amplified by the steep rise of the phase in the interesting region, combined with a 20 MeV binning size for what concerns the CERN-Munich data [119]. ; nevertheless, a comparison of the $P$-wave $\pi\pi$ phase shift extracted from BHLS/BHLS2 with other sources provides an additional cross check.

Both panels of Figure 10 display the experimental data collected by the CERN-Munich Collaboration [119], the phase shift derived using the Roy Equations [120] and the phase shift derived in [43] in the context of a common fit of the $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}$ cross-section and of the dipion spectrum in the $\tau$ decay424242Plots including other data [121] and the ChPT predictions [122] can be found in Figure 5 of [27]; they are not shown here in order to avoid overburdening Figure 10. . The phase shift coming from BHLS2 is displayed in the top panel while the phase shift from the previous BHLS [25] is displayed in the bottom panel. Both HLS based phase shifts displayed correspond to the phase of the pion form factor amplitude amputated from the $\omega$ and $\phi$ terms (see Equation (57)) to meet the usual definition of the $I=1$ $P$-wave $\pi\pi$ phase shift [120, 43]. The full (i.e. non amputated) BHLS/BHLS2 pion form factor phases are visually identical to the ones displayed in Figure 10 except for additional tiny blips in the $\omega$ and $\phi$ energy regions as can also be seen in Figure 5 of [27].

In the case of BHLS2 (red curve in the top panel of Figure 10), the phase shift superimposes almost exactly on the data points [119] and onto the JS11 curve up to the $\phi$ mass region. For this comparison with the external theoretical expectations [120, 43], one should cancel out the $\omega$ and $\phi$ terms from the pion form factor $F_{\pi}^{e}(s)$ given in Equation (57) as, indeed, these expectations do not carry such signals; one should remind that, within the BHLS/BHLS2 frameworks, these signals are generated via the mixing angles $\alpha(s)$ and $\beta(s)$ (see Equations (58)) which are proportional to the charged and neutral kaon loop difference and this vanishes in the isospin symmetry limit of equal kaon masses.

This argument might not fully apply when comparing the BHLS2 phase to the experimental data [119]. However, the CERN–Munich Collaboration reported [123] a tiny $\rho^{0}-\omega$ mixing effect – actually unsignificant (see Table 1 in ths Reference) – in the [119] phase shift data. Nevertheless, let us quickly consider this aspect of the experimental phase shift for completeness.

For a well-grounded comparison between the BHLS2 phase and the CERN–Munich data, one should rather use the BHLS2 phase smeared out over the experimental (20 MeV) bin width. The purpose of Figure 11 is to display the smeared BHLS2 phase together with the CERN-Munich data [119]. In its leftmost panel the smeared BHLS2 phase plotted is those of the full pion form factor $F_{\pi}^{e}(s)$ just as given in Equation (57); in the rightmost panel, the smeared BHLS2 phase is those of the $(\omega/\phi)$-amputated $F_{\pi}^{e}(s)$. In both panels, the CERN-Munich data points are also plotted but their symbols generally vanish underneath those for the smeared BHLS2 data. One clearly observes a good agreement between BHLS2 and the CERN-Munich data up to the $\phi$ mass region. Comparing the phase shift of the full and amputated BHLS2 pion form factors as displayed in Figure 11 clearly shows that the effects of the $\omega$ and $\phi$ signals on $\delta_{11}$ is quite marginal as expected from the [123] analysis.

In a large energy interval (about $\simeq 200$ MeV wide) centered on the $\omega$ mass, the 3 curves for the BHLS/BHLS2, Roy and JS11 phase shifts overlap within the thickness of the curves and, so, it is impossible to distinguish any one from each other. Then, both JS11 and BHLS2 start to depart from the data increasingly with increasing energies. Put together with the remarks already stated about the behavior of the dikaon cross sections slightly above the $\phi$ mass, this gives a hint for an onset of higher mass vector mesons presently not accounted for in the HLS framework.

Thus, the agreement between the BHLS/BHLS2 predictions on the one hand and the theoretical expectations [120, 43] and the CERN-Munich data on the other hand, looks quite satisfactory434343As illustrated in Figure 2 of  [124] the solution of the Roy equation requires a normalization of the energy scale at about 1 GeV. In the BHLS approach, one could adjust $g_{\rho\pi\pi}$ to get a better agreement for the high energy tail of the Roy solution..

One should also note that the bottom panel in Figure 10 exhibits a fair agreement with the datat much above the validity range of the BHLS model (the $\phi$ mass region); this can be traced back to having allowed the pion loop in the $V-\gamma$ transition amplitude $F_{\rho}^{e}$ and in the $\rho$ self-mass amplitude to carry different subtraction polynomials (see Subsection 11 above). This clearly accounts for the onset of the high mass vector mesons as the agreement remains valid up to 1.3 GeV and somewhat above.

20 Lagrangian Parameter Values in BHLS and BHLS2

Table 6 reports on the Lagrangian parameter values444444The full list of parameter values and uncertainties, especially the subtraction polynomial coefficients, are not given; they can be provided upon request. derived when running BHLS with updated kaon data (first data column). Running the basic variant BS of BHLS2 with $\tau$ data discarded gives parameters values displayed in the second data column, while the third data column reports on the numerical results from the BHLS2 Reference Solution (RS) using the largest set of consistent data samples.

Clearly, the Lagrangian parameters common to the various modelings cover a wide range of values; this does not prevent the fit properties (reminded in the last two data lines) to always exhibit a fairly good description of the data.

20.1 The ($z_{V}\div a$) and ($\Sigma_{V}\div a$) Anticorrelations within BHLS2

The value obtained for $z_{V}$ from fitting with the RS variant of BHLS2 may raise questions – see Section 8 – as $z_{V}-1$ is numerically large454545 Even if only a $2.4\epsilon$ effect, referring to the numerical estimate for generic $\epsilon$’s given in Footnote 6., even larger than $z_{A}-1$. For this purpose, it has been worthwhile to explore the RS variant behavior in fits performed discarding the $3\pi$ final state data; this has revealed a numerical anticorrelation between the HLS parameter $a$ and the BKY breaking parameter $z_{V}$, undetectable in the parameter error correlation matrix as $<\!\delta a~{}\delta z_{V}\!>\simeq-0.24$. Likewise, but of less concern, one observes a similar numerical anticorrelation between $a$ and $\Sigma_{V}$ as $<\!\delta a~{}\delta\Sigma_{V}\!>\simeq-0.18$. Several solutions with comparable total $\chi^{2}$ were found and the most striking difference among them just concerns the values for $a$, $z_{V}$ and $\Sigma_{V}$. This is well illustrated by the three following parameter sets :

[TABLE]

which exhibit clear anticorrelations between the numerical values for $a$ and $z_{V}$ on the one hand and for $a$ and $\Sigma_{V}$ on the other hand. Despite the large variations of each of $a$, $z_{V}$ and $\Sigma_{V}$, the ranges covered by the central values for $az_{V}$ and $a(1+\Sigma_{V})$ are much narrower than those for $a$, $z_{V}$ and $\Sigma_{V}$; indeed, one obtains $az_{V}=(2.63\div 2.66)$ and $a(1+\Sigma_{V})=(1.70\div 1.72)$, i.e. both products exhibit spreads at not more than the percent level.

In order to track back the source of these anticorrelations, it is worth considering the $\rho^{0}\pi^{+}\pi^{-}$ and $\phi K\overline{K}$ couplings. These are given by (see Appendix A) :

[TABLE]

which exhibit the dependences upon $\left[az_{V}\right]$ and $\left[a(1+\Sigma_{V})\right]$.

On the other hand, given the measured cross sections for $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}$ and $e^{+}e^{-}\rightarrow K\overline{K}$, the observables $\Gamma(\rho^{0}\rightarrow e^{+}e^{-})\times\Gamma(\rho^{0}\rightarrow\pi^{+}\pi^{-})$ and $\Gamma(\phi\rightarrow e^{+}e^{-})\times\Gamma(\phi\rightarrow K\overline{K})$ should be weakly model/parameter dependent within successful (global) fits. Substantially, these products are the squares of $H_{\rho}=|F^{e}_{\rho\gamma}(m_{\rho^{0}}^{2})~{}g^{\rho^{0}}_{\pi^{+}\pi^{-}}|$ and $H_{\phi}=|F^{e}_{\phi\gamma}(m_{\phi}^{2})~{}g^{\phi}_{K\overline{K}}|$, the $F^{e}_{V\gamma}(s)$’s being given by Equation (61); these quantities stretch ranges of the order $\pm 3\%$ only :$H_{\rho}=0.753\div 0.803$ and $H_{\phi}=0.124\div 0.136$. So, as the fits reported by Equations (103) are successful, they all reproduce successfully the resonance peak values for the $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}$ and $e^{+}e^{-}\rightarrow K\overline{K}$ cross sections, which implies that the products $az_{V}$ and $a(1+\Sigma_{V})$ weakly vary within different fits, even if each of these parameters scans a wider interval464646 The other parameters involved in the expressions for $g^{\phi}_{K\overline{K}}$ and $g^{\rho^{0}}_{\pi^{+}\pi^{-}}$ – i.e. $g$,$\Delta_{A}$, $z_{A}$ and the $\xi$’s – undergo variations at the percent level or better..

One has also run the BHLS2 RS variant by fixing $z_{V}=1$ and found the results shown in the last data column of Table 6. The 6 units difference between the minimum $\chi^{2}$ obtained fixing $z_{V}=1$ and the best fit leaving $z_{V}$ free is almost entirely concentrated in $\chi^{2}(e^{+}e^{-}\rightarrow 3\pi)$ which becomes 146.3 (compared to 140.4 for 158 data points). Running also the BHLS2 RS variant imposing via minuit, a lower bound $a\geq 2$ returns :

[TABLE]

This fit also returns $az_{V}\simeq 2.66$ and $a(1+\Sigma_{V})\simeq 1.75$ and shows fairly good quality (93.6% probability). Moreover, all solutions return consistent information about the physics observables as will be exemplified below. Therefore, the relatively large value for $z_{V}-1$ obtained in the best fit with the (unconstrained) RS variant of BHLS2 is not a real issue. One should also note that the multiplicity of solutions reflected by Equations (103) reduces to a single one once the 3-pion data are accounted for.

20.2 A Remark on the HLS Parameter $a$

Correlated with the remark on $z_{V}$, one should note the value for the standard HLS model parameter $a$ coming out of the best RS fit : $1.590\pm 0.001$. This is much smaller than any of its previous determinations. One should also note the unexpected value for $a$ naturally returned by the BS variant fit : $a=1.970\pm 0.002$, very close to 2 (see Table 6). Therefore, one may conclude that, within the BS variant of BHLS2, breaking parameters supply the needed $\gamma\pi^{+}\pi^{-}$, $\gamma K^{+}K^{-}$ and $\gamma K_{L}K_{S}$ couplings which, otherwise, would vanish when the $a=2$ constraint is imposed.

20.3 The ${\cal O}(p^{4})$ HLS Lagrangian : Effects of the Parameter $z_{3}$

The present work is the first one where a manifest and thorough use of the $z_{3}$ term from the ${\cal O}(p^{4})$ HLS Lagrangian is performed; Appendix A.3 displays its full expression in the BHLS2 context in terms of the ideal combinations of the vector fields. Actually, the existence of such a term was implicitly accounted for in the former BHLS [25] by disconnecting the subtraction polynomial for the $\gamma\rightarrow V$ transitions from the subtraction polynomials accompanying the pion and kaon loops occuring elsewhere – notably inside the $\rho$ propagators.

Within the BHLS2 framework, the loop corrections to the $\gamma\rightarrow V$ transition amplitudes and to the $\rho^{0}$ self-mass are tightly related474747Compare Subsection 11.3 and Appendix A.5. and additional terms of the form $c_{V\gamma}z_{3}s$ come into play – see Equations (63).

Table 6 shows that $z_{3}$ is small but significant. As the $\gamma\rightarrow V$ transition amplitudes, where solely $z_{3}$ occurs, affect all the cross sections considered, the most appropriate way to examine its effects is to compare fits where $z_{3}$ is let free and fits performed by fixing $z_{3}=0$.

The overall picture coming from these fits is reported just below. To be concise, one denotes here by BHLS${}_{2}(z_{3})$ and BHLS${}_{2}(z_{3}=0)$ the RS based fitting framework running with resp. $z_{3}$ free and $z_{3}=0$.

• •

When all data and channels considered in this paper are used (the 3–pion channel is discarded, for simplicity), one gets $\chi^{2}(\mathrm{BHLS}_{2}(z_{3})/N_{data}=982/1079$ (93.1% probability), whereas $\chi^{2}(\mathrm{BHLS}_{2}(z_{3}=0)/N_{data}=1049/1079$ (50.6% probability), a significant degradation.

A careful comparison of the fit results obtained in both cases shows that the partial $\chi^{2}$ for the $e^{+}e^{-}$ annihilation channels ($\pi^{+}\pi^{-}$, $K^{+}K^{-}$, $K_{L}K_{S}$, $\pi^{0}\gamma$, $\eta\gamma$ and even the $\tau$ dipion spectra) are comparable484848 BHLS${}_{2}(z_{3})$ remains, nevertheless, favored against BHLS${}_{2}(z_{3}=0)$ in these channels..

• •

The bulk of the degradation when going from BHLS${}_{2}(z_{3})$ to BHLS${}_{2}(z_{3}=0)$ is located in the account by the latter of the spacelike sector where the NA7 data impose accurate constraints, hardly absorbed by BHLS${}_{2}(z_{3}=0)$.

Indeed, for the spacelike experimental data, the fit performed with BHLS${}_{2}(z_{3}=0)$ returns $\chi^{2}(\mathrm{spacelike}~{}\pi)/N_{data}=78/59$ and $\chi^{2}(\mathrm{spacelike}~{}K)/N_{data}=29/25$, whereas the corresponding numbers in the (standard) BHLS${}_{2}(z_{3})$ framework are resp. $56/59$ and $18/25$, in perfect agreement with the information already reported in Table 4.

To fully compare with the information reported in Table 4, a BHLS${}_{2}(z_{3}=0)$ fit, based on only the timelike data (*i.e. * discarding the spacelike data), is useful as it provides the corresponding BHLS${}_{2}(z_{3}=0)$ predictions for the spacelike sector.

This BHLS${}_{2}(z_{3}=0)$ global fit has been performed and displays an improved fit quality, as $\chi^{2}(\mathrm{BHLS}_{2}(z_{3}=0)/N_{data}=939/995$ (72.6% probability), nevertheless lower than those for BHLS${}_{2}(z_{3})$. The $\chi^{2}$ distances of the $\pi$ and $K$ spacelike data to the corresponding BHLS${}_{2}(z_{3}=0)$ predictions are $\chi^{2}(\rm{spacelike}~{}\pi)/N_{data}=85/59$ and $\chi^{2}(\rm{spacelike}~{}K)/N_{data}=30/25.$

These two numbers should be compared with the first data column in Table 4. So, the nice feature of BHLS${}_{2}(z_{3})$ which predicts accurately the spacelike data is lost with BHLS${}_{2}(z_{3}=0)$.

One should note that, when $z_{3}=0$ is imposed, the effects accounted for by $z_{3}\neq 0$ are absorbed by the vector mesons propagators494949For instance $g$ jumps from $\simeq 7$ to $\simeq 10$ and the pion and kaon loop subtraction polynomials are severely modified. ; this is what degrades the continuation of the form factors to the spacelike region and, fitting the spacelike data, does not solve the issue.

Stated otherwise, a non-zero $z_{3}$ is the best way to treat simultaneously the spacelike and timelike data with high accuracy and comfortable fit probabilities. This may be considered as an evidence pointing toward the relevance of the ${\cal O}(p^{4})$ terms; a range for $z_{3}$ can be proposed $z_{3}=-([0.37\div 0.42]\pm 0.003_{fit})\times 10^{-2}$ for its first experimental determination.

20.4 Other Physics Results

As for other model parameters, focusing on BHLS and the RS variant of BHLS2 :

• •

One may wonder that the value for the vector coupling $g$ increases by more than 20% when going from BHLS to BHLS2, keeping good fit qualities. However, the relevant piece of information, closer to observables, being the coupling $g^{\rho^{0}}_{\pi\pi}=ag(1+\Sigma_{V})(1+\xi_{3})/2$, it is interesting to compare BHLS to the RS variant of BHLS2 : BHLS returns $g^{\rho^{0}}_{\pi\pi}=6.156\pm 0.005$ while the full BHLS2 fits gets $g^{\rho^{0}}_{\pi\pi}=6.06\pm 0.013$ in reasonable correspondence with each other.

• •

The parameters associated with the anomalous sectors show that505050In the present work, one assumes $c_{3}=c_{4}$. $(c_{3}+c_{4})/2$ yields comparable values with BHL and with BHLS2. However, for $c_{1}-c_{2}$, the difference is large; in this case, however, one should remind that the so-called ”model B” variant of BHLS excludes from the fit the $3\pi$ data in the $\phi$ mass region while the whole $3\pi$ spectra are considered within BHLS2 and are well fitted.

• •

The values for $z_{A}$ are also interesting as they are related with $f_{K}$. Indeed, one knows that [25] :

[TABLE]

is still valid within the BHLS2 framework. For this ratio, BHLS returns $r=1.257\pm 0.003_{fit}$, while the RS variant of BHLS2 finds $r=1.247\pm 0.003_{fit}$. The systematic uncertainty has been estimated to $\pm 0.020$ by varying the running conditions of the BHLS2 fitting code515151We made additional runs fixing $z_{V}=1$ (see Table 6), including the CMD-3 dikaon data and, finally, varying also the running of $\alpha_{em}$. .

This result should be compared with the Review of Particle Properties [125] which recommends $r=1.197\pm 0.002\pm 0.006\pm 0.001$ and with the LQCD determination for $N_{f}=2+1+1$ flavors [118] $r=1.195\pm 0.005$.

• •

The additional parameters introduced by the CD breaking and the Primordial Mixing are all found at the few percent level. The difference for $\xi_{0}$ and $\xi_{3}$ between the (full) RS variant and the BS variant should be noted; it may have to be revisited later on, when a complementary mechanism may allow one to absorb satisfactorily the $\tau$ data within the BS global fit procedure. This may also influence the $3\pi$ sector.

21 The muon LO-HVP Evaluation

BHLS2 (or BHLS) fits provide the contributions of the various hadronic channels525252All along this Section, one may use interchangeably $a_{\mu}^{\rm HVP-LO}$ or, preferably, $a_{\mu}$ to lighten writing; the context will be unambiguous. it covers to $a_{\mu}(\sqrt{s}\leq 1.05~{}{\rm GeV})$. An overview of the main results, obtained in different running conditions is given in Table 7. These have been derived by submitting to fit most of the data samples commented in Section 13 which define a set of consistent samples covering the six annihilations considered; altogether, these 6 channels and the 10 radiative meson decay partial widths represent 1237 data points.

21.1 Comparison of BHLS2 Estimates for $a_{\mu}^{\rm HVP-LO}$ with Others

Without going lengthily into detailed comparisons, it looks worthwhile to provide some elements. For instance, one can extract from Figure 3 in [27], values for $G=a_{\mu}^{\pi\pi}(0.630\!<\!\sqrt{s}\!<\!0.958~{}{\rm GeV})$; our updated estimate for $G$ is $G({\rm BHLS})=(355.59\pm 0.58)\times 10^{-10}$, very close to its former evaluation (see Fig. 3 in [27]) while the numerical integration of the NSK+KLOE data was estimated $G({\rm exp})=(356.67\pm 1.69)\times 10^{-10}$; BHLS2 gives535353As the $(0.630<\sqrt{s}<0.958)$ GeV interval is relatively far from both the threshold and the $\phi$ mass regions, one could indeed expect similar results from BHLS and BHLS2. $G({\rm BHLS}_{2})=(356.50\pm 0.75)\times 10^{-10}$.

In a quite different theoretical context, the recent analysis [114], provides the estimate $G=a_{\mu}^{\pi\pi}(\sqrt{s}<0.630~{}{\rm GeV})=(132.5\pm 0.4_{\rm stat}\pm 1.0_{\rm syst})\times 10^{-10}$ while the BHLS gives $G({\rm BHLS})=a_{\mu}^{\pi\pi}(\sqrt{s}<0.630~{}{\rm GeV})=(130.80\pm 0.22_{\rm fit})\times 10^{-10}$. In the same region one also gets $G({\rm BHLS}_{2})=a_{\mu}^{\pi\pi}(\sqrt{s}<0.630~{}{\rm GeV})=(130.03\pm 0.18_{\rm fit}+[~{}^{+0.25}_{-0.15}]_{\rm syst})\times 10^{-10}$.

Reference [114] also provides $H=a_{\mu}^{\pi\pi}(\sqrt{s}<1.0~{}{\rm GeV})=(494.7\pm 1.5_{\rm stat}\pm 2.0_{\rm syst})\times 10^{-10}$ while the updated BHLS gives  $H({\rm BHLS})=(490.65\pm 0.84_{\rm fit})\times 10^{-10}$ and BHLS2 leads to $H({\rm BHLS}_{2})=(490.75\pm 0.90_{\rm fit}+[~{}^{+0.31}_{-0.26}]_{\rm syst}))\times 10^{-10}$. The method in [114] and ours are different, but more importantly for this energy range, the dipion BaBar data included in the study [114] and discarded in ours, certainly contribute to the difference.

One may also compare the BHLS2 results with others for the contribution to $a_{\mu}^{\rm HVP-LO}(\pi^{+}\pi^{-})$ from the energy region $(0.6\leq\sqrt{s}\leq 0.9)$ GeV. KNT18 [71] gets $(369.41\pm 1.32)\times 10^{-10}$ and CHS [114] $(369.8\pm 1.3_{fit}\pm 1.3_{syst})\times 10^{-10}$ using BaBar data at low energies. BHLS2 running with KLOE10/12, CMD-2, SND, BESIII and Cleo-c data in the $\pi^{+}\pi^{-}$ channel yields $(367.81\pm 0.77)\times 10^{-10}$; replacing the KLOE10/12 data by the KLOE85 combination [70], one gets $(367.61\pm 0.64)\times 10^{-10}$. So the observed effect from BaBar in this energy region amounts to $\simeq 2\times 10^{-10}$. One also observes the gain on accuracy by using BHLS2 : A factor of 2 compared to [71], a factor of 3 compared to [114].

21.2 Estimates for $a_{\mu}^{\rm HVP-LO}(\sqrt{s}\leq 1.05~{}{\rm GeV})$

In order to get the full contribution to $a_{\mu}^{\rm HVP-LO}$ of the energy region up to 1.05 GeV, one should complement the contribution from the channels covered by BHLS/BHLS2 by those of the channels presently outside this framework ($4\pi$, $2\pi\eta$, …). The most recent evaluation of this quantity, estimated by direct numerical integration of data is :

[TABLE]

In view of model effect studies, it has been found motivated to update and correct the data samples submitted to the BHLS fit in the way stated in Section 13. One has thus derived the updated BHLS evaluations given in the first data column of Table 7; this supersedes the BHLS results formerly given in Table 4 of [27]. Differences between BHLS and the RS BHLS2 variant for some channel contributions are observed and, in fine, $a_{\mu}({\rm HLS},\sqrt{s}\leq 1.05~{}{\rm GeV})$ only increases by $\simeq 1\times 10^{-10}$ while its uncertainty slightly decreases ($1.03\rightarrow 0.92$, in units of $10^{-10}$). An evaluation of possible additional systematics, specific of the BHLS modeling, has been performed and summarizes by [27] :

[TABLE]

Each piece shown here was found to act by shifting the central value for $a_{\mu}$ rather than enlarging its uncertainty. The first term refers to the uncertainty in the treatment of the $\phi$ mass region in the 3-pion spectra, the second term takes into account the uncertainty on the $\pi\pi$ threshold behavior within BHLS (see Subsection 10.3 above and, especially, Figure 1); finally, the third piece reflects differences observed by running the BHLS fit procedure [25, 26] with the $\tau$ data included and excluded.

As already stated, the BHLS2 modeling has been motivated by the aim to cancel out the first two sources of uncertainties reported in Equation (106) by a better understanding of the $s=0$ region and a full account of the 3-pion data up to the $\phi$ region.

For this purpose, the three data columns in Table 7 referring to BHLS2 show that :

• •

The evaluation for the various contributions to $a_{\mu}(\sqrt{s}\leq 1.05~{}{\rm GeV})$ derived using the BS and both RS variants of BHLS2 carries differences at a few $10^{-12}$ level.

• •

Using the $\tau$ data has some effect. Indeed, Table 7 shows that running the RS variant of BHLS2 including the $\tau$ data improves the uncertainty in the dipion channel by $\simeq 10\%$ while leaving the central value almost unchanged. This illustrates that the $e^{+}e^{-}$ and $\tau$ data are not conflicting within the RS variant of BHLS2.

Another property exhibited in Table 7 should be noted :

• •

Whereas the sum of the various contributions to $a_{\mu}$ derived using several variants of the BHLS2 fit almost coincides with the estimate derived by direct integration of the data, the $\pi^{+}\pi^{-}$ contributions differ by 3.72 in units of $10^{-10}$, a significant difference.

In the direct integration method, the normalization of each sample is the nominal one and all uncertainties participate to the definition of its weight in the combination with the other samples involved. Instead, in the global fit method, the absolute normalization uncertainty is treated specifically and in consistency for all the data samples involved in all the annihilation channels considered ($\pi\pi$, $K\overline{K}$, $\pi\pi\pi$ and $(\pi^{0}/\eta)\gamma$). It has beeen shown [27] that the iterative method which underlies the HLS global fit procedures gives unbiased results 545454This is well reflected, in the present work, by the nice matching observed between the predicted pion form factor in the spacelike region and the NA7 and ETMC data..

Therefore, that values of definite contributions to $a_{\mu}$ may significantly depend on the method used to derive them is not totally unexpected. That the sum of them may almost coincide is, however, accidental555555 In this respect, it is worth noting that a quite similar behavior is observed in the comparison of the KNT18 and DHMZ17 data reported in [71]. Indeed, Table 5 herein shows that the estimated dipion contributions (based on the same data) differ by $3.40\times 10^{-10}$, whereas the total sums giving the HVP-LO’s differ by only $0.2\times 10^{-10}$..

The data update having been performed, BHLS and BHLS2 have been run on the same data and, thus, the differences between their respective evaluations reflect modeling effects. One can also compare these differences to the numbers listed in Equation (106).

As shown by the various topics examined in Section 18, one can legitimately consider the $\pi^{+}\pi^{-}$ threshold region accurately treated. Therefore, the difference $494.08-493.73=0.35$ between the BHLS2 and BHLS estimates for $a_{\mu}(\pi^{+}\pi^{-},\sqrt{s}\leq 1.05~{}{\rm GeV})$ is a fairly well motivated evaluation of the unaccounted for nonet symmetry breaking effects within the BHLS framework; it is therefore justified to replace the former $+1.4$ estimate in Equation (106) by $+0.35$.

Table 3 shows that all variants of BHLS2 take well into account the $3\pi$ annihilation data over the whole energy domain up to the $\phi$ mass region, the fit quality of the RS variant being noticeable. The difference between BHLS2 and BHLS for this channel $\simeq 0.5\times 10^{-10}$ is in the range expected from the estimate [27] reported by the first term in Equation (106). Nevertheless, the $0.06\times 10^{-10}$ difference for the $3\pi$ entry between the RS variant of BHLS2 runnings including and excluding the $\tau$ data might be considered as additional systematics565656There might be a numerical effect reflecting the physics correlation between the $\pi^{+}\pi^{-}\pi^{0}$ annihilation channel and the $\tau$ spectrum due to their common $\rho^{\pm}$ exchanges taking place in their intermediate states. .

It is interesting to compare the outcomes from fits to the values derived by the direct integration of the same data sample set treated in the BHLS2 framework. One may note that the sum of all HLS channel contributions differs from the fit expectations by only $\simeq 0.6\times 10^{-10}$; however, this hides the fact that the dipion channel is fitted $3.7\times 10^{-10}$ smaller and the 3-pion channel is fitted $1.70\times 10^{-10}$ larger.

In summary :

[TABLE]

is the most conservative final BHLS2 answer, where the uncertainty collects statistical and systematic errors as reported by the various experiments. The remarks just stated concerning the 3-pion channel allows to keep an additional uncertainty of $-0.06\times 10^{-10}$, playing as a bias.

The issue is to examine whether more significant additional sources of systematics are not at work. For this purpose, the consequence on $a_{\mu}({\rm HLS})$ of the various fits referred to in the previous Section carry a relevant piece of information, as also for the fits underlying the results given in Equations (103) and (105) or those displayed in the last data column in Table 6. These fits, covering a wide range of parameter values, are relevant for systematics estimates.

The runs leaving aside the 3-pion data briefly reported in Equations (103), show that the dominant $\pi\pi$ contribution can be as low as $(493.80\pm 0.89)\times 10^{-10}$, i.e. $0.30$ units of $10^{-10}$ smaller than the corresponding reference entry in Table 7. Nevertheless, the middle entry line in Equations (103) is very close to Equations (105), where the 3-pion data have been included and exhibits negligible differences with respect to our reference in Table 7. Finally, the fit corresponding to the last data column in Table 6 has also been analyzed and returns a value for $a_{\mu}(\pi^{+}\pi^{-},\sqrt{s}\leq 1.05~{}{\rm GeV})=(494.32\pm 0.90)\times 10^{-10}$. This may indicate another additional systematic able to shift the central value by, at most, $0.22\times 10^{-10}$ affecting our estimate for $a_{\mu}^{\rm HVP-LO}(\sqrt{s}\leq 1.05~{}{\rm GeV})$ given in Equation (107).

The possible sources of additional systematics just emphasized can be considered as a model uncertainty bounded by $(+0.22,-0.36)\times 10^{-10}$. On the other hand, running the fitting code using the available parametrizations for $\alpha_{em}(s)$, returns differences for the estimates for $a_{\mu}({\rm HLS})$ of $\simeq 0.4\times 10^{-10}$, mostly located in in the $\phi$ mass region.

Finally, the question of using the CMD-3 data can be raised. We have run our code using the CMD-3 data [77, 79] for both $K\overline{K}$ decay channels, discarding the kaon data from BaBar, CMD-2 and SND and using only the diagonal part of their error covariance matrices (see Section 16 above). In total, the change for $a_{\mu}(\sqrt{s}\leq 1.05~{}{\rm GeV})$ is noticeable $(572.52\pm 1.12)\times 10^{-10}$, $0.54\times 10^{-10}$ larger than our reference in Equation (107), essentially coming from the $K^{+}K^{-}$ and $K^{0}\overline{K}^{0}$ channels. On the other hand, combining CMD-3, BaBar, (corrected) CMD-2 [75] and SND, leads to $a_{\mu}(\sqrt{s}\leq 1.05~{}{\rm GeV})=(572.52\pm 1.04)\times 10^{-10}$. This leads us to complete our estimate Equation (107) :

[TABLE]

where the additional systematics come from model variations, $\alpha_{em}$ parametrizations, and data sample consistency issues. The model uncertainties include the marginal effect attributable to the $\tau$ data, as this could reflect some (marginal) isospin breaking shortcoming. We have chosen conservatively as reference the RS variant fit excluding the $\tau$ data. One should also note that the tension between dikaon data samples introduces non-negligible systematics which contribute a bias. This may have to be revisited when new data will arise.

Taking into account the data upgrade, the corresponding quantity for BHLS is :

[TABLE]

using the present findings to go beyond the former systematics evaluation reminded in Equation (106). The various additional uncertainties just estimated represent conservative upper bounds.

21.3 Dispersive vs LQCD Methods : Additional Systematics

Experimental $e^{+}e^{-}\to{\rm hadrons}$ data are dressed by photon radiative corrections and disentangling strong interaction from the electromagnetic effects is a non-trivial task. What is ideally needed as input to the dispersion integrals is the one-particle irreducible (1PI) hadronic part. It is given in QCD by the correlator between two hadronic currents, and this is what is primarily calculated as the LO part of order $O(\alpha^{2})$ in lattice QCD. We remind that the bare cross-section undressed from photon radiation effects is not an observable by itself. It requires some theory input, primarily the photon radiative corrections, for its extraction from the physical dressed data. Mandatorily, one has to drop out the photon vacuum polarization effects to get the undressed bare cross-section. By convention, one also includes final state radiation (FSR) in the dipion channel by adopting scalar QED (sQED) for its calculation. Similarly, the radiative decay channels $\pi^{0}\gamma$ and $\eta\gamma$ are included in our evaluation. Like FSR, they are related to light-by-light scattering insertions.

However, there are other important QED effects like the $\rho^{0}-\gamma$ interference which is substantial in the dipion channel because this QED effect is magnified by the resonance enhancement in the $\rho$ region. In the standard dispersive approach, $\rho^{0}-\gamma$ mixing effects – inherent in the $e^{+}e^{-}$ data – are treated as part of the LO-HVP contribution and are included as well in the calculations of the higher-order contributions. An updated $e^{+}e^{-}$ data-based evaluation of the LO and the NLO parts yields the results shown in Table 8575757For the present exercise, the BaBar dipion data have been taken into account.

Here we have to remind that, in contrast to the usual dispersive estimations, LQCD calculations are based on the hadronic current correlator such that $\rho^{0}-\gamma$ mixing is not included in the leading order LQCD results (unless it is added by hand from phenomenological analyses, if not calculated separately).

So the question of how to disentangle the QCD effects from QED ones, has to be reconsidered under this aspect. One possibility is to apply effective field theory methods like the resonance Lagrangian approach, the other is lattice QCD, the only method which can answer the question from first principles. In fact, to disentangle better QCD from QED effects in the $e^{+}e^{-}$ data-based dispersive approach, one has to subtract the $\rho^{0}-\gamma$ interference from the $e^{+}e^{-}\to\pi^{+}\pi^{-}$ data585858The similar $\omega-\gamma$ and $\phi-\gamma$ effects are much smaller and can be neglected., expecting to get the effect back as a part of the higher-order contributions.

Here is the place to remind that the $\rho^{0}-\gamma$ mixing is responsible for the discrepancy between the isospin rotated isospin breaking corrected $\tau$ spectral-function and the $e^{+}e^{-}$ spectral-function [126, 127, 40]; this has been shown in a straightforward “VMD-II+sQED” simple effective resonance Lagrangian model [43] as well as in the BHLS approach [29, 25] as utilized in the present analysis. By convention, to bring $\tau$ data closer to $e^{+}e^{-}$-data, it looked natural to correct the $\tau$-data to conform with the commonly adopted dispersive approach based on the $e^{+}e^{-}$ data.

In fact, the measured $\tau$-decay dipion spectra are devoid of photon interference effects. According to what we just have been arguing, one should actually not apply this mixing correction to the $\tau$-decay data, as it has been done within our BHLS global fit framework or in the much simpler model applied in [43]. One should better subtract the $\rho^{0}-\gamma$ mixing from the $e^{+}e^{-}$ data, to get an object closer to a pure hadronic “blob” and then proceed as usual. As a possible subtlety, this requires to work in the “VMD-I” basis where the $\rho^{0}-\gamma$ system is parametrized in terms of fields where a direct lepton-$\rho$ coupling is absent. This is the basis of the BHLS Lagrangian formulation. Being proportional to the hadronic production cross-section, the hadronic 1PI amplitude and the corresponding $R(s)$ (see Eq. (1)) do not depend on rotations in the $(\gamma,\rho^{0})$ field space.

The field rotation invariance has been directly confirmed by comparing, on the one hand, the “VMD-II + sQED” modeling adopted in [43], which uses the mass eigenstate basis where a direct $\rho$-lepton coupling is induced by the necessary field rotation [128] and, on the other hand, our BHLS modeling where a direct $\rho$-lepton coupling is absent. One thus can apply the standard diagrammatic expansion of Fig. 12 with the modified effective $R(s)$-function and obtains the results presented in Table 9 :

which corresponds to what, in principle, lattice QCD is doing.

The LQCD results displayed in Fig. 13 are full HVP results, which include $ud$ light quark, strange and charm contributions including the connected and disconnected parts as well as isospin breaking contributions either calculated on the lattice or adopted from phenomenological estimates in terms of data. They do not include $\gamma-\rho$ mixing to our knowledge, however. In lattice QCD this effect is counted as a NLO effect, which has to be evaluated separately595959 Contributions from FSR and radiative decays should have been added in a consistent comparison with data-based dispersive results. .

According to our Tables, the commonly adopted dispersive set up gets a shift606060The systematic uncertainties in Tables 8 and 9 are 100% correlated, except for the sub-dominant diagram c. of $\delta a_{\mu}[\pi\pi]_{\rho\gamma}\simeq+(3.10\pm 0.03)\times 10^{-10}$ at NLO, i.e. a 1 $\sigma$ increase. See also the results presented in Table 3 of [129] where a $\rho^{0}-\gamma$ mixing shift of $2.74\times 10^{-10}$ on the I=1 component at LO has been reported as a difference between LQCD data and the phenomenological dispersive approach.The upgraded result for $\delta a_{\mu}[\pi\pi]_{\rho\gamma}$ is based on an actualized data compilation, with mixing correction applied to the full vector form-factor and including NLO corrections.

More or less surprisingly, the sum of LO+NLO differs in the two settings, i.e. the mixing effect subtracted from the LO result does not show up as a higher-order effect of comparable magnitude. It is actually not surprising because the $\rho^{0}-\gamma$ mixing effect, if included in the LO dispersion integral, is larger than if taken into account in higher-order diagrams where the effect is weighted by the higher-order kernel functions. Our result based on correcting the $e^{+}e^{-}$ dipion channel in the range up to 1 GeV, where sQED estimates can be trusted, indirectly compares with the deviation $(9.1\pm 5.0)\times 10^{-10}(1.8\sigma)$ between the $\tau$-based hadronic contribution to $a_{\mu}$ and the $e^{+}e^{-}$-based one reported in [130]. Our result does not include $\tau$ data while the $\tau$ based result of [130] includes $\tau$ data up to 1.8 GeV which are combined with $e^{+}e^{-}$ data which still include $\rho^{0}-\gamma$ mixing. That the results are increasing in both cases is as expected, but there is no reason why the numerical results should agree.

Our suggestion to base the standard bookkeeping such as Figure 12 on a purer QCD $\gamma-\rho$ undressed hadronic blob is motivated by the fact that lattice QCD does not include any QED effects if not taken into account as an extra contribution and also by the fact that $\tau$-data are free of similar mixing effects, and it looks somewhat artificial to add effects to $\tau$-decay spectral data that are absent in the corresponding measurements. Nevertheless, as already said, surprisingly we subtract a non-negligible contribution, which we do not get back in higher orders. It is important to note that the shift $\delta a_{\mu}[\pi\pi]_{\rho\gamma}$ is real as it derives from the manifest discrepancy observed between $\tau$ and $e^{+}e^{-}$ spectral-function data in the corresponding time-like $\pi\pi$ production processes. As mentioned before, the $\gamma-\rho$ interference is substantial because of the manifest resonance enhancement in the time-like regime. Lattice QCD is taking advantage of the fact that, for the calculation of $a_{\mu}^{\rm had}$, only the space-like VP-function is needed. However, there is no $\rho$-resonance peak in the Euclidean region which could enhance a photon-exchange interference effect; the $\rho$-resonance is completely smeared out and a resonance-enhancement cannot be localized there in any obvious way. Because of a model-independent clean separation of QCD is only possible in lattice QCD, one has to wait for the relevant LQCD results to get a better understanding of QED-QCD interference effects.

One may conclude that *all * standard dispersive evaluations should either be upgraded by adding $+(3.10\pm 0.31)\times 10^{-10}$ to the LO+NLO results616161Concerning the model dependence of the prediction for the $\rho^{0}-\gamma$ interference, we remind that this effect depends on only one extra parameter, the well known leptonic width of the $\rho^{0}$ or, equivalently, on the coupling $g_{\rho ee}$, also well known experimentally. Since the effect is small, a few units in $10^{-10}$, we can generously assign a very conservative 10% model uncertainty, without substantially increasing the overall uncertainties coming with the $e^{+}e^{-}\to\pi^{+}\pi^{-}$ data., or accounted for in the systematic errors by adding there the amount of this shift. For the time being, we have adopted this second solution.

21.4 BHLS2 Evaluation of $a_{\mu}^{\rm HVP-LO}$

With the BHLS2 and BHLS evaluations of $a_{\mu}^{\rm HVP-LO}({\rm HLS},\sqrt{s}\leq 1.05$ GeV displayed in Equations (108) and (109) and the information given in the left part of Table 10, one can derive the corresponding values for the full $a_{\mu}^{\rm HVP-LO}$. One thus gets :

[TABLE]

In the case of BHLS, due to the inclusion of $\tau$ channel within the global fit procedure, the corresponding systematics have been added linearly to the others. Indeed, as stated several times for BHLS and BHLS2, the observed systematics rather come as biases and so, should be treated separately from the fit error which fully absorbs the reported statistical and systematics from the various experiments encompassed inside each of our HLS frameworks.

The additional systematics absorb uncertainties which can be attributed to the model workings and to the (marginal) tension observed between some of the selected experimental data samples. For clarity, the error coming from the $\rho^{0}-\gamma$ interference effect discussed just above is given separately.

Our evaluations for $a_{\mu}^{\rm HVP-LO}$, given by Equations (110), are shown in Figure 13 together with some other recent results derived in different phenomenological contexts [4, 130, 71] or via LQCD data samples with $N_{f}=2+1+1$. As for their specific HLS modeling parts, the corresponding fit probabilities shown by the down-most entries in Table 7 are always fairly good.

As for modeling effects, the central values for BHLS and BHLS2 – running on the same set of data samples – are shifted apart by $\simeq 1.5\times 10^{-10}$; however, if one takes into account the magnitude of the additional systematics, Figure 13 clearly indicates that this shift is not significant and so no noticeable modeling effect shows up. Moreover, Figure 13 clearly shows that both HLS based estimates for $a_{\mu}^{\rm HVP-LO}$ are consistent with the evaluation derived by numerical integration of the data626262See the top point in Figure 13, derived by excluding the BaBar dipion sample for a consistent comparison with the BHLS entries. .

Actually, the estimates for the muon HVP-LO derived by the various dispersive methods (the five data points displayed in the bottom part of Figure 13) are quite comparable with each other as they use as input to evaluate $a_{\mu}^{HVP-LO}$ the $e^{+}e^{-}$ annihilation cross sections as derived from the data or as reconstructed from fit. Therefore, the results thus provided do not split out the electromagnetic corrections; in particular the $\gamma-\rho^{0}$ mixing effects is merged in the plotted systematics only for the BHLS/BHLS2 entries, not for the other dispersive results plotted.

21.5 BHLS2 Evaluation of Muon Anomalous Magnetic Moment

Having at hand our estimate Equation (108), one can derive the BHLS2 evaluation for $a_{\mu}=(g-2)/2$, the muon anomalous magnetic moment. The left part of Table 10 collects the pieces to add up with $a_{\mu}({\rm BHLS_{2}})$ to get the full leading order HVP, its value is shown as the top entry in the right-hand part of Table 10 where the other contributions are included. The value for $\Delta a_{\mu}=a_{\mu}({\rm exp})-a_{\mu}({\rm th})$ and its statistical significance are the bottom entries there. Taking into account all additional systematics, one thus gets :

[TABLE]

and, taking into account the shifts produced by all systematics, one gets :

[TABLE]

which means that the statistical significance of $\Delta a_{\mu}$ is predicted greater than $3.82\sigma$. The shift just reported is partly due to the estimates of the various modeling effects involved, partly due to somewhat conflicting aspects of some data samples and parly due to the identified $\gamma-\rho^{0}$ mixing effects which add up linearly.

Our estimate for $a_{\mu}$ given in Equation (111) relies on a variant of the traditional estimate [132] for the Light-by-Light (LBL) contribution ($(10.34\pm 2.88)\times 10^{-10}$). Danilkin, Redmer and Vanderhaegen (DRV) have quite recently published a comprehensive analysis of all ingredients participating to the LBL term and proposed for it the much more precise value [142] $(8.7\pm 1.3)\times 10^{-10}$ which leads to :

[TABLE]

and to $\Delta a_{\mu}^{\rm DRV}=(34.32\pm 7.10)\times 10^{-10}$. This makes the significance for $\Delta a_{\mu}$ jumping from $(3.8\div 4.5)\,\sigma$ to $(4.2\div 4.7)\,\sigma$.

22 Summary and Conclusions

The present paper has addressed several topics covering modeling, phenomenology and the evaluation of the muon anomalous magnetic moment which puts forward, as well known, a hint for a physics effect beyond the Standard Model expectations. For the sake of clarity, each of these topics is considered separately.

From BHLS to BHLS2

The BHLS framework [25], although quite successful [26, 27, 143], has been shown to exhibit a non-optimum behavior in limited and well identified kinematical regions : The threshold region and some aspects of the $\phi$ mass region; the former issue renders delicate the continuation of form factors across the chiral point, the latter has led to discard the $\phi$ mass region for 3-pion final state – the so-called Model B [25].

In order to cure these diseases, one has initiated a new breaking procedure of the general HLS framework [18]. For this purpose, besides the BKY mechanism [28, 34, 35], one has introduced a breaking mechanism taking place at the level of the covariant derivative (CD) itself, i.e. inside the basic ingredient of the HLS Model. This (seemingly novel) CD breaking mechanism allows naturally the $\rho^{0}$, $\omega$ and $\rho^{\pm}$ fields to carry different Higgs-Kibble masses. This takes its importance when going to loop corrections to the vector meson mass matrix [29, 25] which generate a dynamical mixing of the vector mesons calling for an $s$-dependent redefinition of the vector meson fields to recover mass eigenstates.

Thanks to its nonet symmetry breaking piece, the CD breaking scheme generates a non-zero $m_{\rho^{0}}-m_{\omega}$ difference which naturally makes the three $s$-dependent mixing angles to vanish at the chiral point. In this way, the timelike to spacelike analytic continuations become smooth without any ad hoc trick.

Relying on BKY and the CD breaking, the Basic Solution (BS) is defined and represents the first variant of BHLS2. Another mechanism can also be invoked which turns out to state that the neutral vector fields involved in physical processes are not their ideal combinations but mixtures of these, constructed via a rotation. This Primordial Mixing (PM) mechanism636363The Primordial Mixing of neutral vector mesons resembles the mixing scheme for the neutral pseudoscalar mesons developed in [144]. of vector mesons assumes a ${\bf 1}+{\cal O}(\epsilon)$ rotation. Working at first order in breaking parameters, the merging of the BKY and CD breakings together with the PM mechanism define our Reference Solution (RS), the second variant of BHLS2 covering the full HLS scope more easily than the BS variant (as it stands presently).

In both variants of BHLS2, the ${\cal L}_{V}$ piece of the non-anomalous Lagrangian is modified while its ${\cal L}_{A}$ piece remains unchanged646464This might have to be revisited in order to improve the working of BS in the $\tau$ and $\pi^{+}\pi^{-}\pi^{0}$ sectors of the broken HLS Lagrangian. compared to the former BHLS [25].

The $e^{+}e^{-}$ Annihilation Phenomenology

The BKY and CD breaking schemes, together with the Primordial Mixing, allow to construct the non-anomalous and anomalous pieces of BHLS2, a new (broken) HLS Lagrangian. This allows for a simultaneous study of the $e^{+}e^{-}$ annihilation into 6 final states ($\pi^{+}\pi^{-}$, $\pi^{0}\gamma$, $\eta\gamma$, $\pi^{+}\pi^{-}\pi^{0}$, $K^{+}K^{-}$, $K_{L}K_{S}$) within a unified framework where global fits can be worked out.

The fits have been performed by taking into account the fact that most of the data samples carry normalization uncertainties, frequently dominant; these are appropriately considered in order to avoid biased results for physics quantities. The iterative method developed in [27] is thus intensively used in the global fits reported here; the convergence criterion chosen to stop the iteration at order $n+1$ was $\delta\chi^{2}=|\chi^{2}_{\rm tot}(n+1)-\chi^{2}_{\rm tot}(n)|<(0.3\div 0.5)$. As $\chi^{2}_{\rm tot}\simeq 1000\div 1200$ for $\simeq 1150\div 1240$ data points, this convergence criterion is clearly constraining.

The fitting procedure has allowed to detect (and discard) a very few data samples which hardly accommodate the global framework fed with almost all the existing data samples. Beside the already identified [27] two $\pi^{+}\pi^{-}$ samples, one also found that, out of the 7 available $\pi^{+}\pi^{-}\pi^{0}$ data samples, one [87] hardly fits the global context, as also reported by others [92]. Figures 2, 3 and 4 together with Table 3 prove the fairly good quality of the global description for the $\pi^{+}\pi^{-}$, $(\pi^{0}/\eta)\gamma$ and $\pi^{+}\pi^{-}\pi^{0}$ channels. As for the 3-pion channel, one should note that BHLS2 encompasses easily the samples collected in the $\phi$ mass region and so, solves the difficulty encountered there by the former BHLS.

A thorough study of the available data samples has been performed concerning both the $K^{+}K^{-}$ and $K_{L}K_{S}$ channels for which some tension is expected, as reported in [77]. The analysis in Section 16 has shown that the SND, BaBar [76] and the updated [75] CMD-2 data are in fairly good agreement with each other and accommodate easily the HLS context where all other kinds of data are encompassed. Tension has indeed been observed between the CMD-3 data [77, 79] and the other samples. This tension has been traced back to the off-diagonal part of the CMD-3 error covariance matrices for both the $K^{+}K^{-}$ and $K_{L}K_{S}$ spectra; this analysis is summarized by Table 2. In order to include the BaBar (and/or CMD-3 data) within the global framework, energy shifts relative to the CMD-2/SND energy scale (chosen as energy scale reference) are mandatory and have been fitted at values consistent with the reported expectations [76, 77]. To perform this exercise, when using the CMD-3 data, the error covariance matrices of these (only) have been amputated from their off-diagonal part.

One can finally remark that, once the energy scales of the CMD-3 and BaBar [76] data relative to CMD-2/SND are applied, one does not observe a visible tension among them any longer as illustrated by Figure 5, where all data samples are nicely superimposing onto the same fit function.

Meson Form Factors in the Spacelike Region

The BHLS2 form factors supplied with the parameter values derived from global fits to the *timelike region data only * have been shown to nicely predict the behavior of the $\pi^{\pm}$ and $K^{\pm}$ form factors in the close spacelike region, as they can be expected from NA7 and Fermilab data. Moreover, including these (model-independent) spacelike data within the BHLS2 fit procedure does not exhibit any gain in the description of the spacelike region as obvious from Table 4 and from Figure 6.

Our derived $\pi^{\pm}$ and $K^{\pm}$ form factors, which perform fairly well simultaneously in the spacelike and timelike regions, strongly supports the validity of the BHLS2 predictions in the $s$-region located in between. If the fairly good continuation of $F_{\pi^{\pm}}(s)$ is actually performed over a limited gap in $s$ to reach the spacelike region, one may wonder about the so successful continuation of $F_{K^{\pm}}(s)$ over an $s$-gap greater than 1 GeV2. This indicates that our modified Breit-Wigner formulas are performing (unexpectedly) well far beyond the $\phi$ energy region where their parameters are determined.

The LQCD Collaborations HPQCD and ETM have provided $F_{\pi^{\pm}}(s)$, the former as a parametrization, the latter by corresponding spectra, covering the same space-like $s$-region as the NA7 experiment. Figure 8 shows the superposition of the BHLS2 prediction (without any free parameter) together with the parametrization provided by HPQCD and the two ETM spectra; the accord is very good for HPQCD and simply perfect for both the ETM spectra. There is, unfortunately, no reported LQCD data for the kaon form factors to which one could compare.

The success just emphasized for the BHLS2 $\pi^{\pm}$ and $K^{\pm}$ form factors compared to spacelike experimental or LQCD data gives confidence in other physics information involving the chiral point region. Let us limit oneself, in this Subsection, to briefly comment on the neutral kaon electromagnetic form factor. As reminded in the main text, the slope for the $F_{K^{0}}(s<0)$ form factor at $s\rightarrow 0$ is reported negative. This is qualitatively consistent with the BHLS2 finding for $F_{K^{0}}(s)$ shown in Figure 9, as it is found to carry a negative slope at negative $s$ and a positive slope at positive $s$. However, the effects of the neglected ${\bf O}(\epsilon^{2})$ corrections actually prevent to provide confidently accurate slope predictions in the neighborhood of $s=0$ which is a delicate region for our ${\bf O}(\epsilon)$ approximations.

The pion form factor measured by the NA7 and Fermilab Collaborations are model independent and, actually, it cannot be measured in a model-independent way above $|s|\simeq 0.25$ GeV2. To go beyond, one should accept to consider data where the extraction of $F_{\pi^{\pm}}(s)$ is model-dependent. Figure 6 already indicates that the two unfitted desy measurements fall (almost) exactly onto the fitting curve together with the model-independent NA7 and Fermilab measurements. The JLAB data points are more or less closer to the BHLS2 expectation depending on the method used to extract $F_{\pi^{\pm}}(s)$. Nevertheless, our analysis indicates that the BHLS2 prediction could well be valid down to $\simeq-1$ GeV2, possibly slightly beyond.

In summary, our analysis clearly indicates that, inherently, BHLS2 only supplied with consistent $e^{+}e^{-}$ annihilation data, leads precisely to the spacelike expectations supported by the precise NA7 data or the available LQCD input as well. Nevertheless, LQCD information on $F_{K^{\pm}}(s)$ (and $F_{K^{0}}(s)$) would be welcome to better check the open strangeness sector.

Physical Quantities and Model Dependence Effects

BHLS2 has derived evaluations of several physics quantities (see Sections 18, 19 and 20) and some deserve special emphasis. As shown in Table 5, one reaches precise values for the charged pion and kaon charge radii (fm2) :

[TABLE]

our estimate for $<r_{\pi^{\pm}}^{2}>$ is in fairly good agreement with the most recent evaluations [114, 115]. Our evaluation for $<r_{K^{\pm}}^{2}>$ is in accord with the reported values shown in Table 5 and is more precise. The Reference Solution of BHLS2 also leads to :

[TABLE]

where the quoted systematic error reflects different variants of the modeling.

Table 6 collects the values of the BHLS2 model parameters obtained under several fit conditions and discussed in Section 20. An interesting numerical correlation is observed between the fundamental HLS parameter $a$ and the BKY breaking parameter $z_{V}$ within the BHLS2 framework. This can be approximated by $az_{V}\simeq 2.6$; good fit qualities are thus observed with $a$ varying over a large interval, namely from 1.6 to 2.6. The rightmost two data columns in Table 6, which collect the fit values of a large sample of breaking parameters in two extreme cases ($a\simeq 1.6$ and $a\simeq 2.6$), also show correlatedly important changes in some of the other breaking parameter values. Table 6 also shows that the fit qualities reached in all cases are almost identical.

One can herefrom conclude that the various model parameter sets efficiently conspire to provide almost identical cross-sections and form factors which are, actually, the real observables feeding BHLS2. This has been substantiated by comparing the running of BHLS2 fed with our consistent set of data samples under the various conditions underlying the fits reported in Table 6.

A motivated piece of information to detect parametrization effects are the integrals $a_{\mu}(\sqrt{s}\leq 1.05$ GeV$)$ given in Table 7. Running the BS variant (excluding the $\tau$ data) provides $a_{\mu}(BS,\sqrt{s}\leq 1.05$ GeV$)=(571.84\pm 1.06)\times 10^{-10}$, while running the RS variant running on the same data (with the Primordial Mixing angles free) returns $a_{\mu}(RS,\sqrt{s}\leq 1.05$ GeV$)=(571.98\pm 1.20)\times 10^{-10}$. So the difference between the various parametrizations happens indeed to be marginal and quite comparable in magnitude to the uncertainties generated by using somewhat conflicting data samples. Nevertheless, external information constraining or relating some of the model parameters could be useful in order to reduce the parameter freedom in fits and so the parameter correlations.

As for the important muon HVP topic : Complementing the BHLS2 contribution by the whole hadronic contribution not covered by the (BHLS2) model and reported in the left-hand side part of Table 10, one gets :

[TABLE]

Its (statistical + model) uncertainty (3.01) is sharply dominated by those of the whole non-HLS contribution : $2.76\times 10^{-10}$ – reported in Table 10 – versus the BHLS2 uncertainty $(1.06\div 1.20)\times 10^{-10}$ (see Table 7).

So, a significant improvement of the uncertainty for $a_{\mu}^{\rm HVP-LO}$ can only follow from an improved accuracy of the hadronic contributions of higher energies, especially those from the $(1.05-2.00)$ GeV region. Therefore, presently, as global fit methods like BHLS2 constrain the central value for $a_{\mu}^{\rm HVP-LO}$(HLS), they do alike for $a_{\mu}^{\rm HVP-LO}$(HLS+non-HLS). The effect of the $\rho^{0}-\gamma$ mixing vs LQCD is shown separately for clarity.

Our estimate for the full muon anomalous moment is :

[TABLE]

which is, at least, at $3.8\sigma$ from the BNL measurement. The systematics are actually upper bounds to a possible bias affecting $a_{\mu}$.

Appendices

Appendix A The Full HLS Non-Anomalous Lagrangian in the New Scheme

For clarity, the full non-anomalous HLS Lagrangian is written :

[TABLE]

and the various pieces will be given just below.

A.1 The ${\cal L}_{A}+a{\cal L}_{V}$ Lagrangian Piece

First, one displays the part of ${\cal L}_{\rm VMD}={\cal L}_{A}+a{\cal L}_{V}$ relevant for the physics we address, especially $e^{+}e^{-}$ annihilations. This is :

[TABLE]

The pseudoscalar fields occurring there are the renormalized ones defined exactly as in [25]. For simplicity, we use ideal vector fields. The $V-\gamma$ couplings are :

[TABLE]

One has discarded the (irrelevant) photon mass term. The vector meson masses are given in Eq. (26). $\Delta_{V}$ and $\Sigma_{V}$ have been defined in Section 3.

In practical use, one should perform the change of field given in Eq. (27) and collect the terms corresponding to each (neutral) vector meson coupling.

A.2 The ${\cal L}_{\tau}$ Lagrangian Piece

The new ${\cal L}_{\tau}$ is given below (after performing an integration by parts to remove the terms with $\partial W^{\pm}(=0)$ :

[TABLE]

where all fields are renormalized. $\delta_{P}$ is related with the PS mixing in the one angle mixing scheme [31]; its definition is reminded in Appendix B below together with those for the $\epsilon_{i}$.

The appearance of $W^{\pm}(\pi^{0}/\eta/\eta^{\prime})\pi^{\pm}$ couplings should be noted; one should also note, about the $W^{\pm}(\eta/\eta^{\prime})\pi^{\pm}$ couplings, that the contributions from ${\cal L}_{A}$ and ${\cal L}_{V}$ are not proportional in contrast with the $W^{\pm}\pi^{0}\pi^{\pm}$ case. Finally, it is also worth noting that the $\rho^{\pm}(\eta/\eta^{\prime})\pi^{\pm}$ couplings in Equation (116) are needed in the transitions amplitudes for $\eta/\eta^{\prime}\rightarrow\pi^{+}\pi^{-}\gamma$ in order to yield the right behavior at the chiral point.

A.3 The ${\cal L}_{p^{4}}$ Lagrangian Piece ${\cal L}_{z_{3}}$

As for the ${\cal L}_{z_{3}}$ Lagrangian, it writes in terms of ideal fields :

[TABLE]

It clearly allows for different $V\gamma$ and $WV$ couplings, especially for the $\rho$ meson.

A.4 The Renormalized Couplings of Vector Mesons to $K\overline{K}$

Let us first remind the coupling of the fully renormalized vector mesons to $\pi^{+}\pi^{-}$; having defined these as $g^{V}_{\pi^{\pm}}=ag\widetilde{g}_{\pi\pm}/2$, one has :

[TABLE]

Using the Lagrangian in Section A.1 and Eq. (38) one can derive the renormalized coupling of vector mesons to $K\overline{K}$ pairs. Defining:

[TABLE]

the couplings involved in form factors and cross-sections are given by $G^{V}_{K^{+}K^{-}}=ag/[4z_{A}]~{}\widetilde{g}_{K\pm}^{V}$. One also has :

[TABLE]

and the full couplings are given by $G^{V}_{K^{0}\overline{K}^{0}}=ag/[4z_{A}]~{}\widetilde{g}_{K0}^{V}$. One cancels out here and in the following the dependence upon $\Delta_{V}$. One should also note that all these couplings depend on $s$, the square energy flowing through the relevant vector line.

A.5 The $V-\gamma$ Transition Loop Terms

As seen in the main text, the $V-\gamma$ transitions amplitudes write :

[TABLE]

for each of $V=\rho^{0}_{R},~{}\omega_{R},~{}\phi_{R}$. In this Section, one displays the expression for the loop terms $\Pi_{V\gamma}(s)$ in terms of the basic loop functions denoted $\Pi_{\pi}^{e}(s)$, $\Pi_{K\pm}^{e}(s)$ and $\Pi_{K^{0}}^{e}(s)$ and defined in Subsection 11.2. For the sake of conciseness, it is convenient to define :

[TABLE]

neglecting the ${\cal O}(\epsilon^{2})$ terms. These couplings do not depend on the ”angles” $\alpha(s)$, $\beta(s)$ and $\gamma(s)$.

Using these definitions and those displayed in the Subsection just above, on finds :

[TABLE]

Expanding Equations (121) in order to rather parametrize in terms of $\Pi^{e}_{K~{}{\rm sum}}$ and $\Pi^{e}_{K~{}{\rm diff}}$ (see Subsection 11.2) allows to exhibit the terms of order ${\cal O}(\epsilon^{2})$ in breakings which can be dropped out.

Appendix B Pseudoscalar Field Renormalization : A Brief Reminder

As already stated, one goes on using unchanged the breaking procedure in the ${\cal L}_{A}$ sector of the HLS model as defined in [25]. The broken ${\cal L}_{A}$ Lagrangian displayed in Eq. (8) implies a first step field redefinition in order to make the pseudoscalar kinetic energy term canonical. Additionally, as BHLS and BHLS2 also include the ’t Hooft determinant terms [32], the pseudoscalar kinetic energy receives an additional piece :

[TABLE]

where $\eta_{0}$ is the singlet PS field and $\lambda$ is a parameter to be fixed. This term imposes a second step renormalization [31] of the ($\pi^{0},~{}\eta,~{}{\eta^{\prime}}$) sector.

The first step renormalization turns out to define the (step one) renormalized pseudoscalar field matrix $P^{R_{1}}$ in term of the bare one $P$ by [25] :

[TABLE]

Some pseudoscalar (PS) fields happen to be fully renormalized after this first step; they are related to their bare partners by :

[TABLE]

These renormalized fields have already been used to derive the Lagrangians given in Appendix A. As for the other PS fields, the first step renormalization leads to the $R1$ renormalized PS fields :

[TABLE]

The second step renormalization required by the ’t Hooft term gives the $R$ renormalized fields [31] :

[TABLE]

where :

[TABLE]

The parameters affecting the PS fields submitted to fit are thus $\Delta_{A}$, $z_{A}$ and $\lambda$ (or alternatively $v$).

The $\pi^{0}-\eta-\eta^{\prime}$ mixing is supplemented by defining the physically observable PS fields in terms of their $R$ renormalized partners by [144, 25] :

[TABLE]

which exhibits three more parameters to be determined by fit. The singlet-octet mixing angle $\theta_{P}$ occurring here is expected in the range $(-10^{\circ}\div-15^{\circ})$ and is also algebraically related to the usual $\theta_{8}$ [31]. In order to express the forthcoming $VVP$ and $AVP$ coupling constants, the mixing angle $\delta_{P}$ defined by [25] :

[TABLE]

happens to be more appropriate than $\theta_{P}$. On the other hand, it is also useful to define the combinations [25] :

[TABLE]

where $v$ is the nonet symmetry breaking parameter defined in Equations(127). $x$, $x^{\prime}$ and $x^{\prime\prime}$ reflect the various ways by which nonet symmetry breaking occurs in the PS sector of broken HLS Lagrangians.

Appendix C The $VVP$ Lagrangian in Terms of Ideal Vector Fields

The $\pi^{0}$, $\eta$ and ${\eta^{\prime}}$ pieces of the ${\cal L}_{VVP}$ Lagrangian are given separately just below. If the PS fields occurring in this Appendix are renormalized by the procedure reminded in Appendix B, it is simpler to express the various $VVP$ Lagrangian pieces in terms of the ideal vector fields. The derivation of the couplings involving the renormalized fields is given in the main text.

We split up the $VVP$ Lagrangian into the pieces involving the renormalized $\pi^{0}$, $\pi^{\pm}$, $\eta$ and ${\eta^{\prime}}$ pieces, below. As for the $\pi^{0}$, we have :

[TABLE]

where :

[TABLE]

The $VV\eta$ Lagrangian is given by :

[TABLE]

with :

[TABLE]

The $VV{\eta^{\prime}}$ Lagrangian is given by :

[TABLE]

with :

[TABLE]

Finally, it is worth extracting out from Eq. (C) the part involving a charged $\rho$ meson for illustrative purposes. This write :

[TABLE]

with :

[TABLE]

Appendix D The $AAP$ and $APPP$ Lagrangians in the HLS Framework

For the reader’s convenience, we have found appropriate to remind the expressions for the $AAP$ Lagrangian and for the part of the $APPP$ Lagrangian relevant for the annihilation process $e^{+}e^{-}\rightarrow\pi^{0}\pi^{+}\pi^{-}$. These have been derived in [25]. One first have :

[TABLE]

with :

[TABLE]

using the renormalized PS fields defined in Appendix B and their specific parameters.

The pion sector of the $APPP$ Lagrangian is [25] :

[TABLE]

with

[TABLE]

and :

[TABLE]

One may prefer reexpressing this formula by :

[TABLE]

using the angle $\delta_{P}$ defined by Equations (129).

Appendix E The $VPPP$ Lagrangian and its Renormalization

The $VPPP$ anomalous Lagrangian in terms of ideal vector fields can be written :

[TABLE]

where one has limited oneself to display the $P_{0}\pi^{+}\pi^{-}$ sector. The leading terms of these couplings are :

[TABLE]

[TABLE]

[TABLE]

One may note, once again the clarity reached by using $\delta_{P}$ instead of $\theta_{P}$.

The couplings to the renormalized fields can be derived from those to the ideal fields by means of the ${\cal R}(s)$ matrix defined in the body of the text :

[TABLE]

the elements of the $g^{0}_{P_{0}}$ vectors being defined in Equations (145), Equations (146) or Equations (147) for resp. the $P_{0}=~{}\pi^{0}$ or $~{}\eta,~{}{\eta^{\prime}}$ mesons.

For the couplings of the $\phi$ meson, at leading order in breaking, one gets :

[TABLE]

One clearly sees that, at the $\phi$ mass, one obtains a non-vanishing direct coupling of the $\phi$ to 3 pions, additionally $s$-dependent.