Tools / Beam and test / Proton dose and flux calculator

Proton beam flux and dose rate

Beam current → flux (p/cm²/s) → ionizing dose rate (rad/s, krad/min) in a silicon target. Enter the beam energy and current, or the flux you want, and the beam area, and read the time to your target dose and the total dose for a fluence.

Silicon targetBethe–Bloch LET with optional facility-table overrideUC Davis CNL presetruns in your browserHelp & Reference

Beam

Input mode

Beam area

LET data source

Bethe–Bloch with density correction is accurate above ~20 MeV. Below that, paste your facility / NIST PSTAR data here.
Paste LET vs energy table Whitespace, comma, tab, or semicolon separated. First column = MeV, second column = MeV·cm²/g (electronic stopping power; PSTAR's total column differs by under 0.1 percent above 1 MeV). Log-log interpolation between points.

Time to target dose

Fluence → total dose

Results

From the beam to the dose: the current over the beam area is the flux, the flux times the LET in silicon is the dose rate, and the dose rate sets the time to your target dose and the dose a fluence delivers.
Particle Flux
,
p·cm⁻²·s⁻¹
Dose rate
,
krad/min
Required time to the target dose
,
Total dose for the fluence
,
krad
Beam current (derived)
,
p/s
Dose rate
,
rad/s
Dose rate (SI)
,
Gy/s
LET in Si
,
MeV·cm²/g
LET (linear)
,
keV/µm
Source mode
,
Total dose
,
rad
Total dose (SI)
,
Gy
Time at current flux
,
to deliver the fluence

LET vs energy

Solid line: active LET source. Dot marks the current energy.

Quick reference table (active beam & LET source)

E (MeV)LET (MeV·cm²/g)LET (keV/µm) Flux for 1 krad/min
(p·cm⁻²·s⁻¹)
Beam current for 1 krad/min
(across active area)

How it works

Converting beam current to dose rate

01

Current to flux

The beam current divided by the proton charge is the number of protons per second, and spread over the beam area it is the flux. In flux mode the same relation runs backward to the current the beam needs.

02

LET in silicon

The mass stopping power at the beam energy comes from Bethe–Bloch with the Sternheimer density correction, or from the facility or PSTAR table you paste, interpolated log-log.

03

Flux times LET is the dose rate

Flux times LET, times 1.602×10-8 rad per MeV/g, is the dose rate in silicon. It sets the time to your target dose, and a fluence times the LET gives the total dose.

How the numbers are computed

Flux:  φ = (I / e) / A with e = 1.602×10⁻¹⁹ C.

Dose rate (rad/s):  Ḋ = φ · LET · 1.602×10⁻⁸ with LET in MeV·cm²/g (1 MeV/g = 1.602×10⁻⁸ rad).

Linear LET:  (MeV/cm) = LET·ρ with ρSi = 2.33 g/cm³.

Bethe–Bloch (Si):  PDG form with Sternheimer density correction (Si: I=173 eV, x₀=0.2014, x₁=2.8715, a=0.1492, m=3.255, C̄=4.4355).

Disclaimer

This tool is provided as a free community resource. Results should be verified independently. See Help & Reference