# Spot Weight Adaptation for Moving Target in Spot Scanning Proton Therapy

Paul Morel, Xiadong Wu, Guillaume Blin, Stéphane Vialette, Ryan Flynn,

Daniel Hyer, Dongxu Wang

► To cite this version:

Paul Morel, Xiadong Wu, Guillaume Blin, Stéphane Vialette, Ryan Flynn, et al. Spot Weight Adaptation for Moving Target in Spot Scanning Proton Therapy. Frontiers in Oncology, 2015, 5 (119), pp.1-7. 10.3389/fonc.2015.00119. hal-01151923

HAL Id: hal-01151923

https://hal.science/hal-01151923

Submitted on 13 Dec 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Spot weight adaptation for moving target in spot scanning proton therapy

Paul Morel<fnref n="1" />, Xiaodong Wu²,³, Guillaume Blin<fnref n="1" />,⁴, Stéphane Vialette<fnref n="1" />, Ryan Flynn², Daniel Hyer² and Dongxu Wang²*

**Purpose:** This study describes a real-time spot weight adaptation method in spot-scanning proton therapy for moving target or moving patient, so that the resultant dose distribution closely matches the planned dose distribution.

**Materials and methods:** The method proposed in this study adapts the weight (MU) of the delivering pencil beam to that of the target spot; it will actually hit during patient/target motion. The target spot that a certain delivering pencil beam may hit relies on patient monitoring and/or motion modeling using four-dimensional (4D) CT. After the adapted delivery, the required total weight [Monitor Unit (MU)] for this target spot is then subtracted from the planned value. With continuous patient motion and continuous spot scanning, the planned doses to all target spots will eventually be all fulfilled. In a proof-of-principle test, a lung case was presented with realistic temporal and motion parameters; the resultant dose distribution using spot weight adaptation was compared to that without using this method. The impact of the real-time patient/target position tracking or prediction was also investigated.

**Results:** For moderate motion (i.e., mean amplitude 0.5 cm), D95% to the planning target volume (PTV) was only 81.5% of the prescription ($R_X$) dose; with spot weight adaptation PTV D95% achieves 97.7% $R_X$. For large motion amplitude (i.e., 1.5 cm), without spot weight adaptation PTV D95% is only 42.9% of $R_X$; with spot weight adaptation, PTV D95% achieves 97.7% $R_X$. Larger errors in patient/target position tracking or prediction led to worse final target coverage; an error of 3 mm or smaller in patient/target position tracking is preferred.

**Conclusion:** The proposed spot weight adaptation method was able to deliver the planned dose distribution and maintain target coverage when patient motion was involved. The successful implementation of this method would rely on accurate monitoring or prediction of patient/target motion.

**Keywords:** intensity-modulated proton therapy, interplay effect, spot scanning, spot weight adaptation, motion

**OPEN ACCESS**

**Edited by:**

Yi Rong,
The Ohio State University, USA

**Reviewed by:**

Ji-Hong Hong,
Chang Gung Memorial Hospital,
Taiwan
Ahmet Sedat Ayan,
The Ohio State University, USA

**\*Correspondence:**

Dongxu Wang,
Department of Radiation Oncology,
The University of Iowa, 200 Hawkins
Drive, Iowa City, IA 52242, USA
dongxu-wang@uiowa.edu

**Specialty section:**

This article was submitted to
Radiation Oncology, a section of the
journal Frontiers in Oncology

**Received:** 11 March 2015

**Accepted:** 13 May 2015

**Published:** 28 May 2015

**Citation:**

Morel P, Wu X, Blin G, Vialette S, Flynn R, Hyer D and Wang D (2015) Spot weight adaptation for moving target in spot scanning proton therapy.
Front. Oncol. 5:119.
doi: 10.3389/fonc.2015.00119

<footnote n="1">¹Laboratoire Informatique Gaspard Monge (LIGM), UMR CNRS 8049, Université Paris-Est, Paris, France, ²Department of Radiation Oncology, The University of Iowa, Iowa City, IA, USA, ³Department of Electrical and Computer Engineering, The University of Iowa, Iowa City, IA, USA, ⁴UMR 5800, Laboratoire Bordelais de Recherche en Informatique (LaBRI), Université de Bordeaux, Talence, France</footnote>
Introduction

Spot scanning proton therapy has been shown to be a dosimetrically superior radiation delivery modality in cancer treatment (1). The ability to accurately place proton beam's Bragg peak inside a target is crucial to the realization of proton beam's dosimetric benefits. Yet multiple sources of uncertainty exist in the process of treatment planning and treatment delivery (2–4), weakening the potential advantage of spot scanning proton therapy. Among these uncertainties, intra-fraction motion, especially the respiratory motion, is of particular concern (5), due to its less predictable nature.

Motion management and motion mitigation strategies have been proposed for spot scanning proton and particle beam therapy. One of the most common techniques, as used in photon radiation as well, is the beam gating (6, 7), which consists of irradiating the patient only when the patient reaches a certain phase of the breathing cycle. During gating, breathing can also be regulated through coaching by audio or visual feedback in order to limit the cycle's variability (6). Motion may also be suppressed by breathhold or active breathing control techniques (8). Functional apnea or mechanical ventilation may be used under anesthesia (9). These techniques are promising solutions but tend to be constraining for the patient.

Other strategies not constraining the patient have been investigated. Repainting (also called re-scanning) (10) consists of delivering the planned dose into several cycles instead of a single pass. The goal is to obtain a high scanning speed by reducing the time spent at each spot position and reduce the risk of large localized dose discrepancies. While this technique is encouraging, it does not take into account the patient motions during the delivery. Four-dimensional (4D) treatment planning (11) using acquired patient motion information to aid treatment planning is also proposed. Beam tracking method (12, 13) calculates the target displacement and uses this information to adapt the beam position and beam energy to aim the planned target point in the patient. This method is particularly promising; however, it requires the adjustment of both the position and the energy of the delivering pencil beam. A number of reviews and reports (5, 14, 15) provide good summaries of the existing investigative or implemented methods.

In this work, we propose an alternative technique that only adapts the pencil beam (spot) weights to account for the patient intra-fraction motion. Like the beam tracking method, this method too relies on the real-time information of target spot location, but it allows the scanning magnets to change the beam position and energy following a planned and quasi-continuous trajectory, instead of constantly chasing the target spot. The response time in this method is, therefore, potentially shorter; especially considering that energy changes in beam tracking method may be the bottleneck of effectively tracking a target.

Materials and Methods

The Principle of Spot Weight Adaptation

The overall principle of the spot weight adaptation strategy is introduced here. At a given moment, the weight of a pencil beam

at its current energy and direction is adapted to the weight planned for the target spot that is under the irradiation of the current beam at this moment. The weight can be obtained from the patient and tumor displacement at voxel levels, as well as the weight of the pencil beam that was originally planned for this target spot.

More formally, let us consider a treatment plan, for which the beam of one energy layer is represented by the points ($b_i$, $w_i$), where $i$ indicates the beam index, $b_i$ the beam position and $w_i$ the corresponding planned spot weight. This plan induces the irradiation of a set of respective target spot locations ($P_i$) in the static patient, located at the depth of the Bragg peak. For the purpose of implementing spot weight adaptation only, we define ($P_i$) as irradiated by the delivering spot ($b_i$, $w_i$). The total dose to a certain target spot ($P_i$) will include contributions from other delivering spots and other energy layers. In a treatment plan generated for a static patient, the delivering spots ($b_i$, $w_i$) and target spots ($P_i$) have a one-to-one relationship, and therefore often not considered separately (Figure 1).

During the treatment delivery for a static patient, the proton machine changes energy and spot position, and deliver the planned set of pencil beams ($b_i$, $w_i$) one spot at a time. If the target spots are now moving relative to the delivering spots, the delivering spot may no longer irradiate the intended target spot. However, the planned dose distribution may be achieved if the target spot $P_i$ receives the planned beam weight, $w_i$, in terms of Monitor Unit (MU), even if it is from a different delivering spot $b_j$. A simple illustration of this method is shown in Figure 2. The motion is considered in 2D, in the plane orthogonal to the beam direction, as a displacement vector $v(t)$. As an example, at time $t=2$, the beam will be at position $b_2$ and will deliver a weight $w_2$, following the sequence of control points in the treatment control system as defined by the plan. However, at this moment another target point $P_5$ would be irradiated instead of the planned target spot $P_2$, and $P_5$ would receive an irradiation corresponding to $w_2$ instead of its planned value $w_5$. The spot weight adaptation technique relies on the adaptation of the beam weight for delivering spot $b_2$, which is set to $w_5$ instead of $w_2$ (see Figure 3).

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /><figcaption>**FIGURE 2** | Illustration of the delivery for a moving patient without spot weight adaptation. Delivering spot positions are represented in green. Pencil beams are represented in red. Target spots in the patient are represented in blue. The 2D displacement vector is represented in orange.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**FIGURE 3** | Illustration of the delivery for a moving patient with the spot weight adaptation; the weight now is w5. Delivering spot positions are represented in green. Pencil beams are represented in red. Target spots in the patient are represented in blue. The 2D displacement vector is represented in orange.</figcaption></figure>

Implementation of Spot Weight Adaptation

**Algorithm 1** represents the overall process of the delivery with adapted spot weights.

For each field with $l_{\text{pos}}$ is the list of all delivering spots within a certain energy layer, and the list $l_{\text{plan}}$ is built with all the delivering spots with non-zero spots. The pencil beam moves at the planned positions in $l_{\text{plan}}$ in order. At time $t$ the beam is at position $b_i$, and its weight is adapted to value $w_k$ to reflect the new target spot $b_k$. MUs corresponding to entire or part of $w_k$ may be delivered by

**ALGORITHM 1** | **Algorithm of the compensated delivery.** TP is the treatment plan provided by the user. $l_{\text{pos}}$ is the list of all delivering spots within a certain energy layer, including those with zero weights based on the plan. $l_{\text{plan}}$ is the list of the delivering spots with non-zero spots.

begin

foreach Gantry angle do

foreach Energy level do

Build $l_{\text{pos}}$ the list of all the delivering spots;

//Note: if margin is used, include in $l_{\text{pos}}$ the margin beam positions with weight = 0 if margin is used. See description later about margin.

Build $l_{\text{plan}}$ the list of all the delivering spots with weights >0;

while $l_{\text{plan}}$ is not empty do

foreach Delivering spot $b_i$ in $l_{\text{pos}}$ do

Get the patient displacement vector $v(t)$ and target spot ($b_k, w_k$)

if available

//Obtained from the motion monitoring system and/or motion modeling; see Motion Information section below.

//$w_k$ being the weight associated to $b_k$ in treatment plan

Irradiate the target spot with ($b_i, w_k$) or part of $w_k$;

Update $l_{\text{plan}}$: remove $w_k$ or part of $w_k$ from $b_k$;

return;

<figure><img src="image_5.png" /><figcaption>**FIGURE 4** | Illustration of the delivery for a moving patient with spot weight adaptation where a margin is added around the planned delivering spots. Planned delivering spots are represented in green. Beams are represented in red. Target spots in the patient are represented in blue. The patient displacement vector is represented in orange.</figcaption></figure>

$b_i$, and $w_k$ is updated to become its remaining weight or zero for $b_k$. The beam then moves to the next position. To fully cover the target spots, Algorithm 1 needs to be repeated multiple times, i.e., re-scanning is necessary; otherwise the planned target spots may not be fully delivered. The number of re-scans can be reduced by adding a margin around the planned spot positions (see **Figure 4**). The size of the margin can be adjusted by user to its optimal value, which may depend on the extent and speed of patient motion.

Motion Information

Implementing the spot weight adaptation algorithm requires knowing or predicting the patient and target position $P_i(t)$ at each moment of beam delivery. This is the same process as what is performed in the beam tracking method (12). In practice, this can
<figure><img src="image_6.png" /><figcaption>FIGURE 5 | The SBRT lung testing case with the delineation of the tumor and the organ-at-risk.</figcaption></figure>

<table><caption>TABLE 1 | Motion and monitoring settings in different tests for the testing case.</caption><thead><tr><th rowspan="2">Test No.</th><th rowspan="2">Intra-fraction motion</th><th colspan="4">Motion</th><th rowspan="2">Spot weight adaptation?</th><th rowspan="2">Motion monitoring error [σ<sub>x</sub>, σ<sub>y</sub>, σ<sub>z</sub>] in cm</th></tr><tr><th>Mean amplitude [x, y, z] in cm</th><th>Amplitude Gaussian noise [σ<sub>x</sub>, σ<sub>y</sub>, σ<sub>z</sub>] in cm</th><th>Mean cycle period [x, y, z] in s</th><th>Cycle period Gaussian noise [σ<sub>x</sub>, σ<sub>y</sub>, σ<sub>z</sub>] in s</th></tr></thead><tbody><tr><td>1</td><td>No</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>2</td><td>Yes</td><td>[0.5, 0.5, 0]</td><td>[0.1, 0.1, 0]</td><td>[4, 4, 0]</td><td>[0.8, 0.8, 0]</td><td>No</td><td>—</td></tr><tr><td>3</td><td>Yes</td><td>[1.5, 1.5, 0]</td><td>[0.3, 0.3, 0]</td><td>[4, 4, 0]</td><td>[0.8, 0.8, 0]</td><td>No</td><td>—</td></tr><tr><td>4</td><td>Yes</td><td>[0.5, 0.5, 0]</td><td>[0.1, 0.1, 0]</td><td>[4, 4, 0]</td><td>[0.8, 0.8, 0]</td><td>Yes</td><td>[0, 0, 0]</td></tr><tr><td>5</td><td>Yes</td><td>[1.5, 1.5, 0]</td><td>[0.3, 0.3, 0]</td><td>[4, 4, 0]</td><td>[0.8, 0.8, 0]</td><td>Yes</td><td>[0, 0, 0]</td></tr><tr><td>6</td><td>Yes</td><td>[1.5, 1.5, 0]</td><td>[0.3, 0.3, 0]</td><td>[4, 4, 0]</td><td>[0.8, 0.8, 0]</td><td>Yes</td><td>[0.3, 0.3, 0]</td></tr><tr><td>7</td><td>Yes</td><td>[1.5, 1.5, 0]</td><td>[0.3, 0.3, 0]</td><td>[4, 4, 0]</td><td>[0.8, 0.8, 0]</td><td>Yes</td><td>[0.75, 0.75, 0]</td></tr><tr><td>8</td><td>Yes</td><td>[1.5, 1.5, 0]</td><td>[0.3, 0.3, 0]</td><td>[4, 4, 0]</td><td>[0.8, 0.8, 0]</td><td>Yes</td><td>[1.5, 1.5, 0]</td></tr></tbody></table>

It is assumed that patient moves in sinusoidal patterns in X and Y axes independently (X-Y plane is perpendicular to the posterior beam's direction), with Gaussian noise existing for mean amplitudes and cycle periods. Noise is assumed to be 20% of the mean value (16). It is further assumed that a real-time patient monitoring system exists, but with Gaussian error in its detected signals. Test 1 is the static plan; Test 2 is the test with moderate motion without spot weight adaptation, and Test 4 is the test with moderate motion with spot weight adaptation. Test 3 is the test with large motion but no spot weight adaptation; Tests 5-8 are tests with large motion and with spot weight adaptation but different levels of error in motion monitoring.

<figure><img src="image_7.png" /><figcaption>FIGURE 6 | Dose distributions for testing case in static situation [[A), Test 1], in moderate moving situation without spot weight adaptation [[B), Test 2], and in moderate moving situation with spot weight adaptation [[C), Test 4]. Dose color map is relative to max dose in the static situation.</figcaption></figure>

be achieved with reasonable accuracy in two ways: (1) by modeling patient motion through 4DCT acquired before treatment planning and then relating patient motion to respiratory phase observables during treatment; (2) or by monitoring patient motion during treatment in real time through devices such as Calypso® system (Varian Medical Systems, Palo Alto, CA, USA). The spot weight adaptation method itself does not rely on any specific method to supply patient motion information. In this work for the

purpose of proof-of-principle testing, a motion model is created for the testing case patient, and a certain patient monitoring system is assumed to monitor patient motion in real-time, with realistic time delay and measurement noise included.

## Motion Simulation and Dose Computation

Motion simulation and dose computation were performed in in-house software named Motion Simulator for Proton Therapy
(MSPT) (http://code.google.com/p/mspt/). MSPT is an open-source treatment planning software that enables 4D dose computation. MSPT can import patient CT, structure, and spot scanning treatment plan in DICOM RT format. It can simulate user-defined motion pattern on the patient dataset on a voxel-by-voxel basis, and compute delivered dose in each temporal moment considering the temporal model for the beam delivery system, using analytical dose computation algorithm. It then sums the dose for each patient voxel along the time axis and produces the total dose distribution for the patient. A detailed description about MSPT and its validation has been published elsewhere (Morel et al. under review). The weight adaptation algorithm was implemented in MSPT for the test casing described below.

## A Testing Case

To evaluate the compensation method, we perform motion simulation and dose computation on a SBRT lung case shown in **Figure 5**. Motions of moderate amplitude (0.5 cm in any dimension) and large amplitude (1.5 cm in any dimension) are simulated respectively. Different levels of error for motion monitoring are also introduced. The simulation settings are summarized in **Table 1**.

<figure><img src="image_8.png" /></figure>

Dose distribution color map and dose volume histogram (DVH) are used in evaluating the delivered dose distributions in different tests. To evaluate the quality of the delivered dose distribution, we consider the following criterion to evaluate whether the desired tumor coverage is maintained with patient motion: if 95% of the tumor volume receives 95% of the prescribed dose $R_x$ ($D_{95\%} \ge 95\% R_x$) in one fraction, then we consider target coverage is maintained in this fraction. In our situation the dose distribution used as a reference is the planned static delivery. Therefore, we will assume that the dose prescription $R_x$ is the tumor $D_{95\%}$ in the planned static dose distribution.

Patient imaging data used in this study are from consented patients and approved by the Institutional Review Board (IRB) at the University of Iowa (IRB #201211805).

## Results

### Spot Weight Adaptation for Moderate Amplitude Motions

The effect of spot weight adaptation for a moderate motion amplitude (i.e., 0.5 cm) can be seen from the comparison between Tests 1, 2, and 4 (see **Table 1** for test settings). **Figure 6** shows the dose distribution color map of the three tests obtained near the center of the tumor. **Figure 7** shows the DVH for the three tests.

From **Figure 6**, we can observe that with moderate motion (i.e., mean amplitude 0.5 cm) the dose distributions of the static and the spot-weight-adapted delivery are close, whereas dose heterogeneities are visible in the tumor in the case of the non-adapted delivery (Test 2, **Figure 6B**). Based on the DVH in **Figure 7**, PTV coverage is degraded in Test 2 without spot weight adaptation, where PTV $D_{95\%}$ is only 81.5% of $R_x$; but with spot weight adaptation, PTV $D_{95\%}$ achieves 97.7%$R_x$. The hot spot in the non-adapted Test 2 is also higher than the planned value. OAR dose in the non-adapted Test 2 is also higher, though not significantly worse than planned value.

### Spot Weight Adaptation for Large Amplitude Motions

The effect of the spot weight adapted delivery for a large motion amplitude (i.e., 1.5 cm) can be seen from the comparison between Test 1, 3, and 5 (see **Table 1** for the settings). **Figure 8** shows the

<figure><img src="image_9.png" /></figure>
dose distribution color map of the three tests obtained near the center of the tumor. **Figure 9** shows the DVH for the three tests.

From **Figure 8**, we can observe that with large motion (i.e., mean amplitude 1.5 cm) the dose distributions of the static and the spot-weight-adapted delivery are close, whereas dose heterogeneities are very obvious in the tumor in the case of the non-adapted delivery (Test 3, **Figure 8B**). Based on the DVH in **Figure 9**, PTV coverage is degraded in Test 3 without spot weight adaptation, where PTV D<sub>95%</sub> is only 42.9% of R<sub>X</sub>; but with spot weight adaptation, PTV D<sub>95%</sub> achieves 97.7%R<sub>X</sub>. With spot weight adaptation, the delivered dose to OARs significantly deteriorated from the planned values; while with spot weight adaptation,

delivered dose is very close to planned values, as shown from the DVH in **Figure 9**.

## The Impact of the Monitoring System Accuracy

In the evaluations for spot weight adaptation with moderate and large motion (Test 2–5), an ideal patient monitoring system is assumed. In reality, patient motion information as an input to the spot weight adaptation algorithm should always contain errors. While both spatial and temporal errors can be present, only the impact of spatial errors is studied here, because the temporal error in detecting a current position of target spot and the delay of beam delivering can still be reflected in the error of target spot’s spatial location. Gaussian errors of three different levels of standard deviations were added: Tests 5 ($\sigma$ = 0 cm), 6 ($\sigma$ = 0.3 cm), 7 ($\sigma$ = 0.75 cm), and 8 ($\sigma$ = 1.5 cm). **Table 1** has the summary of the settings for these tests. **Figure 10** shows the dose distributions for Tests 5–8, and **Figure 11** shows the DVH for Tests 5–8.

From both **Figures 10** and **11**, we can see that the accuracy of the monitoring system has an impact of the quality of the final dose distribution. PTV D<sub>95%</sub> for Gaussian error of 0.3, 0.75, and 1.5 cm are 91.2% R<sub>X</sub>, 73.1% R<sub>X</sub>, and 55.7% R<sub>X</sub>, respectively.

## Discussion and Conclusion

When moderate (mean amplitude 0.5 cm) or large motion (mean amplitude 1.5 cm) is present in the testing case, delivering the plan without any adaptation deteriorated PTV coverage (Test 2 and Test 3, with PTV D<sub>95%</sub> at only 81.5% R<sub>X</sub> and 42.9% R<sub>X</sub>, respectively). This study proposed a method to adapt spot weight during delivery, with patient motion information as input. With spot weight adaptation, PTV coverage was able to achieve 97.7%R<sub>X</sub> for

<figure><img src="image_10.png" /><figcaption>**FIGURE 9** | Dose volume histogram (DVH) for testing case in static situation (Test 1), in large moving situation without spot weight adaptation (Test 3), and in moving situation with spot weight adaptation (Test 5). 95% of R<sub>X</sub> dose and 95% of planning target volume (PTV) are marked.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**FIGURE 10** | Dose distributions for spot-weighted-adapted delivery with large motion, with different patient monitoring errors: $\sigma$ = 0 cm (A), 0.3 cm (B), 0.75 cm (C), 1.5 cm (D). Dose color map is relative to max dose in the static situation.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>FIGURE 11 | Dose volume histogram for spot-weighted-adapted delivery with large motion, with different levels of patient monitoring errors.</figcaption></figure>

moderate and large motion (Test 3 and Test 5), which is clinically acceptable.

The proof-of-principle example in this study did not include change in radiological depth, which would require spot weight adaptation across energy layers. The spot weight adaptation method in principle can adapt weight across energy layers, i.e., the method will look for the correct weight if the energy layer is different. Similar to the lateral spot margins as shown in **Figure 4**,

## References

1. Lomax AJ, Böhringer T, Bolsi A, Coray D, Emert F, Goitein G, et al. Treatment planning and verification of proton therapy using spot scanning: initial experiences. *Med Phys* (2004) **31**:3150. doi:10.1118/1.1779371

2. Lomax AJ. Intensity modulated proton therapy and its sensitivity to treatment uncertainties 2: the potential effects of inter-fraction and inter-field motions. *Phys Med Biol* (2008) **53**:1043–56. doi:10.1088/0031-9155/53/4/015

3. Lomax AJ. Intensity modulated proton therapy and its sensitivity to treatment uncertainties 1: the potential effects of calculational uncertainties. *Phys Med Biol* (2008) **53**:1027–42. doi:10.1088/0031-9155/53/4/014

4. Paganetti H. Range uncertainties in proton therapy and the role of Monte Carlo simulations. *Phys Med Biol* (2012) **57**:R99. doi:10.1088/0031-9155/57/11/R99

5. Bert C, Durante M. Motion in radiotherapy: particle therapy. *Phys Med Biol* (2011) **56**:R113. doi:10.1088/0031-9155/56/16/R01

6. Keall PJ, Mageras GS, Balter JM, Emery RS, Forster KM, Jiang SB, et al. The management of respiratory motion in radiation oncology report of AAPM Task Group 76a. *Med Phys* (2006) **33**:3874–900. doi:10.1118/1.2240847

7. Minohara S, Kanai T, Endo M, Noda K, Kanazawa M. Respiratory gated irradiation system for heavy-ion radiotherapy. *Int J Radiat Oncol Biol Phys* (2000) **47**:1097–103. doi:10.1016/S0360-3016(00)00524-1

8. Wong JW, Sharpe MB, Jaffray DA, Kini VR, Robertson JM, Stromberg JS, et al. The use of active breathing control (ABC) to reduce margin for breathing motion. *Int J Radiat Oncol Biol Phys* (1999) **44**:911–9. doi:10.1016/S0360-3016(99)00129-7

9. Engelsman M, Bert C. Precision and uncertainties in proton therapy for moving targets. In: Paganetti H, editor. *Proton Therapy Physics. Series in Medical Physics and Biomedical Engineering*. Vol. 1. Boca Raton, FL: CRC Press (2011). p. 435–60.

10. Zenklusen S, Pedroni E, Meer D. A study on repainting strategies for treating moderately moving targets with proton pencil beam scanning at the new Gantry 2 at PSI. *Phys Med Biol* (2010) **55**:5103. doi:10.1088/0031-9155/55/17/014

an energy margin with spots of zero initial weight can be added. The zero-weight spots are to allow the delivering beam to move to that energy/place first, and then choose its MU based on the currently target spot. Re-scanning will be necessary in order to cover all spots and all energy layers. The number of re-scanning is likely smaller when the margin is larger. Since energy layer switch is usually in the order of seconds for spot scanning proton therapy, the number of re-scanning is important in determining the total treatment time. The relation between the number of re-scans and the size of the energy and lateral margins is currently being investigated and will be reported in future.

While the results also show the impact of the patient motion modeling and monitoring accuracy in spot weight adaptation method, just as in beam tracking method (12), this study shows that the spot weight adaptation itself is a promising strategy in mitigating the impact of patient motion during spot scanning proton therapy.

## Acknowledgments

We wish to thank Christina Zacharatou from the Radiotherapy Group of Institut Bergonié (Bordeaux, France) for her valuable advice on proton therapy. This work is partially supported by the French National Research Agency (ANR) project BIRDS JCJC SIMI 2-2010.

11. Richter D, Schwarzkopf A, Trautmann J, Krämer M, Durante M, Jäkel O, et al. Upgrade and benchmarking of a 4D treatment planning system for scanned ion beam therapy. *Med Phys* (2013) **40**:051722. doi:10.1118/1.4800802

12. Bert C, Saito N, Schmidt A, Chaudhri N, Schardt D, Rietzel E. Target motion tracking with a scanned particle beam. *Med Phys* (2007) **34**:4768–71. doi:10.1118/1.2815934

13. Seregni M, Kaderka R, Fattori G, Riboldi M, Pella A, Constantinescu A, et al. Tumor tracking based on correlation models in scanned ion beam therapy: an experimental study. *Phys Med Biol* (2013) **58**:4659. doi:10.1088/0031-9155/58/13/4659

14. Knopf A, Bert C, Heath E, Nill S, Kraus K, Richter D, et al. Special report: workshop on 4D-treatment planning in actively scanned particle therapy – recommendations, technical challenges, and future research directions. *Med Phys* (2010) **37**:4608–14. doi:10.1118/1.3475944

15. Bert C, Graeff C, Riboldi M, Nill S, Baroni G, Knopf A-C. Advances in 4D treatment planning for scanned particle beam therapy – report of dedicated workshops. *Technol Cancer Res Treat* (2014) **13**:485–95. doi:10.7785/tcrtexpress.2013.600274

16. Heath E, Unkelbach J, Oelfke U. Incorporating uncertainties in respiratory motion into 4D treatment plan optimization. *Med Phys* (2009) **36**:3059–71. doi:10.1118/1.3181908

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2015 Morel, Wu, Blin, Viallette, Flynn, Hyer and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.