Redetermination of the strong-interaction width in pionic hydrogen

The hadronic width of the ground state of pionic hydrogen has been redetermined by X-ray spectroscopy to be Γ1sπH=(856±16stat±22sys)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varGamma ^{\pi \mathrm {H}}_{1s}=(856\,\pm \,16_\mathrm{stat}\,\pm \,22_\mathrm{sys})$$\end{document} meV. The experiment was performed at the high-intensity low-energy pion beam of the Paul Scherrer Institute by using the cyclotron trap and a high-resolution Bragg spectrometer with spherically bent crystals. Coulomb de-excitation was studied in detail by comparing its influence on the line shape by measuring the three different transitions Kα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}, Kβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}, and Kγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} at various hydrogen densities. The pion-nucleon scattering lengths and other physical quantities extracted from pionic-atom data are in good agreement with the results obtained from pion-nucleon and nucleon-nucleon scattering experiments and confirm that a consistent picture is achieved for the low-energy pion-nucleon sector with respect to the expectations of chiral perturbation theory.


Introduction
The last decades led to a successful theoretical description of strong-interaction phenomena at low energies within the framework of an effective field theory (EFT), Chiral Perturbation Theory (χ PT). It exploits the chiral symmetry of the QCD Lagrangian in an ideal world of the massless light quarks u, d, and even s (chiral limit) [1]. In the real world, the chiral symmetry is explicitly broken because of the finite quark masses.
The small current quark masses of 2, 5, and 95 MeV/c 2 for the u, d, and s quarks [2], however, are unable to explain the masses of the hadron multiplets. The masses of the pseudoscalar mesons require both finite current quark mass values and a non-vanishing chiral condensate [3]. However, as the current quark masses are small compared to the QCD scale corresponding to about the nucleon mass, the zero mass limit is still closely fulfilled. This especially holds true for pions exemplified by the ratio (m π /m p ) 2 ≈ 2% 1. Thus, χ PT allows a perturbative approach, where underlying symme-try properties of the QCD Lagrangian manifest themselves in observables as scattering lengths by means of low-energy theorems [1,[4][5][6][7].
A chiral expansion ordered by counting the powers of (small) momenta, the quark-mass difference (m d −m u ) when restricting to SU (2), and the fine structure constant α allows one to include on the same footing strong isospin-breaking effects resulting from the u and d quark mass difference and those of electromagnetic origin. The unknown structure of the effective interaction at short distances arising from the QCD Lagrangian is parametrized by so-called low-energy constants (LECs), which must be taken from experiment as long as results from lattice-QCD calculations are not available.
According to its origin, χ PT works best for the lightest quarks u and d as combined in the lightest hadron, the pion, and for the description of the pion-pion interaction [8,9] but requires a demanding experimental approach [10,11]. On the other hand, it has been shown that such an approach can also be successfully applied to the meson-nucleon case, then denoted as Heavy Baryon Chiral Perturbation Theory (HBχ PT) [12][13][14][15]. The theoretical progress in HBχ PT has reached already the 3rd and 4th order in the chiral expansion for the scattering amplitudes [16][17][18].
In view of the theoretical achievements, the extraction of scattering lengths from experiment is desirable at a level of 1% or better. For such studies, hadronic atoms provide an ideal laboratory [19] because normalisation problems that are inherent to partial-wave analyses when using scattering data do not occur. The strong interaction leads to a level shift and a level broadening Γ in atomic systems formed by negatively charged hadrons like pions and the nucleus. As atomic binding energies are in the few keV range, the measurement of and Γ constitutes a scattering experiment at threshold.
The pion-nucleon scattering lengths can be extracted from the hadronic shift π H 1s and broadening Γ π H 1s of the atomic ground state 1s in pionic hydrogen (π H) and, after applying 3-body corrections, also from the hadronic shift π D 1s in pionic deuterium (π D) (see Sect. 2). For recent reviews within the framework of χ PT see [20,21].
Experimentally, 1s and Γ 1s are accessible via the (np − 1s) X-ray transitions. The energy shift 1s is of the order of a few eV and the broadening Γ 1s ≈ 1 eV, whereas the (np − 1s) X-ray energies are around 3 keV (see below and Table 2). In particular, the precise measurement of such a small line broadening requires high-resolution devices like crystal spectrometers, which in turn need strong X-ray sources. For this reason, precision results were not attainable before high-intensity and low-energy pion beams became available. First experiments of this kind achieved Γ π H 1s = (970 ± 112) meV [22,23] and (868 ± 55) meV [24,25], which corresponds to a precision of 12% and 7%, respectively.
The result described here completes a long series of experiments [26] aiming at an improved determination of the strong-interaction effects both in pionic hydrogen and deuterium. The results for the shifts, π H 1s = (7.0858 ± 0.0096) eV and (1) π D 1s = (−2.356 ± 0.031) eV, (2) have been published earlier in [27][28][29]. In this paper, the sign convention for the shift is ≡ E X − E QED , which corresponds to the change of the X-ray transition energy E X compared to the pure electromagnetic value E QED . It is the negative of the change of total energy of the pionic atom due to the strong interaction.
In order to improve the accuracy for Γ π H 1s , a specifically developed cyclotron trap was used together with a Johanntype spectrometer equipped with spherically bent Bragg crystals as well as a large area array of charge-coupled devices (CCDs), and an especially tailored concrete shielding. The experiment benefited from a 4-fold improved beam intensity provided at the Paul Scherrer Institute (Villigen, Switzerland). Statistics and peak-to-background ratio were improved by a factor of about 30 and 50, respectively, in comparison to [24,25]. In addition, a new calibration method was introduced: narrow fluorescence X-rays emitted from highly charged ions were used to determine the spectrometer resolution function with maximum possible accuracy for the given experimental approach (see Sect. 4).
The main difficulty in the extraction of Γ 1s is to quantify the Doppler broadening effect on the line shape of the X-ray transitions. As discussed in Sect. 3, the kinetic energy distribution of the π H atoms changes during the atomic cascade due to collisional processes, which are density dependent. As a result, the influence of the Doppler effect depends on the initial state np and the target density. Therefore, three different ground-state transitions, π H(2 p − 1s), π H(3 p − 1s), and π H(4 p − 1s) (Kα, Kβ, and Kγ ), have been measured at hydrogen densities between 3.9 bar and liquid (LH 2 ), in order to ensure that the Doppler effect is consistently taken into account in the analysis of the line shapes.
The analysis of the spectra as described in Sect. 5 uses a model for the kinetic energy distribution which is derived from a preceding experiment with muonic hydrogen and is introduced in Sect. 5.1. It is applied both in a χ 2 analysis (frequentist analysis Sect. 5.2) and in a Bayesian analysis (Sect. 5.3).
The use of a Bayesian approach was especially motivated by the case of pionic deuterium. There, the occurrence of a bias in a frequentist method [30], influencing the extracted values for Γ π D 1s , has been studied in detail [29]. Such a bias stems from the -principally unknown -difference of the probability distribution of the data itself and the one assumed for the model. An estimate for the bias of the parameter Γ π D 1s was found by applying the same model as used for the analysis to a series of Monte-Carlo simulations generating the X-ray line shape and by comparing input and output parameters. Such a lengthy procedure may be circumvented by using Bayesian methods, which are supposed to be free from bias effects in parameter estimation [31]. However, they suffer from different drawbacks, e. g. the influence of the choice of a prior distribution.
A comparison of the results of the two analysis methods establishes an independent consistency check (Sect. 5.4) and is in line with recent assessments recommending to use both approaches as good practice in evaluating data. Such a comparison of Bayesian and frequentist approaches has been already applied for the study of the twin-system muonic hydrogen [32,33], which completed this series of experiments. The features of the two evaluation methods are outlined in apps. A.1 and A.2 of Ref. [33].
A discussion of the different sources of systematic uncertainties refers to the influence of the different cuts in the treatment of the spectra (Sect. 5.5.1) and investigates the systematic error which is inherent to the model used for the kinetic energy distribution (Sects. 5.5.2 and 5.5.3). In Sect. 6 the result of our experiment and its implications are presented.

Strong-interaction effects
The hadronic broadening of the π H ground state results from the s-wave charge-exchange reaction π − p → π 0 n and radiative capture π − p → γ n. Within a non-relativistic effective field-theoretical approach, the scattering length a cex π − p associated to the charge exchange reaction is related to Γ π H 1s by [35,36] Here, α is the fine structure constant and m π is the mass of the charged pion. μ H = m π /(1 + m π m p ) is the reduced mass of the π − p system, and q 0 is the centre-of-mass momentum of the π 0 . The factor 1 + 1 P takes into account the fraction of radiative capture, where the Panofsky ratio for pion capture at rest was measured to be P = σ (π − p→π 0 n) σ (π − p→γ n) = 1.546 ± 0.009 [37]. The π 0 n elastic scattering length a π 0 n for the virtual intermediate state is well approximated by the isoscalar scattering length a + (see below) [38]. The correc-tion term δ vac = 0.48% accounts for the effect of the vacuum polarization on the wave function [39]. Its uncertainty is assumed to be negligibly small compared to the experimental uncertainties [20].
In the limit of isospin conservation, the scattering lengths of all reactions π N → π N can be expressed by the isoscalar and isovector s-wave scattering lengths a + and a − , which are given in terms of the elastic reactions π − p → π − p and π + p → π + p by In terms of the two possible isospin combinations I = 1/2 and I = 3/2, a + and a − have the form a + = 1 3 a 1/2 + 2a 3/2 and (6) The relations between the scattering lengths in the isospin conservation limit, a ± , and the ones used in Eqs. (3) and (4) are given by The corrections a have been calculated at various levels and approaches during past decades [16,36,[41][42][43][44][45][46][47][48][49]. Numerical values and, in particular, the uncertainties of the various a were subject to significant changes resulting in nonsatisfactory conclusions on the scattering lengths a + and a − in the year 1998 when this experiment was proposed. The most recent determinations for a have been worked out within the framework of χ PT in next-to-leading order (NLO) [21,38,50,51]: These corrections constitute the isospin violating contributions of strong (m u = m d ) and electromagnetic origin and are assumed to be complete within the given order of the chiral expansion. The large errors are due to the poor knowledge of LECs appearing in the chiral expansion. For π D, the modified Deser formula has the form [52] − π D 1s − i where the related pion-deuteron scattering length a π D is complex. Here, μ D is the reduced mass of the π − d system and δ vac D = 0.51% [39]. The real part of the scattering length Re a π D gives the leading contribution to π D 1s and can be expressed as sum of two-and three-body terms a (2) π D and a (3) π D Re a π D = a (2) π D + a The two-body term a (2) π D = a π − p→π − p +a π − n→π − n is essentially 2 · a + and, thus, provides an important constraint on the results for a + and a − as obtained from π H 1s and Γ π H 1s . Though a (2) π D involves substantial isospin-violating corrections and a (3) π D depends in leading order on (a − ) 2 , the corrections are well under control. The theoretical framework is introduced and discussed in detail in Refs. [20,38,[51][52][53][54][55][56][57][58].
The imaginary part Im a π D , which gives the leading contribution to Γ π D 1s , measures the transition strength of swave pions on an isoscalar nucleon-nucleon pair π N N ↔ N N [59] and is an independent quantity not related to the scattering lengths a + and a − [29,60,61].
The following relations hold in the isospin conservation limit: The deviation from the identity Eq. (14) quantifies the isospin-breaking effects. In terms of measurable quantities, it can be described by the ratio and is expected to be of the order of a few % [38]. The isovector scattering length a − and correspondingly Γ π H 1s are connected to other important quantities. These include the pion-nucleon coupling constant, the Goldberger-Treiman discrepancy, and the pion-production amplitude γ n → π − p at threshold. The Goldberger-Miyazawa-Oehme (GMO) sum rule [62], which can be written in the form [51,63] relates the scattering lengths a π − p and a π + p to the charged pion-nucleon coupling constant g 2 c /4π . The GMO sum rule essentially constitutes a forward dispersion relation containing the integral J ranging over all pion momenta which, employing the optical theorem, can be expressed as the difference σ tot π − p − σ tot π + p of the Coulomb-corrected total cross sections. After applying appropriate corrections (see Eqs. (25) and (30)), the scattering length difference a π − p − a π + p can be related to a π − p and a cex π − p , i. e. to the hadronic shift π H 1s and to the level broadening Γ π H 1s in pionic hydrogen [51,63,64].
In the chiral limit, the Goldberger-Treiman relation m N g A = f π g c holds true [65], where g A is the axial-vector coupling constant in β-decay, f π is the pion decay constant, and the average nucleon mass m N = m p +m n 2 . The deviation from the Goldberger-Treiman relation, the Goldberger-Treiman discrepancy [66], can be written as (17) with the pion-nucleon coupling constant g c obtained, e.g. from the GMO sum rule. It measures the explicit chiral symmetry breaking by the quark masses and is expected to be small in the standard chiral symmetry breaking (≈ 1.5%) because of the small u and d quark masses [50,[67][68][69]. Furthermore, the Kroll-Ruderman theorem connects g c with s-wave pion photoproduction γ N → π N [70]. For π − production, the threshold electric dipole amplitude reads Here, the chiral corrections (h.o.) to the leading-order term have been calculated to third order in the framework of χ PT [71][72][73]. The isoscalar π N amplitude a + is related the π N σterm σ π N [74], which encodes the contribution of the lightest quark doublet to the masses of the nucleons. Furthermore, its connection to the scalar coupling of dark matter [75][76][77] or μ → e lepton flavor violation [78,79] is of interest in the context of new physics research.
The long-standing debate about the extraction of σ π N is, among other reasons, due to the difficulties arising from subthreshold extrapolation of amplitudes [80,81]. In the mean-time, substantial progress has been made in the development of concepts for such a procedure. Within the framework of dispersion relations, a direct correlation including chiral corrections between σ π N and the π N s-wave scattering lengths has been established [21,82].
Most analyses of recent pionic-atom [83] and low-energy scattering data lead to σ π N ≈ 60 MeV [84,85]. Selection of one particular data set (σ π N = (44 ± 12) MeV [86]) or the choice of a particular partial-wave analysis [87] lead to non-decisive or scheme-dependent results. In contrast, lattice calculations [88][89][90][91] typically yield such small values being in tension by 3 -4σ with most data-based results [92]. One possibility to support solutions of this puzzle, is to confirm the reliability of experimental information, e. g. the pionatom data base.

Pionic hydrogen atomic cascade
After slowing down to a few eV, a heavy negatively charged particle (x − = μ − , π − , K − , orp), while passing through the atomic or molecular electron cloud and losing further its kinetic energy, is eventually captured in a highly excited bound state where an atomic quantum cascade starts. These initial states are comparable in size with the hydrogen atom and have the principle quantum number n c ≈ √ μ x /m e , where μ x denotes the x − p reduced mass [93,94]. In the case of pionic hydrogen, n c ≈ 16. The initial distributions in n c and the angular momentum quantum number have been calculated to be broad distributions [95][96][97][98].
Before the pions undergo the radiative transitions (np → 1s) that can be observed experimentally, different collisional processes occur in the atomic cascade ( Fig. 1). While these collisional processes are not directly observable, they change the kinetic energy of the π H atom and, therefore, influence the line shape of the measured radiative transitions due to Doppler broadening. For this reason, a good understanding of the atomic cascade is crucial to a proper evaluation of the measured data.
Traditionally, the collisional processes (π − p) nl + H → (π − p) n l + H were classified according to their main effect in the atomic cascade (Table 1). In particular, during the collisions Stark mixing causes transitions with n = n and l = l that can produce ns states where nuclear absorption takes place with high probability [93]. The nuclear absorption via the Stark mixing significantly reduces the intensities of the radiative transitions, especially at high hydrogen density.
The Coulomb transitions [99,100] correspond to the collisions with n < n where the transition energy is shared between colliding particles, which have comparable masses, thus resulting in a significant acceleration of the π H system. Elastic scattering, where the principal quantum number remains unchanged (n = n), plays an important role in the In this experiment, the 3 most intense X-ray lines π H(2 p − 1s), π H(3 p − 1s), and π H(4 p − 1s) have been studied. Because strong-interaction effects in the np states are negligibly small and radiative widths Γ X Γ 1s (see Table 2), the natural line width is practically equal to the 1s level broadening Γ 1s due to the reactions (π − p) 1s → π 0 + n and γ + n. The level shift 1s is caused in leading order by elastic scattering π − p → π − p with the π − H system being in the atomic ground state deceleration of the pionic hydrogen, since the kinetic energy of (π − p) nl is usually much higher than the kinetic energy of the H 2 molecules.
The hydrogen molecular structure can be important for certain collision channels, especially at high n when the size of (π − p) nl is not small in comparison with the one of the hydrogen atom and when the degrees of freedom of H 2 molecules play a role. In the earlier cascade calculations, different approximations were used to model the above mentioned collisional processes [93,101]. The modern calculation of the collisional rates treats them as different final states of the same scattering problem by using a unified framework that improves significantly the accuracy of the theoretical results [102].
The collisional processes dominate the upper and the intermediate part of the atomic cascade, because even at modest densities, the corresponding reaction rates exceed the radiative de-excitation rates. Only for small n, X-ray emission starts to compete and finally dominates in the very last steps of the cascade.
The rates of collisional processes scale with hydrogen density φ. As a result, the cascade time rapidly decreases with the increase of the target density. The cascade times are calculated to be between 10 −9 to 10 −11 s in the range φ = 10 −3 to 1 [103], where φ = 1 corresponds to the density of liquid hydrogen. Therefore, almost all pions reach an atomic ns The nuclear reaction channels and the pion decay terminate the atomic cascade

Fig. 2
Absolute K X-ray yields in pionic hydrogen vs. target density φ (normalized to liquid hydrogen density). The theoretical yields are obtained from a calculation within the framework of ESCM [106]. The experimental data (points with error bars) are from [107] state, where the pion wave function has a large overlap with the nucleus, and the nuclear reactions of charge exchange and radiative capture occur on the femtosecond scale. Hence, the pions reaching s-states disappear from the de-excitation cascade in contrast to the muonic hydrogen case where muon capture from excited atomic states is negligibly small. Consequently, the X-ray yields in pionic hydrogen are about one order of magnitude smaller than the ones in muonic hydrogen. Absolute yields are in the few % range for the densities covered by this experiment (Fig. 2). The competition between the acceleration (Coulomb deexcitation) and the deceleration (elastic and inelastic scattering) leads to the evolution of the kinetic energy distribution during the atomic cascade, and the corresponding energy distributions can be studied experimentally for several particular states. Experimentally, the Doppler broadening was observed for the first time in the time-of-flight distribution of neutrons produced in the charge-exchange reaction π − p → π 0 n after pion capture at rest [104,105] and studied later in detail in connection with a precision determination of the π − −π 0 mass difference [108]. Predominantly, Coulomb de-excitation involves one atom of the H 2 molecule, i.e. the collisions (π − p) nl + H must be considered [109][110][111]. Therefore, the π H atom gains about one half of the de-excitation energy.
In contrast to Coulomb de-excitation, the kinetic energy transfer to the π H atom during the initial capture and from external Auger effect is negligibly small, because the energy difference is carried by the emitted electron. Similarly, the recoil effects from X-ray emission, the Stark collisions (at small scattering angles), and molecular formation need not be considered for a noticeable contribution to Doppler broadening. Furthermore, radiative decay from molecular states causing a line broadening due to satellite transitions is expected to be negligibly small [112][113][114][115]. No evidence for such transitions could be identified from the study of the density dependence of the X-ray energy [27,33].
At densities of φ ≈ 10 −2 , also the contributions of induced Stark width, collisional broadening, and thermal motion add up at most to a negligible value of ≈ 3 meV [33,109]. At φ = 1, an upper limit of 30 meV is assumed for the induced Stark width [109].
Elastic and inelastic scattering are the only efficient deceleration processes, which compete with the acceleration caused by Coulomb de-excitation. This leads to a complex kinetic energy distribution at the time of the Xray emission (Fig. 3). Consequently, the state-dependent development of the kinetic energy is essential for modern cascade calculations (extended standard cascade model (ESCM)) [101,106,[116][117][118][119].

Set-up
The experiment was performed in the π E5 area at the meson factory of the Paul Scherrer Institute (Fig. 4). The X-ray source has been provided by the cyclotron trap. A weakly focussing magnetic field, produced by a superconducting split-coil magnet, guided the pion beam to the centre of the trap [120,121]. In the centre was installed a thin-walled Kap- Table 2 Energies E X , radiative widths Γ X , and Bragg angles Bragg angles for π H lines are taken from the QED energies as calculated by [124], combined with the hadronic shift as given in [27]. The lattice constants are taken to be 2d    [122] allows the simultaneous measurement of an energy interval when using a correspondingly large position-sensitive detector. By using spherically bent crystals a partial vertical focusing is achieved.
The spectrometer consists of a Bragg crystal mounted inside of a vacuum system in order to avoid absorption losses on the X-rays' path [123]. The X-rays were recorded in a 3 × 2 array (ver tical × hori zontal) of charge-coupled devices (CCDs) having altogether a sensitive area of 72 × 48 mm 2 [127,128]. The image area of each of the six CCDs contains 610 × 602 = 367220 pixels of 40 × 40 µm 2 .
The detector plane is oriented perpendicular to the line crystal centre to CCD array centre. The focussing condition (distance crystal-detector) is given by R c · sin Θ B , where Θ B is the Bragg angle. The dispersion in the detector plane reads Table 2). In order to keep X-ray absorption small, the exit window of the target cell was made from a 7.5 μm thick Kapton foil. For a diameter of about 50 mm, adapted to the detector size in the direction of dispersion, such a window safely allows a pressure difference of up to about 1 bar only. Therefore, the hydrogen density was adjusted by temperature variation by means of a cooling finger. In this way, a density range from 4 bar equivalent (at T = 20 • C) to liquid (corresponding to 787 bar) was covered (Table 3). A minimum density equivalent to 10 bar is necessary to stop a sufficiently large fraction of pions inside the target cell in order to achieve the statistics needed for the precise determination of the line shape.
The experimental approach is identical to the one used for the most recent determinations of the hadronic shift in pionic hydrogen, the strong-interaction effects in pionic deuterium, and for the line-shape measurement of muonic hydrogen. Accelerator conditions, experiment set-up and measurement procedure are described in more detail in the corresponding publications [27,29,33].

Spectrometer resolution
For the X-ray diffraction, only silicon and quartz crystals provide the energy resolution needed in the considered energy range. In this experiment, one quartz and one silicon crystal were used cut along the (101) and (111) plane, respectively. Crystal disks of 0.2 and 0.3 mm thickness and of 100 mm in diameter were spherically bent to radii measured to be R c = (2980.6 ± 0.4) and (2982.2 ± 0.3) mm, respectively [129]. The area of the disks was limited twofold: (i) by a spherical aperture to a diameter of 95 mm in order to avoid edge effects and (ii) to ± 30 mm in the direction of dispersion to keep Johann broadening small [130].
The measured total line shape is composed of contributions from the spectrometer response, Doppler broadening due to Coulomb de-excitation, and the natural line width of the X-ray transition. As both the contributions from Doppler broadening and the natural line width are to be determined, the knowledge of the spectrometer response is essential to a level of accuracy that its influence on the results for Γ 1s remains negligible. The spectrometer response itself is composed of the rocking curve reflecting the intrinsic resolution of the crystal material, the imaging properties resulting from the bent-crystal set-up, and possible imperfections from setup geometry and crystal mounting.
The spectrometer response was determined by using Xrays emitted from highly-charged ions (Fig. 5). The ions were produced by means of discharge in a dedicated electron-cyclotron ion resonance trap (ECRIT) installed in the centre of the cyclotron trap [131]. In this way, the count rates up to 10000 per hour could be achieved for calibration lines, which is necessary for a precise mapping of the reflection properties and, in particular, the tails of the response function. The resolution was determined under conditions as close as possible to the π H experiment. The crystal spectrometer was directly attached to the cyclotron trap with the embedded ECRIT in the same geometrical set-up as realized for the experiment [132]. The plasma cloud provided an extended X-ray source similar to the dimensions of the hydrogen target cell.
The M1 X-ray lines from helium-like argon, chlorine, and sulphur have energies of 3104, 2765, and 2430 eV, respectively, which are very close to the energies of the three π H lines measured in this experiment ( Table 2). M1 transitions 1s2s 3 S 1 → 1s 2 1 S 0 originate from metastable initial states and have natural line widths of about 10 neV [133]. The Doppler effect on the transitions has been measured to be 15−80 meV for Ar, which corresponds to a kinetic energy of the ions of 1 − 30 eV [134][135][136]. Such a Doppler broadening contributes at maximum with a few meV to the spectrometer resolution (Table 2). Therefore, the measured line shape is taken as the direct response of the apparatus.
As a main result, the crystal response function has been measured for the three different energies with an accuracy of better than 1% in terms of FWHM. In particular, the precise knowledge of the asymmetric line shape assures that the contribution of spectrometer response to the uncertainty of Γ 1s is negligibly small within the statistics achieved for the π H spectra.
The measurements were compared to a modelling of the resolution and found to be close to the theoretical limit. Here, the resolution function at a given energy was constructed from the rocking curve as calculated from the dynamical theory of diffraction for a perfect flat crystal (Fig. 6) convoluted with the imaging properties of the whole experimental set-up by means of a Monte-Carlo ray-tracing code. The result constitutes the ideal response with a corresponding line width E ideal (FWHM) for a given set-up [123]. The intrinsic resolution (rocking curve) was calculated with the code XOP [139,140].
Evaluation of the results showed that the line shape measured with the ECRIT set-up is already very close to the ideal response. It is sufficient to introduce an additional Gaussian contribution Γ G to take into account imperfections of the crystal material and mounting for a good description of the line shape. The values for Γ G were found from fits to the measured spectra of the sulphur, chlorine, and argon M1 transitions [137,138]. The corresponding values and uncertainties The H 2 target densities are given as equivalent pressure p eq (in bar corresponding to a temperature of 20 • C (NTP)) and as fraction φ of liquid hydrogen density. Counts are results of fits to the line shape. The accelerator current varied between 1.2 and 2 mA, which corresponds to pion fluxes arriving in the π E5 area of about (3 − 5) · 10 9 π − /s  Entries are given for 1/10 of the detector's pixel size or 4 μm for 10 6 Monte-Carlo events each of the total width E exp are given in Table 2. The response functions including the Gaussian contribution are shown in Fig. 7.
As described in Refs. [137,138,141], the data taken with the ECRIT set-up were also used to determine precisely the crystal's bending radii and resolutions for various apertures, to identify possible miscuts and their orientations, as well as to validate the X-ray tracking code which includes the aberrations to all orders for the assumed geometry [123].

Data processing
The Bragg reflection of a mono-energetic X-ray line generates a cone-section like hit pattern in the detector plane. The position resolution is given by the granularity of the CCD array. The low-energetic X-rays create charge only in one or two adjacent pixels, because of the small photo-electron range of a few μm. A substantial background reduction has been achieved by a pattern analysis followed by the rejection of larger clusters stemming from beam-induced background, which has been further improved by applying an appropriate window in the pixel charge spectrum (Fig. 8).
The detector resolution in terms of the collected charge has been determined by a Gaussian fit to the various π H peaks and found to be between 120 and 300 eV (FWHM) for the individual CCDs. Of the events attributed to the π H lines, about 75% are single and about 25% two-pixel clusters. Defect (horizontal) CCD lines (0−4 per CCD) and additional individual hot pixels (≈ 10 per CCD) were masked by software.
The curvature of the reflection has been corrected by means of a parabolic fit before projecting onto the axis of dispersion (Fig. 9). The projection constitutes the energy spectra to be analysed for the Doppler and the hadronic broadening (Fig. 10). The uncertainty of the curvature correction influences the line width by less than 0.1 pixels, which corresponds to 3 − 7 meV for the various transitions. Such an

Data analysis
The measured line profile L P (see examples for the three transitions in Fig. 10) is formed by a convolution of the natural line width (given by a Lorentzian L with width Γ π H 1s ), the spectrometer response R, and a sum of the Doppler profiles D i which correspond to the various kinetic-energy components originating from the Coulomb de-excitation. By exploiting distributivity, it can be written as For the analysis of the spectra, both a frequentist and a Bayesian approach were used. Apart from trivial factors like intensity, position, and background level (assumed to be constant), the important parameters for the analysis are -the hadronic width Γ π H 1s and -the kinetic energy distribution given by the model chosen for i D i .

Modelling of the kinetic energy distribution
The model used for the kinetic energy distribution is based on the results obtained from the muonic hydrogen experiment, where the line-shape model reduces to [32,33]. The distribution is built up by a number of narrow and discrete energy regions, i. e. energy intervals of a few eV width ("components" (Fig. 11)). The evaluation of muonic hydrogen showed that the quality of the description of the line shape does not depend on the width of the intervals representing the components of the kinetic-energy distribution, when the interval widths were varied between 1 and 8 eV.
Therefore, the low-energy component was fixed to an energy interval of 0 − 1 or 0 − 2 eV and for the width of the high-energetic components always a few eV was used. The central energies of the high-energy components and the relative intensities of all components were free parameters. The relative intensities were normalized to add up to 100%. The construction of the corresponding Doppler spectra is described in detail in Ref. [137].

Frequentist data analysis
The spectra were fitted using the MINUIT package [142]. In a first Poisson likelihood χ 2 optimization [143,144] it turned out that at least a low-energy component and in addition one Fig. 10 The spectra of the π H(2 − 1) (A1), π H(3 − 1) (B2), and π H(4 − 1) transitions (C1) measured at 10 bar equivalent target density  The relative intensities of low-and high-energy component at minimal χ 2 are shown in Table 4 for the most relevant transitions measured.
Such two-component analyses provided starting parameters for the 3-component fits including the second highenergy component. However, even in the case of the (2 − 1) transition, a 3-component fit does not significantly improve the value of χ 2 . The weighted averages for Γ π H 1s as obtained from the 2-and 3-component fits are almost identical ( Table 7).
The 2-component analyses of two low-statistics spectra at the extreme densities 3.9 bar and LH 2 (measurements B1 and B4 (Table 3)) yielded Γ π H 1s = (890 ± 150) meV and (815 ± 109) meV in good agreement with the result obtained for the other spectra. As the upper limit for the induced Stark width of 30 meV at LH 2 density [109] is much smaller than the statistical error, the corresponding correction is omitted for the measurement B4.
For the kinetic energy models, the results corresponding to the minimal χ 2 were used. The errors reflect statistical uncertainties only As discussed in the earlier analysis of pionic hydrogen and deuterium data from this series of experiments [29,61,145], the possible occurrence of a bias in the parameter determination must be considered [30], in particular, because the kinetic energy distribution can only roughly be approximated by a simple few-component model. Based on these results, a series of Monte-Carlo simulations has been performed for each case (400-800 per spectrum) in order to estimate the bias correction. The values are shown in Table 5.

Bayesian data analysis
The Bayesian analysis, though the evidence for a second high-energy component is weak, always used the 3component model for the (2 p −1s), (3 p −1s), and (4 p −1s) transitions.
The evaluation procedure follows a suggestion by R. Trotta [146]. It uses a number of models inside the range of physically reasonable assumptions. Each of the models is subject to a Bayesian evaluation which leads to a value for the so-called Bayesian evidence, which is proportional to the probability of the model itself, as well as for a posterior distribution.
In the present case, 190 different 3-component models have been used for the (2 − 1) transition. For the study of the (3 − 1) and the (4 − 1) transitions, a number of 78 and 36 models, respectively, turned out to be sufficient because Table 6 The schematics of kinetic energy distributions assumed for the (2 − 1) transition with one low-energy (LE) and two high-energy components (HE) having a width of 1 eV each In total 190 parameter combinations ("models") are treated in a Bayesian evaluation of the smaller energy range. The distribution in energy of the three components was systematically varied over the physical range permitted, with the exception that the low-energy component was always fixed to the interval 0 − 1 eV. The scheme for the (2 − 1) analysis is given in Table 6. As a test, for the (2 − 1) transition the 190 models were extended from the upper energy of 200 eV up to 300 eV. This resulted in a change in the value of Γ π H 1s of 2 meV only. Here, the advantage of this method comes into effect because for each of the models the value of the Bayesian evidence provides the corresponding weight of the model itself (proportional to its assigned probability) for averaging. In order to do this, the model with maximum evidence was determined and normalized to a value of 1, leading then to lower values for the transformed evidences of all the other models. These lower values are the measure by which the number of points in the posterior distribution for the parameter (here Γ π H 1s ) is reduced for each model. This was done by a Monte Carlo routine. In a final step, the different weighted posterior distributions were summed up. From this final distribution, the values for the mean of Γ π H 1s and its uncertainty as given in Table 7 were found.

Comparison of frequentist and Bayesian results
The results for the width Γ π H 1s from the individual spectra are compared in Table 7. The weighted average values are in good agreement with each other where the associated uncertainty can be regarded as the statistical error of the data available. An estimate for the systematic contribution to the error of the final average value is derived in Sect. 5.5.3.
Though we can conclude from the results given in Table 7 that the assumptions about the bias correction are reasonable, the frequentist fit can suffer from the particular assumptions for the kinetic energy distribution. In view of this uncertainty, the Bayesian analysis with its conservative error estimate is chosen for the final result. In order to estimate systematic effects in the evaluation of the data, the influence of the following parameters has been studied: 1. Window applied to the CCD charge spectrum (see Fig. 8): A variation of the window between 2/3 and 4/3 of FWHM changes the value of Γ π H 1s by less than 50% of its statistical error for the corresponding spectrum. The average variation of Γ π H 1s was found to be The value of Γ π H 1s stabilizes when the interval containing the peak is 15 eV or larger. Smaller intervals resulted in too large values for Γ π H 1s because the interval boundaries then cut into the tails of the line. This could be proven by an evaluation of the (3 − 1) and (4 − 1) transitions by systematically shortening the fit intervals to values smaller than 36 or 44 eV, respectively. 3. Fluctuation of data subsets: The total measuring time for an individual transition spanned from 3 to 25 days. Therefore, subsets of the larger periods were evaluated separately. No time dependence of the result for Γ π H 1s could be found.

Influence of the kinetic energy distribution
It is tempting to use the theoretical prediction for the kinetic energy distributions, like the ones shown in Fig. 3, for i D i in Eq. (19). The corresponding results for Γ π H 1s are in the range 781 − 963 meV. The statistical errors are only ±7 meV for the measurement with the highest statistics (measurement A2) and ±(17−20) meV else (A1, B2, B3, C1, C2). Such small uncertainties are due to the fact that this analysis assumes no error for the kinetic energy distributions. Therefore, we emphasize that these values for Γ π H 1s and their errors do not constitute experimental result.
The above-mentioned range of the results for Γ π H 1s , which is much larger than the individual errors, yields an uncertainty of about 40 meV for the average of the six measurements. We interpret this as an incomplete matching of theoretical description and data.
To quantify this mismatch, results from the muonichydrogen analysis are used, where it was possible to reconstruct the coarse structure of the kinetic energy distribution by de-convolution. This revealed that experiment and cascade theory agree within about 3σ and, in particular, that the intensity of the low-energy component is reproduced with a value of about 20% lower than theoretically predicted. Such a de-convolution turned out to be not feasible for pionic hydrogen because the additional parameter Lorentzian broadening, Γ π H 1s (Eq. (19)), hides any detail.
To estimate the order of magnitude of the mismatch on Γ π H 1s , the strength of the low-energetic component was reduced for the model distributions to values of 80% compared to the values of the theoretical distributions as shown in Fig. 3. For this study, the six spectra (2 for each of the transitions (np − 1s)) with the highest accumulated statistics were used (measurements A1, A2, B2, B3, C1, C2).
In the case of the (2 p − 1s), (3 p − 1s), (4 p − 1s) transition, the change of the strength of the low-energy component reduces the value for Γ π H 1s by 20, 26, and 12 meV, respectively. The differences are consistent with the systematic error derived from the scattering of the results from the cascade-model-free analyses (Table 7).

Estimate of systematic error
From the discussion in Sect. 5.5.2 it can be concluded that details of the kinetic energy distribution, which are not accessible to our analysis, are important for the scattering of the results for Γ π H 1s ( Table 7). The method recommended by the Particle-data-Group is to quantify such a scattering by the x i and σ i are the individual results and their uncertainties, andx the weighted average of the n measurements [2]. The total uncertainty is then given by the average statistical error multiplied by the enhancement factor √ S. One can get an estimate for the systematic error by subtracting quadratically the average statistical error from the total error (Table 7 -systematic error). It turns out, that the range of values obtained in this way is close to the scattering found from the variation of the fraction of the low-energy component of the kinetic energy distribution.

Results
The final result for Γ π H 1s of the experiment series described here, taken from the Bayesian analysis, is The total experimental error obtained by quadratic addition amounts to Γ π H 1s,exp = ± 27 meV or ± 3.1%.
This improved value in comparison with [28] is due to a newer result for the LEC c 1 [21].

Fig. 12
Constraints on the scattering lengthsã + and a − by the hadronic 1s level shift π H 1s in pionic hydrogen, π D 1s in pionic deuterium, and the 1s level broadening Γ π H 1s in pionic hydrogen as obtained within the framework of [21,51,57,147]. The ellipse displays the combined constraints as given in Table 8 (3 rd line) As already mentioned, three experimental results, π H 1s , Γ π H 1s , and π D 1s , are available to determine the two quantities a + and a − . Table 8 compares the results obtained from different data sets: only the pionic hydrogen (1 st line), only the shift data with a − dependent corrections (2 nd line), and after applying a constraint calculation with all three results (3 rd line).
The errors of a + and with it a + + a − (Eqs. (22) and (24)) are dominated by the uncertainty of the chiral corrections, which is mainly due to the poor knowledge of the LEC f 1 that appears in combination with the LEC c 1 and a + itself. Therefore, in theoretical analyses the corresponding combination denotedã + is more convenient, which has a much smaller uncertainty [51,57] (Table 8). The most recent χ PT calculation provides for the differenceã + − a + = (− 6.1 ± 2.5) · 10 −3 m −1 π [21]. In contrast, the uncertainty of a − is dominated by the experimental error.
Given the final values of the three experiments, in view of the larger uncertainty of Γ π H 1s , the values of a + and a − are essentially determined by the results derived from the energy shifts π H 1s and π D 1s . Therefore, discarding the constraint from Γ π H 1s leads to marginal changes only. It is worthwhile to note that a + remains positive by 3 standard deviations. Figure 12 shows the constraints imposed on the scattering lengthsã + and a − by the three pionic-atom results. It constitutes the update of the analyses of [21,51,57,147] for the new result for Γ π H 1s . The π N scattering lengths directly accessible from pionic-atom experiments via Γ π H 1s and π H 1s are a cex π − p and a π − p (Eqs. (3) and (4)). In order to quantify the isospinbreaking corrections via the triangle relation (Eq. (15)), the scattering-length difference a π − p − a π + p is needed. In an approach chosen by [51,63], emphasising the experimental information for π H 1s together with corrections depending amongst others on π D 1s and a − as derived from Γ π H 1s , the scattering-length difference can be written as a π − p − a π + p = 2 · a π − p + a( π D 1s , a − ). (25) Numerically, the scattering length difference is given to 97% by the experimental results for a π − p , i. e. the result is very robust against the input for a( π D 1s , a − ). The correction term worked out in the framework of the chiral expansion has the form [38,51,57] a π D 1s , a − = −2 · ã + + ã + + a π − n − a π + p = (2.1 ± 1.7) · 10 −3 m −1 π .

Related quantities
In order to determine the π N coupling via the GMO sum rule (Eq. (16)), the scattering-length difference a π − p − a π + p must be corrected for virtual-photon contributions [51,85]. Starting from Eq. (25) one gets Alternatively, the scattering length difference can be expressed as Here, the experimental result for Γ π H 1s dominates the numerical value by 99%. With the appropriate correction to include also virtual-photon effects, which depends among others on Using the constraint given by π D 1s , the full information available from pionic hydrogen and deuterium is exploited and, therefore, constitutes the final result of this experiment series (3 rd line) Table 9 Results for pion-nucleon coupling constant g 2 c 4π , the Goldberger-Treiman discrepancy GT , and the threshold amplitude for π − photoproduction E thr 0+ (π − p) The value for the dispersion integral J = 1.073 ± 0.034 mb is adopted from the analysis of [51] and contributes (by quadratic addition) with about 50% to the total error b Using f π = (92.28 ± 0.09) MeV [2] and g A = 1.27641 ± 0.00056 [156] c Errors reflect mainly statistical uncertainties The two sets of values derived from pionic atoms, combining π H 1s with π D 1s and from Γ π H 1s , are compared to the ones obtained from pion-nucleon (π N ) dispersion-relation (DR) applications and more recent analyses of nucleon-nucleon (N N) scattering data. The results for E thr 0+ (π − p) include the chiral correction (− 1.0 ± 0.6) · 10 −3 m −1 π as given by [71] the chiral prediction for R [147], the result is a γ / π − p − a γ / π + p = − √ 2 · a cex π − p + (0.7 ± 1.0) · 10 −3 m −1 π = (176.6 ± 3.0) · 10 −3 m −1 π .
Equations (29) and (31) represent the same physical quantity derived from the two practically independent pionic-atom results: on the one hand from π H 1s together with π D 1s and on the other from Γ π H 1s . They agree with each other within 2σ . This comparison constitutes an important check of consistency for the pionic-atom data. Table 9 shows the corresponding two sets of values for the pion-nucleon coupling constant g 2 c 4π via the GMO sum rule, the Goldberger-Treiman discrepancy GT (Eq. (17)), and the threshold electric dipole amplitude E thr 0+ (π − p) (Eq. (18)). Within the experimental and theoretical uncertainties, the slight tension between the two sets of results does not constitute any contradiction. The results show also good agreement with more recent values obtained from the analyses of π N , N N, and photoproduction data.

Summary
The strong-interaction width Γ π H 1s of the ground state of pionic hydrogen has been determined from the line shape of the X-ray transitions (2 p − 1s), (3 p − 1s), and (4 p − 1s) for H 2 densities ranging from 4 bar (effective) to liquid. Using different transitions, a consistency check can be done on the influence of the density-dependent Doppler broadening due to Coulomb de-excitation, which is the main obstacle for a high-precision determination of Γ π H 1s . Compared to the previous experiment [24,25], which measured the (3 p − 1s) transition at one target density only, the uncertainty for Γ π H 1s was reduced to about 50% of the one quoted there. The scattering of our results, however, demonstrates that the measurement of one transition at one density cannot avoid the noticeable errors arising from the complexity of the corrections for the Doppler broadening, even though the background levels have been reduced by factors up to 50 in the present series of measurements.
Both the frequentist approach as well as the Bayesian analysis clearly show the presence of a distinctive lowenergy component of the kinetic energy distribution containing about 2/3 of the total intensity. The complex structure of the kinetic energy spectrum at higher energies can be effectively modelled by one or two high-energy components for the present experimental statistics. The details of these components are not critical for the numerical result, however, the accuracy for Γ π H 1s is limited to about 3%. To summarize the current situation one can assert that the theoretical prediction for the kinetic energy distributions and the experimental results do not show any drastic inconsistencies.
Taking into account the difficulties related to the lineshape corrections due to Coulomb-de-excitation, future experiments should preferably use X-ray transitions from higher np initial states where the Doppler broadening is expected to be smaller. Since the accelerator currents improve continuously with time, future measurements with much improved statistics for the (4 − 1) or even (5 − 1) transition will become feasible. In this respect it is important to emphasize that for measurements using Johann-type setups, the coverage of a sufficiently large energy interval, i. e. a sufficiently large detector being equivalent to choose a well adapted dispersion, is essential to fix the background on both sides of the line. In this way, an accuracy of 1% or better can be achieved for the hadronic width of pionic hydrogen.