Esra Neufeld

The Virtual Family—development of surface-based anatomical models of two adults and two children for dosimetric simulations

The objective of this study was to develop anatomically correct whole body human models of an adult male (34 years old), an adult female (26 years old) and two children (an 11-year-old girl and a six-year-old boy) for the optimized evaluation of electromagnetic exposure. These four models are referred to as the Virtual Family. They are based on high resolution magnetic resonance (MR) images of healthy volunteers. More than 80 different tissue types were distinguished during the segmentation. To improve the accuracy and the effectiveness of the segmentation, a novel semi-automated tool was used to analyze and segment the data. All tissues and organs were reconstructed as three-dimensional (3D) unstructured triangulated surface objects, yielding high precision images of individual features of the body. This greatly enhances the meshing flexibility and the accuracy with respect to thin tissue layers and small organs in comparison with the traditional voxel-based representation of anatomical models. Conformal computational techniques were also applied. The techniques and tools developed in this study can be used to more effectively develop future models and further improve the accuracy of the models for various applications. For research purposes, the four models are provided for free to the scientific community.

Simulation techniques in hyperthermia treatment planning

Abstract Clinical trials have shown that hyperthermia (HT), i.e. an increase of tissue temperature to 39–44 °C, significantly enhance radiotherapy and chemotherapy effectiveness [1]. Driven by the developments in computational techniques and computing power, personalised hyperthermia treatment planning (HTP) has matured and has become a powerful tool for optimising treatment quality. Electromagnetic, ultrasound, and thermal simulations using realistic clinical set-ups are now being performed to achieve patient-specific treatment optimisation. In addition, extensive studies aimed to properly implement novel HT tools and techniques, and to assess the quality of HT, are becoming more common. In this paper, we review the simulation tools and techniques developed for clinical hyperthermia, and evaluate their current status on the path from ‘model’ to ‘clinic’. In addition, we illustrate the major techniques employed for validation and optimisation. HTP has become an essential tool for improvement, control, and assessment of HT treatment quality. As such, it plays a pivotal role in the quest to establish HT as an efficacious addition to multi-modality treatment of cancer.

The clinical feasibility of deep hyperthermia treatment in the head and neck: new challenges for positioning and temperature measurement.

To apply high-quality hyperthermia treatment to tumours at deep locations in the head and neck (H&N), we have designed and built a site-specific phased-array applicator. Earlier, we demonstrated its features in parameter studies, validated those by phantom measurements and clinically introduced the system. In this paper we will critically review our first clinical experiences and demonstrate the pivotal role of hyperthermia treatment planning (HTP). Three representative patient cases (thyroid, oropharynx and nasal cavity) are selected and discussed. Treatment planning, the treatment, interstitially measured temperatures and their interrelation are analysed from a physics point of view. Treatments lasting 1 h were feasible and well tolerated and no acute treatment-related toxicity has been observed. Maximum temperatures measured are in the range of those obtained during deep hyperthermia treatments in the pelvic region but mean temperatures are still to be improved. Further, we found that simulated power absorption correlated well with measured temperatures illustrating the validity of our treatment approach of using energy profile optimizations to arrive at higher temperatures. This is the first data proving that focussed heating of tumours in the H&N is feasible. Further, HTP proved a valuable tool in treatment optimization. Items to improve are (1) the transfer of HTP settings into the clinic and (2) the registration of the thermal dose, i.e. dosimetry.

Noninvasive Deep Brain Stimulation via Temporally Interfering Electric Fields

Summary We report a noninvasive strategy for electrically stimulating neurons at depth. By delivering to the brain multiple electric fields at frequencies too high to recruit neural firing, but which differ by a frequency within the dynamic range of neural firing, we can electrically stimulate neurons throughout a region where interference between the multiple fields results in a prominent electric field envelope modulated at the difference frequency. We validated this temporal interference (TI) concept via modeling and physics experiments, and verified that neurons in the living mouse brain could follow the electric field envelope. We demonstrate the utility of TI stimulation by stimulating neurons in the hippocampus of living mice without recruiting neurons of the overlying cortex. Finally, we show that by altering the currents delivered to a set of immobile electrodes, we can steerably evoke different motor patterns in living mice.

Measurement, simulation and uncertainty assessment of implant heating during MRI

The heating of tissues around implants during MRI can pose severe health risks, and careful evaluation is required for leads to be labeled as MR conditionally safe. A recent interlaboratory comparison study has shown that different groups can produce widely varying results (sometimes with more than a factor of 5 difference) when performing measurements according to current guidelines. To determine the related difficulties and to derive optimized procedures, two different generic lead structures have been investigated in this study by using state-of-the-art temperature and dosimetric probes, as well as simulations for which detailed uncertainty budgets have been determined. The agreement between simulations and measurements is well within the combined uncertainty. The study revealed that the uncertainty can be kept below 17% if appropriate instrumentation and procedures are applied. Optimized experimental assessment techniques can be derived from the findings presented herein.

MIDA: A Multimodal Imaging-Based Detailed Anatomical Model of the Human Head and Neck

Computational modeling and simulations are increasingly being used to complement experimental testing for analysis of safety and efficacy of medical devices. Multiple voxel- and surface-based whole- and partial-body models have been proposed in the literature, typically with spatial resolution in the range of 1–2 mm and with 10–50 different tissue types resolved. We have developed a multimodal imaging-based detailed anatomical model of the human head and neck, named “MIDA”. The model was obtained by integrating three different magnetic resonance imaging (MRI) modalities, the parameters of which were tailored to enhance the signals of specific tissues: i) structural T1- and T2-weighted MRIs; a specific heavily T2-weighted MRI slab with high nerve contrast optimized to enhance the structures of the ear and eye; ii) magnetic resonance angiography (MRA) data to image the vasculature, and iii) diffusion tensor imaging (DTI) to obtain information on anisotropy and fiber orientation. The unique multimodal high-resolution approach allowed resolving 153 structures, including several distinct muscles, bones and skull layers, arteries and veins, nerves, as well as salivary glands. The model offers also a detailed characterization of eyes, ears, and deep brain structures. A special automatic atlas-based segmentation procedure was adopted to include a detailed map of the nuclei of the thalamus and midbrain into the head model. The suitability of the model to simulations involving different numerical methods, discretization approaches, as well as DTI-based tensorial electrical conductivity, was examined in a case-study, in which the electric field was generated by transcranial alternating current stimulation. The voxel- and the surface-based versions of the models are freely available to the scientific community.

Novel conformal technique to reduce staircasing artifacts at material boundaries for FDTD modeling of the bioheat equation

The modeling of thermal effects, often based on the Pennes Bioheat Equation, is becoming increasingly popular. The FDTD technique commonly used in this context suffers considerably from staircasing errors at boundaries. A new conformal technique is proposed that can easily be integrated into existing implementations without requiring a special update scheme. It scales fluxes at interfaces with factors derived from the local surface normal. The new scheme is validated using an analytical solution, and an error analysis is performed to understand its behavior. The new scheme behaves considerably better than the standard scheme. Furthermore, in contrast to the standard scheme, it is possible to obtain with it more accurate solutions by increasing the grid resolution.

Thermal Tissue Damage Model Analyzed for Different Whole‐Body SAR and Scan Durations for Standard MR Body Coils

This article investigates the safety of radiofrequency induced local thermal hotspots within a 1.5T body coil by assessing the transient local peak temperatures as a function of exposure level and local thermoregulation in four anatomical human models in different Z‐positions.

Children and adults exposed to low-frequency magnetic fields at the ICNIRP reference levels: theoretical assessment of the induced electric fields

To avoid potentially adverse health effects, the International Commission on Non-Ionizing Radiation Protection (ICNIRP) has defined reference levels for time varying magnetic fields. Restrictions on the electric fields induced in the human body are provided based on biological response data for peripheral nerve stimulation and the induction of phosphenes. Numerical modeling is commonly used to assess the induced electric fields for various exposure configurations. The objective of this study was to assess the variations of the electric fields induced in children and adults and to compare the exposure at reference levels with the basic restrictions as function of anatomy. We used the scalar potential finite element method to calculate the induced electric fields in six children and two adults when exposed to uniform magnetic fields polarized in three orthogonal directions. We found that the induced electric fields are within the ICNIRP basic restrictions in nearly all cases. In PNS tissues, we found electric fields up to 95% (upper uncertainty limit due to discretization errors, k = 2) of the ICNIRP basic restrictions for exposures at the general public reference levels. For occupational reference levels, we found an over-exposure of maximum 79% (k = 2) in PNS tissues. We further found that the ICNIRP recommendations on spatial averaging in 2 × 2 × 2 mm³ contiguous tissue volumes and removal of peak values by the 99th percentile cause the results to depend strongly on the grid discretization step (i.e. an uncertainty of more than 50% at 2 mm) and the number of distinguished tissues in the anatomical models. The computational results obtained by various research institutes should be robust for different discretization settings and various anatomical models. Therefore, we recommend considering alternative routines for small anatomical structures such as non-contiguous averaging without taking the 99th percentile in future guidelines leading to consistent suppression of peak values amongst different simulation settings and anatomical models. The peak electric fields depend on the local tissue distribution in the various anatomical models, and we could not find a correlation with the size of the anatomy. Therefore, we recommend extending the evaluation using a sufficient set of anatomies including other than standing postures to assess the worst-case exposure setting and correspondence to the basic restrictions.

Children and adults exposed to electromagnetic fields at the ICNIRP reference levels: theoretical assessment of the induced peak temperature increase

To avoid potentially adverse health effects of electromagnetic fields (EMF), the International Commission on Non-Ionizing Radiation Protection (ICNIRP) has defined EMF reference levels. Restrictions on induced whole-body-averaged specific absorption rate (SAR(wb)) are provided to keep the whole-body temperature increase (T(body, incr)) under 1 °C during 30 min. Additional restrictions on the peak 10 g spatial-averaged SAR (SAR(10g)) are provided to prevent excessive localized tissue heating. The objective of this study is to assess the localized peak temperature increase (T(incr, max)) in children upon exposure at the reference levels. Finite-difference time-domain modeling was used to calculate T(incr, max) in six children and two adults exposed to orthogonal plane-wave configurations. We performed a sensitivity study and Monte Carlo analysis to assess the uncertainty of the results. Considering the uncertainties in the model parameters, we found that a peak temperature increase as high as 1 °C can occur for worst-case scenarios at the ICNIRP reference levels. Since the guidelines are deduced from temperature increase, we used T(incr, max) as being a better metric to prevent excessive localized tissue heating instead of localized peak SAR. However, we note that the exposure time should also be considered in future guidelines. Hence, we advise defining limits on T(incr, max) for specified durations of exposure.