Python code for making plots relevant to the implantation of thorium ions, including density profile plots.

The goal of density_profile_plots_gaussians.ipynb is to plot the density of Th ions (number of ions per cm^3 per beam current in pps) as a function of perpendicular distance from the front surface of a CaF2 crystal, for only one given beam energy. To do this, we define three normalized shape functions that describe the density of Th ions in space after the implantation process. So far, we are assuming that the density rho(x,y,z) can be factorized into the product of 3 independent Gaussians. Calculations in this notebook rely on values derived from LISE++ and SRIM calculations, and useful references (Google Sheet with all calculated values, MJT logbook pages, elog entries such as elog:Pa229/167, ...) are linked in the first cell.

The goal of Th_impl_plots_20251111.ipynb is to make several plots which were relevant to the MSU/UCLA Th implantation meeting on Nov 11, 2025 (see meeting slides, as well as Aiden's notes from this meeting at elog:Pa229/149). All of these plots are based on mean and stdev implantation values calculated in LISE++, with beam energies from 0.0001 MeV/u to 4 MeV/u. Those calculated values are read from a csv file, lise_th_impl_data_20251111.csv. The comparison of different beam energies in the third plot shows that for a given beam energy within the ReA range, the straggling is low and the implantation is very nonuniform -- most ions will end up more or less at the same depth. List of plots created in this Python notebook:

• Average Th ion stopping range vs beam energy per nucleon
• Standard deviation of Th ion stopping vs beam energy per nucleon
• Yield vs implantation depth in CaF2 for several beam energies, ranging from 0.1 MeV/u to 4 MeV/u (uses rate of Th ions from beam as calculated by FRIB beam rate calculator)
• Yield vs implantation depth in CaF2 for a single beam energy of 0.0001 MeV/u (total beam energy 0.023 MeV, lower than ReA could do)