Absorbed and effective dose estimates in HR-pQCT of the distal radius and tibia: virtual dosimetry with a GPU-accelerated Monte Carlo code.

Abstract

Objective. To obtain maps of absorbed dose and an estimate of the effective dose to an adult patient for a high resolution peripheral quantitative computed tomography (HR-pQCT) (XtremeCT II) examination of the distal tibia and radius, using a graphical processing unit (GPU)-based Monte Carlo (MC) code.Approach. We adapted the validated code gCTD (GPU-based CT Dose calculator), to replicate the HR-pQCT configuration. MC simulations were performed on digital phantoms of the tibia and radius obtained from bilateral scans of the ankle and wrist of a 25 year-old female volunteer. Scans were segmented using an ad hoc algorithm in Fiji. Simulations run on an NVIDIA GeForce RTX 3090 GPU board. MC dose estimates were validated via computed tomography dose index measurements.Main Results.We obtained the absorbed dose distribution in the skin, bone, bone marrow, fat, and muscle tissues. The effective dose for the HR-pQCT examination were 2.05μSv and 2.13μSv for the right and left tibia, and 1.48μSv and 1.49μSv for the right and left radius, respectively, with a Type A statistical uncertainty of 0.06% (k= 3) with 4.66 × 10¹¹photon histories. Corresponding effective dose conversion coefficients (k-factors) were 0.185μSv ∙ mGy^-1· cm^-1(tibia), and 0.133μSv mGy^-1· cm^-1(radius).Significance.We reported the first independent estimate of the effective dose for standard HR-pQCT clinical scans of the distal tibia and radius with the XtremeCT II scanner. Effective dose estimates (considering a total relative uncertainty of less than 40%) were lower than those indicated by the manufacturer and commonly reported for these scans. With 4.66 × 10⁹photon histories, the gCTD MC code can produce 3D dose maps from segmented HR-pQCT images in less than 12 s (GPU time), with 0.9% (k= 3) statistical uncertainty, making real-time personalized dose estimate feasible.