Trainee question answered by Gemini Pro.
Not bad, but I had to drive it into the depths.
Photons per mu
Let us determine the number of photons incident on a water surface from a 6 MV clinical beam to deliver 100 cGy (1 Gy) to \(d_{max}\) (the depth of maximum dose).
To get a concrete number, we can use the fundamental relationship between dose, fluence, and energy. A 6 MV beam is polyenergetic, so we must make a few standard approximations.
- Define the Knowns:
Target Dose (\(D\)): 100 cGy = 1 Gray = \(1 \text{ J/kg}\)
Depth of max dose (\(d_{max}\)): For a 6 MV beam, \(d_{max} \approx 1.5 \text{ cm}\).
Average Energy (\(\bar{E}\)): A 6 MV bremsstrahlung spectrum has a maximum energy of 6 MeV, but the average photon energy is roughly one-third of the maximum: \(\approx 2 \text{ MeV}\). \(2 \text{ MeV} = 2 \times 10^6 \text{ eV} \times (1.602 \times 10^{-19} \text{ J/eV}) \approx 3.20 \times 10^{-13} \text{ J/photon}\)
- Calculate Fluence at \(d_{max}\): The dose deposited at a specific point is related to the photon fluence (\(\Phi\), photons per unit area) by the mass energy-absorption coefficient (\(\mu_{en}/\rho\)). \[D = \Phi_{dmax} \cdot \bar{E} \cdot \left(\frac{\mu_{en}}{\rho}\right)\]
For water at 2 MeV, the mass energy-absorption coefficient is \(\approx 0.026 \text{ cm}^2/\text{g}\), which is \(0.0026 \text{ m}^2/\text{kg}\).
Rearranging for Fluence (\(\Phi_{dmax}\)): \[\Phi_{dmax} = \frac{D}{\bar{E} \cdot (\mu_{en}/\rho)}\] \[\Phi_{dmax} = \frac{1 \text{ J/kg}}{(3.20 \times 10^{-13} \text{ J/photon}) \cdot (0.0026 \text{ m}^2/\text{kg})}\] \[\Phi_{dmax} \approx 1.20 \times 10^{15} \text{ photons/m}^2 = 1.20 \times 10^{11} \text{ photons/cm}^2\]
- Relate \(d_{max}\) Fluence back to Surface Fluence: The fluence we just calculated is at a depth of 1.5 cm. To find the fluence incident on the surface, we must account for the attenuation of the primary beam as it travels through that 1.5 cm of water. \[\Phi_{surface} = \Phi_{dmax} \cdot e^{+\mu \cdot d}\]
For water at 2 MeV, the linear attenuation coefficient (\(\mu\)) is \(\approx 0.049 \text{ cm}^{-1}\). \[\Phi_{surface} = (1.20 \times 10^{11}) \cdot e^{(0.049 \cdot 1.5)}\] \[\Phi_{surface} = (1.20 \times 10^{11}) \cdot e^{0.0735} \approx (1.20 \times 10^{11}) \cdot 1.076\] \[\Phi_{surface} \approx 1.29 \times 10^{11} \text{ photons/cm}^2\]
- Determine Total Number of Incident Photons: To get the total absolute number of photons, simply multiply the surface fluence by the field size area on the water surface. For a standard \(10 \text{ cm} \times 10 \text{ cm}\) reference field (Area = \(100 \text{ cm}^2\)): \[\text{Total Photons} = 1.29 \times 10^{11} \text{ photons/cm}^2 \times 100 \text{ cm}^2 \approx 1.29 \times 10^{13} \text{ total photons}\]
The monitor chamber is located upstream (above) the secondary collimating jaws.
Because the monitor chamber sits above the X and Y collimating jaws, it is entirely unaware that you have set a 10x10 cm field! (there is some backscatter though)
The primary collimator (which sits above the chamber) restricts the beam to its maximum possible size. For a standard Linac this corresponds to roughly a \(40 \text{ cm} \times 40 \text{ cm}\) square projected at isocenter (or a circle tightly encompassing it).
The flattening filter (also above the chamber) ensures the fluence is roughly uniform across this entire massive cone. Therefore, to find the total number of photons traversing the gas volume of the monitor chamber, you must project the fluence across the maximum field size:
Maximum Field Area (at isocenter): \(40 \text{ cm} \times 40 \text{ cm} = 1600 \text{ cm}^2\).
Total Photons inside Chamber: \((1.29 \times 10^{11} \text{ photons/cm}^2) \times 1600 \text{ cm}^2 = \mathbf{2.06 \times 10^{14} \text{ photons}}\).
Dose Rate
To find the dose rate during the ultra-short pulse (the instantaneous dose rate), you can break the calculation down into a few simple steps.
Medical linear accelerators don’t produce a continuous beam; they fire in extremely brief, intense bursts. Here is the math:
- Calculate the Average (Mean) Dose Rate
Total Dose: \(100 \text{ cGy}\)
Total Time: \(10 \text{ seconds}\)
Average Dose Rate: \(100 / 10 = \mathbf{10 \text{ cGy/sec}}\) (which is \(600 \text{ cGy/min}\), a standard clinical delivery rate).
- Calculate the Dose Deposited per Pulse
Total Time: \(10 \text{ seconds}\)
Pulse Rate: \(400 \text{ pulses/second}\)
Total Pulses Fired: \(10 \text{ s} \times 400 \text{ pulses/s} = \mathbf{4,000 \text{ pulses}}\)
Dose Per Pulse: \(100 \text{ cGy} / 4,000 \text{ pulses} = \mathbf{0.025 \text{ cGy/pulse}}\)
- Calculate the Instantaneous Dose Rate (During the Pulse)
Now that we know exactly how much dose is squeezed into a single pulse, we divide it by the length of time that pulse actually lasts.
Pulse Width: \(3 \mu\text{sec} = 3 \times 10^{-6} \text{ seconds}\)
Dose Rate During Pulse: \(0.025 \text{ cGy} / (3 \times 10^{-6} \text{ seconds})\)
Result: \(\mathbf{\approx 8,333 \text{ cGy/sec}}\)
But this isn’t quite right!
Within that \(3 \mu\text{s}\) window, the linac doesn’t produce a steady stream of electrons; it produces train of electron “bunches” separated by the frequency of the microwaves pushing them. This is the Micropulse structure.
Here is exactly how you apply the 3 GHz (\(3 \times 10^9\) Hz) microwave frequency to find the true peak instantaneous dose rate:
- Calculate the Number of Bunches (Micropulses)
Because the RF wave oscillates at 3 GHz, it completes three billion cycles every second. During one \(3 \mu\text{s}\) macropulse, the number of RF cycles is: \(3 \times 10^9 \text{ cycles/second} \times (3 \times 10^{-6} \text{ seconds}) = \mathbf{9,000 \text{ cycles}}\)
Since electrons are bunched and pushed by every wave cycle, there are exactly 9,000 individual electron bunches (micropulses) squeezed into a single \(3 \mu\text{s}\) macropulse.
- Calculate the Dose per Micropulse
From our previous math, one \(3 \mu\text{s}\) macropulse deposits \(0.025 \text{ cGy}\). If that dose is divided evenly among the 9,000 bunches:
\(0.025 \text{ cGy} / 9,000 = \mathbf{0.00000277 \text{ cGy per bunch (micropulse)}}\)
- Calculate the True “Beam On” Time (Bunch Length)
To find the ultra-peak instantaneous dose rate, we need to know the exact duration of a single micropulse.
The RF Period: One full 3 GHz cycle takes \(1 / (3 \times 10^9 \text{ Hz}) = \mathbf{333.3 \text{ picoseconds}}\).
The Phase Angle (Bunching): Electrons are not accelerated during the entire 333 ps wave. They only ride the “crest” of the positive half-cycle. The width of this electron bunch depends on the linac’s electron gun and waveguide tuning, but a typical phase acceptance angle (\(\Delta \phi\)) is around 15° to 30° out of the full \(360^{\circ}\) RF cycle. Let’s assume a typical bunch width of roughly 20° of the RF cycle:
Fraction of the cycle = $ $ of the cycle.
Actual Bunch Duration: \(333.3 \text{ ps} \times 0.055 = \mathbf{18.5 \text{ picoseconds}}\). So, those electrons are only hitting the patient for 18.5 picoseconds at a time!
- The True Peak Instantaneous Dose Rate Finally, to get the dose rate during that 18.5 picosecond impact, we divide the dose per bunch by the bunch duration:
\(0.00000277 \text{ cGy} / (18.5 \times 10^{-12} \text{ seconds}) \approx \mathbf{150,000 \text{ cGy/sec} \text{!!!}}\)
Photons per sec
Here is the “Russian Doll” of dose rates reimagined completely in terms of total photons per second hitting the patient (across the full 10x10 cm field).
We will use our original base number: It takes roughly \(1.29 \times 10^{13}\) total photons (12.9 trillion) to deliver 100 cGy, and the delivery takes 10 seconds.
- Clinical Mean Rate (The whole 10-second treatment)
This is the simple average over the whole treatment time.
Total Photons: \(1.29 \times 10^{13}\)
Total Time: \(10\) seconds
Mean Photon Rate: \(1.29 \times 10^{13} / 10 = \mathbf{1.29 \times 10^{12} \text{ photons/second}}\) (1.29 Trillion photons per second)
- Macropulse Rate (During the \(3\mu\text{s}\) burst)
The treatment occurs in \(3\mu\text{s}\) bursts at 400 Hz. Over 10 seconds, there are \(4,000\) bursts.
Photons per Burst (Macropulse): \(1.29 \times 10^{13} / 4,000 \text{ pulses} = \mathbf{3.23 \times 10^9 \text{ photons/pulse}}\)
Rate During the \(3\mu\text{s}\) Burst: \(3.23 \times 10^9 \text{ photons} / (3 \times 10^{-6} \text{ seconds})\)
Macropulse Photon Rate: \(\mathbf{1.07 \times 10^{15} \text{ photons/second}}\) (1.07 Quadrillion photons per second)
- Micropulse Rate (During the 18.5 ps 3 GHz bunch impact)
Inside each \(3\mu\text{s}\) burst, the 3 GHz microwave frequency creates 9,000 individual electron bunches (producing 9,000 photon bursts), with each burst lasting just 18.5 picoseconds.
Photons per Bunch (Micropulse): \(3.23 \times 10^9 \text{ photons} / 9,000 \text{ bunches} \approx \mathbf{358,000 \text{ photons/bunch}}\)
Rate During the 18.5 ps Bunch: \(358,000 \text{ photons} / (18.5 \times 10^{-12} \text{ seconds})\)
Micropulse Photon Rate: \(\mathbf{1.93 \times 10^{16} \text{ photons/second}}\)