Scholarly article on topic 'Measurement of the charged pion mass using X-ray spectroscopy of exotic atoms'

Measurement of the charged pion mass using X-ray spectroscopy of exotic atoms Academic research paper on "Physical sciences"

CC BY
0
0
Share paper
Academic journal
Physics Letters B
OECD Field of science
Keywords
{"Charged pion mass" / "Exotic atoms" / "X-ray spectroscopy"}

Abstract of research paper on Physical sciences, author of scientific article — M. Trassinelli, D.F. Anagnostopoulos, G. Borchert, A. Dax, J.-P. Egger, et al.

Abstract The 5 g − 4 f transitions in pionic nitrogen and muonic oxygen were measured simultaneously by using a gaseous nitrogen–oxygen mixture at 1.4 bar. Due to the precise knowledge of the muon mass the muonic line provides the energy calibration for the pionic transition. A value of (139.57077 ± 0.00018) MeV/c2 (± 1.3 ppm) is derived for the mass of the negatively charged pion, which is 4.2 ppm larger than the present world average.

Academic research paper on topic "Measurement of the charged pion mass using X-ray spectroscopy of exotic atoms"

Contents lists available at ScienceDirect

Physics Letters B

www.elsevier.com/locate/physletb

Measurement of the charged pion mass using X-ray spectroscopy of exotic atoms

M. Trassinellia *. D.F. Anagnostopoulosb, G. Borchertc1, A. Daxd, J.-P. Eggere, D. Gottac M. Hennebachc 2, P. Indelicatof, Y.-W. Liud'3, B. Manilf'4, N. Nelmsg'5, L.M. Simonsd, A. Wellsg

a Institut des NanoSciences de Paris, CNRS-UMR 7588, Sorbonne Universités, UPMC Univ. Paris 06, 75005, Paris, France b Dept. of Materials Science and Engineering, University of ¡oannina, GR-45110 ¡oannina, Greece c Institut für Kernphysik, Forschungszentrum Jülich GmbH, D-52425Jülich, Germany d Laboratory for Particle Physics, Paul Scherrer Institut, CH 5232-Villigen PSI, Switzerland e Institut de Physique de l'Université de Neuchâtel, CH-2000 Neuchâtel, Switzerland

f Laboratoire Kastler Brossel, UPMC-Sorbonne Universités, CNRS, ENS-PSL Research University, College de France, 4, place Jussieu, 75005 Paris, France g Dept. of Physics and Astronomy, University of Leicester, Leicester LEI7RH, England, United Kingdom

CrossMark

A R T I C L E I N F 0

A B S T R A C T

Article history:

Received 10 May 2016

Received in revised form 10 June 2016

Accepted 12 June 2016

Available online 15 June 2016

Editor: V. Metag

Keywords: Charged pion mass Exotic atoms X-ray spectroscopy

The 5g — 4f transitions in pionic nitrogen and muonic oxygen were measured simultaneously by using a gaseous nitrogen-oxygen mixture at 1.4 bar. Due to the precise knowledge of the muon mass the muonic line provides the energy calibration for the pionic transition. A value of (139.57077 ± 0.00018)MeV/c2 (± 1.3 ppm) is derived for the mass of the negatively charged pion, which is 4.2 ppm larger than the present world average.

© 2016 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license

(http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.

X-ray spectroscopy of exotic atoms allows the determination of the mass of captured negatively charged particle like muons, pions, and antiprotons from the energies of the characteristic X-radiation. X-ray transitions occur during the de-excitation cascade of the exotic atom which is formed at principal quantum numbers of n & 16 in the case of pions [1,2]. The precise determination of the pion mass requires the use of X-ray lines which are not affected either by strong-interaction effects nor by collisions with surrounding atoms. Such conditions are found in the intermediate part of the cascade for exotic atoms formed in gases.

The most recent X-ray measurements were performed at the Paul Scherrer Institute (PSI) and used either a DuMond [3-5] or

* Corresponding author.

E-mail address: martino.trassinelli@insp.jussieu.fr (M. Trassinelli).

1 Present address: TU Munich, D-85747 Garching, Germany.

2 Present address: DAHER NUCLEAR TECHNOLOGIES GmbH, D-63457 Hanau, Germany.

3 Present address: Phys. Depart., National Tsing Hua Univ., Hsinchu 300, Taiwan.

4 Present address: Lab. de Physique des Lasers, Université Paris 13, Sorbonne Paris Cité, CNRS, France.

5 Present address: ESA-ESTEC, PO Box 299, 2200 AG, Noordwijk, The Netherlands.

a Johann-type crystal spectrometer [6]. In the case of the Du-Mond spectrometer, the energy calibration for the pionic magnesium (4f — 3d) transition was performed with a nuclear y-ray, while for the Johann set-up Ka fluorescence radiation from copper was used to determine the energy of the pionic nitrogen (5g — 4 f) transition.

In the n Mg experiment, electron refilling is unavoidable due to the use of a solid state target. Different assumptions on the K electron population lead to differences in the pion mass up to 16 ppm [5]. The previous n N experiment, as well as the present one, used a nitrogen gas target at pressures around 1 bar, where electron refilling is unlikely [7,8], i.e. the de-excitation cascade is decoupled from the environment. The absence of refilling of the electrons ejected already during the upper part of the cascade by internal Auger effect manifests in the appearance of X-ray lines at n > 5, which otherwise would be converted into Auger transitions [9-11]. Furthermore, a large Doppler broadening was measured for (5 — 4) transitions [12]. It originates from Coulomb explosion during the formation process of the exotic atom with molecules and indicates that the velocity at the time of X-ray emission is essentially unchanged since the breakup of the molecule. Thus, the

http://dx.doi.org/10.1016/j.physletb.2016.06.025

0370-2693/© 2016 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.

Table 1

Calculated contributions to the total QED transition energy of /O and nN(5g — 4f) lines (in eV) [17]. For the pionic transition, the world average pion mass value as given in [14] is used. The /O line constitutes a triplet due to the muon spin. The total uncertainty of the QED calculation (excluding the uncertainty of the pion mass) is ±1meV.

Transition n 14N

(5g9/2 - 4/7/2) (5g7/2 - 4/7/2) (5g7/2 - 4/5/2) (5g - 4 / )

Coulomb 4022.8625 4022.6188 4023.4124 4054.1180

self energy -0.0028 -0.0013 -0.0013 -0.0001

vac. pol. (Uehling) 0.8800 0.8800 0.8807 1.2485

vac. pol. Wichman-Kroll -0.0007 -0.0007 -0.0007 -0.0007

vac. pol. two-loop Uehling 0.0003 0.0004 0.0004 0.0008

vac. pol. Kallén-Sabry 0.0084 0.0084 0.0084 0.0116

relativistic recoil 0.0025 0.0025 0.0025 0.0028

hyperfine structure - - - -0.0008

Total 4023.7502 4023.5079 4024.2983 4055.3801

absence of screening effects from remaining electrons in the intermediate part of the atomic cascade leads to a unique solution for the mass [6]. In addition, in dilute targets the line intensity is already mostly collected in the circular transitions (n, I = n — 1) ^ (n — 1, I = n — 2), where corrections owing to the hadronic potential are still tiny.

From the nN experiment mn — = (139.57071 ± 0.00053) MeV/c2 [6] is obtained which suggests that both K electrons are present when the nMg(4f — 3d) transition occurs (solution B: mn — = (139.56995 ± 0.00035) MeV/c2 [5]). This is corroborated by the fact that the result, assuming 1 K electron only (solution A: mn — = (139.56782 ± 0.00037) MeV/c2), is in conflict with the measurement of the muon momentum for charged pion decay at rest n + ^ i+V/ [13]. For solution A, the mass squared of the muon neutrino becomes negative by six standard deviations, whereas the average of solution B and the result of the n N(5g — 4f) measurement (mn — = (139.57018 ± 0.00035) MeV/c2 [14]) yields the upper limit mlv < 190keV/c2 (90% c.l.).

The experiment described here resumes the strategy of the gas target, but exploits (i) the high precision of 0.033 ppm for the mass of the positively charged muon being = (105.6583715 ± 0.0000035) MeV/c2 [14] and (ii) the unique feature that in nN and /O transition energies almost coincide (Table 1). Using a N2/O2 gas mixture in the target allows the simultaneous measurement of n N and /O lines, with the muonic transition serving as an on-line calibration. Hence, systematic shifts during the unavoidably long measuring periods are minimised.

In the case of nitrogen and oxygen, (6h — 5g), (5 g — 4 f), and (4f — 3d) transitions meet the operating conditions of the crystal spectrometer. Finally, the (5g — 4f) transition was chosen because: (i) for the (6h — 5g) lines (2.2 keV) absorption in the target gas itself and windows significantly reduces the count rate and (ii) the 3d-level energy in n N requires a substantial correction because of the strong interaction. Electromagnetic transition energies (Tables 1 and 2) were calculated using a multi-configuration Dirac-Fock approach [15,16] to a precision of ±1 meV and include relativistic and quantum electrodynamics contribution (relativistic recoil, self-energy, vacuum polarisation) as well as the hyperfine structure of pionic nitrogen [17].

Energy shifts due to nuclear finite size are found to be as small as 4aeV and 2peV for the 5g and 4f levels in n N. Values for nuclear masses, radii, and moments were taken from recent compilations [18-20]. The strong-interaction shifts of the n N levels were estimated from interpolating the measured hadronic 2p-level shifts in nC and nO [21] and by using scaling relations based on the overlap of nucleus and a hydrogen-like wave function for the pion orbit (see Table 3). Details on the calculation of the transition energies may be found elsewhere [22].

The measurement was performed at the high-intensity pion beam line n E5 of the Paul Scherrer Institute (PSI) using a set-up

Table 2

Transition energies £qED [17] and Bragg angles 6B of the /O and n N lines used in the fit to the spectrum. The relative intensities within the fine structure multiplets of /O (FS int.) have been fixed in the fit to the statistical weight. The Bragg angle includes the index of refraction shift calculated with the code XOP [32]. For twice the lattice distance 2d = 0.768 062 286 (13) nm is assumed at a temperature of 22.5°C [34]. The conversion constant used is hc = 1.239 841 930 (28) nmkeV [14]. The nN(5g — 4f) and nN(5f — 4d) transition energies include the strong-interaction shift (see Table 3).

Transition FS int. EQED/eV ®e

M16O(5g7/2 - 4/7/2 ) 1 4023.5079 53c 21'51.48''

M16O(5g9/2 - 4/7/2 ) 35 4023.7503 53c 21'34.77''

M16O(5g7/2 - 4 /5/2) 27 4024.2984 53c 20'57.01''

M16O(5 /5/2 - 4d5/2) 1 4025.3956 53c 19'41.47''

M16O(5 /7/2 - 4d5/2) 20 4025.8031 53c 19'13.44''

M16O(5 /5/2 - 4d3/2) 14 4026.9922 53c 17'51.70''

M16O(5d5/2 - 4p3/2) 9 4028.5625 53c 16'3.90''

M16O(5d3/2 - 4pi/2) 5 4033.5273 53c 10'24.10''

M18O(5g7/2 - 4/7/2 ) 1 4026.6692 53c '18'13.90''

M18O(5g9/2 - 4 /7/2) 35 4026.9132 53c 17'57.13''

M18O(5g7/2 - 4 /5/2) 27 4027.4642 53c 17'19.28''

n 14N(5g - 4 / ) 4055.3802 52°45'46.76"

n 14N(5/ - 4d) 4057.6984 52°43'11.81"

n 14N(5d - 4p) QED only 4061.9460 52°38'28.76"

n 15N(5g - 4 / ) 4058.2394 52°42'35.67''

n 15N(5/ - 4d) 4060.5605 52°40' 0.95''

similar to the one used by Lenz et al. [6]. Major improvements are: (i) The use of cyclotron trap II [23] having a larger gap between the magnet coils yielding a substantially increased muon stop rate, (ii) a Bragg crystal of superior quality and (iii) a large-area X-ray detector in order to simultaneously cover the reflections of the muonic and pionic transitions (see Fig. 1). In addition, the average proton current of the accelerator was about 1.4 mA, which is 40% higher than in the previous experiment.

The N2/O2 gas mixture was enclosed in a cylindrical target cell placed at the centre of the cyclotron trap. The cell wall was made of a 50jam thick Kapton® foil. Towards the crystal spectrometer a circular 7.5 jam Mylar® window was used supported by a stainless steel honeycomb structure with a free area of 90%. The target was operated at 1.4 bar and room temperature.

The muons used originate from the decay of slow pions inside the cyclotron trap, because the stop density for muons at the high-intensity pion beam is still superior to the one at a dedicated muon channel. For the simultaneous measurement comparable count rates are required for the nN and the ¡iO line. This was achieved with a N2/O2 mixture of 10%/90% by adapting the set of polyethylene degraders inside the magnet gap and optimised by

Table 3

Corrections to the measured angle difference between the n 14N(5g — 4f) and the /160(5g9/2 — 4f7/2) transitions and associated uncertainties. A 1ppm change in the pion mass corresponds to 4.055 meV in transition energy, to 0.27 arcsec in diffraction angle, or to a displacement of 3.2 pm in the detector plane. Contributions to the mass uncertainty from lattice and conversion constant cancel in leading order because the measurement principle is based on the angular difference. For more details see text.

Type of uncertainty

/ arcsec

/ arcsec

Total / arcsec

Uncertainty / ppb

index of refraction shift silicon lattice constant bending correction penetration depth correction

14.01 —0.07

13.71 —0.07

0.30 0

±20 ±2 ±20 ±4

focal length temperature correction CCD alignment pixel distance

alignment of detector normal detector height offset shape of target window

— 30 + 0

— 35

shape of reflection

individual curvature correction

response function and Doppler broadening

line pattern modelling

fit interval

±225 ±150

— 350 + 190

— 290

muon mass QED energy conversion constant hc 4f strong interaction 45|ieV 5g strong interaction 0.2 |ieV K electron screening

total systematic error statistical error

0.003 0.000

—0.003 0.000

+ 960 — 1000

means of an X-ray measurement using a Si(Li) semiconductor detector.

The crystal spectrometer is set up in Johann geometry [24] using a spherically bent Bragg crystal and optimised to the needs of exotic-atom X-ray spectroscopy [25]. Such a configuration allows the simultaneous measurement of two different energies within an energy interval, the limits of which are given by the extension of the target in the direction of dispersion and correspondingly by the size the detector. Spherical bending leads to a partial vertical focusing [26] which increases the count rate.

The Bragg crystal was made from a silicon crystal disk of 290]am thickness and of a diameter of 100 mm. The disk is attached to a high-quality polished glass lens defining a spherical segment. The average radius of curvature of the crystal surface was measured to Rc = (2981.31 ± 0.33) mm by sampling 500 points at the surface with a mechanical precision sensor (performed by Carl Zeiss AG, D-73447 Oberkochen, Germany). An upper limit for the cut angle (angle between crystal surface and reflecting lattice planes) was determined in a dedicated measurement to be 120 seconds of arc [27]. Hence, the focal condition corresponds to the symmetric Bragg case being Rc • sin ©B. The measurement uses the second order reflection at the (110) planes. An aluminium aperture of 90 mm diameter covered the boundary region of the Si disk in order to avoid edge effects. For source geometry as given here, the overall efficiency of the crystal set-up is & 5 • 10—8. About 85% of the reflected intensity is covered by the sensitive area of the detector.

The detector with a total sensitive area of about 48x72 mm2 (width x height) was built up by a 2x3 array of charge-coupled

devices (CCDs) of 24 mmx 24 mm (600x 600 pixels) with frame storage option [28]. Having a depletion depth of about 30 ^m these CCDs reach their maximum in detection efficiency of almost 90% at 4keV. The detector surface is oriented perpendicular to the direction of the incoming X-rays. Excellent background conditions are achieved (i) by using an especially tailored concrete shielding of at least 1m thickness between the X-ray detector and the target region and (ii) by exploiting the different pixel topology of low-energy X-rays and background events, which are mainly caused by neutron induced high energetic nuclear y rays [2,6].

The Bragg angle for the nN(5g — 4f) transition and thereby its energy is determined from the position difference to the ¡0(5 g — 4f) line. The positions are determined from the projection of the pattern on the CCD onto the direction of dispersion after correction for curvature by means of a parabola fit (Fig. 1). The main transitions ¡i,O(5g — 4f) and n N(5g — 4f) are separated by about 25 mm.

About 9000 events for each element were collected in each of the (5g — 4f) transitions during 5 weeks of data taking. The count rates for the n N and ¡0 transitions were about 15 events per hour each. Only a common small drift was observed for the line positions of less than one pixel in total. Because of the simultaneous measurement the position difference is not affected. Bragg angle dependent corrections are small because the leading order cancels in such a difference measurement performed in the same order of reflection.

In fourth order, the Bragg angles of the Cu Ka lines are very close to the ones of the ¡0(5g — 4f) transitions. Therefore, in addition Cu X-rays were repeatedly recorded as a stability monitor

Fig. 1. Simultaneously measured (5g — 4f) transitions in muonic oxygen (calibration) and pionic nitrogen. Top: Distribution of the Bragg reflections on the surface of the 2 x 3 CCD array. The binning corresponds to the pixel size of the CCDs (note the different scales vertically and horizontally). Straight dashed lines indicate CCD boundaries. Middle: Projection on the axis of dispersion after correction for curvature (see text). Bottom: Details of the fit to line patterns.

corroborating the amount of the small common drift observed for the fiO/n N pair.

Various parameters of the analysis and of the set-up enter in the determination of the line positions and their difference. These contributions and their uncertainties are summarised in Table 3 and are discussed in detail below.

Index of refraction shift. The systematic uncertainty of the index shift correction is assumed to be about 5% [29,30], i.e. the uncertainty of the difference is negligibly small.

Silicon lattice constant and wavelength conversion. Both the silicon lattice constant 2d and the conversion constant hc are known to an accuracy of ^ 10—8.

Bending and penetration depth corrections. The energy dependent penetration depths of the X-rays lead to different corrections for the lattice constant of the Bragg crystal due to its curvature. The difference of the shift due to the average penetration depths itself turns out to be negligible. The primary extinction lengths including absorption were calculated both with the codes XOP [32] and DIXI [33], where results were found to coincide perfectly. We assume that the crystal behaves like an ideal one for such large bending radii [35]. The corrections for the Bragg angle were calculated following the approach of [36,37] using for the Poisson number the value v = 0.208 obtained from [38,39].

Focal length. Because of the different focal lengths for the n N and /O lines of 18.4 mm, the detector was placed in an intermediate position, which was determined by a survey measurement to be (2388.27 ± 0.20) mm. The uncertainty of the distance crystal-to-detector represents the largest contribution to the systematic error.

Temperature correction. The temperature during the measurement varied between 19° C and 21° C during the measurement. All periods were rescaled to 22.5° C by using the appropriate thermal expansion coefficient. The main correction comes from the change of the lattice constant. A smaller contribution arises from the variation of the distance crystal detector.

CCD alignment and pixel distance. In the CCD array small gaps of the order of 0.3 mm emerge between the individual devices. Secondly, the nominal pixel size of the CCDs, reported to be 40 jam x 40 jam at room temperature, changes for the operating temperature of — 100°C. Both the relative orientations of the six CCD devices and the average pixel distance have been measured precisely in a separate experiment using a nanometric quartz mask [31]. The average pixel distance was found to be (39.9775 ± 0.0006) jm, substantially different from the nominal value.

Alignment of detector normal. The surface of the CCD array was set-up perpendicular to the direction crystal-detector to better than ±0.14°. The uncertainty also includes the imperfectness of the vacuum tubes, of their connections, and of the support structures of the CCDs.

Detector height offset. A possible offset in height of the detector from the ideal geometry defined by the plane through the centres of X-ray source, crystal, and detector leads to a distortion of the reflections. The size of such an effect was quantified by means of a Monte-Carlo simulation.

Shape of the target window. The circular shape of the target window leads to boundaries of different inclination for the n N and /O reflections. The corresponding possible uncertainty for the position difference was determined from a Monte-Carlo simulation.

Shape of reflection. The curvature of the n N and /O reflections is determined from a parabola fit to the hit pattern of the circular transitions. The assumption of a parabolic shape for the curvature is valid only close to the above-mentioned central plane. In addition, the curvature fit assumes a constant width of the reflection. A possible effect on the position difference over the height of the CCD array, which principally increases with increasing distance from the central plane, was studied by restricting the detector surface in height. The deviations are found to be far below the statistical error of the line positions.

Individual curvature correction. The parabola parameters for the n N and /O reflections are slightly different because of different

focal lengths. No difference could be verified from the fits which, however, is expected within the available statistics. The uncertainty is therefore given by the error of the fit to the curvature. For curvature correction, the average values were taken of the n N and ¡ O reflection.

Response function and Doppler broadening. The response is found by a convolution of the intrinsic crystal response with the aberration caused by the imaging properties of a spherically bent crystal. The crystal response was calculated with the code XOP [32], and the geometry was taken into account by means of Monte-Carlo ray-tracing [25]. The resulting response shows a significant asymmetry having a width of 450 meV (FWHM).

Measured line widths of n N and ¡O transitions, however, are dominated by Doppler broadening due to Coulomb explosion [12], which was underestimated in the analysis reported by Lenz et al. [6] because of an inferior quality of the Bragg crystal. The line shapes are almost symmetric having a width of about 750 meV (FWHM). The Doppler broadening was accounted for best by folding in an additional Gaussian of about 40 seconds of arc. The Gaussian was determined from the analysis of a dedicated measurement optimised for pion stops, where in total 60000 events were accumulated in the nN(5g — 4f) transition.

The defocusing due to the different focal lengths is included in the Monte-Carlo based response, which is calculated for the appropriate distance in each case. In addition, it was verified that the parameters found in the curvature fit to the data are reproduced for the Monte-Carlo result.

Line pattern. The total line pattern to be considered is a superposition of the circular (5g — 4f) and the inner transitions (5f — 4d) and (5d — 4p) together with the corresponding contributions from the other isotopes (Table 2). The isotope abundances are fixed as tabulated (16O/18O: 99.76%/0.21%, 14N/15N: 99.64%/0.36%). The relative intensities of the inner transitions are due to the cascade dynamics and, therefore, free parameters of the fit.

The line positions within the nN and ¡O(5g — 4f) patterns were fixed according to the QED energies. In the case of ¡ O, all fine structure components were included in the fit. For a proper description of the background, the two strong components of the ¡16O(5d — 4p) triplet and the nN(5d — 4p) transition were included in the fit. For the pionic line, position and width were free parameters, because it is shifted and broadened by about 1eV compared to the electromagnetic value by the strong interaction [6].

Fit interval. Changing the interval used in the fit of the line positions does affect the result insignificantly.

K electron screening. From the analysis of the high-statistics n N(5g — 4f) data, we exclude the influence of satellites lines due to remaining K electrons. The energy shift of the pionic transition is calculated to be —456 (—814) meV in the case of one (two) K electron(s). Two hypothesis (presence of satellites or not) are compared via the Bayes factor [40-43] yielding an upper limit of less than 3 • 10—6 for the relative intensity of possible satellites.

The measured energy of the nN(5g — 4f) transition was found to be (4055.3970 ± 0.0033stat ± 0.0038sys) eV. Basically two facts limit the accuracy of the method of a simultaneous measurement as described here: (i) The low rate obtainable from the muonic transitions hinders to accumulate as high statistics as would be achievable when using a set-up optimised for pionic atoms. For pi-onic transitions, count rates being a factor of 20 larger than for

Fig. 2. Results for the mass of the charged pion. Also shown are previous exotic-atom results (Jeckelmann et al. (86B [3,4]), Lu et al. [44], Carter et al. [45], Marushenko et al. [46]) and n + decay at rest (Daum et al. [47]). The shaded region indicates the world average before this experiment [14].

muonic X-rays can be achieved. (ii) The large Doppler broadening induced by Coulomb explosion when using diatomic gases, which approximately doubles the line width as expected from the spectrometer response.

To summarise, the mass of the negatively charged pion has been measured by means of equivalent X-ray transitions in hydrogen-like pionic nitrogen and muonic oxygen, where the muonic line serves as energy calibration. The value of (139.57077 ± 0.00018) MeV/c2 is 4.2 ppm larger than the present world average [14]. Repeating the procedure as described in ref. [6] by using the Cu Ka1 line for calibration, yields a value of mn = (139.57090 ± 0.00056) MeV/c2. The accuracy of ±4.0ppm represents the limit for a calibration with broad X-ray fluorescence lines. Both results are in good agreement with the mass obtained by [6], but 5.4 ppm and 6.8 ppm, respectively, above the result of the pionic magnesium experiment (solution B [5]) using a nuclear y ray for calibration (Fig. 2).

The analysis shows no evidence for any satellite lines from remaining electrons at the time of X-ray emission of the (5 g — 4 f) transition. This corroborates strongly our assumption for a complete depletion of the electron shell during the preceding steps of the atomic cascade.

In conclusion, the present study demonstrates the potential of crystal spectroscopy with bent crystals in the field of exotic atoms. Its limits are given, on one hand, by statistics for the present beam and detector technologies. On the other hand, the systematic uncertainties discussed at length above illustrate the level of sophistication which must be applied.

Facing the fact that pion beams at PSI provide a flux of about 109/s, the use of double-flat crystal spectrometers may be considered allowing for absolute angle calibrations without a (muonic or X-ray) reference line. Choosing pionic transitions not affected by Coulomb explosion, e.g. from pionic neon, a precision for the pion mass determination of the order of 0.5 ppm is feasible which, however, may be regarded to be the maximum achievable by means of exotic-atom X-ray spectroscopy.

As a result, X-rays of hydrogen-like pionic atoms are useful to provide calibration standards in the few keV range, where suitable radioactive sources are not available [22,49]. At present, the accuracy is given by the uncertainty of the pion mass [50]. The quality of such standards may benefit substantially from laser spec-

troscopy of metastable high-lying pionic states which is proposed [13

to be performed in pionic helium also at PSI [51]. [14

Combined with the measurement of the muon momentum af- [15

ter pion decay at rest [13], a non-zero value for the muon neutrino ^ mass is obtained of mVl = 183 + 63 keV/c2 (c.l. 90%) when using [17

the statistical approach of [48]. The result is far above the cosmo- [18

logical limit of at least 11 eV/c2 for the sum of all neutrino flavours [19

[14]. However, extending the error limits to 3a either for the pion [:?!

with zero Acknowledgements

References

mass or the muon momentum yields values for mVl consistent [22

We are grateful to N. Dolfus, H. Labus, B. Leoni and K.-P. Wieder

for solving numerous technical problems. We thank the PSI staff [26

for providing excellent beam conditions and appreciate the support [27

by the Carl Zeiss AG, Oberkochen, Germany, which fabricated the [29

Bragg crystals. We thank Prof. Dr. E. Förster and his collaborators [30

at the University of Jena, and A. Freund and his group at ESRF, for [31

the help in characterising the crystal material as well as A. Blech- [32

mann for a careful study of the CCD performance. We are indebted [33

to PSI for supporting the stay during the run periods (D.F.A.). This [34 work is part of the PhD thesis of B.M. (Université Pierre et Marie

Curie, 2001), N.N. (University of Leicester, 2002) and M.T. (Univer- [35

sité Pierre et Marie Curie, 2005). [36

[37 [38

[1] L.M. Simons, D. Horvâth, G. Torelli (Eds.), Proc. of the Fifth Course of the International School of Physics of Exotic Atoms, May 14-20, 1989, Erice, Italy, [41 Plenum Press, New York, 1990, and references therein.

[2] D. Gotta, Prog. Part. Nucl. Phys. 52 (2004) 133. [42

[3] B. Jeckelmann, et al., Phys. Rev. Lett. 56 (1986) 1444. [43

[4] B. Jeckelmann, et al., Nucl. Phys. A 457 (1986) 709. [44

[5] B. Jeckelmann, P.F.A. Goudsmit, H.J. Leisi, Phys. Lett. B 335 (1994) 326. [45

[6] S. Lenz, et al., Phys. Lett. B 416 (1998) 50. [46

[7] R. Bacher, et al., Phys. Rev. A 39 (1989) 1610. [47

[8] K. Kirch, et al., Phys. Rev. A 59 (1999) 3375. [48

[9] G.R. Burbidge, A.H. de Borde, Phys. Rev. 89 (1953) 189. [49

[10] P. Vogel, Phys. Rev. A 22 (1980) 1600. [50

[11] R. Bacher, et al., Phys. Rev. Lett. 54 (1985) 2087. [51

[12] T. Siems, et al., Phys. Rev. Lett. 84 (2000) 4573.

K. Assamagan, et al., Phys. Rev. D 53 (1996) 6065.

K.A. Olive, et al., Particle Data Group, Chin. Phys. C 38 (9) (2014) 090001.

J. Desclaux, et al., Computational Approaches of Relativistic Models in Quantum

Chemistry, vol. 10, Elsevier, 2003.

J.P. Santos, et al., Phys. Rev. A 71 (2005) 032501.

M. Trassinelli, P. Indelicato, Phys. Rev. A 76 (2007) 012510.

G. Audi, A. Wapstra, C. Thibault, Nucl. Phys. A 729 (2003) 337.

I. Angeli, At. Data Nucl. Data Tables 87 (2004) 185.

P. Raghavan, At. Data Nucl. Data Tables 42 (1989) 189.

G. de Chambrier, et al., Nucl. Phys. A 442 (1985) 637.

M. Trassinelli, PhD thesis, Univ. Piere et Marie Curie, Paris, 2005, http://tel.ccsd.

cnrs.fr/tel-00067768.

L.M. Simons, Phys. Scr. T 22 (1988) 90;

L.M. Simons, Hyperflne Interact. 81 (1993 ) 253.

H.H. Johann, Z. Phys. 69 (1931) 185.

D.E. Gotta, L.M. Simons, Spectrochim. Acta, Part B 120 (2016) 9, http://dx.doi. org/10.1016/j.sab.2016.03.006.

J. Eggs, K. Ulmer, Z. Angew. Phys. 20 Band (Heft 2) (1965) 118.

D.S. Covita, et al., Rev. Sci. Instrum. 79 (2008) 033102.

N. Nelms, et al., Nucl. Instrum. Methods, Sect. A 484 (2002) 419.

B.L. Henke, E.M. Gullikson, J.C. Davies, At. Data Nucl. Data Tables 54 (1993) 181.

C.T. Chantler, J. Phys. Chem. Ref. Data 24 (1995) 71.

P. Indelicato, et al., Rev. Sci. Instrum. 77 (2006) 043107.

M. Sanchez del Rio, R.J. Dejus, XOP version 2.4: recent developments of the

X-ray optics toolkit, SPIE Proc. 8141 (2011) 814115.

G. Hölzer, O. Wehrhan, E. Förster, Cryst. Res. Technol. 33 (1998) 555.

P.J. Mohr, B.N. Taylor, D.B. Newell, C0DATA2010, Rev. Mod. Phys. 84 (2012)

I. Uschmann, et al., J. Appl. Crystallogr. 26 (1993) 405. F. Cembali, et al., J. Appl. Crystallogr. 25 (1992) 424.

C.T. Chantler, R.D. Deslattes, Rev. Sci. Instrum. 66 (1995) 5123. J.J. Wortman, R.A. Evans, J. Appl. Phys. 36 (1965) 153.

F.N. Chukhovskii, et al., J. Appl. Crystallogr. 29 (1996) 438.

E.T. Jaynes, G.L. Bretthorst, Probability Theory: The Logic of Science, Cambridge University Press, 2003.

D.S. Sivia, J. Skilling, Data Analysis: A Bayesian Tutorial, second ed., Oxford University Press, 2006.

R.E. Kass, A.E. Raftery, J. Am. Stat. Assoc. 90 (1995 ) 773.

C. Gordon, R. Trotta, Mon. Not. R. Astron. Soc. 382 (2007) 1859.

D.C. Lu, et al., Phys. Rev. Lett. 45 (1980) 1066. A.L. Carter, et al., Phys. Rev. Lett. 37 (1976) 1360. V.N. Marushenko, et al., JETP Lett. 23 (1976) 72. M. Daum, et al., Phys. Lett. B 265 (1991) 425.

G.J. Feldman, R.D. Cousins, Phys. Rev. D 57 (1998) 3873.

D.F. Anagnostopoulos, et al., Phys. Rev. Lett. 91 (2003) 240801. S. Schlesser, et al., Phys. Rev. C 84 (2011) 015211.

M. Hori, V.I. Korobov, Anna Soter, Phys. Rev. A 89 (2014) 042515, http:// dx.doi.org/10.1103/PhysRevA.89.042515.