We have developed a novel 4D dynamic liver blood flow model, capable of accurate dose estimation to circulating blood cells during liver-directed external beam radiotherapy, accounting for blood recirculation and radiation delivery time structure. Adult male and adult female liver computational phantoms with detailed vascular trees were developed to include the hepatic arterial, hepatic portal venous, and hepatic venous trees. A discrete time Markov Chain approach was applied to determine the spatiotemporal distribution of 105 blood particles (BP) in the human body based on reference values for cardiac output and organ blood volumes. For BPs entering the liver, an explicit Monte Carlo simulation was implemented to track their propagation along ~2000 distinct vascular pathways through the liver. The model tracks accumulated absorbed dose from time-dependent radiation fields with a 0.1s time resolution. The computational model was then evaluated for 3 male and 3 female patients receiving photon (VMAT and IMRT) and proton (passive SOBP and active PBS) treatments. The dosimetric impact of treatment modality, delivery time, and fractionation on circulating blood cells was investigated and quantified using the mean dose (μdose,b), V>0Gy, V>0.125Gy, and D2%. Average reductions in μdose,b, V>0Gy, V>0.125Gy, and D2% of 45%, 6%, 53%, 19% respectively, were observed for proton treatments as compared to photon treatments. Our simulation also showed that V>0Gy, V>0.125Gy, and D2% were highly sensitive to the beam-on time. Both V>0Gy and V>0.125Gy increased with beam-on time, whereas D2% decreased with increasing beam-on time, demonstrating the tradeoff between low dose to a large fraction of blood cells and high dose to a small fraction of blood cells. Consequently, proton treatments are not necessarily advantageous in terms of dose to the blood simply based on integral dose considerations. Instead, both integral dose and beam-on time can substantially impact relevant dosimetric indices.

Purpose To develop a model of the internal vasculature of the adult liver and demonstrate its application to the differentiation of radiopharmaceutical decay sites within liver parenchyma from those within organ blood. Method Computer-generated models of hepatic arterial (HA), hepatic venous (HV), and hepatic portal venous (HPV) vascular trees were algorithmically created within individual lobes of the ICRP adult female and male livers (AFL/AML). For each iteration of the algorithm, pressure, blood flow, and vessel radii within each tree were updated as each new vessel was created and connected to a viable bifurcation site. The vascular networks created inside the AFL/AML were then tetrahedralized for coupling to the PHITS radiation transport code. Specific absorbed fractions (SAF) were computed for monoenergetic alpha particles, electrons, positrons, and photons. Dual-region liver models of the AFL/AML were proposed, and particle-specific SAF values were computed assuming radionuclide decays in blood within two locations: (1) sites within explicitly modeled hepatic vessels, and (2) sites within the hepatic blood pool residing outside these vessels to include the capillaries and blood sinuses. S values for 22 and 10 radionuclides commonly used in radiopharmaceutical therapy and imaging, respectively, were computed using the dual-region liver models and compared to those obtained in the existing single-region liver model. Results Liver models with virtual vasculatures of ~ 6000 non-intersecting straight cylinders representing the HA, HPV, and HV circulations were created for the ICRP reference. For alpha emitters and for beta and auger-electron emitters, S values using the single-region models were approximately 11% (AML) to 14% (AFL) and 11% (AML) to 13% (AFL) higher than the S values obtained using the dual-region models, respectively. Conclusions The methodology employed in this study has shown improvements in organ parenchymal dosimetry through explicit consideration of blood self-dose for alpha particles (all energies) and for electrons at energies below ~ 100 keV.

No abstract

Objective: Phantoms of the International Commission on Radiological Protection provide a framework for standardized dosimetry. The modeling of internal blood vessels – essential to tracking circulating blood cells exposed during external beam radiotherapy and to account for radiopharmaceutical decays while still in blood circulation – is, however, limited to the major inter-organ arteries and veins. Intra-organ blood is accounted for only through the assignment of a homogeneous mixture of parenchyma and blood [single-region (SR) organs]. Our goal was to develop explicit dual-region (DR) models of intra-organ blood vasculature of the adult male brain (AMB) and adult female brain (AFB). Approach: A total of 4000 vessels were created amongst 26 vascular trees. The AMB and AFB models were then tetrahedralized for coupling to the PHITS radiation transport code. Absorbed fractions were computed for monoenergetic alpha particles, electrons, positrons, and photons for both decay sites within the blood vessels and for tissues outside these vessels. Radionuclide S-values were computed for 22 and 10 radionuclides commonly employed in radiopharmaceutical therapy and nuclear medicine diagnostic imaging, respectively. Main Results: For radionuclide decays, values of S(brain tissue ← brain blood) assessed in the traditional manner (SR) were higher than those computed using our DR models by factors of 1.92, 1.49, and 1.57 for therapeutic alpha-emitters, beta-emitters, and auger electron-emitters, respectively in the AFB and by factors of 1.65, 1.37, and 1.42 for these same radionuclide categories in the AMB. Corresponding ratios of SR and DR values of S(brain tissue ← brain blood) were 1.34 (AFB) and 1.26 (AMB) for four SPECT radionuclides, and were 1.32 (AFB) and 1.24 (AMB) for six common PET radionuclides. Significance: The methodology employed in this study can be explored in other organs of the body for proper accounting of blood self-dose for that fraction of the radiopharmaceutical still in general circulation.

The dosimetric dependence of ocular structures on eye size and shape was investigated within the standard ICRP Publication 116 irradiation geometries. A realistic transport geometry was constructed by inserting a scalable and deformable stylised eye model developed in our previous study within the head of the ICRP Publication 110 adult male reference computational phantom. Beam irradiations of external electrons, photons, and neutrons on this phantom were simulated using the Monte Carlo radiation transport code PHITS in the geometries of AP, RLAT, PA and ROT. Absorbed doses in ocular structures such as ciliary body, retina, and optic nerves were computed as well as that in lens. A clear dosimetric dependence of ocular structures on eye size and shape was observed for external electrons while only a small dependence was seen for external photons and neutrons. Difference of the tendency was attributed to their depth-dose distributions where spread dose distributions were created by photons and neutrons while more concentrated distributions were created by external electrons.