obtained, provides a uniquely powerful tool to get insight into the intrinsic properties of the
material. However, to apply hydrostatic pressure (e.g. at GPD, which will be the spectrometer
used as the example throughout this paper), the samle needs to be compressed into a pellet
and inserted into the pressure cell made of either MP45N or Cu-Be. Before starting any
measurements, the muon beam and its momentum need to be adjusted to guarantee the maximal
signal from the sample. With surface muon beamlines (e.g. GPS at PSI), the energy of the muons
is usually 5 - 40 MeV/c. This is enough to penetrate a thin layer of mylar foil or capton tape
which is used to mount the sample to the sample holder. However, this energy is not enough to
penetrate the thick layers of metal that make up a pressure cell. For this, decay muon beamlines
are needed, which can provide muons in the energy range of 40 - 125 MeV/c. For the surface
muon case, alignment is simple as there is essentially only the sample signal to detect. When the
pressure cell is introduced, it gets slightly more difficult as there is now another material in very
close proximity to the sample that will give a significant signal. This is not an issue when the
sample gives a very strong response, as the sample signal will be much greater than the pressure
cell signal. Consequently, determining the optimal muon momentum for alignment is not very
difficult. However, if the sample response is weak, then it may be difficult to distinguish the
sample contributions from the pressure cell (as mentioned in Ref [8]). Additionally, there are
different types of pressure cells (MP35N, CuBe and a combination of the two) with different
signals. Further, each individual pressure cell made from same material are in fact also slightly
different. As a result, there are several cases where it is incredibly difficult to extract the sample
response from the total signal.
One way around this is to first employ an indirect alignment using another sample with a
strong response (but same pressure cell), and then using this alignment result also for the sample
of interest (with weak signal). However, this can be a problem as slight differences in the sample
densities can drastically change the stopping fraction of muons in the sample (a higher density
means less muons pass through the sample). This then translates into issues in the fitting process
as you need to have a correct estimate for the signal fraction from the sample. Consequently, if
your assumed stopping fraction is incorrect, the obtain fitting results do not completely reflect
the intrinsic physical properties and behaviour of the sample.
In this paper, we present an efficient and user friendly simulation method based using TRIM
[9] combined with the pySRIM python module [1]. In such approach, it is possible to simulate,
with higher level of accuracy, the number of muons that will be stopping in the pressure cell
and the sample, respectively, for any given muon momentum input. This is useful as a tool for
user to employ both before and after experiments. This is because it will both give the user
an idea on feasibility of their proposed experiment (in the respect of how easy it will be to see
their sample signal) and also for supporting the data analysis when fitting the percentages of
the contributions from each fit function (pressure cell and sample).
2. Basis of Simulation
2.1. Beam setup
The GUI software is a python package that determines the stopping fractions in a beamline
setup. The software runs TRIM calculations in the background and presents the results as early
interpretable figures and stopping percentages.
There are three key aspects that make this software more accurate for modelling stopping
fractions compared to running the standard TRIM simulations:
(i) The software utilises a Gaussian distribution of momentums/energies.
(ii) The input beam can be collimated to any area.
(iii) There is a random small angular divergence on each muon as each muon will not be
completely parallel to the beam.