Joy-Luo Stopping Power of Kapton

% Define parameters
E = logspace(1, 8, 1000); % ion energy (eV)
v = sqrt(2*E*1.602e-19/1.5e-25); % ion velocity (m/s)
N = 9.09e22; % number density of Kapton atoms (atoms/cm^3)
z = 1; % ion charge (proton)
m = 1.67e-27; % ion mass (kg)
I = 56; % mean excitation energy of Kapton atoms (eV)
gamma = 0.149;
delta = 0.0014;
Delta = 0.00039;

% Compute stopping power
S = 4 * pi * N * (1.602e-19)^2 * z^2 ./ (m*v.^2) .* ...
    (log(2*m*v.^2/I) + gamma - delta - Delta);

% Convert stopping power to MeV/(mg/cm^2)
rho = 1.42; % Kapton density (g/cm^3)
S_mev = S * rho * 10^-3;

% Plot stopping power vs. energy
loglog(E, S_mev);
xlabel('Ion Energy (eV)');
ylabel('Stopping Power (MeV/(mg/cm^2))');
title('Joy-Luo Stopping Power of Kapton');
grid on;