The Essential Role of Monte Carlo Simulations for Lung Dosimetry in Liver Radioembolization with 90Y Microspheres

This study compares various methodologies for lung dosimetry in radioembolization using Monte Carlo (MC) simulations. A voxelized anthropomorphic phantom, created from a real patient’s CT scan, preserved the actual density distribution of the lungs. Lung dosimetry was evaluated for five lung-shunt (LS) cases using traditional methods: the mono-compartmental organ-level approach (MIRD), local energy deposition (LED), and convolution with voxel S-values, either with local density corrections (SVOX_L) or without (SVOX_ST). Additionally, a novel voxel S-value (VSV) kernel for lung tissue with an ICRU density of 0.296 g/cm3 was developed. Calculations were performed using either the ICRU lung density (Lung_296), the average lung density of the phantom (Lung_221), or the local density (Lung_L). The comparison revealed significant underestimations in the mean absorbed dose (AD) for the classical approaches: approximately −40% for MIRD, −27% for LED, −28% for SVOX_L, and −88% for SVOX_ST. Similarly, calculations with the lung VSV kernel showed underestimations of about −62% for Lung_296, −50% for Lung_221, and −35% for Lung_L. Given the high heterogeneity of lung tissue, traditional dosimetric methods fail to provide accurate estimates of the mean AD for the lungs. Therefore, MC dosimetry based on patient images is recommended as the preferred method for precise assessment of lung AD during radioembolization.