Closed-form solution for the thermal dose delivered during single pulse thermal therapies

Numerical Model Study of In Vivo Magnetic Nanoparticle Tumor Heating

Iron oxide nanoparticles are currently under investigation as heating agents for hyperthermic treatment of tumors. Major determinants of effective heating include the biodistribution and minimum iron oxide loading required to achieve adequate heating at practically achievable magnetic field strengths. These inter-related criteria ultimately determine the practicality of this approach to tumor treatment. Further, in our experience the currently used treatment assessment criterion for hyperthermia treatment-cumulative equivalent minutes at 43 °C, CEM₄₃ -provides an inadequate description of the expected treatment effectiveness.

OBJECTIVES

Couple numerical models to experimental measurements to study the relative heating effectiveness described by cell death predictions.

METHODS

FEM numerical models were applied to increase the understanding of a carefully calibrated series of experiments in mouse mammary adenocarcinoma.

RESULTS

The numerical model results indicate that minimum tumor loadings between approximately 1.3 to 1.8 mg of Fe per cm³ of tumor tissue are required to achieve the experimentally observed temperatures in magnetic field strengths of 32 kA/m (rms) at 162 kHz.

CONCLUSION

We show that including multiple cell death processes operating in parallel within the numerical models provides valuable perspective on the likelihood of successful treatment.

SIGNIFICANCE

We show and believe that these assessment methods are more accurate than a single assessment figure of merit based only on the comparison of thermal histories, such as the CEM method.

Mathematical formulation and analysis of the nonlinear system reconstruction of the online image-guided adaptive control of hyperthermia

PURPOSE

A nonlinear system reconstruction can theoretically provide timely system reconstruction when designing a real-time image-guided adaptive control for multisource heating for hyperthermia. This clinical need motivates an analysis of the essential mathematical characteristics and constraints of such an approach.

METHODS

The implicit function theorem (IFT), the Karush-Kuhn-Tucker (KKT) necessary condition of optimality, and the Tikhonov-Phillips regularization (TPR) were used to analyze and determine the requirements of the optimal system reconstruction. Two mutually exclusive generic approaches were analyzed to reconstruct the physical system: The traditional full reconstruction and the recently suggested partial reconstruction. Rigorous mathematical analysis based on IFT, KKT, and TPR was provided for all four possible nonlinear reconstructions: (1) Nonlinear noiseless full reconstruction, (2) nonlinear noisy full reconstruction, (3) nonlinear noiseless partial reconstruction, and (4) nonlinear noisy partial reconstruction, when a class of nonlinear formulations of system reconstruction is employed.

RESULTS

Effective numerical algorithms for solving each of the aforementioned four nonlinear reconstructions were introduced and formal derivations and analyses were provided. The analyses revealed the necessity of adding regularization when partial reconstruction is used. Regularization provides the theoretical support for one to uniquely reconstruct the optimal system. It also helps alleviate the negative influences of unavoidable measurement noise. Both theoretical analysis and numerical examples showed the importance of having a good initial guess for accomplishing nonlinear system reconstruction.

CONCLUSIONS

Regularization is mandatory for partial reconstruction to make it well posed. The Tikhonov-Phillips regularized Gauss-Newton algorithm has nice theoretical performance for partial reconstruction of systems with and without noise. The Levenberg-Marquardt algorithm is a more robust algorithmic option compared to the Gauss-Newton algorithm for nonlinear full reconstruction. A severe limitation of nonlinear reconstruction is the time consuming calculations required for the derivatives of temperatures to unknowns. Developing a method of model reduction or implementing a parallel algorithm can resolve this. The results provided herein are applicable to hyperthermia with blood perfusion nonlinearly depending on temperature and in the presence of thermally significant blood vessels.

Effective learning strategies for real-time image-guided adaptive control of multiple-source hyperthermia applicators

PURPOSE

This paper investigates overall theoretical requirements for reducing the times required for the iterative learning of a real-time image-guided adaptive control routine for multiple-source heat applicators, as used in hyperthermia and thermal ablative therapy for cancer.

METHODS

Methods for partial reconstruction of the physical system with and without model reduction to find solutions within a clinically practical timeframe were analyzed. A mathematical analysis based on the Fredholm alternative theorem (FAT) was used to compactly analyze the existence and uniqueness of the optimal heating vector under two fundamental situations: (1) noiseless partial reconstruction and (2) noisy partial reconstruction. These results were coupled with a method for further acceleration of the solution using virtual source (VS) model reduction. The matrix approximation theorem (MAT) was used to choose the optimal vectors spanning the reduced-order subspace to reduce the time for system reconstruction and to determine the associated approximation error. Numerical simulations of the adaptive control of hyperthermia using VS were also performed to test the predictions derived from the theoretical analysis. A thigh sarcoma patient model surrounded by a ten-antenna phased-array applicator was retained for this purpose. The impacts of the convective cooling from blood flow and the presence of sudden increase of perfusion in muscle and tumor were also simulated.

RESULTS

By FAT, partial system reconstruction directly conducted in the full space of the physical variables such as phases and magnitudes of the heat sources cannot guarantee reconstructing the optimal system to determine the global optimal setting of the heat sources. A remedy for this limitation is to conduct the partial reconstruction within a reduced-order subspace spanned by the first few maximum eigenvectors of the true system matrix. By MAT, this VS subspace is the optimal one when the goal is to maximize the average tumor temperature. When more than 6 sources present, the steps required for a nonlinear learning scheme is theoretically fewer than that of a linear one, however, finite number of iterative corrections is necessary for a single learning step of a nonlinear algorithm. Thus, the actual computational workload for a nonlinear algorithm is not necessarily less than that required by a linear algorithm.

CONCLUSIONS

Based on the analysis presented herein, obtaining a unique global optimal heating vector for a multiple-source applicator within the constraints of real-time clinical hyperthermia treatments and thermal ablative therapies appears attainable using partial reconstruction with minimum norm least-squares method with supplemental equations. One way to supplement equations is the inclusion of a method of model reduction.

Online feedback focusing algorithm for hyperthermia cancer treatment

PURPOSE

Magnetic resonance (MR) imaging is increasingly being utilized to visualize the 3D temperature distribution in patients during treatment with hyperthermia or thermal ablation therapy. The goal of this work is to lay the foundation for improving the localization of heat in tumors with an online focusing algorithm that uses MR images as feedback to iteratively steer and focus heat into the target.

METHODS

The algorithm iteratively updates the model that quantifies the relationship between the source (antenna) settings and resulting tissue temperature distribution. At each step in the iterative process, optimal settings of power and relative phase of each antenna are computed to maximize averaged tumor temperature in the model. The MR-measured thermal distribution is then used to update/correct the model. This iterative procedure is repeated until convergence, i.e. until the model prediction and MR thermal image are in agreement. A human thigh tumor model heated in a 140 MHz four-antenna cylindrical mini-annular phased array is used for numerical validation of the proposed algorithm. Numerically simulated temperatures are used during the iterative process as surrogates for MR thermal images. Gaussian white noise with a standard deviation of 0.3 degrees C and zero mean is added to simulate MRI measurement uncertainty. The algorithm is validated for cases where the source settings for the first iteration are based on erroneous models: (1) tissue property variability, (2) patient position mismatch, (3) a simple idealized patient model built from CT-based actual geometry, and (4) antenna excitation uncertainty due to load dependent impedance mismatch and antenna cross-coupling. Choices of starting heating vector are also validated.

RESULTS

The algorithm successfully steers and focuses a tumor when there is no antenna excitation uncertainty. Temperature is raised to > or = 43 degrees C for more than about 90% of tumor volume, accompanied by less than about 20% of normal tissue volume being raised to a temperature > or = 41 degrees C. However, when there is antenna excitation uncertainty, about 40% to 80% of normal tissue volume is raised to a temperature > or = 41 degrees C. No significant tumor heating improvement is observed in all simulations after about 25 iteration steps.

CONCLUSIONS

A feedback control algorithm is presented and shown to be successful in iteratively improving the focus of tissue heating within a four-antenna cylindrical phased array hyperthermia applicator. This algorithm appears to be robust in the presence of errors in assumed tissue properties, including realistic deviations of tissue properties and patient position in applicator. Only moderate robustness was achieved in the presence of misaligned applicator/tumor positioning and antenna excitation errors resulting from load mismatch or antenna cross coupling.