425
Views
24
CrossRef citations to date
0
Altmetric
RESEARCH ARTICLE

3D versus 2D steering in patient anatomies: A comparison using hyperthermia treatment planning

, , &
Pages 74-85 | Received 29 Mar 2010, Accepted 11 Oct 2010, Published online: 04 Jan 2011

Abstract

Purpose: In this study hyperthermia treatment planning is used to investigate whether the target temperature during hyperthermia treatment can be increased using the 3D AMC-8 instead of the 2D AMC-4 system (AMC: Academic Medical Center).

Methods and materials: The heating ability of the AMC-4 and AMC-8 system was analysed for five patients with cervix uteri carcinoma. Dielectric and thermal models were generated, based on a hyperthermia planning computerised tomography (CT), at a resolution of 2.5 × 2.5 × 5.0 mm3. Calculation of the electric fields with the finite-difference time-domain method was followed by SAR- and temperature-based optimisation. The ability to correct for axial shifts of the patient by phase/amplitude steering was investigated for both systems. Finally, it was investigated whether adjusting the ring-to-ring distance of the AMC-8 system can be used for further optimisation.

Results: An average increase in T90 of ∼0.5°C (0.2°–0.8°C) was found for the AMC-8 system compared to the AMC-4 system. The gain in T50 and T10 was also 0.5°C on average. The additional power required to achieve this gain was 36% to 71% of the power required for the AMC-4 system. The AMC-8 system has the capability of correcting changes in axial position (−8 cm, +8 cm), contrary to the AMC-4 system. For both systems the axial position should be known within 1–2 cm.

Conclusions: Hyperthermia treatment with the AMC-8 system can lead to a clinically relevant increase of the target temperature compared to treatment with the AMC-4 system. The AMC-8 system provides large freedom in the axial positioning of the patient.

Introduction

Hyperthermia has proven its added value when combined with radiotherapy and/or chemotherapy in various randomised phase III trials for different tumour sites Citation[1–3].

These good clinical results were achieved in spite of the fact that the temperature aimed for, 43°C during one hour of treatment, was seldom achieved Citation[4]. The thermal dose administered to the target is generally lower because the potential local overheating of normal tissue prevents applying sufficient power. Increasing the thermal dose is expected to lead to improved clinical results, given that the fraction of tumour cells surviving combined radiotherapy and hyperthermia was found to decrease significantly with increasing target temperature and treatment duration Citation[5].

In most institutes, deep-seated tumours such as cervix uteri carcinoma are heated with loco-regional phased-array systems operating at frequencies in the range of 70–150 MHz Citation[6], Citation[7]. The distribution of the absorbed power can be steered by varying the amplitude and the phase of the electric fields imposed by the individual antennas. In this way the thermal dose delivered to the tumour can be maximised within the limits of patient tolerance.

An important strategy to improve thermal dose delivery during hyperthermia treatment is the use and design of systems with a higher level of control on the power distribution in the patient. The level of control of a loco-regional hyperthermia device is determined by a number of system parameters. Frequency determines the wavelength and thus the scale of the interference pattern. With increasing frequency, electromagnetic waves penetrate with less depth into the patient; at the same time, constructive interference between multiple sources takes place on a smaller spatial scale resulting in a smaller focus as a result of the decreased wavelength. In addition, if the frequency becomes too large, i.e. if the wavelength is much smaller than the typical dimension of the patient, it may no longer be possible to create a single focus.

The number of antennas and their geometric configuration is another important design parameter. Different previous studies Citation[8–11] investigated the effect of antenna number, configuration and applied frequency for different tumour sites on the achieved target temperature or power absorption. In these studies treatment quality was found to improve with an increasing number of antennas organised in multiple rings. The advantage of having a configuration with multiple rings of antennas is that the power or temperature distribution can be steered in 3D whereas a single ring system only allows for steering within a plane.

At our department, loco-regional hyperthermia treatments were carried out using a single ring, four waveguide system operating at 70 MHz: the 2D AMC-4 system Citation[12]. Recently, a new version of the system operating at the same frequency of 70 MHz with an extra ring of four waveguides was introduced in the clinic: the AMC-8 system () Citation[6].

Figure 1. The AMC-8 phased-array 70 MHz waveguide system. The system consists of two rings of four waveguides. Every waveguide has a separate water bolus that provides superficial cooling of the patient and coupling of the incident electromagnetic field into the patient. The two bottom waveguides share the same water bolus. The distance between the two rings is adjustable.

Figure 1. The AMC-8 phased-array 70 MHz waveguide system. The system consists of two rings of four waveguides. Every waveguide has a separate water bolus that provides superficial cooling of the patient and coupling of the incident electromagnetic field into the patient. The two bottom waveguides share the same water bolus. The distance between the two rings is adjustable.

Simulations and measurements in phantoms showed that the full-width half-maximum of the heating pattern was found to be prolonged with at least 50% in the axial direction, depending on the aperture size and the ring-to-ring distance Citation[13]. At equal superficial power deposition, the power in a target located at the centre of the phantom was found to increase when moving from four to eight waveguides with up to 30% Citation[6].

All these pre-implementation studies involved homogeneous phantoms; in this study hyperthermia treatment planning will be used to address the question of whether the target temperature during loco-regional hyperthermia treatment can be increased using the AMC-8 instead of the AMC-4 system. To this end, for both systems, the heating efficacy under standard positioning is investigated as well as its sensitivity with respect to changes in patient position. Furthermore, the ability of both systems to correct for those changes in position by adjustment of phases and amplitudes is investigated. Finally, it is investigated whether the variable ring-to-ring distance provides a useful additional degree of freedom to optimise heating with the AMC-8 system.

Materials and methods

Definition of the patient model

For a group of five patients with a cervix uteri carcinoma, a data-set was acquired with computerised tomography (0.9 × 0.9 × 5.0 mm3) and was automatically segmented into inner-air, fat, bone and muscle tissue, based on the Hounsfield units Citation[14]. The axial extent of the scans was approximately 60 cm. This is not sufficient for simulation of heating with the AMC-8 system and for this reason the scans were extended by copying the top and bottom slices. Appropriate dielectric and thermal parameters were assigned as reported in Citation[15], Citation[16]. A radiation oncologist delineated the gross target volume (GTV) as the target volume for optimisation. gives an overview of a number of patient characteristics.

Table I.  Tissue properties in the electromagnetic and thermal simulations.

Table II.  Overview of patient characteristics for the selected patients.

To reduce the computational costs, the data-set was re-sampled on a 2.5 × 2.5 × 5.0 mm3 grid. The technique used for re-sampling is referred to as the winner-takes-all technique: a low-resolution voxel is assigned the properties of the tissue type that takes up to largest fraction of its volume.

Definition of the phased-array wave guide system and patient positioning

The AMC-4 and AMC-8 system consist of four and eight 70 MHz waveguides (aperture size 20.2 × 34.3 cm2) organised in one and two rings, respectively (). Each waveguide can be translated independently in the direction normal to its aperture. In the clinic, the position of the waveguides is chosen such that a gap of 5 cm exists between the waveguides and the patient. Water boluses are placed between these gaps to ensure adequate coupling of the incident electromagnetic field. Furthermore, these water boluses, circulated with distilled water kept at approximately 12°C, provide superficial cooling that is essential to prevent overheating. For the AMC-8 system the distance between the two rings can be varied. Unless mentioned otherwise, this distance was set to 1.5 cm ().

Figure 2. Coronal cross-section of the waveguide configuration of the AMC-4 (left) and AMC-8 (right) phased-array system and the water boluses. Waveguide aperture size equals 20.2 × 34.3 cm2. The distance between the two rings of the AMC-8 system (ring-to-ring distance (RRD)) can be varied and by default in this study is 1.5 cm.

Figure 2. Coronal cross-section of the waveguide configuration of the AMC-4 (left) and AMC-8 (right) phased-array system and the water boluses. Waveguide aperture size equals 20.2 × 34.3 cm2. The distance between the two rings of the AMC-8 system (ring-to-ring distance (RRD)) can be varied and by default in this study is 1.5 cm.

By default, the patient model was positioned in the AMC-4 model such that the centre of gravity (COG) of the target volume coincides with centre of the applicator ring. For the AMC-8 system, the COG of the target volume coincided with the centre of the two rings. This means that the axial position of the COG is in between the two rings.

Calculation of the electric field distribution

The electric field distribution in the patient was determined by superposition of the fields imposed by the different antennas. Therefore, in the simulations, all waveguides were excited separately with unit amplitude and phase zero to calculate the individual fields.

To compute the electric field distribution, the finite difference time domain (FDTD) method was used Citation[17]. The computational domain was truncated using a perfectly matched layer Citation[18]. By discrete Fourier transformation the complex electric field was stored after 35 000 steps in time (Δt = 0.95 ΔtCourant). To accelerate the FDTD algorithm, a graphical processing unit (GPU) implementation was used Citation[19]. The computation time for a single field was 30–45 minutes for a simulation volume of typically 15 million cells.

Optimisation

To evaluate the performance of the two systems for a series of five patients, two optimisation procedures were followed: SAR- and temperature-based optimisation.

SAR-based optimisation

An optimisation procedure based on the specific absorption rate (SAR) was used to investigate the ability to focus the power absorption on the target region for varying axial positions of the patient for the AMC-8 and AMC-4 system. The patient was shifted in the caudal–cranial direction from −8 cm to +8 cm with 2-cm steps. For each position, the ratio of SAR delivered to the target and delivered to the patient in total was optimised Citation[20]. The optimisation procedure was constrained by requiring that the power delivered by a single antenna is within 10% and 40% and 5% and 25% of the total power for the AMC-4 and the AMC-8 system, respectively.

A drawback of this SAR-based optimisation procedure is that phase and amplitude settings that give optimal coverage of the target region may lead to unacceptably high SAR levels elsewhere. This problem could potentially be solved if there was a rational basis to constrain SAR to specific levels. However, in an (off-line) optimisation procedure, such constraints are hard to define since the corresponding temperature levels are largely dependent on variables such as tissue perfusion and other forms of cooling, e.g. water bolus cooling. Although SAR-based methods can be tailored to suppress overheating with certain efficiency Citation[21], temperature-based optimisation suppresses overheating naturally by setting the maximum tolerable normal tissue temperature as a constraint Citation[22].

Temperature-based optimisation

To investigate the ability of the two systems to heat the target while respecting normal tissue constraints, temperature-based optimisation was performed for varying axial patient positions. As for the SAR-based optimisation, the patient was shifted in a caudal-cranial direction from −8 to +8 cm in 2-cm steps. This shows which patient position leads to optimal heating and, by shifting the patient without changing the phase/amplitude settings, how sensitive the thermal dose is for positioning errors. Das et al. showed, by using Pennes’ bio-heat equation Citation[23], how phase and amplitude can be related to the local temperature T(x, y, z) using the previously calculated electric fields Citation[24]. This method was used in this study. The relevant steady-state temperature fields Citation[24] were calculated using a GPU accelerated finite-difference method. The relevant parameters for the temperature calculations are presented in . We defined the following object function for optimisation Citation[25]:where Vtarget is the target volume and is a complex valued vector holding the amplitudes and phases of the different sources. This function is minimised under the constraint that the temperature in normal tissue does not exceed a constraint temperature here chosen to be 45°C Citation[26]. The value of Tthreshold was set to 43°C. Additional constraints were specified for the contribution to the total power of a single antenna. For the AMC-4 system this contribution should be within 10% and 40%, for the AMC-8 system within 5% and 25%. Solutions to this non-linear constrained optimisation problem were computed using the CFSQP package (eps = 1e-5, epseqn = 1e-5, miter = 35/50 (AMC-4/AMC-8)) Citation[27]. Since the optimisation process might halt on a local minimum, 20 runs with random starting points in phase/amplitude space were performed for the AMC-4 system and 10 runs for the AMC-8 system (to limit the computation time).

In order to determine the potential gain in tumour temperature when treating with the eight-antenna system instead of treating with four antennas, T90, (for 90% of the tumour volume T ≥ T90), T50 and T10 of the target volume were compared. Furthermore the total power absorbed in the patient was calculated. Cumulative temperature–volume histograms (TVH) were made to analyse the differences in the overall temperature distribution.

Variation of the ring-to-ring distance.

Since the distance between the two rings of the AMC-8 system can be varied, it was investigated whether this is a relevant degree of freedom in the optimisation procedure. Hence, temperature-based optimisation was performed setting the ring-to-ring distance to 1.5, 3.5, 5.5 and 7.5 cm. In addition, the sensitivity of the thermal dose to the ring-to-ring distance was investigated. This was done by applying the settings resulting from optimisation at 1.5 cm ring-to-ring distance and evaluating the thermal dose at larger distances.

Results

SAR-based optimisation

shows the relative SAR ratio (SAR in the target divided by SAR in the total volume) as a function of axial position for both systems for all five patients after optimisation of this ratio. As a scaling factor, the SAR ratio under optimal positioning for the AMC-4 system is used. This figure shows that for the AMC-8 system the SAR ratio increases moving the patient to off-centre positions, whereas for the AMC-4 system the SAR ratio decreases.

Figure 3. The normalised ratio of target- and normal tissue SAR for the AMC-4 and the AMC-8 system for different axial positions after optimisation for patients 1 to 5. The SAR ratios are normalised by the SAR ratio found for the AMC-4 system under optimal positioning. Antennas were not allowed to contribute more to the total delivered power than 40%, or less than 10% for the AMC-4 system, and no more than 25% or less than 5% for the AMC-8 system.

Figure 3. The normalised ratio of target- and normal tissue SAR for the AMC-4 and the AMC-8 system for different axial positions after optimisation for patients 1 to 5. The SAR ratios are normalised by the SAR ratio found for the AMC-4 system under optimal positioning. Antennas were not allowed to contribute more to the total delivered power than 40%, or less than 10% for the AMC-4 system, and no more than 25% or less than 5% for the AMC-8 system.

Temperature-based optimisation

shows the increase in target T10, T50 and T90 for the AMC-8 compared to the AMC-4 system after temperature-based optimisation applying standard positioning. T90 shows an average increase of ∼0.5°C (range 0.2°–0.8°C). T50 and T10 show a similar increase.

Figure 4. The temperature gain for T10, T50 and T90 when moving from the AMC-4 to the AMC-8 system for patients 1 to 5. The applied settings resulted from temperature-based optimisation. For the AMC-8 system a 1.5-cm ring-to-ring distance was applied.

Figure 4. The temperature gain for T10, T50 and T90 when moving from the AMC-4 to the AMC-8 system for patients 1 to 5. The applied settings resulted from temperature-based optimisation. For the AMC-8 system a 1.5-cm ring-to-ring distance was applied.

shows the additional power needed to achieve the higher target temperatures for the AMC-8 system. The increase in total power is found to be between 36% and 71% of the power that was required for heating with the AMC-4 system.

Figure 5. Total absorbed power for the AMC-4 and the AMC-8 system for patients 1 to 5. The applied settings result from temperature-based optimisation. For the AMC-8 system a 1.5-cm ring-to-ring distance was applied.

Figure 5. Total absorbed power for the AMC-4 and the AMC-8 system for patients 1 to 5. The applied settings result from temperature-based optimisation. For the AMC-8 system a 1.5-cm ring-to-ring distance was applied.

As illustrated in , a gain in target temperature is associated with an increase of the volume with a temperature close to the constraint level. As an example, for patient 1 the volume with a temperature of 44°C is ten times larger when treating with the AMC-8 compared to the volume when treating with the AMC-4 system.

Figure 6. Cumulative temperature−volume histograms of the temperature distribution resulting from temperature-based optimisation for the five patients. The inserts zoom in to the 43−45°C range of the histogram. Indicated in the lower left corner of each plot is the patient number.

Figure 6. Cumulative temperature−volume histograms of the temperature distribution resulting from temperature-based optimisation for the five patients. The inserts zoom in to the 43−45°C range of the histogram. Indicated in the lower left corner of each plot is the patient number.

An example of the temperature distribution after temperature-based optimisation for the AMC-4 and the AMC-8 system for one of the five patients (patient 1) is shown in . It can be seen that the temperature distribution extends over a larger length in the axial direction (here particularly clear for the legs). Higher tumour temperatures are achieved, however, for the same maximum normal tissue temperature.

Figure 7. An example of the resulting temperature distribution (right) after temperature-based optimisation with the corresponding coronal slices (left) for the AMC-4 system (top) and the AMC-8 system (bottom).

Figure 7. An example of the resulting temperature distribution (right) after temperature-based optimisation with the corresponding coronal slices (left) for the AMC-4 system (top) and the AMC-8 system (bottom).

Patient positioning

compares the T90 and T50 found using settings from temperature-based optimisation for varying axial patient positions. For every shift T90 and T50 are evaluated using settings that result from I) optimisation with the patient model at the correct/actual position (unconnected markers in ) and II) optimisation without a shift of the patient (connected markers). It shows the ability of the two systems to correct for axial shifts of the patient. For the AMC-8 system, both the T90 and T50 are found to be relatively stable for all patients over a range of −8 cm to +8 cm. Without correction for axial patient shifts, the thermal dose is found to be equally sensitive for both systems (power was scaled to the highest tolerable level, i.e. such that the maximum temperature equals the constraint temperature). The ability of the AMC-4 system to correct for patient position variations is clearly very limited.

Figure 8. T90 and T50 for different axial shifts (z-axis) of patients 1 to 5. For every shift the T90 and T50 were evaluated for two sets of amplitudes and phases. The first is the result of optimisation with the patient modelled at the correct/actual position (unconnected markers). The second results from optimisation without shifting the patient (connected markers). The figure shows that the AMC-8 system has the capability of correcting for changes in the axial position of the patient, contrary to the AMC-4 system.

Figure 8. T90 and T50 for different axial shifts (z-axis) of patients 1 to 5. For every shift the T90 and T50 were evaluated for two sets of amplitudes and phases. The first is the result of optimisation with the patient modelled at the correct/actual position (unconnected markers). The second results from optimisation without shifting the patient (connected markers). The figure shows that the AMC-8 system has the capability of correcting for changes in the axial position of the patient, contrary to the AMC-4 system.

Ring-to-ring distance

shows the T90 heating with the AMC-8 system using different distances between the two rings relative to the results found using a default 1.5 cm ring-to-ring distance. The left figure shows the effect of increasing the distance without adjustment of phase/amplitude settings (total power is set to the highest tolerable level). Minor deviations (<2 cm) in the ring-to-ring distance led to a decrease of the T90 below 0.3°C (with exception of patient 3). The figure on the right shows the difference in T90 after optimising for 3.5, 5.5 and 7.5 cm distance compared to the T90 found using the ‘standard’ distance of 1.5 cm. The optimal ring-to-ring distance is 5.5 cm for patient 1, 3.5 cm for patient 2, 1.5 cm for patient 3, 7.5 cm for patient 4 and 1.5 cm for patient 5, but in general adjusting the ring-to-ring distance does not lead to a profound improvement in thermal dose delivery.

Figure 9. Differences in T90 after temperature optimisation for the AMC-8 system using varying ring-to-ring distances. Results are relative to the results for a ring-to-ring distance of 1.5 cm. The figure on the left shows the change in T90 applying optimised settings for a 1.5 cm ring-to-ring distance at larger distances without adjustment of settings. The figure on the right shows the T90 optimising for the actual distances relative to the optimal T90 for a ring-to-ring distance of 1.5 cm.

Figure 9. Differences in T90 after temperature optimisation for the AMC-8 system using varying ring-to-ring distances. Results are relative to the results for a ring-to-ring distance of 1.5 cm. The figure on the left shows the change in T90 applying optimised settings for a 1.5 cm ring-to-ring distance at larger distances without adjustment of settings. The figure on the right shows the T90 optimising for the actual distances relative to the optimal T90 for a ring-to-ring distance of 1.5 cm.

Discussion

This study applied hyperthermia treatment planning to investigate whether the use of the eight antenna, double ring, 3D AMC-8 system is expected to result in higher target temperatures relative to treatment with the single ring, 2D AMC-4 system. A similar comparison is relevant for the 2D BSD Sigma 60 and the 3D BSD Sigma Eye applicator.

showed that the position of the patient in either the AMC-4 or the AMC-8 system plays an important role in the fraction of the total power that can be delivered to the target. This analysis did not take hot-spots into account. Consequently, no statements can be made about the tolerable power level which is to be determined by the maximum temperature in the patient. For the AMC-4 system, moving the patient away from the mid-plane reduces the fraction of power that can be delivered to the target. For the AMC-8 system, moving the patient away from the transversal mid-plane in between the two rings is beneficial for the power delivered to the target. However, superficial levels of power absorption will increase leading to a lower tolerable power level in clinical practice. Different behaviour is expected to occur when comparing the BSD Sigma 60 (2D, 1 ring) and the BSD Sigma Eye applicator (3D, 3 rings). Under standard positioning, in a three ring set-up the target is in the centre of the middle ring, while for the two ring AMC-8 system the target is in between the rings.

Temperature-based optimisation can be constrained more straightforwardly by setting a maximum normal tissue temperature and is therefore preferred over SAR-based optimisation. However, a disadvantage of temperature-based optimisation is the uncertainty in the perfusion distribution under hyperthermic conditions. With the present-day lack of perfusion data, accurate clinical patient-specific temperature-based optimisation is not feasible. The gain found in this study is therefore expected to be dependent on thermal modelling parameters. Our estimates are assumed to represent the correct average but will deviate for the individual patients. However, temperature-based optimisation is clinically the most relevant method to evaluate the heating of the same patient with different apparatus by means of simulation.

Temperature-based optimisation of the target temperature showed an average increase of the target T90 with 0.5°C with large inter-patient variability (0.2°–0.8°C) for the AMC-8 system compared to the AMC-4 system. T50 and T10 increased typically with 0.5°C as well (). This increase is clinically relevant since the cumulative minutes at 43°C, a parameter commonly used to describe thermal dose, approximately doubles with every 0.5°C increase in temperature (for temperatures below 42.5°C) Citation[28].

As shown in , for the AMC-8 system the tumour temperature expressed as T90 and T50 is relatively insensitive to changes in the axial position of the patient when the optimisation takes the actual position into account. These results are comparable to those of Gellermann et al. for the 3D BSD Sigma Eye applicator Citation[29]. Altering the axial position of the patient might prove to be a strategy to relieve complaints reported in earlier treatment sessions. The absorption at a specific hot-spot can be diminished while the tumour temperature is maintained at the desired level. The impact on the temperature distribution of shifting the patient in the axial direction differs from changing phases and amplitudes. By shifting the patient the electric field distribution of the individual fields is changed while those under phase/amplitude steering obviously remain unchanged. In addition, shifting the patient in the axial direction can be a way to treat patients that under standard positioning have insufficient length to be treated with two rings. The results of and show that positioning accuracy in the axial direction should be within 1–2 cm for both the AMC-4 and AMC-8 system. If this accuracy is achieved, differences in target temperature remain within 0.25°–0.5°C. These are substantial differences compared to the increase in target temperature heating with the AMC-8 instead of the AMC-4 system. Canters et al. Citation[30] found that for the 2D BSD-2000 Sigma 60 system, shifting the patient in the positive z-direction (equivalent to the shifting of the applicator in the negative z-direction) and applying optimised settings for standard positioning reduced the hot-spot–target quotient (HTQ) for which they presented a strong correlation with the target T50. Similar behaviour is not observed in our study since, by optimising directly for temperature, changing the position while not adjusting the settings is expected to make tumour heating less effective. However, it is observed for the AMC-4 system that if the temperature distribution is optimised for varying axial positions, for patients 1, 2, 4 and 5 () a shift in caudal direction of 2 cm (4 cm is optimal for patient 5) causes an increase in T50. This shift corresponds to an applicator shift of 2 cm as discussed in Citation[30].

addresses the question whether the variable ring-to-ring distance (RRD) provides a useful extra degree of freedom in the optimisation of the thermal dose. As the figure illustrates small, patient-specific, differences are found when changing this distance. For patients 1 and 4 relevant improvements were found, while for patients 2, 3 and 5 the standard RRD was either the best choice or improvements were very small (<0.1°C). Given the computational costs associated with including this parameter in the optimisation procedure (those increase with a factor of 4), a standard ring-to-ring distance of 1.5 cm was considered to be acceptable. Moreover, the uncertainties currently present in HTP render it questionable whether relatively small improvements can be predicted accurately for individual patients.

The variable ring-to-ring distance is a unique feature of the AMC-8 system. However, when designing multi-ring systems a choice has to be made for the axial spacing of the antennas as well. Extrapolating our results, in general the impact on heating efficacy will be minor and other requirements such as limiting applicator coupling might have priority.

In our study the patient model is based on the automated segmentation of a CT scan resulting in only a limited number of tissue types. The accuracy of the segmentation in terms of dielectric and thermally alike tissues/tissue types is, however, expected to influence the predicted temperature distribution. However, hot-spots are dominantly located at interfaces of bone, fat, and muscle tissue i.e. high dielectric contrast, in the superficial fat layers (sufficiently deep for cooling to be ineffective), and in the proximity of bony structures due to resonance effects. These hot-spots are expected to dominate the optimisation result. The purpose of our study was to compare two different systems based on patient models that are representative for the population and hence the level of detail used in our study is considered appropriate.

The volume of the patient that is heated with both the AMC-4 as well the AMC-8 system is relatively large compared to the target volume, since these systems are loco-regional hyperthermia devices. The consequences for thermal modelling when large volumes are heated are presently unclear. All thermal modelling done was based on Pennes’ bio-heat equation that is particularly suited for the modelling of small-scale heating problems where arterial pre-heating is expected to be absent. Whether heating of such large volumes can still modelled accurately is unclear. Van den Berg et al. Citation[31] studied the modelling of vascular heat transfer for prostate hyperthermia. They concluded that pre-heating in the arteries feeding the prostate vasculature can be neglected whereas accurate thermal modelling requires including the prostate vasculature. Based on the equilibration length, increasing the heated volume by heating with eight versus four antennas is not expected to cause important differences.

Another effect, not covered in the simulations is the increase in systemic temperature Citation[32]. As shown, a major increase of power is required, and moreover feasible, when heating with eight antennas. As shown in , the total power absorbed in the patient to realise the optimised temperature distribution for the AMC-8 system was in the range of 752–927 W as compared to 468–607 W for the AMC-4 system. This is also reflected in and , showing that the heated volume is significantly larger for the AMC-8 system. This may result in an extra increase of the systemic temperature that has an important impact on the heat transfer via the vasculature and for example lead to an additional increase in T90. To relate the increase of power to an increase in systemic temperature, further clinical data has to be acquired.

Conclusions

Loco-regional hyperthermia treatment with the 3D AMC-8 system can lead to a clinically relevant increase of the target temperature compared to treatment with the 2D AMC-4 system. The T90 and T50 of the target increased with 0.5°C but with large patient variability. The additional power required to realise this increase ranged from 36% to 71% of the power required for the AMC-4 system. Ring-to-ring distance is of minor importance, but should be known for proper optimisation within 2 cm. For the AMC-8 system, the achievable thermal dose is stable over a range of −8 cm to +8 cm shifting the patient in the caudal–cranial direction. Axial position should be known within 1–2 cm both for the AMC-4 and the AMC-8 system.

Declaration of interest: The work described in this article was financially supported by the Dutch Cancer Society (grant UVA – 2004 – 3834). The authors alone are responsible for the content and writing of the paper.

References

  • Van der Zee J, Gonzalez GD, Van Rhoon GC, van Dijk JDP, van Putten WLJ, Hart AA. Comparison of radiotherapy alone with radiotherapy plus hyperthermia in locally advanced pelvic tumours: A prospective, randomised, multicentre trial. Dutch Deep Hyperthermia Group. Lancet 2000; 355: 1119–1125
  • Overgaard J, Gonzalez GD, Hulshof MCCM, Arcangeli G, Dahl O, Mella O, Bentzen SM. Randomised trial of hyperthermia as adjuvant to radiotherapy for recurrent or metastatic malignant melanoma. European Society for Hyperthermic Oncology. Lancet 1995; 345: 540–543
  • Issels RD. Hyperthermia adds to chemotherapy. Eur J Cancer 2008; 44: 2546–2554
  • Franckena M, Fatehi D, de Bruijne M, Canters RA, van NY, Mens JW, Van Rhoon GC, Van der Zee J. Hyperthermia dose-effect relationship in 420 patients with cervical cancer treated with combined radiotherapy and hyperthermia. Eur J Cancer 2009; 45: 1969–1978
  • Sapareto SA, Dewey WC. Thermal dose determination in cancer therapy. Int J Radiat Oncol Biol Phys 1984; 10: 787–800
  • Crezee J, Van Haaren PMA, Westendorp H, De Greef M, Kok HP, Wiersma J, Van Stam G, Sijbrands J, Zum Vörde Sive Vörding PJ, Van Dijk JDP, et al. Improving locoregional hyperthermia delivery using the 3-D controlled AMC-8 phased array hyperthermia system: A preclinical study. Int J Hyperthermia 2009; 25: 581–592
  • Van Rhoon GC, Van Der Heuvel DJ, Ameziane A, Rietveld PJ, Volenec K, van der ZJ. Characterization of the SAR-distribution of the Sigma-60 applicator for regional hyperthermia using a Schottky diode sheet. Int J Hyperthermia 2003; 19: 642–654
  • Seebass M, Beck R, Gellermann J, Nadobny J, Wust P. Electromagnetic phased arrays for regional hyperthermia: Optimal frequency and antenna arrangement. Int J Hyperthermia 2001; 17: 321–336
  • Paulsen KD, Geimer S, Tang J, Boyse WE. Optimization of pelvic heating rate distributions with electromagnetic phased arrays. Int J Hyperthermia 1999; 15: 157–186
  • Kroeze H, van de Kamer JB, De Leeuw AAC, Lagendijk JJW. Regional hyperthermia applicator design using FDTD modelling. Phys Med Biol 2001; 46: 1919–1935
  • Paulides MM, Bakker JF, Zwamborn AP, Van Rhoon GC. A head and neck hyperthermia applicator: Theoretical antenna array design. Int J Hyperthermia 2007; 23: 59–67
  • van Dijk JDP, Schneider CJ, van Os RM, Blank LE, Gonzalez DG. Results of deep body hyperthermia with large waveguide radiators. Adv Exp Med Biol 1990; 267: 315–319
  • Wiersma J. Hyperthermia treatment planning. PhD thesis (chapter 6), University of Amsterdam, 2005.
  • James BJ, Sullivan DM. Creation of 3-dimensional patient models for hyperthermia treatment planning. IEEE Trans Biomed Eng 1992; 39: 238–242
  • Gabriel C, Gabriel S, Corthout E. The dielectric properties of biological tissues: I. Literature survey. Phys Med Biol 1996; 41: 2231–2249
  • ESHO Taskgroup Committee. Treatment Planning and Modelling in Hyperthermia, a Task Group Report of the European Society for Hyperthermic Oncology. Rome: Tor Vergata.
  • Yee KS. Numerical solutions of initial boundary value problems involving Maxwell's equations in isotropic media. IEEE Trans Antennas Propagat 1966; 14: 302–307
  • Berenger JP. A perfectly matched layer for the absorption of electromagnetic-waves. J Comput Physics 1994; 114: 185–200
  • Correia D, De Greef M, Kok HP, Bel A, Crezee J. Fast FDTD computation using graphics card: The orientation of the yee-cell. In: ESHO-2009: Abstract Book, 2009, p. 32. Verona, Italy: European Society for Hyperthermic Oncology.
  • Bardati F, Borrani A, Gerardino A, Lovisolo GA. SAR optimization in a phased array radiofrequency hyperthermia system. Specific absorption rate. IEEE Trans Biomed Eng 1995; 42: 1201–1207
  • Canters RA, Wust P, Bakker JF, Van Rhoon GC. A literature survey on indicators for characterisation and optimisation of SAR distributions in deep hyperthermia, a plea for standardisation. Int J Hyperthermia 2009; 25: 593–608
  • De Greef M, Kok HP, Bel A, Crezee J. Optimization in hyperthermia treatment planning: The impact of tissue perfusion uncertainty. Med Phys 2010; 37: 4540–4550
  • Pennes HH. Analysis of tissue and arterial blood temperatures in the resting human forearm. 1948. J Appl Physiol 1948; 1: 93–122
  • Das SK, Clegg ST, Samulski TV. Computational techniques for fast hyperthermia temperature optimization. Med Phys 1999; 26: 319–328
  • Kok HP, Van Haaren PMA, van de Kamer JB, Wiersma J, Van Dijk JDP, Crezee J. High-resolution temperature-based optimization for hyperthermia treatment planning. Phys Med Biol 2005; 50: 3127–3141
  • Lagendijk JJ, Van Rhoon GC, Hornsleth SN, Wust P, De Leeuw AC, Schneider CJ, van Dijk JD, van der ZJ, Van Heek-Romanowski R, Rahman SA, et al. ESHO quality assurance guidelines for regional hyperthermia. Int J Hyperthermia 1998; 14: 125–133
  • Lawrence C, Zhou JL, Tits AL. User's Guide for CFSQP version 2.5d: A C Code for Solving (Large Scale) Nonlinear (Minimax) Optimization Problems, Generating Iterates Satisfying All Inequality Constraints. Electrical Engineering Department and Institute for Systems Research, University of Maryland.
  • Hand JW, Lagendijk JJ, Bach AJ, Bolomey JC. Quality assurance guidelines for ESHO protocols. Int J Hyperthermia 1989; 5: 421–428
  • Gellermann J, Goke J, Figiel R, Weihrauch M, Cho CH, Budach V, Felix R, Wust P. Simulation of different applicator positions for treatment of a presacral tumour. Int J Hyperthermia 2007; 23: 37–47
  • Canters RA, Franckena M, Paulides MM, Van Rhoon GC. Patient positioning in deep hyperthermia: Influences of inaccuracies, signal correction possibilities and optimization potential. Phys Med Biol 2009; 54: 3923–3936
  • Van den Berg CAT, van de Kamer JB, De Leeuw AAC, Jeukens CRLPN, Raaymakers BW, Van Vulpen M, Lagendijk JJW. Towards patient-specific thermal modelling of the prostate. Phys Med Biol 2006; 51: 809–825
  • Van Haaren PM, Hulshof MC, Kok HP, Oldenborg S, Geijsen ED, van Lanschot JJ, Crezee J. Relation between body size and temperatures during locoregional hyperthermia of oesophageal cancer patients. Int J Hyperthermia 2008; 24: 663–674

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.