Lateral Spot Research Articles

TEAM: Triangular-mEsh Adaptive and Multiscale proton spot generation method.

Proton therapy is preferred for its dose conformality to spare normal tissues and organs-at-risk (OAR) via Bragg peaks with negligible exit dose. However, proton dose conformality can be further optimized: (1) the spot placement is based on the structured (e.g., Cartesian) grid, which may not offer conformal shaping to complex tumor targets; (2) the spot sampling pattern is uniform, which may be insufficient at the tumor boundary to provide the sharp dose falloff, and at the same time may be redundant at the tumor interior to provide the uniform dose coverage, for example, due to multiple Coulomb scattering (MCS); and (3) the lateral spot penumbra increases with respect to the depth due to MCS, which blurs the lateral dose falloff. On the other hand, while (1) the deliverable spots are subject to the minimum-monitor-unit (MMU) constraint, and (2) the dose rate is proportional to the MMU threshold, the current spot sampling method is sensitive to the MMU threshold and can fail to provide satisfactory plan quality for a large MMU threshold (i.e., high-dose-rate delivery). This work will develop a novel Triangular-mEsh-based Adaptive and Multiscale (TEAM) proton spot generation method to address these issues for optimizing proton dose conformality and plan delivery efficiency. Compared to the standard clinically-used spot placement method, three key elements of TEAM are as follows: (1) a triangular mesh instead of a structured grid: the triangular mesh is geometrically more conformal to complex target shapes and therefore more efficient and accurate for dose shaping inside and around the target; (2) adaptive sampling instead of uniform sampling: the adaptive sampling consists of relatively dense sampling at the tumor boundary to create the sharp dose falloff, which is more accurate, and coarse sampling at the tumor interior to uniformly cover the target, which is more efficient; and (3) depth-dependent sampling instead of depth-independent sampling: the depth-dependent sampling is used to compensate for MCS, that is, with increasingly dense sampling at the tumor boundary to improve dose shaping accuracy, and increasingly coarse sampling at the tumor interior to improve dose shaping efficiency, as the depth increases. In the TEAM method the spot locations are generated for each energy layer and layer-by-layer in the multiscale fashion; and then the spot weights are derived by solving the IMPT problem of dose-volume planning objectives, MMU constraints, and robustness optimization with respect to range and setup uncertainties. Compared to the standard clinically-used spot placement method UNIFORM, TEAM achieved (1) better plan quality using <60% number of spots of UNIFORM; (2) better robustness to the number of spots; (3) better robustness to a large MMU threshold. Furthermore, TEAM provided better plan quality with fewer spots than other adaptive methods (Cartesian-grid or triangular-mesh). A novel triangular-mesh-based proton spot placement method called TEAM is proposed, and it is demonstrated to improve plan quality, robustness to the number of spots, and robustness to the MMU threshold, compared to the clinically-used spot placement method and other adaptive methods.

Mixed-size spot scanning with a compact large momentum acceptance superconducting (LMA-SC) gantry beamline for proton therapy

Objective. Lowering treatment costs and improving treatment quality are two primary goals for next-generation proton therapy (PT) facilities. This work will design a compact large momentum acceptance superconducting (LMA-SC) gantry beamline to reduce the footprint and expense of the PT facilities, with a novel mixed-size spot scanning method to improve the sparing of organs at risk (OAR). Approach. For the LMA-SC gantry beamline, the movable energy slit is placed in the middle of the last achromatic bending section, and the beam momentum spread of delivered spots can be easily changed during the treatment. Simultaneously, changing the collimator size can provide spots with various lateral spot sizes. Based on the provided large-size and small-size spot models, the treatment planning with mixed spot scanning is optimized: the interior of the target is irradiated with large-size spots (to cover the uniform-dose interior efficiently), while the peripheral of the target is irradiated with small-size spots (to shape the sharp dose falloff at the peripheral accurately). Main results. The treatment plan with mixed-size spot scanning was evaluated and compared with small and large-size spot scanning for thirteen clinical prostate cases. The mixed-size spot plan had superior target dose homogeneities, better protection of OAR, and better plan robustness than the large-size spot plan. Compared to the small-size spot plan, the mixed-size spot plan had comparable plan quality, better plan robustness, and reduced plan delivery time from 65.9 to 40.0 s. Significance. The compact LMA-SC gantry beamline is proposed with mixed-size spot scanning, with demonstrated footprint reduction and improved plan quality compared to the conventional spot scanning method.

Extension of matRad with a modified microdosimetric kinetic model for carbon ion treatment planning: Comparison with Monte Carlo calculation.

Treatment planning is essential for in silico particle therapy studies. matRad is an open-source research treatment planning system (TPS) based on the local effect model, which is a type of relative biological effectiveness (RBE) model. This study aims to implement a microdosimetric kinetic model (MKM) in matRad and develop an automation algorithm for Monte Carlo (MC) dose recalculation using the TOPAS code. In addition, we provide the developed MKM extension as open-source tool for users. Carbon beam data were generated using TOPAS MC pencil beam irradiation. We parameterized the TOPAS MC beam data with a double-Gaussian fit and modeled the integral depth doses and lateral spot profiles in the range of 100-430 MeV/u. To implement the MKM, the specific energy data table for Z=1-6 and integrated depth-specific energy data were acquired based on the Kiefer-Chatterjee track structure and TOPAS MC simulation, respectively. Generic data were integrated into matRad, and treatment planning was performed based on these data. The optimized plan parameters were automatically converted into MC simulation input. Finally, the matRad TPS and TOPAS MC simulations were compared using the RBE-weighted dose calculation results. A comparison was made for three geometries: homogeneous water phantom, inhomogeneous phantom, and patient. The RBE-weighted dose (DRBE ) distribution agreed with TOPAS MC within 1.8% for all target sizes for the homogeneous phantom. For the inhomogeneous phantom, the relative difference in the range of 80% of the prescription dose in the distal fall-off region (R80) between the matRad TPS and TOPAS MC was 0.6% (1.1mm). DRBE between the TPS and the MC was within 4.0%. In the patient case, the difference in the dose-volume histogram parameters for the target volume between the TPS and the MC was less than 2.7%. The relative difference in R80 was 0.7% (1.2mm). The MKM was successfully implemented in matRad TPS, and the RBE-weighted dose was comparable to that of TOPAS MC. The MKM-implemented matRad was released as an open-source tool. Further investigations with MC simulations can be conducted using this tool, providing a good option for carbon ion research.

Intensity modulated proton arc therapy via geometry-based energy selection for ependymoma.

We developed and tested a novel method of creating intensity modulated proton arc therapy (IMPAT) plans that uses computing resources similar to those for regular intensity-modulated proton therapy (IMPT) plans and may offer a dosimetric benefit for patients with ependymoma or similar tumor geometries. Our IMPAT planning method consists of a geometry-based energy selection step with major scanning spot contributions as inputs computed using ray-tracing and single-Gaussian approximation of lateral spot profiles. Based on the geometric relation of scanning spots and dose voxels, our energy selection module selects a minimum set of energy layers at each gantry angle such that each target voxel is covered by sufficient scanning spots as specified by the planner, with dose contributions above the specified threshold. Finally, IMPAT plans are generated by robustly optimizing scanning spots of the selected energy layers using a commercial proton treatment planning system (TPS). The IMPAT plan quality was assessed for four ependymoma patients. Reference three-field IMPT plans were created with similar planning objective functions and compared with the IMPAT plans. In all plans, the prescribed dose covered 95% of the clinical target volume (CTV) while maintaining similar maximum doses for the brainstem. While IMPAT and IMPT achieved comparable plan robustness, the IMPAT plans achieved better homogeneity and conformity than the IMPT plans. The IMPAT plans also exhibited higher relative biological effectiveness (RBE) enhancement than did the corresponding reference IMPT plans for the CTV in all four patients and brainstem in three of them. The proposed method demonstrated potential as an efficient technique for IMPAT planning and may offer a dosimetric benefit for patients with ependymoma or tumors in close proximity to critical organs. IMPAT plans created using this method had elevated RBE enhancement associated with increased linear energy transfer (LET) in both targets and abutting critical organs.

Intrafraction tumor motion monitoring and dose reconstruction for liver pencil beam scanning proton therapy

Pencil beam scanning (PBS) proton therapy can provide highly conformal target dose distributions and healthy tissue sparing. However, proton therapy of hepatocellular carcinoma (HCC) is prone to dosimetrical uncertainties induced by respiratory motion. This study aims to develop intra-treatment tumor motion monitoring during respiratory gated proton therapy and combine it with motion-including dose reconstruction to estimate the delivered tumor doses for individual HCC treatment fractions. Three HCC-patients were planned to receive 58 GyRBE (n=2) or 67.5 GyRBE (n=1) of exhale respiratory gated PBS proton therapy in 15 fractions. The treatment planning was based on the exhale phase of a 4-dimensional CT scan. Daily setup was based on cone-beam CT (CBCT) imaging of three implanted fiducial markers. An external marker block (RPM) on the patient's abdomen was used for exhale gating in free breathing. This study was based on 5 fractions (patient 1), 1 fraction (patient 2) and 6 fractions (patient 3) where a post-treatment control CBCT was available. After treatment, segmented 2D marker positions in the post-treatment CBCT projections provided the estimated 3D motion trajectory during the CBCT by a probability-based method. An external-internal correlation model (ECM) that estimated the tumor motion from the RPM motion was built from the synchronized RPM signal and marker motion in the CBCT. The ECM was then used to estimate intra-treatment tumor motion. Finally, the motion-including CTV dose was estimated using a dose reconstruction method that emulates tumor motion in beam's eye view as lateral spot shifts and in-depth motion as changes in the proton beam energy. The CTV homogeneity index (HI) The CTV homogeneity index (HI) was calculated as . The tumor position during spot delivery had a root-mean-square error of 1.3mm in left-right, 2.8mm in cranio-caudal and 1.7mm in anterior-posterior directions compared to the planned position. On average, the CTV HI was larger than planned by 3.7%-points (range: 1.0-6.6%-points) for individual fractions and by 0.7%-points (range: 0.3-1.1%-points) for the average dose of 5 or 6 fractions. A method to estimate internal tumor motion and reconstruct the motion-including fraction dose for PBS proton therapy of HCC was developed and demonstrated successfully clinically.

The “True Colours” of Golden Loaches (Teleostei: Cobitidae)

Loaches of the genus Sabanejewia are model organisms for many ichthyologic studies focusing on morphological and colour variation. The bright lateral spots—a part of the mimicry of this fish—cause a taxonomic puzzle. Common practice is that morphometric and molecular studies lack habitat assessments of species originating from different environments. In this study, we analysed body shape and colour variation within two species of Sabanejewia in order to reveal whether they are more affected by phylogeny or environmental conditions. The environmental parameters were obtained using satellite mapping tools, body shape using geometric morphometry, and the molecular identification of specimens using an analysis of mitochondrial DNA (sequence of the cytochrome b gene). Generalized linear mixed models were used for the analysis of lateral spot numbers, and principal components analysis to describe the changes in landmark position defined by the permutation tests. Results show a significant impact of the locality elevation and specimen’s sex on the change of the body shape. We identified significant differences in the number of lateral spots, which are considered mimicry variation across an elevation gradient. Morphological differences of Sabanejewia spp. vary with ecological factors along the river continuum.

Ventral and lateral spot patterns differentiation between three smooth newt species (Amphibia: Salamandridae: Lissotriton)

The location, shape and area of dark spots on the belly, throat and sides of the body of three closely related species of smooth newts (Lissotriton kosswigi, L. schmidtleri and L. vulgaris) were analysed. Differences were found between species and sexes in the location and shape of dark spots on the belly, throat and lateral surfaces of the body. Discriminant analysis of these spot characteristics allows to identificate of males of all three species with a high degree of confidence (85–91%). In females, only L. vulgaris was correctly distinguished from the other two species (accuracy 81–94%). Anatolian and Thracian populations of L. schmidtleri have very similar patterns of dark spots, which confirm their conspecificity. The differences in the location and size of dark spots make possible successful identification of the species during field research and study of museum specimens. The developed method could be useful for comparative studies of other animal species which have spotted camouflage or aposematic colouration.

Ultrahigh-aspect-ratio beam generation with super-resolution spot

Focusing to a small spot along a large depth of focus (DOF) beyond the limits imposed by diffraction has attracted extensive theoretical as well as practical interest. The aspect ratio, defined as the DOF divided by the diameter of spot, has been extended by several schemes with reasonable success, but there is no report that enables one to modulate the aspect ratio as long as possible with sidelobe suppression. In this Letter, we propose a method to generate ultrahigh-aspect-ratio beams (UARBs) by introducing an analytic procedure using Zernike polynomials in the modified iterative Fourier-transform algorithm, and a super-resolution spot is also realized to further extend the aspect ratio. We design the diffractive optical elements that generate UARBs with sidelobe suppression and further experimentally demonstrate a stretched 1602:1 aspect ratio of UARBs with a lateral super-resolution spot, suggesting potential applications for further super-resolution imaging or precision fabrication.

Variable-focus optoacoustic lens with wide dynamic range and long focal length by using a flexible polymer nano-composite membrane

We demonstrate a variable-focus optoacoustic lens (VFOL) by pneumatically controlling a flexible polymer nano-composite membrane, which can produce laser-generated focused ultrasound (LGFU) with a high peak amplitude (>30 MPa) and a tight focal dimension (several hundred μm) over a wide dynamic range of focus variation (>20 mm) together with a long focal length up to 60 mm, each of which is widest and longest among optoacoustic lenses developed so far. Two different designs in lens dimension have been fabricated and characterized: VFOL-L with a 40-mm diameter and VFOL-S with 10 mm. VFOL-L exhibits a wide dynamic range of focal length variation from 38.52 to 60.39 mm with a center frequency near ~ 10 MHz, which is proper for practical long-range applications with several-cm depth. In comparison, VFOL-S covers a focal variation range from 6.75 to 11.1 mm with ~ 14 MHz, producing a relatively higher-pressure amplitude, which allows the inception of acoustic cavitation at an impedance-mismatched boundary. The nano-composite membrane of VFOL is actuated from a planar to deeply curved shape by externally injecting liquid into the VFOL, resulting in a focal gain up to 255 and a positive peak pressure of > 30 MPa in the VFOL-L case. The minimum-geometrical f-number as low as 0.963 is achieved at the shortest focal length (38.52 mm) with 6-dB lateral and axial spot dimensions of 304 μm and 2.86 mm, respectively. We expect that the proposed VFOL-based LGFU with a high peak pressure, a wide dynamic axial range, and a tight focal dimension are suitably applied for depth-dependent characterization of tissues and shockwave treatment, taking advantages of optoacoustic pulses as input with inherent broadband high-frequency characteristics.

An amyotrophic lateral sclerosis hot spot in the French Alps associated with genotoxic fungi

Between 1990 and 2018, 14 cases of amyotrophic lateral sclerosis (ALS) were diagnosed in residents of, and in visitors with second homes to, a mountainous hamlet in the French Alps. Systematic investigation revealed a socio-professional network that connected ALS cases. Genetic risk factors for ALS were excluded. Several known environmental factors were scrutinized and eliminated, notably lead and other chemical contaminants in soil, water or home-grown vegetation used for food, radon and electromagnetic fields. Some lifestyle-related behavioral risk factors were identified: Prior to clinical onset of motor neuron disease, some patients had a high degree of athleticism and smoked tobacco. Recent investigations on site, based on a new hypothesis, showed that all patients had ingested wild mushrooms, notably poisonous False Morels. Half of the ALS cohort reported acute illness following Gyromitra gigas mushroom consumption. This finding supports the hypothesis that genotoxins of fungal origin may induce motor neuron degeneration.

High spatial resolution dosimetry with uncertainty analysis using Raman micro-spectroscopy readout of radiochromic films.

The purpose of this work is to develop a new approach for high spatial resolution dosimetry based on Raman micro-spectroscopy scanning of radiochromic film (RCF). The goal is to generate dose calibration curves over an extended dose range from 0 to 50Gy and with improved sensitivity to low (<2Gy) doses, in addition to evaluating the uncertainties in dose estimation associated with the calibration curves. Samples of RCF (EBT3) were irradiated at a broad dose range of 0.03-50Gy using an Elekta Synergy clinical linear accelerator. Raman spectra were acquired with a custom-built Raman micro-spectroscopy setup involving a 500 mW, multimode 785 nm laser focused to a lateral spot diameter of 30µm on the RCF. The depth of focus of 34µm enabled the concurrent collection of Raman spectra from the RCF active layer and the polyester laminate. The preprocessed Raman spectra were normalized to the intensity of the 1614cm-1 Raman peak from the polyester laminate that was unaltered by radiation. The mean intensities and the corresponding standard deviation of the active layer Raman peaks at 696, 1445, and 2060cm-1 were determined for the 150×100µm2 scan area per dose value. This was used to generate three calibration curves that enabled the conversion of the measured Raman intensity to dose values. The experimental, fitting, and total dose uncertainty was determined across the entire dose range for the dosimetry system of Raman micro-spectroscopy and RCF. In contrast to previous work that investigated the Raman response of RCFs using different methods, high resolution in the dose response of the RCF, even down to 0.03Gy, was obtained in this study. The dynamic range of the calibration curves based on all three Raman peaks in the RCF extended up to 50Gy with no saturation. At a spatial resolution of 30×30µm2 , the total uncertainty in estimating dose in the 0.5-50Gy dose range was [6-9]% for all three Raman calibration curves. This consisted of the experimental uncertainty of [5-8]%, and the fitting uncertainty of [2.5-4.5]%. The main contribution to the experimental uncertainty was determined to be from the scan area inhomogeneity which can be readily reduced in future experiments. The fitting uncertainty could be reduced by performing Raman measurements on RCF samples at further intermediate dose values in the high and low dose range. The high spatial resolution experimental dosimetry technique based on Raman micro-spectroscopy and RCF presented here, could become potentially useful for applications in microdosimetry to produce meaningful dose estimates in cellular targets, as well as for applications based on small field dosimetry that involve high dose gradients.

Spatially Resolved Optical Characterization of Functional Materials Using Coherence Scanning Interferometry

Herein, some of the work in adapting white light interference microscopy to perform local spectroscopy for measuring the optical properties of microscopic structures is reviewed. Theoretical and experimental results are shown in which Fourier transform processing of the polychromatic fringe signal combined with careful calibration of the optical system is used to make measurements of local reflectance spectra. This approach captures the spectral information within the entire field of view in a single scan, allowing rapid spectral mapping of spatially extended surfaces. Results are shown of local reflectance spectra measured on different materials and buried under transparent layers with a lateral spot size of between 0.5 μm and several micrometer over a field of view up to 650 × 650 μm. The technique is extended to the measurement of local refractive index and thickness of transparent layers as well as to the size determination of small spherical beads buried in transparent and scattering layers with a priori information. The significance of the work is the potential for local materials characterization in complex, hybrid, and functional materials and even disease detection through the study of living cells.

Laser processing over a large area by wavefront-controlled scanning

In the processing of laser materials, the scanning area is limited by the extent to which the laser focus can be maintained within a plane. In order to increase the processing area, one has to enlarge the size of the laser focal spot, and, therefore, reduce the peak intensity and spatial resolution. In this paper, the authors report on a method that extends the processing area by controlling the wavefront of a laser beam using adaptive optics. By adding Zernike modes to the laser beam wavefront, the axial (along the beam propagation direction) location of the focal spot can be tuned within a range larger than the confocal length (two times the Rayleigh length), while the lateral spot size (and, thus, peak intensity and resolution) is maintained throughout the tuning range. A wave propagation model that simulates the propagation of a wavefront-controlled beam is used to compare with the experimental results, which are also confirmed by Zemax calculation. Laser trepanning of circular holes with consistent size on a wide polypropylene sheet is demonstrated. This work shows the potential of tailoring the phase of a laser beam for improving the accuracy and throughput in large-area laser processing of materials.

Experimental and Monte Carlo characterization of a dynamic collimation system prototype for pencil beam scanning proton therapy

There has been a growing interest in the development of energy-specific collimators for low-energy pencil beam scanning (PBS) to reduce the lateral penumbra. One particular device that has been the focus of several recent published works is the dynamic collimation system (DCS), which provides energy-specific collimation by intercepting the scanned proton beam as it nears to target edge with a set of orthogonal trimmer blades. While several computational studies have shown that this dynamic collimator can provide additional healthy tissue sparing, there has not been any rigorous experimental work to benchmark the theoretical models used in these initial studies. Therefore, it was the purpose of this work to demonstrate an experimental method that could integrate an experimental prototype with a clinical PBS system and benchmark the Monte Carlo methods that have been used to model the DCS. An experimental DCS prototype was designed and built in house to actively collimate individual proton beamlets during PBS within a well-characterized experimental setup. Monte Carlo methods were initially used to assess construction tolerances and later benchmarked against measurements, including integral depth dose and lateral asymmetric beamlet profiles. The experimental apparatus and measurement geometry were modeled using MCNP6 benchmarked from measurements performed at the Northwestern Chicago Proton Center. Gamma analysis tests were used to evaluate the agreement between the measured and simulated profiles with a strict 1mm/1% criteria and 5% dose threshold. Excellent agreement was observed between the simulated and measured profiles, which included 1mm/1% gamma analysis pass rates of at least 100% and 95% for the integral depth dose (IDD) profiles and lateral profiles, respectively. Differences in the relative profile shape were observed experimentally between beamlets collimated on- and off-axis, which was attributed to the partial transmission of the beam through an unfocused collimator. Exposure rates resulting from the activation of the device were monitored with survey meter measurements and were found to agree with Monte Carlo estimates of the exposure rate to within 20%. A DCS prototype was constructed and integrated into a clinical dose delivery system. While the results of this work are not exhaustive, they demonstrate the effects of beam source divergence, device activation, and beamlet deflection during scanning, which were found to be successfully modeled using Monte Carlo methods and experimentally benchmarked. Excellent agreement was achieved between the simulated and measured lateral spot profiles of collimated beamlets delivered on- and off-axis in PBS. The Monte Carlo models adequately predicted the measured elevated plateau region in the integral depth-dose profiles from the low-energy scatter off the collimators.

Electron Beam Aberration Correction Using Optical Near Fields.

The interaction between free electrons and optical near fields is attracting increasing attention as a way to manipulate the electron wave function in space, time, and energy. Relying on currently attainable experimental capabilities, we design optical near-field plates to imprint a lateral phase on the electron wave function that can largely correct spherical aberration without the involvement of electric or magnetic lenses in the electron optics, and further generate on-demand lateral focal spot profiles. Our work introduces a disruptive and powerful approach toward aberration correction based on light-electron interactions that could lead to compact and versatile time-resolved free-electron microscopy and spectroscopy.

A vector control technology of waterjet propelled crafts using synchronous and mirror-equiangular thrust control strategy

Vector control technology with single joystick for ship maneuvering is investigated on waterjet propelled crafts for excellent maneuverability and easy manipulation. In this paper, a vector control strategy of ship motion based on synchronous thrust adjustment of two waterjets with mirror-equiangular magnitude which takes the hull cross section as the benchmark was selected firstly. Then the relationship among three manipulation parameters of joystick, six operation parameters of two sets of waterjet propulsion systems, and three motion parameters of hull was established to obtain the vector correspondence between joystick operation and hull motion. On this basis, a vector controller prototype was constructed and mounted on a ship model propelled by two sets of waterjets. The vector control functions of the lateral motion, spot turn and other complex motions such as sectorial motion were realized successfully. The technology of single joystick vector control can be applied to many operation conditions of waterjet propelled crafts, such as berthing, entering and leaving port, position operations at sea, and tactical maneuver to improve the maneuverability and safety of ships. In addition, the vector control function can also be applied to the motion control module for low-speed maneuvering condition of smart ship.

Automatic phase space generation for Monte Carlo calculations of intensity modulated particle therapy

Monte Carlo (MC) is generally considered as the most accurate dose calculation tool for particle therapy. However, a proper description of the beam particle kinematics is a necessary input for a realistic simulation. Such a description can be stored in phase space (PS) files for different beam energies. A PS file contains kinetic information such as energies, positions and travelling directions for particles traversing a plane perpendicular to the beam direction. The accuracy of PS files plays a critical role in the performance of the MC method for dose calculations. A PS file can be generated with a set of parameters describing analytically the beam kinematics. However, determining such parameters can be tedious and time consuming. Thus, we have developed an algorithm to obtain those parameters automatically and efficiently. In this paper, we presented such an algorithm and compared dose calculations using PS automatically generated for the Shanghai Proton and Heavy Ion Center (SPHIC) with measurements. The gamma-index for comparing calculated depth dose distributions (DDD) with measurements are above 96.0% with criterion 0.6%/0.6 mm. For each single energy, the mean difference percentage between calculated lateral spot sizes at 5 different locations along beam direction and measurements are below 3.5%.

A standardized commissioning framework of Monte Carlo dose calculation algorithms for proton pencil beam scanning treatment planning systems.

Treatment planning systems (TPSs) from different vendors can involve different implementations of Monte Carlo dose calculation (MCDC) algorithms for pencil beam scanning (PBS) proton therapy. There are currently no guidelines for validating non-water materials in TPSs. Furthermore, PBS-specific parameters can vary by 1-2 orders of magnitude among different treatment delivery systems (TDSs). This paper proposes a standardized framework on the use of commissioning data and steps to validate TDS-specific parameters and TPS-specific heterogeneity modeling to potentially reduce these uncertainties. A standardized commissioning framework was developed to commission the MCDC algorithms of RayStation 8A and Eclipse AcurosPT v13.7.20 using water and non-water materials. Measurements included Bragg peak depth-dose and lateral spot profiles and scanning field outputs for Varian ProBeam. The phase-space parameters were obtained from in-air measurements and the number of protons per MU from output measurements of 10×10cm2 square fields at a 2cm depth. Spot profiles and various PBS field measurements at additional depths were used to validate TPS. Human tissues in TPS, Gammex phantom materials, and artificial materials were used for the TPS benchmark and validation. The maximum differences of phase parameters, spot sigma, and divergence between MCDC algorithms are below 4.5µm and 0.26mrad in air, respectively. Comparing TPS to measurements at depths, both MC algorithms predict the spot sigma within 0.5mm uncertainty intervals, the resolution of the measurement device. Beam Configuration in AcurosPT is found to underestimate number of protons per MU by ~2.5% and requires user adjustment to match measured data, while RayStation is within 1% of measurements using Auto model. A solid water phantom was used to validate the range accuracy of non-water materials within 1% in AcurosPT. The proposed standardized commissioning framework can detect potential issues during PBS TPS MCDC commissioning processes, and potentially can shorten commissioning time and improve dosimetric accuracies. Secondary MCDC can be used to identify the root sources of disagreement between primary MCDC and measurement.

Benchmarking a GATE/Geant4 Monte Carlo model for proton beams in magnetic fields.

PurposeMagnetic resonance guidance in proton therapy (MRPT) is expected to improve its current performance. The combination of magnetic fields with clinical proton beam lines poses several challenges for dosimetry, treatment planning and dose delivery. Proton beams are deflected by magnetic fields causing considerable changes in beam trajectories and also a retraction of the Bragg peak positions. A proper prediction and compensation of these effects is essential to ensure accurate dose calculations. This work aims to develop and benchmark a Monte Carlo (MC) beam model for dose calculation of MRPT for static magnetic fields up to 1 T.MethodsProton beam interactions with magnetic fields were simulated using the GATE/Geant4 toolkit. The transport of charged particle in custom 3D magnetic field maps was implemented for the first time in GATE. Validation experiments were done using a horizontal proton pencil beam scanning system with energies between 62.4 and 252.7 MeV and a large gap dipole magnet (B = 0–1 T), positioned at the isocenter and creating magnetic fields transverse to the beam direction. Dose was measured with Gafchromic EBT3 films within a homogeneous PMMA phantom without and with bone and tissue equivalent material slab inserts. Linear energy transfer (LET) quenching of EBT3 films was corrected using a linear model on dose‐averaged LET method to ensure a realistic dosimetric comparison between simulations and experiments. Planar dose distributions were measured with the films in two different configurations: parallel and transverse to the beam direction using single energy fields and spread‐out Bragg peaks. The MC model was benchmarked against lateral deflections and spot sizes in air of single beams measured with a Lynx PT detector, as well as dose distributions using EBT3 films. Experimental and calculated dose distributions were compared to test the accuracy of the model.ResultsMeasured proton beam deflections in air at distances of 465, 665, and 1155 mm behind the isocenter after passing the magnetic field region agreed with MC‐predicted values within 4 mm. Differences between calculated and measured beam full width at half maximum (FWHM) were lower than 2 mm. For the homogeneous phantom, measured and simulated in‐depth dose profiles showed range and average dose differences below 0.2 mm and 1.2%, respectively. Simulated central beam positions and widths differed <1 mm to the measurements with films. For both heterogenous phantoms, differences within 1 mm between measured and simulated central beam positions and widths were obtained, confirming a good agreement of the MC model.ConclusionsA GATE/Geant4 beam model for protons interacting with magnetic fields up to 1 T was developed and benchmarked to experimental data. For the first time, the GATE/Geant4 model was successfully validated not only for single energy beams, but for SOBP, in homogeneous and heterogeneous phantoms. EBT3 film dosimetry demonstrated to be a powerful dosimetric tool, once the film response function is LET corrected, for measurements in‐line and transverse to the beam direction in magnetic fields. The proposed MC beam model is foreseen to support treatment planning and quality assurance (QA) activities toward MRPT.

Extension of longitudinal measuring range of a confocal surface profilometer using a gradient-intensity probe beam

A confocal surface profilometer with a long vertical range was developed by using a gradient-intensity probe beam. In conventional confocal microscopy, vertical scanning of a tested surface by either stage height movement or the shifting of objective focus is time-consuming, thus making unacceptable measurement efficiency for in situ inspection. To overcome this, in this research, a quasi-Bessel beams (QBB) with the properties of narrow lateral spot and gradient-intensity axial distribution is adopted as a probe beam to generate an accurate intensity-to-height conversion for 3-D profile measurement. The influences of the incident Gaussian beam and the real axicon lens with an oblate-tip shape on the axial probe beam are calculated and analyzed. By switching the axicon lenses with different base angles, the measurable height of the surface profilometer can be ranged from a few millimeters to several hundred millimeters. According to the measured results, it was found that by using the axicon lens with a base angle of 20°, the measurable height reached at least 2 mm. By integrating a high-speed acquisition unit, the proposed profilometer can further increase the measurement efficiency, especially fitting in situ automatic optical inspection (AOI).