Next Article in Journal
The Analysis of Fractional-Order Kersten–Krasil Shchik Coupled KdV System, via a New Integral Transform
Previous Article in Journal
Different Drop Heights in Bilateral Asymmetry and Interjoint Coordination during Repetitive Drop-Jumps
 
 
symmetry-logo
Article Menu

Article Menu

Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Search for μ+ → e+γ with 10–14 Sensitivity: The Upgrade of the MEG Experiment

by
Alessandro M. Baldini
1,
Vladimir Baranov
2,
Michele Biasotti
3,4,
Gianluigi Boca
5,6,
Paolo W. Cattaneo
5,
Gianluca Cavoto
7,8,
Fabrizio Cei
1,9,
Marco Chiappini
1,9,
Gianluigi Chiarello
7,8,
Alessandro Corvaglia
10,11,
Federica Cuna
10,11,
Giovanni dal Maso
1,9,12,
Antonio de Bari
5,6,
Matteo De Gerone
3,
Marco Francesconi
1,9,
Luca Galli
1,
Giovanni Gallucci
3,
Flavio Gatti
3,4,
Francesco Grancagnolo
10,
Marco Grassi
1,
Dmitry N. Grigoriev
13,14,
Malte Hildebrandt
12,
Kei Ieki
15,
Fedor Ignatov
13,
Toshiyuki Iwamoto
15,
Peter-Raymond Kettle
12,
Nikolay Khomutov
2,
Satoru Kobayashi
15,
Alexander Kolesnikov
2,
Nikolay Kravchuk
2,
Victor Krylov
2,
Nikolay Kuchinskiy
2,
William Kyle
16,
Terence Libeiro
16,
Vladimir Malyshev
2,
Manuel Meucci
7,8,
Satoshi Mihara
17,
William Molzon
16,
Toshinori Mori
15,
Alexander Mtchedlishvili
12,
Mitsutaka Nakao
15,
Donato Nicolò
1,9,
Hajime Nishiguchi
17,
Shinji Ogawa
15,
Rina Onda
15,
Wataru Ootani
15,
Atsushi Oya
15,
Dylan Palo
16,
Marco Panareo
10,11,
Angela Papa
1,9,12,
Valerio Pettinacci
7,
Alexander Popov
13,
Francesco Renga
7,
Stefan Ritt
12,
Massimo Rossella
5,
Aleksander Rozhdestvensky
2,
Patrick Schwendimann
12,18,
Kohei Shimada
15,
Giovanni Signorelli
1,
Alexey Stoykov
12,
Giovanni F. Tassielli
10,11,
Kazuki Toyoda
15,
Yusuke Uchiyama
15,*,
Masashi Usami
15,
Cecilia Voena
7,*,
Kosuke Yanai
15,
Kensuke Yamamoto
15,
Taku Yonemoto
15 and
Yury V. Yudin
13
add Show full author list remove Hide full author list
1
INFN Sezione di Pisa, Largo B. Pontecorvo 3, 56127 Pisa, Italy
2
Joint Institute for Nuclear Research, 141980 Dubna, Russia
3
INFN Sezione di Genova, Via Dodecaneso 33, 16146 Genoa, Italy
4
Dipartimento di Fisica, dell’Università di Genova, Via Dodecaneso 33, 16146 Genoa, Italy
5
INFN Sezione di Pavia, Via Bassi 6, 27100 Pavia, Italy
6
Dipartimento di Fisica, dell’Università di Pavia, Via Bassi 6, 27100 Pavia, Italy
7
INFN Sezione di Roma, Piazzale A. Moro 2, 00185 Rome, Italy
8
Dipartimento di Fisica, dell’Università “Sapienza” di Roma, Piazzale A. Moro 2, 00185 Rome, Italy
9
Dipartimento di Fisica, dell’Università di Pisa, Largo B. Pontecorvo 3, 56127 Pisa, Italy
10
INFN Sezione di Lecce, Via per Arnesano, 73100 Lecce, Italy
11
Dipartimento di Matematica e Fisica, dell’Università del Salento, Via per Arnesano, 73100 Lecce, Italy
12
Paul Scherrer Institut PSI, 5232 Villigen, Switzerland
13
Budker Institute of Nuclear Physics of Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
14
Physico-Technical Department, Novosibirsk State Technical University, 630092 Novosibirsk, Russia
15
ICEPP, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
16
Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA
17
KEK, High Energy Accelerator Research Organization 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan
18
Department of Physics, Swiss Federal Institute of Technology ETH, 8093 Zurich, Switzerland
*
Authors to whom correspondence should be addressed.
Symmetry 2021, 13(9), 1591; https://doi.org/10.3390/sym13091591
Submission received: 23 July 2021 / Revised: 15 August 2021 / Accepted: 17 August 2021 / Published: 29 August 2021
(This article belongs to the Section Physics)

Abstract

:
The MEG experiment took data at the Paul Scherrer Institute in the years 2009–2013 to test the violation of the lepton flavor conservation law, which originates from an accidental symmetry that the Standard Model of elementary particle physics has, and published the most stringent limit on the charged lepton flavor violating decay μ + e + γ : BR( μ + e + γ ) < 4.2 × 10 13 at 90 % confidence level. The MEG detector has been upgraded in order to reach a sensitivity of 6 × 10 14 . The basic principle of MEG II is to achieve the highest possible sensitivity using the full muon beam intensity at the Paul Scherrer Institute ( 7 × 10 7 muons/s) with an upgraded detector. The main improvements are better rate capability of all sub-detectors and improved resolutions while keeping the same detector concept. In this paper, we present the current status of the preparation, integration and commissioning of the MEG II detector in the recent engineering runs.

1. Introduction

Lepton flavor conservation is associated with an approximate symmetry in the Standard Model (SM) of elementary particle physics owing to the tiny neutrino masses. It is not protected by an explicit gauge symmetry and could be broken in New Physics (NP) beyond the SM. Because of the approximate symmetry, the processes of charged lepton flavor violation (CLFV) are allowed only at extremely small branching ratios (≪10 50 ) in the SM. Therefore, they are free of SM background and they are ideal probes for NP searches. CLFV transitions have been searched for in a variety of channels, but no evidence has been found so far. Nonetheless, they are predicted at measurable rates, not far from present experimental limits, by many SM extensions [1] and in particular in light of the recent muon g 2 precision measurement [2] at Fermilab, for example in [3] (flavor symmetries) and in [4] (several extensions of SM including MSSM).
Among the CLFV processes, the decay μ + e + γ is very sensitive to NP. The MEG experiment at the Paul Scherrer Institute (PSI) searched for this decay in the period 2009–2013 and set the current best world limit to its branching ratio: B R ( μ + e + γ ) < 4.2 × 10 13 at 90% confidence level [5].
Other CLFV channels currently investigated are μ N e N , μ 3 e , τ γ and τ ( = e , μ ). The μ N e N conversion will be searched for by the DeeMe [6] and COMET [7] experiments in preparation at J-PARC. COMET will reach a sensitivity of O ( 10 15 ) in its first phase, to be compared with the extant limit of 7 × 10 13 [8] while in a second phase, it aims for a sensitivity of O ( 10 17 ) . Another experiment, Mu2e [9] under construction at Fermilab, will search for μ N e N aiming to a sensitivity of 3 × 10 17 . The μ 3 e is being searched for in a new experiment currently in the commissioning phase at PSI: Mu3e [10], that in a staged approach plans to reach a sensitivity of 10 16 , while the extant limit is 1 × 10 12 [11]. CLFV decays involving τ leptons will be studied at the Belle II [12] experiment at Super KEKB with a sensitivity goal of O ( 10 9 ) . Experiments at a super Charm-Tau factory [13], approved by the Russian government and presently in the R&D phase, were proposed with a similar sensitivity reach. The sensitivity to NP in the different μ and τ channels depends on the specific NP models assumed, and all the above modes can be considered powerful probes to explore NP, complementary to direct searches at LHC and HL-LHC.
The MEG II experiment, under commissioning at the PSI, aims for a sensitivity enhancement compared to the MEG final result, O ( 6 × 10 14 ) . MEG II has the same experimental concept as MEG, which is discussed in the next section, but faces the challenge of a more intense muon beam with the need for keeping high efficiency and detector resolutions. The full MEG II proposal can be found in [14]. In this paper, we review the main aspects of the design and we report the results obtained in the engineering runs in the years 2017–2020, in view of the 2021 run that will be the first with the whole detector instrumented. An update of the expected sensitivity is also presented.

2. Experimental Components and Methods

2.1. The Experimental Approach

In μ + e + γ search experiments, positive muons are stopped in a target and the signature of a two-body decay at rest is exploited: an e + and a γ emitted simultaneously, moving back-to-back with their energies equal to half of the muon mass ( m μ / 2 = 52.8 MeV). The signal of the μ + e + γ decay can be reconstructed by measuring the photon energy E γ , the positron momentum p e + , their relative angle Θ e + γ and time difference t e + γ .
The background comes from two sources: radiative muon decays (RMD, μ + e + ν ν ¯ γ ) with the neutrinos carrying away a small amount of energy, or from an accidental coincidence of a positron from a Michel decay ( μ + e + ν ν ¯ ) with a photon coming from RMD, bremsstrahlung or positron annihilation-in-flight (AIF, e + e γ γ ). Since the rate of accidental background events varies quadratically with the muon stopping rate R μ + , it must be chosen appropriately in order to limit the signal to background ratio and optimize the discovery potential.
The MEG II detector, shown schematically in Figure 1, is based on the same experimental concept as MEG and uses the same beam, but at the increased intensity of 7 × 10 7 μ + /s.
The positron spectrometer uses a gradient magnetic field and detects positron tracks with a low-mass single-volume cylindrical drift chamber (CDCH) with high rate capability. The positron time is measured by a pixelated timing counter (pTC). The acceptance of the spectrometer was more than doubled from MEG. The photon is measured by an upgraded liquid xenon calorimeter (LXe) where the collection of scintillation light is more uniform. This was obtained by replacing the MEG photomultiplier tubes (PMTs) on the photon entrance face with smaller vacuum-ultraviolet (VUV) sensitive silicon photomultipliers (SiPMs). A novel device for active suppression of the accidental background was introduced, the radiative decay counter (RDC). The trigger and data-acquisition system was also upgraded to deal with an increased number of readout channels and to cope with the required bandwidth. Moreover, several auxiliary devices have been developed to monitor the beam and calibrate the detectors. The MEG II subsystems are described in the following subsections focusing on the changes from MEG, while the details of the MEG system can be found in [15].

2.2. The Muon Beam and the Target

Various combinations of beam momentum and target have been investigated to increase the muon stopping rate from R μ + = 3 × 10 7   μ + /s to 7 × 10 7   μ + /s while preserving a low material budget along the beam line and in the target region. The configuration chosen is a surface muon beam of 28 MeV/c (as in MEG) with a scintillating target of 140 μ m thickness, placed at an angle of 15 with respect to the beam axis, which allows a simultaneous non-destructive beam intensity and profile measurement. A photogrammetric method to monitor the target position, orientation and shape was implemented, thus addressing one of the dominant systematic effects in MEG. The system includes two CCD cameras and LEDs placed at one of the CDCH endplates, which optically monitors printed patterns (dots) on the target [16,17].
Three new detectors were built to measure the beam profile and rate, to complement the well-established methods used in MEG, i.e., a scanning pill-counter at the collimator and a scanning APD detector at the center of the spectrometer. The first is made by scintillating fibers and mounted at the entrance to the spectrometer. The second is a luminophore foil detector (CsI on a Mylar support) with a CCD camera installed at the intermediate focus collimator system. The third one is a matrix detector consisting of a 9 × 9 matrix of small scintillation counters, which is used at the center of the spectrometer.

2.3. The Constant Bending Radius Magnet

The COBRA (COnstant Bending RAdius) solenoid, inherited from MEG is a superconducting magnet with a graded magnetic field in the axial direction, ranging from 1.27 T at the center to 0.49 T at both ends of the magnet cryostat. Compared to a uniform solenoidal field, COBRA has the advantage that particles produced with small longitudinal momentum have a shorter latency time in the spectrometer, allowing operations at a high rate. Moreover, the field was designed so that positrons emitted from the target followed a trajectory with an almost constant projected bending radius, weakly dependent on the emission polar angle. The magnet is equipped with a pair of compensation coils to reduce the stray field to a suitable level for LXe PMT operations.

2.4. The Detector

2.4.1. The Cylindrical Drift Chamber (CDCH)

The new MEG II positron tracker is a cylindrical drift chamber, 191 cm long, with the radial extent ranging from 17 to 29 cm. Nine layers of 192 drift cells, each a few mm wide, allows a high granularity measurement of the positron trajectories. The wires form an angle with the CDCH axis, varying from 6 in the innermost layer to 8.5 in the outermost one. The stereo angle depends on the layer with an alternating sign, allowing reconstruction of the longitudinal hit coordinate. The drift cell shape is quasi-square, with a 20 μ m Au-plated W sense-wire surrounded by 40 or 50 μ m Ag-plated Al field wires, with a 5:1 field-to-sense wire ratio. The sensitive volume is filled with a low-mass He:iC 4 H 10 (90:10) gas mixture (plus some small quantities of additives to improve the operational stability), which allows a good compromise between low scattering (radiation length 1.5 × 10 3 X 0 ) and single-hit resolution (<120 μ m, measured on prototypes [18]).
Wire breaking problems arose during the CDCH assembly and commissioning, even if the operations were performed inside clean rooms under controlled environmental conditions. The problem affected aluminum wires and was investigated by performing optical inspections with microscopes, chromatography, practical tests and SEM/EDS (Scanning Electron Microscopy with Energy Dispersive Spectroscopy) analyses. The origin of the breaking phenomenon turned out to be chemical corrosion of the Al core in presence of water condensation on wires from ambient humidity. The development of corrosion was stopped by keeping the wire volume in an absolutely dry atmosphere with a continuous flow of inert gas (nitrogen or helium). The total number of broken cathode wires amounted to 97 (10) for diameters of 40 (50) μ m. The effect of a missing cathode wire on the reconstruction was evaluated by means of Garfield [19] and ANSYS [20] simulations and found to be negligible.
The front-end electronics are constituted by custom boards, which also provide HV to the wires and read out signals at both wire ends. The active electronics are cooled by a chiller system. The HV is supplied by a commercial HV [21] system made of high precision modules, while a custom-made gas system [22] flushes and monitors the gas mixture to the chamber and helium to COBRA.

2.4.2. The Pixelated Timing Counter (pTC)

The e + time must be measured with a resolution of about 40 ps, at a hit rate of 5 MHz. This is needed to allow for a precise measurement of the relative time of the e + γ pair. The e + time is measured by the pTC which also generates trigger signals by providing prompt timing and direction information. The pTC is composed of 512 scintillator counters organized in two sectors, placed on the upstream side and downstream side of the target, each consisting of 256 counters. The signal e + , curled by the COBRA magnetic field, hits on average nine counters, significantly improving the total time resolution compared to MEG. The background hit rate on a single counter is less than 100 kHz even though the beam intensity is twice as high as in MEG.
A single counter is composed of a fast plastic scintillator plate, BC-422 ( 120 × 40 × 5 mm 3 or 120 × 50 × 5 mm 3 depending on the counter location) and six SiPMs connected in series at each end. The resolution of all counters was measured with a 90 Sr source and found to be below 100 ps. The average resolution was 72 (81) ps for 120 × 40 ( 50 ) × 5 mm 3 counters. The variations of the counter resolutions come from the variations of the photon detection efficiency amongst SiPMs and variations of the light yield between the scintillators.
The pTC is equipped with a cooling system to reduce dark count rates in the SiPMs [23] and a laser calibration system [24] to calibrate the timing between the counters.

2.4.3. The Liquid Xenon Calorimeter (LXe)

The liquid xenon detector allows the measurement of the photon’s energy, time and position. The MEG LXe detector was one of the world’s largest xenon based detectors with 900 L of liquid xenon. It was surrounded by 846 2-inch PMTs to detect the scintillation light emitted in the VUV range at a temperature of 165 K. Its performance was limited due to the non-uniform coverage of the PMT sensitive area, especially on the photon entrance face. While we reuse the liquid xenon and the cryostat in MEG II, we have upgraded it by replacing the PMTs on the photon entrance face with smaller photosensors, in order to reduce the non-uniformity of the detector response. In addition, the layout of the PMTs on the lateral faces was rearranged to improve light collection for events near the lateral walls. We inherited several methods to carefully calibrate the detector from MEG, with some modifications to match the upgraded configuration. The designed energy and position resolutions are ∼0.5 MeV and ∼2.5 mm, respectively.
A new type of VUV sensitive SiPM, the VUV-MPPC [25], was developed. It has high immunity to magnetic fields and is sensitive to single photons, which enables an easier and more reliable calibration of the detector. Moreover, a finer readout granularity allows more precise reconstruction of events where the initial high energy photons convert near the entrance face (shallow events) and to reduce pile-up of these photons in the same acquisition window. Thus, 216 PMTs at the entrance face were replaced by 4092 SiPMs. A VUV-MPPC consists of four 6 × 6 mm 2 sensor chips and the signals are read out by one channel per SiPM, connecting the four chips in series inside a PCB.

2.4.4. The Radiative Decay Counter (RDC)

Photons contributing to accidental background come from either RMD, bremsstrahlung, or positron AIF. The AIF background decreases in MEG II because of the reduced mass of the CDCH compared to the MEG drift chambers. The yield of the AIF background photons above 48 MeV per muon decay expected from MC simulations is 1.4 ×10 6 , to be compared with 2.3 × 10 6 in MEG. On the other hand, the RMD photon background is unchanged. The RDC, a detector newly installed in MEG II, is capable of identifying a fraction of low-energy positrons from RMD decays with photon energies close to the kinematic limit.
The detector consists of 12 plastic scintillators BC-418 for time measurement and 76 LYSO crystals for energy measurement. It is placed 140 cm downstream of the muon stopping target and covers the region within 10 cm from the beam axis. It identifies low-energy positrons (1–5 MeV) in time coincidence with the detection of a high energy photon in the LXe detector. According to simulations, it can detect about ∼40% of RMD background with E γ > 48 MeV improving the sensitivity of the μ + e + γ search by ∼10%.
If we place another module of RDC upstream of the target, we can identify the other half of the RMD positrons emitted towards the upstream side. However, the upstream RDC is technically challenging because it must be placed in the beam path. Several detector techniques have been examined; currently, the best candidate is a resistive plate chamber with extremely low-mass diamond-like carbon based electrodes. Though the upstream RDC is not included in the baseline design of MEG II at this moment, intensive R&D work is underway.

2.5. The Trigger and Data Acquisition System

The readout electronics was re-designed to deal with a factor of two increase in muon stopping rate from MEG and an almost tripling in the number of electronic channels. Since the physical space for the electronics could not be expanded, the only option was to combine the previous DAQ, trigger and high-voltage systems into a combined single system. The result is the new WaveDAQ system [26], consisting of four distinguished subsystems. First, the system is housed in a full custom crate with an integrated remote controlled power supply (24 V/360 W) and a custom backplane featuring gigabit serial links in a dual star topology. The second system is the WaveDREAM (Drs4 based REAdout Module) board (WDB), which on one hand digitizes all input signals with a sampling frequency up to 5 GSPS/12 bit and on the other hand performs continuous trigger operations at 80 MHz such as summing all input channels and comparing the result to a predefined threshold. In addition, the WDB can house a high voltage generator for SiPM biasing. The third system is a dedicated Trigger Concentrator Board (TCB), which receives all digital trigger data from one crate, combines and processes them, then sends the result to a central trigger system also consisting of the same board type, where the final global trigger decision is made every 12.5 ns. The fourth system is the Data Concentrator Board (DCB), which receives the DAQ data stream from each WDB and sends it to the central computer via Gigabit Ethernet lines.

3. Results

In this section the results of the engineering runs carried out in the years 2017–2020 are reported.
During the engineering runs, the muon beam was delivered at different intensities, up to the MEG II nominal intensity ( 7 × 10 7   μ + /s) and stopped in the target. Beam tuning and measurements were performed taking advantage of the new beam monitoring tools, which allowed continuous measurement of the beam conditions. The scintillating fiber and matrix detectors were frequently used during the runs, providing the beam profile and rate consistent with those measured by the pill and APD counter at the beginning of the runs. The target camera system was extensively tested and showed the capability to detect target displacements at the required precision level (100 μ m for movements along the direction transverse to the target plane). In 2020, a π p charge-exchange run was performed, in which the beam was changed from positive muons to 70.5 MeV/c negative pions and the muon stopping target was replaced with a liquid hydrogen target to induce the reaction π p π 0 n . Photons from the subsequent decay π 0 γ γ were used for LXe calibration purposes.
The DAQ and the trigger system operated successfully enabling stable data taking during the runs. A limited number of DAQ channels was available because the mass production of the full electronics was not yet started, awaiting the detailed results from the prototype test. Various problems were identified and solved; in particular common-mode noise was reduced to a level that does not affect the photon energy resolution. The trigger was successfully commissioned and performed well during the runs. The full system has been deployed in spring 2021 with about 9000 channels reading out the whole MEG II detectors. A first test showed that the system could run at an event rate of 50 Hz resulting in a data stream of about 8 GBit/s.
The pTC was the first sub-detector that was fully commissioned already during the first engineering run, in 2017. It operated in the nominal MEG II muon beam. The background e + hit rate was confirmed to be less than 100 kHz as expected. Radiation damage to SiPMs was carefully studied. The dark current in the SiPMs gradually increased during the use of the beam while the impact on the time resolution was kept under control by cooling the detector to 10 C. The e + time was reconstructed by combining hit times after subtracting the time-of-flight between counters. It is expected to be calculated from the e + trajectory given by the CDCH but the highly segmented design of pTC allows the tracking of e + using the hit pattern. This pTC-alone tracking was used to evaluate the pTC time resolution; Figure 2 shows the resolution measured with the 2017 data as a function of the number of hits. The measured resolution was worse than the expectation from the tests with a 90 Sr source made before the installation because of the larger noise level in the MEG II environment. By weighting with the distribution of the number of the hits for the signal e + s obtained by a Monte Carlo (MC) simulation and correcting for a bias in the evaluation method, the overall time resolution was estimated to be 35 ps, a factor of two better than that in MEG (76 ps). After 3 years of physics run, the resolution is expected to deteriorate to 41 ps due to the radiation damage, still satisfying the requirement of about 40 ps. In the 2018–2020 runs, only one sector, either upstream or downstream, was installed to provide a trigger to the DAQ for taking Michel e + data together with the CDCH.
The CDCH has been integrated into the experiment since 2018. Due to the limited DAQ channels available during the engineering runs a complete particle tracking test was not possible. Nonetheless, the chamber was tested with cosmic rays and different beam intensities up to the MEG II nominal one; tests with different gas mixtures were also performed. A few wire breaks, due to corrosion that started before the chamber was put in an inert atmosphere, occurred during the runs, requiring electrical insulation of part of the chamber to continue data-taking. The broken wires had to be removed at the end of the beam periods after re-opening the detector. Anomalous high currents (up to 300 μ A) were observed in parts of the CDCH during the runs, that were cured after optimizing the gas mixture and operating the chamber with up to 2% oxygen. As a result, the chamber was operated in stable conditions at the full MEG II beam intensity with He/isobutane (90/10) + isopropyl alcohol (1%) + O 2 (0.5%) for a period of about one week during the 2020 run. Front-end (FE) electronics with three different amplifications were used to investigate the optimum electronics gain value. The total gain (gas gain × FE gain) was measured using cosmic ray data by comparing amplitude and charge spectra to MC simulations which contained a detailed description of ionization patterns, drift and diffusion from Garfield, single-electron pulse shapes from prototype measurements under UV laser light as well as the noise spectra measured on data. The configuration with the highest FE gain was chosen and is currently under implementation (a modification of all the FE boards was required) since it has been demonstrated that it provided the best signal to noise ratio. The gas gain was extracted from the currents of the HV power supplies, taking into account the hit rate that was also directly measured from data. The measured gain was in the range of 4–7 × 10 5 in the mixture He/isobutane 90/10 plus 1% isopropyl alcohol. Studies with simulations and data were carried out to understand the impact of oxygen to optimize the gas mixture. We do not expect a significant deterioration of the chamber performances caused by the use of oxygen concentrations at the level 0.5% or below considering also the expected improvements deriving from the larger signal to noise ratio of the modified electronics.
The limited number of hits on reconstructed tracks, due to the small number of readout channels, did not allow an accurate measurement of the spatial resolution and the CDCH performances. Thus, the expected performances of the fully instrumented detector were evaluated by MC simulations, updated with noise, gain and electronics response observed in data. Positron efficiency ϵ e + and angular resolution σ θ e + are currently worse than the design (65% and 6.7 mrad vs. 70% and about 5 mrad, respectively), the vertex resolutions are already at the design level (about 1.7 mm and 0.8 mm along the beam and vertical directions, respectively) while the momentum resolution σ p e + is better than the design (100 keV/c vs. 130 keV/c). The efficiency loss is due to the higher than expected noise levels. We observed both low-frequency coherent noise over sets of wires on the same WaveDREAM board as well as incoherent noise extending to higher frequencies. Digital filters were developed to cope with it, but an intensive campaign for the identification and suppression of hardware noise sources is foreseen for the 2021 engineering run. The improved momentum resolution is due to an improved track fitting method for tracks with multiple turns in the chamber.
The upgraded LXe detector was put into operation in 2017 and calibrated with sources, LED, low-energy photons from a nuclear reaction using a Cockcroft–Walton accelerator [27] and finally with high-energy photons in the 2020 π p charge-exchange run. The RMD and background spectra were acquired with the muon beam at different intensities. The long-term stability of the PMTs and MPPCs was investigated and degradation of the photosensor performance was observed under muon beam conditions. A PMT gain decrease at the MEG II nominal intensity was first observed to be faster than the expected based on the MEG experience in 2018, but the deterioration speed moderated year by year. We decided to halve the gain to further mitigate the gain decrease and the current gain decrease rate measured in 2020 (0.16%/day) was now consistent with the expectation, which enabled the PMTs to operate over the full MEG II data taking period.
The MPPC VUV light sensitivity (PDE) was measured to be 12% in the experiment, lower than that measured in laboratory tests, which was >15%. Moreover, the PDE further decreased during the runs. In the 2020 run, the deterioration speed was lower than in 2019 (0.03 %/h instead of 0.06 %/h, where 0.03 %/h means that the PDE will go to zero in 100/(0.03 %/h), about 139 days.). Possible causes considered are the hole trapping at the sensor surface induced by VUV light or irradiation at low temperature. A possible solution was identified by annealing at high temperatures which removes accumulated charges. The effect was demonstrated with several MPPCs in the LXe detector; the PDEs recovered to values >15%. The optimum annealing strategy for all MPPCs is currently under investigation. Since the timing resolution would be dominated by PMT timing resolution in the case of the small MPPC PDE, MC simulations showed that almost no degradation of the MEG II sensitivity occurs down to 5% PDE, which is good enough to guarantee the sensitivity reach. If the degradation is saturated at this value, annealing may not be necessary.
A preliminary estimate of the detector performance was made. The position resolution σ x γ was measured using photons from the 7 Li ( p , γ ) 8 Be reaction excited using the Cockcroft–Walton proton accelerator with a 5 mm wide slit collimator. Improvements in both horizontal and vertical resolution were observed ( σ x γ 2.5 mm), twice better than MEG. The energy resolution σ E γ was evaluated by fitting background photon spectra from muon decays using the MC distributions as shown in Figure 3 and from monochromatic 55 MeV photons from the π p charge-exchange reaction, which provide consistent results of σ E γ / E γ = 1.7 % in the signal region. The improved granularity of the photon incident face almost eliminated the depth dependence of the energy resolution. The resolution seems at the level of the worst-case scenario considered in [14]; however, the data analysis is still preliminary and the results were evaluated using a limited number of electronics channels. The intrinsic time resolution, which was estimated by the difference of times reconstructed by even and odd channels, is 39 ps, consistent with the MC expectation.
The downstream RDC successfully demonstrated the RMD detection capability under muon beam conditions during the engineering runs, by observing positrons in time coincidence with high energy photons in LXe. The peak in Figure 4 represents coincident hits by a γ and a positron from RMD. The fraction of events at the time coincident peak was found to be consistent with the expectation from simulation. An upstream RDC prototype based on resistive plate chambers with diamond-like-carbon based electrodes is being developed and showed a good efficiency and time resolution in a laboratory. The first tests under muon beam conditions were performed during the 2020 run.

4. Discussion

All the detector components were integrated into the MEG II experiment and tested in the engineering runs, which enables us to identify and solve issues prior to the start of physics data collection. Based on the results from the engineering runs, we re-evaluated the projected sensitivity of MEG II. The updated detector performance is summarized in Table 1 (MEG II updated) in comparison with that in MEG and MEG II design. We input these updated values to the sensitivity calculation; however, note that some of the parameters were evaluated on more realistic MC simulations, updated on the basis of observation on data, and others on the data but the results are still preliminary. The following is assumed here:
  • The DAQ time is 20 weeks per year with 84% live time for three years;
  • The muon stopping rate is R μ + = 7 × 10 7   μ + /s;
  • The PDE of VUV-MPPCs is at a constant value of 6%;
  • The upstream RDC is not included.
The sensitivity ( S 90 ) is defined as the median of 90% confidence branching-ratio upper limits for the background-only hypothesis. In three years, we will reach a sensitivity of S 90 = 6 × 10 14 , which is roughly one order of magnitude better than MEG’s sensitivity of 5.3 × 10 13 . The expected single event sensitivity is 1 × 10 14 . In the realistic evaluation of the detector performance, some parameters improved while others became worse than the assumed design parameters; however, the overall projected sensitivity stays almost the same as that reported in Ref. [14]. The data statistics is 5.6 times larger than MEG, and the remaining improvement of a factor of 1.6 comes from the lower background due to the improved detector resolutions. Actually, the sensitivity improvement is not yet saturated and an extension of DAQ time can improve the sensitivity reach further.
If the degradation of VUV-MPPC PDE is not saturated at a value above 5%, we have to recover it by the annealing process. Because annealing all the MPPCs takes a couple of months, it can be done only during the annual accelerator shutdown period. In this case, the obtainable statistics is limited by the accumulated beam intensity up to which the MPPC is operational. For the given statistics, it is better to run for the full beam time at a reduced beam intensity since the number of accidental background events has a quadratic dependence on the beam intensity. Moreover, the detector performance is better with a lower beam intensity due to a low pileup probability. For example, the positron efficiency improves from 65% at 7 × 10 7   μ + /s to 74% at 3.5 × 10 7   μ + /s with the current reconstruction algorithms. As a result, it was shown by a preliminary study that a comparable sensitivity is achievable at a stopping rate down to 3.5 × 10 7   μ + /s, with which the VUV-MPPCs can operate for 20 weeks of DAQ based on the degradation speed observed in the engineering runs.
Cathode and field wires of the CDCH suffer from corrosion caused by exposure to humidity during construction. The risk of possible wire breakings during data taking imposed the construction of a new backup chamber with different cathode and field wires. The new chamber will use more robust 50 μ m diameter Ag-coated aluminum wires, instead of the old 40 μ m ones, for which the final drawing process, which was verified to weaken the wires and make them more prone to corrosion, will be avoided. Construction of the new chamber will take approximately 18 months; in the meantime, the existing one will be used in the experiment.
At the time of writing, preparation for the 2021 engineering run is underway. In this run, all the detector signals will be collected for the first time with the full readout electronics. The stability and performance of the existing CDCH and the long-term behavior of the PDE degradation of the VUV-MPPCs will be extensively studied. We plan to collect data with a μ + e + γ trigger at the end of the run, which will be the first physics data of MEG II.

Author Contributions

Investigation, A.M.B., V.B., M.B., G.B., P.W.C., G.C. (Gianluca Cavoto), F.C. (Fabrizio Cei), M.C., G.C. (Gianluigi Chiarello), A.C., F.C. (Federica Cuna), G.d.M., A.d.B., M.D.G., M.F., L.G., G.G., F.G. (Flavio Gatti), F.G. (Francesco Grancagnolo), M.G., D.N.G., M.H., K.I., F.I., T.I., P.-R.K., N.K. (Nikolay Khomutov), S.K., A.K., N.K. (Nikolay Kravchuk), V.K., N.K. (Nikolay Kuchinskiy), W.K., T.L., V.M., M.M., S.M., W.M., T.M., A.M., M.N., D.N., H.N., S.O., R.O., W.O., A.O., D.P., M.P., A.P. (Angela Papa), V.P., A.P. (Alexander Popov), F.R., S.R., M.R., A.R., P.S., K.S., G.S., A.S., G.F.T., K.T., Y.U., M.U., C.V., K.Y. (Kosuke Yanai), K.Y. (Kensuke Yamamoto), T.Y. and Y.V.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by DOE DEFG02-91ER40679 (USA), INFN, MIUR Montalcini D.M. 2014 n. 975 (Italy), JSPS Core-to-Core Program, A. Advanced Research Networks JPJSCCA20180004, JSPS KAKENHI Grant Numbers JP26000004, JP18K13557, JP19J13635, JP19J21535, JP19J21730, JP20H00154, JP21H00065 (Japan), the Russian Federation Ministry of Science and Higher Education (Russia), and Schweizerischer Nationalfonds (SNF) Grant 200020_172706 and Grant 206021_177038 (Switzerland).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study can be provided by the authors upon reasonable request.

Acknowledgments

We are grateful for the support and co-operation provided by PSI as the host laboratory and to the technical and engineering staff of our institutes.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Calibbi, L.; Signorelli, G. Charged Lepton Flavor Violation: An experimental and theoretical introduction. Riv. Del Nuovo C. 2018, 41, 71–174. [Google Scholar] [CrossRef]
  2. Abi, B.; Albahri, T.; Al-Kilani, S.; Allspach, D.; Alonzi, L.P.; Anastasi, A.; Anisenkov, A.; Azfar, F.; Badgley, K.; Baeßler, S.; et al. Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm. Phys. Rev. Lett. 2021, 126, 141801. [Google Scholar] [CrossRef] [PubMed]
  3. Calibbi, L.; López-Ibáñez, M.L.; Melis, A.; Vives, O. Implications of the Muon g-2 result on the flavour structure of the lepton mass matrix. arXiv 2021, arXiv:2104.03296. [Google Scholar]
  4. Lindner, M.; Platscher, M.; Queiroz, F. A call for new physics: The muon anomalous magnetic moment and lepton flavor violation. Phys. Rep. 2018, 731, 1–82. [Google Scholar] [CrossRef] [Green Version]
  5. Baldini, A.M.; Bao, Y.; Baracchini, E.; Bemporad, C.; Berg, F.; Biasotti, M.; Boca, G.; Cascella, M.; Cattaneo, P.W.; Cavoto, G.; et al. Search for the lepton flavour violating decay μ+→e+γ with the full dataset of the MEG experiment-MEG Collaboration. Eur. Phys. J. C 2016, 76, 434. [Google Scholar] [CrossRef] [Green Version]
  6. Kinsho, M.; Ikegami, M.; Kawamura, N.; Kobayashi, H.; Matsumoto, H.; Mibe, T.; Miyake, Y.; Mihara, S.; Nishiguchi, H.; Ohomori, C.; et al. Proposal of an Experimental Search for μ→e Conversion in Nuclear Field at Sensitivity of 10−14 with Pulsed Proton Beam from RCS. 2010. Available online: https://j-parc.jp/researcher/Hadron/en/pac_1101/pdf/KEK_J-PARC-PAC2010-13.pdf (accessed on 17 August 2021).
  7. Abramishvili, R.; Adamov, G.; Akhmetshin, R.R.; Allin, A.; Angélique, J.C.; Anishchik, V.; Aoki, M.; Aznabayev, D.; Bagaturia, I.; Ban, G.; et al. (COMET collaboration) COMET Phase-I Technical Design Report. Prog. Theor. Exp. Phys. 2020, 2020, 033C01. [Google Scholar] [CrossRef] [Green Version]
  8. Bertl, W.; Engfer, R.; Hermes, E.A.; Kurz, G.; Kozlowski, T.; Kuth, J.; Otter, G.; Rosenbaum, F.; Ryskulov, N.M.; van der Schaaf, A.; et al. (SINDRUM II Collaboration). A search for μ→e conversion in muonic gold. Eur. Phys. J. C 2006, 47, 337–346. [Google Scholar] [CrossRef]
  9. Abrams, R.; Alezander, D.; Ambrosio, G.; Andreev, N.; Ankenbrandt, C.M.; Asner, D.M.; Arnold, D.; Artikov, A.; Barnes, E.; Bartoszek, L.; et al. (Mu2e Collaboration). Mu2e Conceptual Design Report. arXiv 2012, arXiv:1211.7019. [Google Scholar]
  10. Arndt, K.; Augustin, H.; Baesso, P.; Berger, N.; Berg, F.; Betancourt, C.; Bortoletto, D.; Bravar, A.; Briggl, K.; vom Bruch, D.; et al. Technical design of the phase I Mu3e experiment. arXiv 2021, arXiv:2009.11690. [Google Scholar]
  11. Bellgardt, U.; Otter, G.; Eichler, R.; Felawka, L.; Niebuhr, C.; Walter, H.K.; Bertl, W.; Lordong, N.; Martino, J.; Egli, S.; et al. Search for the decay μ+→e+e+e. Nucl. Phys. B 1988, 299, 1–6. [Google Scholar] [CrossRef]
  12. Kou, E.; Urquijo, P.; Altmannshofer, W.; Beaujean, F.; Bell, G.; Beneke, M.; Bigi, I.I.; Bishara, F.; Blanke, M.; Bobeth, C.; et al. The Belle II Physics Book. Prog. Theor. Exp. Phys. 2019, 2019, 123C01. [Google Scholar] [CrossRef]
  13. Zhou, H.; Zhang, R.Y.; Han, L.; Ma, W.G.; Guo, L.; Chen, C. Searching for τ→μγ lepton-flavor-violating decay at super Charm-Tau factory. Eur. Phys. J. C 2016, 76, 421. [Google Scholar] [CrossRef] [Green Version]
  14. Baldini, A.M.; Baracchini, E.; Bemporad, C.; Berg, F.; Biasotti, M.; Boca, G.; Cattaneo, P.W.; Cavoto, G.; Cei, F.; Chiappini, M.; et al. (MEG II Collaboration). The design of the MEG II experiment. Eur. Phys. J. C 2018, 78, 380. [Google Scholar] [CrossRef]
  15. Adam, J.; Bai, X.; Baldini, A.M.; Baracchini, E.; Bemporad, C.; Boca, G.; Cattaneo, P.W.; Cavoto, G.; Cei, F.; Cerri, C.; et al. The MEG detector for μ+→e+γ decay search. Eur. Phys. J. C 2013, 73, 2365. [Google Scholar] [CrossRef] [Green Version]
  16. Palo, D.; Hildebrandt, M.; Hofer, A.; Kyle, W.; Lad, D.; Libeiro, T.; Molzon, W. Precise Photographic Monitoring of MEG II Thin-film Muon Stopping Target Position and Shape. Nucl. Instrum. Methods Phys. Res. Sect. A 2019, 944, 162511. [Google Scholar] [CrossRef] [Green Version]
  17. Cavoto, G.; Chiarello, G.; Hildebrandt, M.; Hofer, A.; Ieki, K.; Meucci, M.; Milana, S.; Pettinacci, V.; Renga, F.; Voena, C. A photogrammetric method for target monitoring inside the MEG II detector. Rev. Sci. Instrum. 2021, 92, 043707. [Google Scholar] [CrossRef] [PubMed]
  18. Baldini, A.M.; Baracchini, E.; Cavoto, G.; Cascella, M.; Cei, F.; Chiappini, M.; Chiarello, G.; Chiri, C.; Dussoni, S.; Galli, L.; et al. Single-hit resolution measurement with MEG II drift chamber prototypes. J. Instrum. 2016, 11, P07011. [Google Scholar] [CrossRef]
  19. Garfield++ – Simulation of Gaseous Detectors. Available online: https://garfieldpp.web.cern.ch/garfieldpp/ (accessed on 17 August 2021).
  20. Ansys® Electronic Desktop 2019. Available online: https://www.ansys.com/ (accessed on 17 August 2021).
  21. ISEG EHS 8630p-305F. Available online: http://www.iseg-hv.com (accessed on 17 August 2021).
  22. Baldini, A.M.; Baracchini, E.; Cavoto, G.; Cei, F.; Chiappini, M.; Chiarello, G.; Chiri, C.; Francesconi, M.; Galli, L.; Grancagnolo, F.; et al. Gas Distribution and Monitoring for the Drift Chamber of the MEG-II Experiment. J. Instrum. 2018, 13, P06018. [Google Scholar] [CrossRef] [Green Version]
  23. Boca, G.; Cattaneo, P.W.; De Gerone, M.; Gatti, F.; Nakao, M.; Nishimura, M.; Ootani, W.; Rossella, M.; Uchiyama, Y.; Usami, M.; et al. Timing resolution of a plastic scintillator counter read out by radiation damaged SiPMs connected in series. Nucl. Instrum. Methods Phys. Res. Sect. A 2021, 999, 165173. [Google Scholar] [CrossRef]
  24. Boca, G.; Cattaneo, P.W.; De Gerone, M.; Francesconi, M.; Galli, L.; Gatti, F.; Koga, J.; Nakao, M.; Nishimura, M.; Ootani, W.; et al. The laser-based time calibration system for the MEG II pixelated Timing Counter. Nucl. Instrum. Methods Phys. Res. Sect. A 2019, 947, 162672. [Google Scholar] [CrossRef] [Green Version]
  25. Ieki, K.; Iwamoto, T.; Kaneko, D.; Kobayashi, S.; Matsuzawa, N.; Mori, T.; Ogawa, S.; Onda, R.; Ootani, W.; Sawada, R.; et al. Large-area MPPC with enhanced VUV sensitivity for liquid xenon scintillation detector. Nucl. Instrum. Methods Phys. Res. Sect. A 2019, 925, 148–155. [Google Scholar] [CrossRef] [Green Version]
  26. Galli, L.; Baldini, A.M.; Cei, F.; Chiappini, M.; Francesconi, M.; Grassi, M.; Hartmann, U.; Meucci, M.; Morsani, F.; Nicolò, D.; et al. WaveDAQ: An highly integrated trigger and data acquisition system. Nucl. Instrum. Methods Phys. Res. Sect. A 2019, 936, 399–400. [Google Scholar] [CrossRef]
  27. Adam, J.; Bai, X.; Baldini, A.; Baracchini, E.; Bemporad, C.; Boca, G.; Cattaneo, P.W.; Cavoto, G.; Cei, F.; Cerri, C.; et al. Calibration and monitoring of the MEG experiment by a proton beam from a Cockcroft–Walton accelerator. Nucl. Instrum. Methods Phys. Res. Sect. A 2011, 641, 19–32. [Google Scholar] [CrossRef]
Figure 1. A scheme of the MEG II experiment. Photons and positrons from the decay of muons stopped in a target placed at the center of the detector system are measured by the liquid xenon photon detector and the positron spectrometer consisting of a superconducting magnet, a cylindrical drift chamber, and two sets of pixelated timing counter, respectively. A new detector identifying background photons (RDC) is installed on the beam axis.
Figure 1. A scheme of the MEG II experiment. Photons and positrons from the decay of muons stopped in a target placed at the center of the detector system are measured by the liquid xenon photon detector and the positron spectrometer consisting of a superconducting magnet, a cylindrical drift chamber, and two sets of pixelated timing counter, respectively. A new detector identifying background photons (RDC) is installed on the beam axis.
Symmetry 13 01591 g001
Figure 2. The pTC time resolution as a function of the number of hit counters obtained from the 2017 run. The signal e + hits nine counters on average.
Figure 2. The pTC time resolution as a function of the number of hit counters obtained from the 2017 run. The signal e + hits nine counters on average.
Symmetry 13 01591 g002
Figure 3. A background photon spectrum measured in the 2019 run at the MEG II beam intensity (black points) as well as MC simulation spectra convoluted with different energy resolution assumptions. The red histogram shows the 1.7% resolution case (best fit), and green/blue show 1.0%/2.3%, respectively.
Figure 3. A background photon spectrum measured in the 2019 run at the MEG II beam intensity (black points) as well as MC simulation spectra convoluted with different energy resolution assumptions. The red histogram shows the 1.7% resolution case (best fit), and green/blue show 1.0%/2.3%, respectively.
Symmetry 13 01591 g003
Figure 4. Distribution of time difference between an e + measured in RDC and a γ in LXe detector for the events triggered by a high energy γ at the LXe detector (black). The peak corresponds to the RMD events. The flat time distribution corresponds to accidental hits of high energy e + from Michel decays, which is reduced by applying a cut on the e + energy measured with the LYSO crystals of RDC ( E e + LYSO < 8 MeV, red histogram).
Figure 4. Distribution of time difference between an e + measured in RDC and a γ in LXe detector for the events triggered by a high energy γ at the LXe detector (black). The peak corresponds to the RMD events. The flat time distribution corresponds to accidental hits of high energy e + from Michel decays, which is reduced by applying a cut on the e + energy measured with the LYSO crystals of RDC ( E e + LYSO < 8 MeV, red histogram).
Symmetry 13 01591 g004
Table 1. Summary of detector performance and sensitivity. σ X denotes the resolution of variable X and ϵ Y denotes the detection efficiency of particle Y. Two values for σ E γ are values for shallow (<2 cm)/deep (>2 cm) events. σ t e + γ is the coincidence time resolution for an e + γ pair and evaluated by combining the resolutions of the positron spectrometer and the LXe photon detector. S 90 denotes the sensitivity of the experiment; see the text for the definition.
Table 1. Summary of detector performance and sensitivity. σ X denotes the resolution of variable X and ϵ Y denotes the detection efficiency of particle Y. Two values for σ E γ are values for shallow (<2 cm)/deep (>2 cm) events. σ t e + γ is the coincidence time resolution for an e + γ pair and evaluated by combining the resolutions of the positron spectrometer and the LXe photon detector. S 90 denotes the sensitivity of the experiment; see the text for the definition.
R μ + σ p e + σ θ e + σ E γ σ x γ σ t e + γ ϵ e + ϵ γ S 90
MEG 3 × 10 7 s 1 380 keV/c9.4 mrad2.4%/1.7%5 mm122 ps30%63% 5.3 × 10 13
MEG II design 7 × 10 7 s 1 130 keV/c5.3 mrad1.1%/1.0%2.4 mm84 ps70%69% 6 × 10 14
MEG II updated 7 × 10 7 s 1 100 keV/c6.7 mrad1.7%/1.7%2.4 mm70 ps65%69% 6 × 10 14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baldini, A.M.; Baranov, V.; Biasotti, M.; Boca, G.; Cattaneo, P.W.; Cavoto, G.; Cei, F.; Chiappini, M.; Chiarello, G.; Corvaglia, A.; et al. The Search for μ+ → e+γ with 10–14 Sensitivity: The Upgrade of the MEG Experiment. Symmetry 2021, 13, 1591. https://doi.org/10.3390/sym13091591

AMA Style

Baldini AM, Baranov V, Biasotti M, Boca G, Cattaneo PW, Cavoto G, Cei F, Chiappini M, Chiarello G, Corvaglia A, et al. The Search for μ+ → e+γ with 10–14 Sensitivity: The Upgrade of the MEG Experiment. Symmetry. 2021; 13(9):1591. https://doi.org/10.3390/sym13091591

Chicago/Turabian Style

Baldini, Alessandro M., Vladimir Baranov, Michele Biasotti, Gianluigi Boca, Paolo W. Cattaneo, Gianluca Cavoto, Fabrizio Cei, Marco Chiappini, Gianluigi Chiarello, Alessandro Corvaglia, and et al. 2021. "The Search for μ+ → e+γ with 10–14 Sensitivity: The Upgrade of the MEG Experiment" Symmetry 13, no. 9: 1591. https://doi.org/10.3390/sym13091591

APA Style

Baldini, A. M., Baranov, V., Biasotti, M., Boca, G., Cattaneo, P. W., Cavoto, G., Cei, F., Chiappini, M., Chiarello, G., Corvaglia, A., Cuna, F., dal Maso, G., de Bari, A., De Gerone, M., Francesconi, M., Galli, L., Gallucci, G., Gatti, F., Grancagnolo, F., ... Yudin, Y. V. (2021). The Search for μ+ → e+γ with 10–14 Sensitivity: The Upgrade of the MEG Experiment. Symmetry, 13(9), 1591. https://doi.org/10.3390/sym13091591

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop