A comparative study of target fabrication strategies for microgram muonic atom spectroscopy
We use essential cookies to make sure the site can function. We also use optional cookies for advertising, personalisation of content, usage analysis, and social media.
By accepting optional cookies, you consent to the processing of your personal data – including transfers to third parties. Some third parties are outside of the European Economic Area, with varying standards of data protection.
See our privacy policy for more information on the use of your personal data.
Manage preferences for further information and to change your choices.
Advertisement
Scientific Reports
15, Article number: 6939 (2025) Cite this article
Metrics
Muonic atom spectroscopy is a method that can determine absolute nuclear charge radii with typical relative precision of
. Recent developments have enabled to extend muonic atom spectroscopy to microscopic target quantities as low as
μ
. This substantial reduction from the traditional limit of the order of
is based on a transfer mechanism in a high-pressure hydrogen gas cell, which transports the muon to the surface of the target material rather than stopping it over a broad depth range. This approach enables the measurement of absolute nuclear charge radii of long-lived radioactive isotopes (half-life above
20 years), but the production of appropriate targets for the technique has presented some major challenges, such as the formation of organic layers on the substrate. This study presents a systematic investigation of the stopping efficiency for different target preparation methods: ion implantation, drop-on-demand printing, and molecular plating. Notable differences between the three methods were discovered in terms of their performance allowing to further fine tune the method of choice for future target preparations. Our findings show that implantation provides appropriate targets for our method with negligible losses. This achievement opens the landscape of potential measurements to isotopes where high mass separation is required not achievable with other methods. Furthermore, molecular plated targets performed substantially better than those prepared using drop-on-demand printing.
Muonic atom spectroscopy exploits the distinctive properties of the muon in order to study the atomic nucleus. The muon is a lepton with a mass approximately 207 times larger than that of the electron. When negative muons are brought close to a nucleus, the Coulomb attraction can lead to the formation of a muonic atom. One can understand the main differences between regular atoms and muonic atoms by replacing the mass of the electron by that of the muon in the Bohr model. Due to the large muon-to-electron mass ratio, the atomic binding energy is much higher than for electrons, while the Bohr radius is severely reduced compared to its electronic counterparts (both approximately by a factor
). Because of their reduced distance to the nucleus, muonic atoms are much more sensitive to nuclear finite size effects, by a factor
. This enhanced sensitivity forms the cornerstone of the muonic atom spectroscopy method. Unlike other atomic spectroscopy methods, which primarily probe changes in mean square radii via the isotope shift1, muonic atom spectroscopy allows to extract the absolute nuclear mean square charge radii directly from transition energies of the corresponding x rays. These absolute radii provide invaluable input for benchmarking radii attained from other techniques, such as laser spectroscopy, thus enhancing their precision2. An external input that is still required is the shape of the nucleus. This input can be extracted from elastic electron scattering, recent advancements of which show the feasibility of measurements on long-lived radioactive species3.
A central challenge in muonic atom spectroscopy lies in the formation of these exotic atoms. To capture a muon in the Coulomb potential of a nucleus, the muon’s energy must be at most a few electronvolt. Traditional methods rely on the direct stopping in the target material of interest. However, given the typical momentum of negative muon beams at accelerator facilities, the stopping range is of the order of a millimeter. Hence, the required amount of target material is typically at least
, corresponding to a target thickness of about 100 μm − 1 mm, depending on the target material. In the case of isotopically pure or radioactive samples, this is often not achievable.
Recent developments introduced a high-pressure hydrogen gas cell (100 bar) with a small deuterium admixture (
0.25%) in order to slow down muons sufficiently for atomic capture on the surface of a substrate4. First, the muon is stopped in a gas target of 100 bar H2 with a small deuterium admixture, forming muonic hydrogen (
) in the process. The muonic hydrogen atom will diffuse in the gas cell with a very short mean free path, until it encounters a deuterium atom. As the binding energy for a muon to the deuteron is slightly higher than that of the proton, the muon is transferred, such that muonic deuterium (
) is formed. Due to the Ramsauer-Townsend minimum in the scattering cross section of
in H25, the
observes an almost transparent gas cell, such that a large fraction of them can diffuse to the end of the gas cell, transferring the muon to the nucleus of interest. A schematic representation of the transfer process is shown in Fig. 1. This method reduces the required target material to approximately
μ
, opening the door to studies with long-lived radioactive isotopes (
20 years), as well as rare isotopes that are not available in larger quantities. The half-life limitation is at this time primarily limited by radioprotection limits and gamma background in the detector array emanating from the target.
Schematic representation of the transfer process explained above. (a) Muon enters the gas cell, (b) muon is stopped and forms
, (c)
diffuses until it gets close to D2, (d) muon is transferred to deuterium, (e)
diffuses to the back of the gas cell, (f) muonic atom of interest is formed and corresponding x rays are emitted6.
This work aims to build on previous developments by investigating the muonic x-ray yields for different target production methods. First, the attenuation of the diffusing
atoms through different layer thicknesses of graphite was investigated. The goal of these measurements was to show the feasibility of using implanted targets, and to probe the optimal implantation depth for measurements. This optimal depth still has to be combined with implantation limitations, which are specific to the implanted species of interest. Furthermore, the range of
atoms provides an indication of the importance of thin layers of organic material being present on top of the substrate. These organic layers may often form on top of substrates during chemical preparation, originating from the solutions used for the target preparation or contact with air. Next, samples implanted with gold or potassium were investigated and compared to surface deposited targets of similar masses. Gold is used as a benchmark for heavy elements, as it is stable and mono-isotopic. 39K, on the other hand, was measured as a test for a future measurement on 40K. For potassium, metallic surface deposition is not possible. Hence, chromium surface targets were used as a proxy. Finally, a comparison was made between barium targets prepared using molecular plating (MP)7,8,9 and drop-on-demand printing (DoD)10, both of which stand at the forefront of radioactive target development.
The systematic study was performed at the
E1 beamline of the Paul Scherrer Institute’s High Intensity Proton Accelerator (HIPA) facility11. Figure 2a shows the experimental apparatus used for this study. A set of scintillators is used for coincidence and veto logic. The first type is a muon veto, which collimates the beam and vetoes background outside of the acceptance of the gas cell. Downstream of this veto detector, a muon entrance detector is positioned. The main purpose of this detector is to relate the timing of measured x rays to the muon’s arrival time, such that the background can be substantially reduced using coincidence logic. Finally, a set of electron veto scintillator detectors is placed around the gas cell in order to detect Michel electrons (electrons originating from muon decay). By using these as veto signals, direct hits of these electrons on the germanium detectors and Bremsstrahlung generated by them can be substantially reduced.
The targets are placed at the back of a high-pressure gas cell made of aluminium, which is in turn closed off towards the muon beamline using a
μ
-thick carbon fiber window. To withstand a pressure of
, a titanium and carbon fiber support grid was added. This entrance window of the gas cell has a transmission of about 55% for muons, as shown in Ref.4. The gas injection is handled in a similar way as in some of our previous work4. A CAD drawing of the gas cell is shown in Fig. 2b.
Schematic of the setup4 (left) and CAD drawing of the gas cell6 (right).
An array of germanium detectors of various types is placed around the gas cell for the detection of the muonic x rays. The array was constructed together with the MIXE collaboration12. The detectors in the frame consisted of a Miniball cluster detector13; a telescope (coax + planar crystals); reverse electrode coaxial germanium (REGe) detectors with relative efficiencies of 95% (x1), 70% (x2), and 60% (x1); standard electrode coaxial germanium (SEGe) detectors with relative efficiencies of 50% (x3) and 75% (x1), and broad energy germanium (BEGe) detectors (x5).
The data acquisition was performed using SIS3316 digitizer modules from Struck Innovative Systeme. The modules have 16 channels each, with a sampling rate of
. The outputs of the plastic scintillators and germanium detectors were saved without hardware triggers, such that coincidence and veto logic could be applied offline after optimization. The energy determination of the germanium detectors was extracted from a trapezoidal filter on field programmable gate arrays (FPGAs) onboard of the digitizer modules. Besides the trigger time and energy, the first
μ
(300 samples) are also saved in order to calculate improved timestamps, as explained in the analysis methods section.
As a backing material, glassy carbon (SIGRADUR K) disks with a thickness of
and a diameter of
were purchased from Hochtemperatur-Werkstoffe GmbH. This material is a high-purity porous form of carbon with greater strength and surface quality than graphite. A low-Z material is preferred for the backing material, such that no muonic x-ray background is generated in the energy regions of interest. Furthermore, for the implanted targets, such materials should display less sputtering than higher proton number alternatives. Over the years, gold has become the reference for optimization of the setup. The primary reason for this is that it is mono-isotopic while being available in many different forms. Consequently, it was the isotope of interest chosen to investigate the attenuation of the x-ray signal through graphite and as a first proof of principle for implanted targets. Additionally, implanted targets of 39K were prepared as a test for a future campaign evaluating 39, 40, 41K, where 40K can only be attained at a sufficient purity through magnetic mass separation. A measurement of such a triplet could be used to benchmark existing laser spectroscopy experiments14. Finally, drop-on-demand and molecular plating targets of barium were prepared. Barium is a homologue of radium, for which the charge radius measurement is a key goal of our program15.
In the interest of having reference targets to compare to the implanted targets, we prepared several targets using physical vapor deposition. The first set of targets was prepared by evaporating the material under vacuum in a tungsten crucible heated electrically. The resulting targets had 50, 10, and 3
thick gold layers, amounting to 170, 34, and
μ
and a 25 and
chromium target (32 and
μ
). The
target served for optimization purposes, while the others acted as a reference for the gold and the potassium (chromium being reasonably close in Z to potassium) implanted targets, respectively. Additionally, we prepared targets via DC magnetron sputtering in an Ar atmosphere. They consisted of a
gold layer (thick enough such that all the
atoms that reach it transfer to gold) and a layer of graphite on top (0, 10, 25, 50, 100,
).
To evaluate implanted targets produced at different implantation energies, glassy carbon disks were implanted with stable 197Au at 4.5, 27, and 90
at the UK National Ion Beam Centre in Surrey. These energies were initially chosen to correspond to projected ranges of 10, 25, and 50
as predicted by SRIM16. However, due to sputtering effects, this range can be substantially be reduced17,18. Similarly, a
implantation of the stable 39K on glassy carbon was performed at Helmholtz-Zentrum Dresden-Rossendorf (HZDR). In order to investigate self-sputtering, where earlier implanted atoms are removed from the sample by implanting new ions, TRIDYN19 simulations were performed. Rutherford backscattering spectrometry (RBS) was performed as verification of the simulations. These measurements did not fully match our observations, likely due to a combination of the porosity of glassy carbon and off-diagonal elements in the surface binding energy matrix. A more detailed discussion of these simulations can be found in Ref.17. The RBS profile for the implanted potassium sample is shown in Fig. 3. From this spectrum, one can see that there is a substantial low-energy tail for potassium (deeper in the substrate). This tail is most likely caused by the high reactivity of potassium and/or temperature-induced diffusion in the sample during implantation. Consequently, only about 60% of the particles are located near the surface of the substrate. Furthermore, there is a large peak in the spectrum originating from oxygen (not present before implantation). We suspect that this is primarily caused by the strong reactivity of elemental potassium.
RBS spectrum for and empty glassy carbon sample (grey) and implanted 39K target (color). The data was obtained using an incident
4He+ beam and a detector placed at a scattering angle of
.
In order to evaluate the different chemical target preparation techniques, drop-on-demand (DoD) and molecular plating (MP) samples were prepared on glassy carbon backing with deposition masses of 5, 10, 25, and
μ
. Due to a tight experimental schedule and a clear global conclusion, the
μ
DoD target and the
μ
MP target were not measured. The molecular plated samples were prepared with the setup described in Ref.7. The sample preparation was performed in galvanostatic mode. MP parameters, such as voltage, current density and deposition time were defined to maximize the deposition yield while producing thin, uniform and homogeneous films. By using trace amounts of the radioactive 133Ba, the homogeneity and yield could be monitored. Radiographic imaging, as shown in Fig. 4 (left), indicates a homogeneous spatial distribution over the entire surface, except for the edges. The DoD samples were prepared using the method described in Ref.10. A commercial Ba ICP-MS-standard solution with a concentration of
μ
μ
was used as stock solution. This stock solution was evaporated and subsequently redissolved in
HNO3 and diluted with acetonitrile. The following printing sequence was used: 261 drops,
each, in a circular pattern with a diameter of
This results in a deposition of
μ
of the solution, which corresponds to
μ
of Ba. For the PipeJet printer by BioFluidix, the tip size “200-L” was used, at a stroke velocity of
μ
. In order to produce 4 targets with different amounts of Ba, this sequence was printed once (
μ
Ba), twice (
μ
Ba), five times (
μ
Ba), and ten times (
μ
Ba) on the individual glassy carbon backings. These backings were cleaned with deionized water, isopropanol and acetone before usage. Similar to the MP samples, radiography was performed to check the sample homogeneity. This time, stable 130Ba was irradiated with neutrons before printing to artificially induce radioactivity in the sample. The radiographic image, as shown in Fig. 4 (right), indicated that the activity is distributed along the substrate, but there are localized clusters of higher concentrations. These clusters were further investigated using scanning electron microscopy (SEM), which showed micocrystallites formation around the radioactivity hotspots, which is not favored for microgram muonic atom experiments. An overview of the targets used for this study is given in Table 1.
Radiography images of one of the MP samples (left) and one of the DoD samples (right). Apart from the sample edges, the MP sample shows an excellent homogeneity. The DoD sample, on the other hand, shows some clustering.
In this study, we investigated the quality of targets by evaluating the observed intensity of the
–
and
–
transitions in the muonic atoms of interest (39K, 52Cr, 138Ba, and 197Au). These lines are the most important, as they have a high intensity and display the strongest effect due to finite size effects. In order to transform count rates to atomic capture rates, the corresponding x-ray intensities are needed. Calculations for these intensities were performed using the muon cascade program developed by Akylas and Vogel20. The input parameters for this code include the atomic mass and the electronic binding energies for the K, L and M orbits, as referenced in Ref.21,22,23, as well as the effective nuclear charge for the same electronic shells24,25. Additionally, the transition energies for 2p–1s and 2s–2p in the muonic atoms were calculated using Mudirac26, and used as an input for the cascade code. Mudirac typically reproduces the experimental 2p–1s transition energies with a reliability of 0.1%. Therefore, the extracted transition energies are deemed sufficiently precise for calculations on the x-ray intensity.
The cascade code also requires information about the electronic filling fraction of the K, L, and M shells, along with the principal quantum number n, which indicates the starting point of the muon cascade. Another important input is the distribution of muons over the orbital angular momentum states. The initial distribution of muons in the atomic levels after their capture via transfer, denoted as
and
, deviates from the typically assumed statistical distribution after the muon’s direct atomic capture where
and
. The
and
after transfer are calculated based on the formulas presented in Ref.27.
The aforementioned parameters are believed to be well under control. The main unknown in the cascade simulation is the refilling rate of the electronic K-shell after Auger emissions. Following the findings in28,29, these widths are estimated as
. The significance of the choice of the refilling parameter in the cascade simulation was studied by halving or doubling its value or assuming no refilling of the K-shell at all. The findings show a relatively small change on the 2p–1s yield (
relative deviation at most). The computed results for the
–
and
–
transition intensities are presented in Table 2. The summed 2p–1s intensity of 39K (42.74%) corresponds well to the literature value of 40Ar (43.4(7)%)30, indicating that the calculations provide reasonable results.
Simulations in earlier works have shown that the maximum
atomic capture rate occurs for a deuterium concentration of around 0.25%. In this case, the distribution of events in time extends up to
after the muon enters the gas cell4. Applying a time cut with respect to the incoming muon, as well as some scintillator veto cuts, is critical to retrieve a good spectrum. Figure 5 shows a comparison between the total energy spectrum in the germanium detector, and the spectrum after cuts.
Comparison of data before and after applying cuts. After cuts, the 2p–1s signal of interest (
,
, and
) becomes substantially more prominent, while the continuous in time
–
of lead (at
) is suppressed.
With regular fixed-threshold timing on germanium detectors, it becomes difficult to suppress background from uncorrelated and longer timescale events. In order to deal with this, extrapolated leading edge timing (ELET), as described in Ref.30, was employed, which allows for an improvement of the timing resolution by an order of magnitude. This method uses the leading edge trace of the waveform, to calculate a more accurate event time. The algorithm introduces two parameters, the threshold Th and the factor f. The time of the pulse is calculated using the crossings of the pulse waveform at Th and
above its baseline. The timestamps of these crossings are denoted as
and
and the time difference is defined as
. In this work, two ELET algorithms were used, which will be referred to as method A and method B. For these methods, the timing of a given event is calculated using Eqs. (1) and (2), respectively.
The former algorithm assumes that the time difference between the calculated ELET time and the lower threshold is equal to that of the lower threshold and the upper threshold. The latter algorithm performs an exact extrapolation to the baseline. Most detectors performed better under method A, which can be understood as a partial cancellation of non-linear behaviour between the two thresholds. Some detectors, notably BEGe detectors, performed substantially better using the exact extrapolation from method B in combination with lower thresholds. We believe this originates from the broader range of rise times, characteristic of the slow charge collection in these detectors. A graphical representation of both methods is shown in Fig. 6.
Graphical representation of the two ELET algorithms applied in this work. The red dotted lines are the preset thresholds, which can be extrapolated using the cyan dashed lines to the predicted timestamp (green circle).
To test the timing resolution for each set of parameters, the time difference spectrum of the germanium detector of interest is taken with respect to the muon entrance detector. Prompt muonic x rays are produced due to direct capture in parts of the setup, e.g., the titanium support grid of the gas cell. Hence, an energy cut can be taken around such a prompt line, which in this case was taken to be the
2p–1s transition of titanium. After the optimization, the detectors were time aligned and the spectra were added, showing a FWHM of
.
It is well known that high-purity germanium detectors show a changing energy calibration over time due to fluctuations in temperature, humidity, and sudden changes to the setup such as vibrations. To counter these effects, calibration data was taken continuously during the measurement by placing sources in the setup. As conventional calibration sources do not reach up to the 2p–1s muonic x-ray lines of gold, isotopically pure 208Pb was placed in front of the muon veto detector. This way, lead muonic x rays of known energy are emitted (described in Ref.31), such that they can be used for calibration purposes. To complement these, lines originating from some conventional calibration sources were used as well (22Na, 60Co, and 137Cs). The obtained data was used to recalibrate the detectors offline approximately every 8 hours, and whenever a new sample was placed in the setup.
The efficiency response function was investigated using similar techniques to those described in Ref.4,6,30. A
μ
-thick lead target was placed in the target position and was subsequently irradiated with muons. Assuming no scattering, all muons are captured such that the x-ray intensities are in principle known. Accordingly, one can calculate how many photons are emitted in each of the muonic x-ray transitions. Hence, the efficiency of the detector array can be evaluated. Due to previously observed discrepancies in the measured lead x-ray intensities, the literature x-ray intensities for gold were used32, which in previous work showed more consistent results4. Similar to the energy calibration, these were complemented with efficiency measurements using conventional sources (22Na, 88Y, 133Ba, and 152Eu). The main difference, however, was that these could not be placed precisely in the target position. Hence, the efficiencies of these sources had to be rescaled to be in agreement with those calculated from the muonic lead lines. The statistical uncertainty on the efficiency fit indicated a statistical uncertainty on the order of a percent. A relative uncertainty of 10% was assumed for the efficiency, deemed conservative compared to systematic uncertainties originating from x-ray intensities and calibration source activities. The detection efficiencies in the regions of interest are given in Table 3.
Besides these photons detection efficiencies, the muon transmission efficiency through the entrance window was measured to be 55.1(11)%. This value is in good agreement with those observed in previous works4.
To quantify the performance of the different samples, the atomic muon capture rate
, and atomic capture rate per unit mass
are calculated according to
where
is the number of photons of a certain transition that were detected,
is the number of incoming muons,
is the photon detection efficiency,
is the window transmission,
is the relative yield (intensity) of the x ray of interest, and m the mass of the investigated target. The primary reason for quoting
is to normalize for different target masses when making comparisons between different targets. The studied lines were those originating from the 2p–1s transitions.
The range of
atoms can be probed through the attenuation of the x-ray yield. The resulting measurements for this attenuation as a function of graphite layer thickness are shown in Fig. 7 and given in Table 4. For the measurements on gold, only the
–
transition was considered, as the
–
displays an only partially resolved hyperfine structure, and has significant overlap with peaks originating from the continuous 208Pb calibration. The resulting trend, shown in Fig. 8, is near-exponential with a half-value thickness of about
This is in reasonable agreement with estimates made using the value obtained from the measured transfer rate for carbon of
at liquid hydrogen number density33,34 after adjusting for the density of graphite (transfer rate of
).
Spectra obtained for gold samples with different graphite thicknesses. The 2p–1s signal of gold (
,
, and
) reduces with increasing graphite thickness (
), while the
–
of lead (at
) remains constant.
Attenuation of the atomic capture rate as a function of graphite thickness. The error bars represent the statistical error and the grey band represents the systematic error on the measurements, which is the same relative factor for all data points (
).
Using similar calculations, the gold and potassium implanted targets were compared to surface deposition targets. For potassium and chromium, the fine structure could not be fully resolved. Hence, the 2p–1s peak was fitted with a double Gaussian. The results are shown in Table 5. As the outer parts of the target disks are covered by the sample holder, quoted target masses only include the inner
diameter. As mentioned before, only about 60% of the particles are located near the surface of the substrate for the implanted 39K target. Accordingly, both the total mass and the part near the surface are quoted, marked by an asterisk (*) and a dagger (
) in Table 5, respectively.
In order to compare different samples, the atomic capture rate is used. Targets with a higher mass typically have a higher
, as more particles are present, but this does not provide a good indication to the quality of a sample. In contrast, if masses are similar,
is a good comparison when trying to determine which target preparation method is better. Another important note is that the atomic capture rate depends on the spacial extent of the Coulomb potential. Hence, the ratio of atomic capture rates between medium-Z and high-Z targets are listed separately. If we thus compare the implanted targets of similar atomic number and mass, we see that the difference in capture rate is at most
. The implanted targets performed equally well as surface targets, showing there are only limited losses to the substrate material. This result indicates that implantation offers a good alternative target fabrication strategy for muonic atom experiments with microgram targets.
Finally, a comparison was made between barium targets produced via MP and DoD printing. For barium, the fine structure is fully resolved, such that its components can be fitted individually. The samples produced through MP showed a significantly larger signal for the same target mass than those prepared through DoD printing. Even the
μ
DoD printed target did not show a clear signal. By fixing the energy positions and relative intensities of the two fine structure peaks, it was possible to extract a non-zero
for the
μ
MP sample. For the other two DoD printed samples, only an upper limit could be extracted. The results are given in Table 6 and a comparison of the spectra is shown in Fig. 9. The lower atomic capture rates in the DoD printed samples are suspected to be caused by organic surface contamination originating from the solutions used in the target making process/contact with air, or by clustering of barium on the substrate, thus greatly reducing the active surface.
It is rather difficult to make a quantitative comparison between the barium targets and the implanted gold and potassium targets, as transfer rates depend on the atomic number and the needed scattering cross sections for
atoms in these materials are not well known. Given that barium has atomic number
, one would expect that its atomic capture rate per unit mass would be substantially closer to that of gold (
) than that of potassium (
). However, even for the best
(
μ
MP sample), it is only comparable to that of the implanted potassium and thus about half of the one for gold. This is not unexpected as the barium is deposited in molecular form and some organic contamination is codeposited during the process thus leading to additional muon losses. However, these results also show that good x-ray yields can be achieved from such chemical target preparation methods (given a uniform deposition and low amount of organic contamination) when implantation or surface deposition is not an option.
Comparison of spectra obtained by measuring different barium targets showing the
–
(
) and
–
(
) transitions. The molecular plated (MP) targets had a substantially higher signal than the drop-on-demand (DoD) targets.
A systematic study was performed to investigate the most appropriate target preparation techniques for microgram muonic atom spectroscopy. The range of
atoms in carbon was probed by measuring the x-ray yield with different thicknesses of graphite coating on gold. The resulting trend shows a near-exponential behavior with a half-value thickness of approximately
, which is in reasonably good agreement with first order estimates based on the transfer rate of carbon. This shows the importance of minimizing organic surface layers originating from the solutions used for chemical preparation or contact with air, which may reduce the quality of reactive species, such as radium. The range of
atoms also indicates the feasibility of using implanted targets, which was in turn tested successfully by measurement of implanted gold and potassium targets. An additional benefit of implantation is the enhanced protection of chemically reactive species. However, an important factor remains the effect of self-sputtering for larger implanted quantities, which needs to be controlled. These sputtering simulations are discussed in more detail in our earlier work17. Finally, our findings show that targets with good x-ray yield can also be prepared by chemical deposition methods, such as molecular plating, when a uniform deposition with low organic contamination can be achieved.
This work has shown the viability of using implanted targets for microgram muonic atom spectroscopy, and paves the way to future measurements using this approach, e.g. the measurement of the 39, 40, 41K triplet, where magnetic mass separation using an ion beam is a must.
The datasets generated for and/or analysed during the current study are available from the corresponding author on reasonable request.
Yang, X., Wang, S., Wilkins, S. & Ruiz, R. G. Laser spectroscopy for the study of exotic nuclei. Prog. Part. Nucl. Phys. 129, 104005 (2023).
Article
CAS
MATH
Google Scholar
Fricke, G. et al. Nuclear ground state charge radii from electromagnetic interactions. At. Data Nucl. Data Tables 60, 177–285 (1995).
Article
ADS
CAS
MATH
Google Scholar
Tsukada, K. et al. First observation of electron scattering from online-produced radioactive target. Phys. Rev. Lett. 131, 092502 (2023).
Article
ADS
CAS
PubMed
MATH
Google Scholar
Adamczak, A. et al. Muonic atom spectroscopy with microgram target material. Eur. Phys. J. A 59, 15 (2023).
Article
ADS
CAS
PubMed
MATH
Google Scholar
Mulhauser, F. et al. Ramsauer–Townsend effect in muonic atom scattering. Phys. Rev. A 73, 034501 (2006).
Article
ADS
MATH
Google Scholar
Vogiatzi, S. M. Studies of muonic
,
, and
for the extraction of nuclear charge radii. Ph.D. thesis, ETH Zurich (2023).
Maugeri, E. et al. Preparation of 7Be targets for nuclear astrophysics research. J. Instrum. 12, P02016 (2017).
Article
MATH
Google Scholar
Parker, W. & Falk, R. Molecular plating: A method for the electrolytic formation of thin inorganic films. Nucl. Instrum. Methods 16, 355–357 (1962).
Vascon, A. et al. Elucidation of constant current density molecular plating. Nucl. Instrum. Methods Phys. Res. Sect. A 696, 180–191 (2012).
Article
ADS
CAS
Google Scholar
Renisch, D. et al. Actinide and lanthanide thin-layer developments using a drop-on-demand printing system. In EPJ Web of Conferences, Vol. 285, 04001 (2023).
Grillenberger, J., Baumgarten, C. & Seidel, M. The high intensity proton accelerator facility. SciPost Phys. Proc. 5, 002 (2021).
Article
CAS
MATH
Google Scholar
Gerchow, L. et al. Germanium array for non-destructive testing (GIANT) setup for muon-induced x-ray emission (MIXE) at the Paul Scherrer Institute. Rev. Sci. Instrum. 94, 045106 (2023).
Article
ADS
CAS
PubMed
MATH
Google Scholar
Warr, N. et al. The miniball spectrometer. Eur. Phys. J. A 49, 1–32 (2013).
Article
MATH
Google Scholar
Koszorús, Á. et al. Charge radii of exotic potassium isotopes challenge nuclear theory and the magic character of N=32. Nat. Phys. 17, 439–443 (2021).
Article
Google Scholar
Wauters, F. & Knecht, A. The muX project. SciPost Phys. Proc. 5, 022 (2021).
Article
CAS
MATH
Google Scholar
Ziegler, J. F., Ziegler, M. D. & Biersack, J. P. SRIM-the stopping and range of ions in matter (2010). Nucl. Instrum. Methods Phys. Res. Sect. B 268, 1818–1823 (2010).
Article
ADS
CAS
MATH
Google Scholar
Heines, M. et al. Muonic x-ray spectroscopy on implanted targets. Nucl. Instrum. Methods Phys. Res. Sect. B 541, 173–175 (2023).
Article
ADS
CAS
Google Scholar
Gamer, L. et al. Simulation and optimization of the implantation of holmium atoms into metallic magnetic microcalorimeters for neutrino mass determination experiments. Nucl. Instrum. Methods Phys. Res. Sect. A 854, 139–148 (2017).
Article
ADS
CAS
MATH
Google Scholar
Möller, W. & Eckstein, W. Tridyn—a TRIM simulation code including dynamic composition changes. Nucl. Instrum. Methods Phys. Res. Sect. B 2, 814–818 (1984).
Article
ADS
MATH
Google Scholar
Akylas, V. & Vogel, P. Muonic atom cascade program. Comput. Phys. Commun. 15, 291–302 (1978).
Article
ADS
CAS
Google Scholar
Bearden, J. A. & Burr, A. F. Reevaluation of x-ray atomic energy levels. Rev. Mod. Phys. 39, 125–142 (1967).
Article
ADS
CAS
MATH
Google Scholar
Cardona, M. & Ley, L. (eds) Photoemission in Solids I: General Principles (Springer, Berlin, Heidelberg, 1978).
MATH
Google Scholar
Fuggle, J. C. & Mårtensson, N. Core-level binding energies in metals. J. Electron Spectrosc. Relat. Phenom. 21, 275–281 (1980).
Article
ADS
CAS
MATH
Google Scholar
Clementi, E. & Raimondi, D. L. Atomic screening constants from SCF functions. J. Chem. Phys. 38, 2686–2689 (1963).
Article
ADS
CAS
MATH
Google Scholar
Clementi, E., Raimondi, D. L. & Reinhardt, W. P. Atomic screening constants from SCF functions. II. Atoms with 37 to 86 electrons. J. Chem. Phys. 47, 1300–1307 (1967).
Article
ADS
CAS
Google Scholar
Sturniolo, S. & Hillier, A. Mudirac: A Dirac equation solver for elemental analysis with muonic X-rays. X-Ray Spectrom. 50, 180–196 (2021).
Article
ADS
CAS
Google Scholar
Haff, P., Rodrigo, E. & Tombrello, T. Muon transfer in gas targets. Ann. Phys. 104, 363–379 (1977).
Article
ADS
CAS
Google Scholar
Bambynek, W. et al. X-ray fluorescence yields, Auger, and Coster-Kronig transition probabilities. Rev. Mod. Phys. 44, 716–813 (1972).
Article
ADS
CAS
MATH
Google Scholar
Vogel, P. Atomic aftereffects and the line shape of muonic X rays. Phys. Rev. A 8, 2292–2297 (1973).
Article
ADS
CAS
MATH
Google Scholar
Skawran, A. Development of a new method to perform muonic atom spectroscopy with microgram targets. Ph.D. thesis, ETH Zurich (2021).
Bergem, P. et al. Nuclear polarization and charge moments of 208Pb from muonic x rays. Phys. Rev. C 37, 2821 (1988).
Article
ADS
CAS
MATH
Google Scholar
Hartmann, F. et al. Measurement of the muonic x-ray cascade in Mg, AI, In, Ho, and Au. Zeitschrift für Physik A Atoms and Nuclei 305, 189–204 (1982).
Article
ADS
CAS
MATH
Google Scholar
Dupays, A. Muon transfer from muonic hydrogen to carbon. Phys. Rev. A 72, 054501 (2005).
Article
ADS
MATH
Google Scholar
Basiladze, S., Ermolov, P. & Oganesyan, K. Measurement of the rate of transfer of a muon from a p
atom to nuclei of other elements. Soviet J. Exp. Theor. Phys. 22, 725 (1966).
ADS
MATH
Google Scholar
Download references
The experiments were performed at the
E1 beam line of PSI. We would like to thank the accelerator and support groups for the excellent conditions. Additional technical support by F. Barchetti, F. Burri, M. Hildebrandt, M. Meier, L. Noorda, and A. Stoykov from PSI is gratefully acknowledged. We thank the PSI laboratories LNQ and LIN – especially A. Weber and H. Kovacikova – for the preparation of the sputtered and evaporated targets. The germanium detector setup is shared with the MIXE project at PSI (https://www.psi.ch/en/smus/muon-induced-x-ray-emission-mixe-project), who have greatly contributed to its construction, providing a fantastic platform for several muonic atom experiments taking place at PSI. The authors thank Ernst Artes and Carl-Cristian Meyer for their help in the production and characterization of the DoD targets.
The authors want acknowledge the following funding institutions. The Swiss National Science Foundation, Sinergia project “Deep
μ
”, Grant: 193691 (MIXE); FWO Vlaanderen, through proposal numbers G0G3121N (NSHAPE), 1167324N (S. Bara) and 11P6V24N (M. Deseyn); KU Leuven BOF under contract number C14/22/104; the European Research Council (ERC) through proposal number 101088504 (NSHAPE); The RADIATE programme under the Grant Agreement 824096 from the EU research and Innovation programme HORIZON 2020 (proposal 22002733-ST and 22002894-ST); The EPSRC under contract number EP/X015491/1; the ETH Research Grant 22-2 ETH-023 (K. von Schoeler); Center of Excellence in High energy Physics and Astrophysics, Suranaree University of Technology.
Ion Beam Centre, University of Surrey, Guildford, UK
L. Antwis & J. Fletcher
Institute for Nuclear and Radiation Physics, KU Leuven, Leuven, Belgium
S. Bara, T. E. Cocolios, M. Deseyn, M. Heines, W. W. M. M. Phyo & S. M. Vogiatzi
Department of Chemistry, Johannes Gutenberg University Mainz, Mainz, Germany
C. Bruhn, Ch. E. Düllmann & D. Renisch
Institute for Particle Physics and Astrophysics, ETH Zürich, Zurich, Switzerland
A. Doinaki, K. Kirch, K. von Schoeler, A. Zendour & S. Zweidler
PSI Center for Neutron and Muon Sciences, Villigen, Switzerland
A. Doinaki, K. Kirch, A. Knecht & A. Zendour
Helmholtz Institute Mainz, Mainz, Germany
Ch. E. Düllmann, T. Kieck & D. Renisch
GSI Helmholtzzentrum für Schwerionenforschung GmbH, Darmstadt, Germany
Ch. E. Düllmann & T. Kieck
Helmholtz-Zentrum Dresden-Rossendorf, Dresden, Germany
R. Heller & U. Kentsch
Laboratory Kastler Brossel, Sorbonne University, CNRS, ENS-PSL Research University, Collège de France, Paris, France
P. Indelicato
PSI Center for Nuclear Engineering and Sciences, Villigen, Switzerland
E. A. Maugeri
RIKEN Nishina Center, Wako, Saitama, Japan
M. Niikura
Institute of Physics, QUANTUM, Johannes Gutenberg University Mainz, Mainz, Germany
A. Ouf & R. Pohl
Quantum Solid-State Physics, KU Leuven, Leuven, Belgium
L. M. C. Pereira & A. Vantomme
PRISMA+ Cluster of Excellence, Johannes Gutenberg University Mainz, Mainz, Germany
R. Pohl & F. Wauters
Suranaree University of Technology, Nakhon Ratchasima, Thailand
N. Ritjoho
Insitute of Nuclear Physics, Johannes Gutenberg University Mainz, Mainz, Germany
F. Wauters
L.A., T.E.C., J.F., M.H., R.H., U.K., T.K., L.P., and A.V. were involved in the making and characterizing of the implanted samples. C.B., C.E.D., and D.R. made the DoD samples. E.A.M. made the MP samples. S.B., T.E.C., M.H., P.I., K.K., A.K., M.N., A.O., R.P., N.R., K.v.S., S.M.V., F.W., and S.Z. contributed to the data collection. T.E.C., M.H., A.K., S.M.V., and K.v.S. were involved in the analysis. T.E.C., M.D., A.D., M.H., P.I., K.K., A.K., E.A.M., W.W.M.M.P., R.P., K.v.S., S.M.V., F.W., and A.Z. were involved in discussing the results. M.H. wrote the manuscript, and all co-authors reviewed it.
Correspondence to M. Heines.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Reprints and permissions
Antwis, L., Bara, S., Bruhn, C. et al. A comparative study of target fabrication strategies for microgram muonic atom spectroscopy. Sci Rep 15, 6939 (2025). https://doi.org/10.1038/s41598-025-90958-5
Download citation
Received
09 August 2024
Accepted
17 February 2025
Published
26 February 2025
DOI
https://doi.org/10.1038/s41598-025-90958-5
Anyone you share the following link with will be able to read this content:
Provided by the Springer Nature SharedIt content-sharing initiative
Advertisement
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
© 2025 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.
