On-line range verification for proton beam therapy using spherical ionoacoustic waves with resonant frequency

In contrast to conventional X-ray therapy, proton beam therapy (PBT) can confine radiation doses to tumours because of the presence of the Bragg peak. However, the precision of the treatment is currently limited by the uncertainty in the beam range. Recently, a unique range verification methodology has been proposed based on simulation studies that exploit spherical ionoacoustic waves with resonant frequency (SPIREs). SPIREs are emitted from spherical gold markers in tumours initially introduced for accurate patient positioning when the proton beam is injected. These waves have a remarkable property: their amplitude is linearly correlated with the residual beam range at the marker position. Here, we present proof-of-principle experiments using short-pulsed proton beams at the clinical dose to demonstrate the feasibility of using SPIREs for beam-range verification with submillimetre accuracy. These results should substantially contribute to reducing the range uncertainty in future PBT applications.

A novel range-verification method using ionoacoustic wave generated from spherical gold markers for particle-beam therapy: a simulation study

This study proposes a novel alternative range-verification method for proton beam with acoustic waves generated from spherical metal markers. When proton beam is incident on metal markers, most of the resulting pressure waves are confined in the markers because of the large difference in acoustic impedance between the metal and tissue. However, acoustic waves with frequency equal to marker’s resonant frequency escape this confinement; the marker briefly acts as an acoustic transmitter. Herein, this phenomenon is exploited to measure the range of the proton beam. We test the proposed strategy in 3-D simulations, combining the dose calculations with modelling of acoustic-wave propagation. A spherical gold marker of 2.0 mm diameter was placed in water with a 60 MeV proton beam incident on it. We investigated the dependence of pressure waves on the width of beam pulse and marker position. At short beam pulse, specific high-frequency acoustic waves of 1.62 MHz originating from the marker were observed in wave simulations, whose amplitude correlated with the distance between the marker and Bragg peak. Results indicate that the Bragg peak position can be estimated by measuring the acoustic wave amplitudes from the marker, using a single detector properly designed for the resonance frequency.

An analysis of beam parameters on proton-acoustic waves through an analytic approach

It has been reported that acoustic waves are generated when a high-energy pulsed proton beam is deposited in a small volume within tissue. One possible application of proton-induced acoustics is to get real-time feedback for intra-treatment adjustments by monitoring such acoustic waves. A high spatial resolution in ultrasound imaging may reduce proton range uncertainty. Thus, it is crucial to understand the dependence of the acoustic waves on the proton beam characteristics. In this manuscript, firstly, an analytic solution for the proton-induced acoustic wave is presented to reveal the dependence of the signal on the beam parameters; then it is combined with an analytic approximation of the Bragg curve. The influence of the beam energy, pulse duration and beam diameter variation on the acoustic waveform are investigated. Further analysis is performed regarding the Fourier decomposition of the proton-acoustic signals. Our results show that the smaller spill time of the proton beam upsurges the amplitude of the acoustic wave for a constant number of protons, which is hence beneficial for dose monitoring. The increase in the energy of each individual proton in the beam leads to the spatial broadening of the Bragg curve, which also yields acoustic waves of greater amplitude. The pulse duration and the beam width of the proton beam do not affect the central frequency of the acoustic wave, but they change the amplitude of the spectral components.

Experimental observation of acoustic emissions generated by a pulsed proton beam from a hospital-based clinical cyclotron

PURPOSE

To measure the acoustic signal generated by a pulsed proton spill from a hospital-based clinical cyclotron.

METHODS

An electronic function generator modulated the IBA C230 isochronous cyclotron to create a pulsed proton beam. The acoustic emissions generated by the proton beam were measured in water using a hydrophone. The acoustic measurements were repeated with increasing proton current and increasing distance between detector and beam.

RESULTS

The cyclotron generated proton spills with rise times of 18 μs and a maximum measured instantaneous proton current of 790 nA. Acoustic emissions generated by the proton energy deposition were measured to be on the order of mPa. The origin of the acoustic wave was identified as the proton beam based on the correlation between acoustic emission arrival time and distance between the hydrophone and proton beam. The acoustic frequency spectrum peaked at 10 kHz, and the acoustic pressure amplitude increased monotonically with increasing proton current.

CONCLUSIONS

The authors report the first observation of acoustic emissions generated by a proton beam from a hospital-based clinical cyclotron. When modulated by an electronic function generator, the cyclotron is capable of creating proton spills with fast rise times (18 μs) and high instantaneous currents (790 nA). Measurements of the proton-generated acoustic emissions in a clinical setting may provide a method for in vivo proton range verification and patient monitoring.

Ionoacoustic tomography of the proton Bragg peak in combination with ultrasound and optoacoustic imaging

Ions provide a more advantageous dose distribution than photons for external beam radiotherapy, due to their so-called inverse depth dose deposition and, in particular a characteristic dose maximum at their end-of-range (Bragg peak). The favorable physical interaction properties enable selective treatment of tumors while sparing surrounding healthy tissue, but optimal clinical use requires accurate monitoring of Bragg peak positioning inside tissue. We introduce ionoacoustic tomography based on detection of ion induced ultrasound waves as a technique to provide feedback on the ion beam profile. We demonstrate for 20 MeV protons that ion range imaging is possible with submillimeter accuracy and can be combined with clinical ultrasound and optoacoustic tomography of similar precision. Our results indicate a simple and direct possibility to correlate, in-vivo and in real-time, the conventional ultrasound echo of the tumor region with ionoacoustic tomography. Combined with optoacoustic tomography it offers a well suited pre-clinical imaging system.

How proton pulse characteristics influence protoacoustic determination of proton-beam range: simulation studies

The unique dose deposition of proton beams generates a distinctive thermoacoustic (protoacoustic) signal, which can be used to calculate the proton range. To identify the expected protoacoustic amplitude, frequency, and arrival time for different proton pulse characteristics encountered at hospital-based proton sources, the protoacoustic pressure emissions generated by 150 MeV, pencil-beam proton pulses were simulated in a homogeneous water medium. Proton pulses with Gaussian widths ranging up to 200 μs were considered. The protoacoustic amplitude, frequency, and time-of-flight (TOF) range accuracy were assessed. For TOF calculations, the acoustic pulse arrival time was determined based on multiple features of the wave. Based on the simulations, Gaussian proton pulses can be categorized as Dirac-delta-function-like (FWHM < 4 μs) and longer. For the δ-function-like irradiation, the protoacoustic spectrum peaks at 44.5 kHz and the systematic error in determining the Bragg peak range is <2.6 mm. For longer proton pulses, the spectrum shifts to lower frequencies, and the range calculation systematic error increases (⩽ 23 mm for FWHM of 56 μs). By mapping the protoacoustic peak arrival time to range with simulations, the residual error can be reduced. Using a proton pulse with FWHM = 2 μs results in a maximum signal-to-noise ratio per total dose. Simulations predict that a 300 nA, 150 MeV, FWHM = 4 μs Gaussian proton pulse (8.0 × 10(6) protons, 3.1 cGy dose at the Bragg peak) will generate a 146 mPa pressure wave at 5 cm beyond the Bragg peak. There is an angle dependent systematic error in the protoacoustic TOF range calculations. Placing detectors along the proton beam axis and beyond the Bragg peak minimizes this error. For clinical proton beams, protoacoustic detectors should be sensitive to <400 kHz (for -20 dB). Hospital-based synchrocyclotrons and cyclotrons are promising sources of proton pulses for generating clinically measurable protoacoustic emissions.

Ionoacoustic characterization of the proton Bragg peak with submillimeter accuracy

PURPOSE

Range verification in ion beam therapy relies to date on nuclear imaging techniques which require complex and costly detector systems. A different approach is the detection of thermoacoustic signals that are generated due to localized energy loss of ion beams in tissue (ionoacoustics). Aim of this work was to study experimentally the achievable position resolution of ionoacoustics under idealized conditions using high frequency ultrasonic transducers and a specifically selected probing beam.

METHODS

A water phantom was irradiated by a pulsed 20 MeV proton beam with varying pulse intensity and length. The acoustic signal of single proton pulses was measured by different PZT-based ultrasound detectors (3.5 and 10 MHz central frequencies). The proton dose distribution in water was calculated by Geant4 and used as input for simulation of the generated acoustic wave by the matlab toolbox k-WAVE.

RESULTS

In measurements from this study, a clear signal of the Bragg peak was observed for an energy deposition as low as 10(12) eV. The signal amplitude showed a linear increase with particle number per pulse and thus, dose. Bragg peak position measurements were reproducible within ±30 μm and agreed with Geant4 simulations to better than 100 μm. The ionoacoustic signal pattern allowed for a detailed analysis of the Bragg peak and could be well reproduced by k-WAVE simulations.

CONCLUSIONS

The authors have studied the ionoacoustic signal of the Bragg peak in experiments using a 20 MeV proton beam with its correspondingly localized energy deposition, demonstrating submillimeter position resolution and providing a deep insight in the correlation between the acoustic signal and Bragg peak shape. These results, together with earlier experiments and new simulations (including the results in this study) at higher energies, suggest ionoacoustics as a technique for range verification in particle therapy at locations, where the tumor can be localized by ultrasound imaging. This acoustic range verification approach could offer the possibility of combining anatomical ultrasound and Bragg peak imaging, but further studies are required for translation of these findings to clinical application.

Review of ultrasound image guidance in external beam radiotherapy: I. Treatment planning and inter-fraction motion management

In modern radiotherapy, verification of the treatment to ensure the target receives the prescribed dose and normal tissues are optimally spared has become essential. Several forms of image guidance are available for this purpose. The most commonly used forms of image guidance are based on kilovolt or megavolt x-ray imaging. Image guidance can also be performed with non-harmful ultrasound (US) waves. This increasingly used technique has the potential to offer both anatomical and functional information.This review presents an overview of the historical and current use of two-dimensional and three-dimensional US imaging for treatment verification in radiotherapy. The US technology and the implementation in the radiotherapy workflow are described. The use of US guidance in the treatment planning process is discussed. The role of US technology in inter-fraction motion monitoring and management is explained, and clinical studies of applications in areas such as the pelvis, abdomen and breast are reviewed. A companion review paper (O'Shea et al 2015 Phys. Med. Biol. submitted) will extensively discuss the use of US imaging for intra-fraction motion quantification and novel applications of US technology to RT.

Proton beam characterization by proton-induced acoustic emission: simulation studies

Due to their Bragg peak, proton beams are capable of delivering a targeted dose of radiation to a narrow volume, but range uncertainties currently limit their accuracy. One promising beam characterization technique, protoacoustic range verification, measures the acoustic emission generated by the proton beam. We simulated the pressure waves generated by proton radiation passing through water. We observed that the proton-induced acoustic signal consists of two peaks, labeled α and γ, with two originating sources. The α acoustic peak is generated by the pre-Bragg peak heated region whereas the source of the γ acoustic peak is the proton Bragg peak. The arrival time of the α and γ peaks at a transducer reveals the distance from the beam propagation axis and Bragg peak center, respectively. The maximum pressure is not observed directly above the Bragg peak due to interference of the acoustic signals. Range verification based on the arrival times is shown to be more effective than determining the Bragg peak position based on pressure amplitudes. The temporal width of the α and γ peaks are linearly proportional to the beam diameter and Bragg peak width, respectively. The temporal separation between compression and rarefaction peaks is proportional to the spill time width. The pressure wave expected from a spread out Bragg peak dose is characterized. The simulations also show that acoustic monitoring can verify the proton beam dose distribution and range by characterizing the Bragg peak position to within ~1 mm.