Transient simulation of laser ablation based on Monte Carlo light transport with dynamic optical properties model

Laser ablation is a minimally invasive therapeutic technique to denature tumors through coagulation and/or vaporization. Computational simulations of laser ablation can evaluate treatment outcomes quantitatively and provide numerical indices to determine treatment conditions, thus accelerating the technique's clinical application. These simulations involve calculations of light transport, thermal diffusion, and the extent of thermal damage. The optical properties of tissue, which govern light transport through the tissue, vary during heating, and this affects the treatment outcomes. Nevertheless, the optical properties in conventional simulations of coagulation and vaporization remain constant. Here, we propose a laser ablation simulation based on Monte Carlo light transport with a dynamic optical properties (DOP) model. The proposed simulation is validated by performing optical properties measurements and laser irradiation experiments on porcine liver tissue. The DOP model showed the replicability of the changes in tissue optical properties during heating. Furthermore, the proposed simulation estimated coagulation areas that were comparable to experimental results at low-power irradiation settings and provided more than 2.5 times higher accuracy when calculating coagulation and vaporization areas than simulations using static optical properties at high-power irradiation settings. Our results demonstrate the proposed simulation's applicability to coagulation and vaporization region calculations in tissue for retrospectively evaluating the treatment effects of laser ablation.