SciPost Submission Page
Muonium reaction in MgO: A showcase for the final steps of ion implantation
by Rui C. Vilão, Ali Roonkiani, Apostolos G. Marinopoulos, Helena V. Alberto, João M. Gil, Ricardo B. L. Vieira, Robert Scheuermann and Alois Weidinger
This is not the latest submitted version.
Submission summary
| Authors (as registered SciPost users): | Apostolos Marinopoulos · Rui Vilão |
| Submission information | |
|---|---|
| Preprint Link: | scipost_202504_00035v2 (pdf) |
| Date submitted: | June 17, 2025, 4:42 p.m. |
| Submitted by: | Rui Vilão |
| Submitted to: | SciPost Physics |
| Ontological classification | |
|---|---|
| Academic field: | Physics |
| Specialties: |
|
| Approaches: | Theoretical, Experimental |
Abstract
We present an in-depth investigation of the implantation of positive muons in magnesium oxide (MgO). Muonium, the positive muon plus an electron is an analogue of the hydrogen atom. This study describes the final stage of the implantation process, from muon diffusion over the potential barrier and the stopping by an inelastic reaction to the final embedding of the muon into the lattice structure. A special aspect is a relatively long-lived intermediate configuration which lasts for several hundred nanoseconds or more and is accessible to muon spin spectroscopy. The model presented here provides a framework for the analysis of the general case of ion implantation.
Author indications on fulfilling journal expectations
- Provide a novel and synergetic link between different research areas.
- Open a new pathway in an existing or a new research direction, with clear potential for multi-pronged follow-up work
- Detail a groundbreaking theoretical/experimental/computational discovery
- Present a breakthrough on a previously-identified and long-standing research stumbling block
Author comments upon resubmission
List of changes
As a reply to the comments of Referee 1:
-
“Doorway model” expresses more clearly than “transition state model” the entrance character of the initial reactions: The passage through the “doorway” by the inelastic reaction and the formation of an initial configuration through which all other reactions proceed. We added a corresponding sentence in the paper.
-
We added in the description of Fig. 4 in the main manuscript: “About 20 % fraction is missing. It is attributed to muon spin polarization loss due to rapid fluctuations of the hyperfine interaction in the initial hot phase after the muon stopping.”
-
We added the following in the manuscript when introducing Fig. 4, before the last paragraph of section 2: “As mentioned in the experimental section, the observed muonium amplitudes show a peaking slightly above 100 K. We note that the increase is parallel to the increase of diamagnetic-like fraction, suggesting that both states are formed in the same process. The decrease of the muonium fraction at temperatures above about 120 K may either be due to a decrease of the formation probability or to dephasing effects which are rather strong at these high frequencies.”
4 – C is a constant fitting parameter. It is related to the coupling strength of the muonium electron with the phonons, but no details can be given here. We added this in the main text, together with the remaining parameters of the equation.
5 – We corrected this and apologize for the wrong citation.
6 – We replaced this reference by the PSI webpage.
Current status:
Reports on this Submission
Report #3 by Anonymous (Referee 4) on 2025-7-19 (Invited Report)
- Cite as: Anonymous, Report on arXiv:scipost_202504_00035v2, delivered 2025-07-19, doi: 10.21468/SciPost.Report.11601
Report
On the basis of the requested experiments (see below) the major revision should be done.
I do not recommend publication of the current version of the manuscript.
Requested changes
The authors must set a solid experimental grounds before suggesting a model. For that: 1. The authors must carry out electric field experiments at various magnetic fields. 2. The authors must measure magnetic field dependences of the muonium amplitude and the initial phase of the muonium precession.
After that is done: 3. The suggested model should be set against different models of Mu formation, in particular, the model of delayed Mu formation via an electron capture by the muon.
Recommendation
Ask for major revision
Report #2 by Jess Brewer (Referee 3) on 2025-7-4 (Invited Report)
- Cite as: Jess Brewer, Report on arXiv:scipost_202504_00035v2, delivered 2025-07-04, doi: 10.21468/SciPost.Report.11515
Strengths
Weaknesses
Report
Requested changes
See attached PDF.
Recommendation
Publish (surpasses expectations and criteria for this Journal; among top 10%)
Referee 3 (Prof. Jess Brewer)
We thank Prof. Jess Brewer for this fair, thoughtful, critical and inspiring review of our paper. We appreciate his comments and questions and answer them as follows
Referee 3
The paper begins with a list of assorted applications of energetic ions (especially protons) in which (the authors declare) the final stages of the stopping process are poorly understood, despite being essential to effective utilization. Since no description is given of what is unknown in the several cases, the reader is left to either read all the references or just take the authors’ word for it.
Our answer
We have rephrased this part as: "Although the stopping process is rather well understood for the high-energy region, very little is known about the final thermalization process [8–12]. The reason for this is that conventional experimental methods are not able to investigate this process at the atomic level. Muon spectroscopy offers the unique possibility that the final steps can be observed directly on the implanted particle and thus detailed information about these final steps can be obtained.”
Referee 3
Footnote 1: For some reason they refer to the 1947 Fermi and Teller paper concerning µ− capture.
Our answer
We wanted to express that stopping of muons was already on the agenda at the beginning of muon research. But this mentioning is not necessary here and we have skipped this reference.
Referee 3
Footnote 2: At the bottom of p. 11 they say, “The inelastic reaction is caused by the force exerted by the squeezed muonium on the surrounding atoms.” If the Mu0 is “squeezed” rather than “expanded”, then why is the hyperfine interaction weaker than in vacuum?
Our answer
We skipped the word “squeezed” here and the word “compressed” in a similar context in the conclusion section.
Referee 3
I have some problems with this picture.
First, the “doorway state” is proposed to last for “hundreds of nanoseconds” before losing its electron completely, during which time it has a slightly shifted precession frequency due to the loosely associated, thermally-polarized electron, after which it precesses at the free muon Larmor frequency ωD. If so, then a Fourier transform omitting the first few hundred ns should yield a sharp line at ωD. Was this done? If not, it should be. If so, and it yielded the same broadened and shifted line as the full time spectrum, then this picture cannot be correct.
Our answer
Analyzing the time spectra separately in the earlier and later times regime gives similar information as a two-component fit with a fast and a slow relaxation. We have tried such fits and obtained of course a fraction with slow relaxation. But the data were not good enough to allow a reliable extraction of the lifetime. We discuss that in section IV.2, third paragraph.
Referee 3
Second, the energy scales seem inconsistent. A fully-bound, compact Mu0 state is bound by some substantial fraction of a Rydberg unless its electron is already shared with the lattice. (Weakly bound Mu0 states are frequently observed in solids, but they always involve the lattice intimately.) Therefore the picture of a epithermal weakly-bound Mu0 state diffusing rapidly through the lattice is self-contradictory. Such diffusion is only possible for a compact, strongly-bound Mu0 state, which therefore arrives at the crucial inelastic collision with a significant fraction of a Rydberg of energy. This energy is supposedly released in said inelastic process, contributing to the local “thermal spike” that raises the effective temperature of the local lattice (which must dramatically reduce the thermal equilibrium polarization of the weakly-bound electron in the “doorway state”). How is it then possible that the resultant frequency shift is so strongly temperature dependent?
Our answer
The energetics discussion is based on the DFT calculations and the data shown in Fig. 6. The diffusing muonium is a compact atom similar to ground state muonium at the interstitial site. At the top of the barrier, the “forward directed” kinetic energy is only a few meV or less and is fully transferred to the oscillating atoms after the inelastic reaction. The inelastic reaction and the subsequent relaxation of the lattice lead then to the long-lived doorway state. The total energy transfer to heat is about 1.2 eV (difference of the barrier heights in Fig. 6).
Referee 3
It may be that my difficulty in following the logic of this paper is due to the somewhat artificial separation of experimental observations from the theory that explains their interpretation; this is probably structurally mandatory, but I get confused when flipping back and forth between assertions and their separate justifications.
Our answer
We had the same difficulty when writing the manuscript, but found that the current structure, showing very clearly the experimental data first and then building the model, is the one which gives the reader the most clear separation of experimental data and theoretical model.
Report #1 by Anonymous (Referee 1) on 2025-6-30 (Invited Report)
- Cite as: Anonymous, Report on arXiv:scipost_202504_00035v2, delivered 2025-06-30, doi: 10.21468/SciPost.Report.11483
Report
The authors addressed all questions well, and I recommend publication after the authors corrected the following problem with the references:
- In line 534, ref. [76] should be replaced by the ref. [58]. There is no relation of ref. [76] (this was ref. [58] in version 1 of the manuscript) to the paper and it can be removed.
Requested changes
- In line 534, ref. [76] should be replaced by the ref. [58]. There is no relation of ref. [76] (this was ref. [58] in version 1 of the manuscript) to the paper and it can be removed.
Recommendation
Ask for minor revision

Author: Rui Vilão on 2025-07-23 [id 5672]
(in reply to Report 3 on 2025-07-19)Referee 4
We do not agree with the assessment of the paper of Referee 4. We think the model presented in the paper has a solid foundation by the theoretical and experimental data presented in the manuscript.
Referee 4
The authors must set a solid experimental grounds before suggesting a model. For that: 1. The authors must carry out electric field experiments at various magnetic fields. 2. The authors must measure magnetic field dependences of the muonium amplitude and the initial phase of the muonium precession.
Our answer
In the present case, electric field measurements are not necessary. Such measurements cannot distinguish between the different models because the behavior of the weakly bound electron is the same in both models. The same applies to magnetic field measurements. The justification for our interpretation is based on the broadened diamagnetic line and the frequency shift.
Referee 4
After that is done: 3. The suggested model should be set against different models of Mu formation, in particular, the model of delayed Mu formation via an electron capture by the muon.
Our answer
In the introduction, we describe the competing models and mention the relevant literature. We also mention the controversy about these models and cite the references about it. In the other sections of the paper, a comparison of the models is implicitly addressed when we talk about the broadened diamagnetic-like line and the frequency shift, which shows that the paramagnetic interaction is present from the beginning.
In our opinion, the model is sufficiently well-founded experimentally and theoretically to be presented to the scientific community for discussion.