Equivalent Path Length Research Articles

Investigation of four-dimensional (4D) Monte Carlo dose calculation in real-time tumor tracking stereotatic body radiotherapy for lung cancers.

To investigate the dosimetric variations and radiobiological impacts as a consequence of delivering treatment plans of 3D nature in 4D manner based on the 4D Monte Carlo treatment planning framework implemented on Cyberknife. Dose distributions were optimized on reference 3D images at end of exhale phase of a 4DCT dataset for twenty-five lung cancer patients treated with 60 Gy / 3Fx or 48 Gy / 4Fx. Deformable image registrations (DIR) between individual 3DCT images to the reference 3DCT image in the 4DCT study were performed to interpolate doses calculated on multiple anatomical geometries back on to the reference geometry to compose a 4D dose distribution that included the tracking beam motion and organ deformation. The 3D and 4D dose distributions that were initially calculated with the equivalent path-length (EPL) algorithm (3DEPL dose and 4DEPL dose) were recalculated with the Monte Carlo algorithm (3DMC dose and 4DMC dose). Dosimetric variations of V60Gy / 48Gy and D99 of GTV, mean doses to the lung and the heart and maximum dose (D1 ) of the spinal cord as a consequence of tracking beam motion in deforming anatomy, dose calculation algorithm, and both were quantified by the relative change from 4DMC to 3DMC doses, from 4DMC to 4DEPL doses, and from 4DMC to 3DEPL doses, respectively. Comparing 4DMC to 3DEPL plans, V60Gy / 48Gy and D99 of GTV decreased considerably by 13 ± 22% (mean ± 1SD) and 9.2 ± 5.5 Gy but changes of normal tissue doses were not more than 0.5 Gy on average. The generalized equivalent uniform dose (gEUD) and tumor control probability (TCP) were reduced by 14.3 ± 8.8 Gy and 7.5 ± 5.2%, and normal tissue complication probability (NTCP) for myelopathy and pericarditis were close to zero and NTCP for radiation pneumonitis was reduced by 2.5 ± 4.1%. Comparing 4DMC to 4DEPL plans found decreased V60Gy / 48Gy and D99 by 12.3 ± 21.6% and 7.3 ± 5.3 Gy, the normal tissues doses by 0.5 Gy on average, gEUD and TCP by 13.0 ± 8.6 Gy and 7.1 ± 5.1%. Comparing 4DMC to 3DMC doses, V60Gy / 48Gy and D99 of GTV was reduced by 5.2 ± 8.8 %and 2.6 ± 3.3 Gy, and normal tissues hardly changed from 4DMC to 3DMC doses. The corresponding decreases of gEUD and TCP were 2.8 ± 4.0 Gy and 1.6 ± 2.4%. The large discrepancy between original 3DEPL plan and benchmarking 4DMC plan is predominately due to dose calculation algorithms as the tracking beam motion and organ deformation hardly influenced doses of normal tissues and moderately decreased V60Gy / 48Gy and D99 of GTV. It is worth to make a thoughtful weight of the benefits of full 4D MC dose calculation and consider 3D MC dose calculation as a compromise of 4D MC dose calculation considering the multifold computation time. This article is protected by copyright. All rights reserved.

Circular multireflection cell for midinfrared laser spectroscopic trace gas analysis

Carbon dioxide (CO2) is an atmospheric trace gas; therefore its accurate sensing is of great interest. Optical sensors exploiting the midinfrared (mid-IR) light absorption of CO2 provide high sensitivity and are widely used in medical diagnostics, atmospheric monitoring, remote sensing, and industrial applications. In this work, we demonstrated accurate CO2 gas sensing at a 4.2-μm wavelength. In addition, detecting the weak mid-IR molecular absorption bands of gases at low concentrations requires using increased optical path lengths. The most obvious method that can expand the potential beam path in a spectroscopic system is to use a longer linear gas cell, which in some situations may be adequate; however, space and volume requirements need to be considered. A circular multireflection (CMR) cell was used to reflect the radiation back and forth through the sample medium multiple times, thus greatly reducing the size footprint compared with a linear cell of equivalent optical path length. A CMR cell was designed and constructed to allow for multireflections within the cell. The optical alignment of the cell and the convenience of changing the optical path length by adjusting its position with respect to the incident light beam were also used to maximize the advantages of the device. The cell has an inner diameter of 6.5 cm, and the path length can reach up to 123.5 cm, which is equivalent to 19 beam passes of reflections. This work will be used as the groundwork for designing an instrument for the high-resolution measurement of CO2 gas in planetary atmospheres.

Relative Humidity Measurement of Air in Low-Temperature Ranges Using Low-Frequency Acoustic Waves and Correlation Signal Processing Techniques.

Air relative humidity (RH) is an important control parameter in many industrial processes. The acoustic method is a novel technique to measure air humidity non-intrusively. Relevant research is limited. Existing methods use ultrasonic waves as a sound source and air humidity is measured by measuring the sound attenuation. In this paper, a novel air humidity measurement system using low-frequency sound waves as a sound source and two acoustic sensors is proposed. Air humidity is acquired by measuring sound speed in the air. Sound speed mainly depends on air temperature, humidity, atmospheric pressure, and air composition. The influence of air temperature, atmospheric pressure, and air constituent concentrations on the RH measurement is analyzed theoretically. A 0.1 s linear chirp signal in the frequency range of 200–500 Hz is selected as the sound source. Sound travel time is calculated by cross-correlating the sound signals received by the two acoustic sensors. To improve the accuracy of the sound speed measurement, sound speed under different RH points is obtained through reference RH experiments and substituted into the calibration equation. Then, equivalent sound path length and systematic delay are estimated using the least squares method. After obtaining these two parameter values, the sound speed measured by the system is closer to the theoretical value at the same RH point. In validation experiments using RH measured by a thermo-hygrometer as a comparison, the relative errors of the acoustically measured RH are within 9.9% in the RH range of 40.7–87.1%, and the standard deviation is within 4.8%.

Proton radiography using discrete range modulation method – A Monte Carlo study

The advantage of proton therapy is the Bragg peak, but its range uncertainties can also lead to a certain dose deviation. This uncertainty can be reduced by using proton radiography to determine the water equivalent path length (WEPL). In this study, a discrete range modulation (DRM) method is introduced and its image quality compared with that of the previously reported continuous range modulation (CRM) method. MCNPX2.7.0 was used to simulate parallel proton beams with energies from 70 to 230 MeV at 2.5 MeV/step. In the DRM method, the WEPL can be related to R80 using the relationship between the proton energy and energy deposition at 80% of the fall-off (E80). For the CRM method, a dose gradient plan was created to estimate the relationship between WEPL and accumulated dose. DRM and CRM were compared using four different phantoms. For the slab, CIRS 062M head and cake phantoms, DRM had a good agreement with the theoretic WEPL: within 2 mm or 1%, and with a standard deviation of less than 0.5 mm. For the sphere phantom, DRM can achieve a GammapRG passing rate of 0.999, passing the criteria of a WEPL deviation of less than 2 mm and a distance-to-agreement of less than 3 mm. Proton radiographs made by the DRM method are more accurate and precise than those made by CRM. Two algorithms are proposed to solve the range mixing problems. DRM can also save unnecessary doses that do not contribute to imaging.

A reconstruction approach for proton computed tomography by modeling the integral depth dose of the scanning proton pencil beam.

To present a proton computed tomography (pCT) reconstruction approach that models the integral depth dose (IDD) of the clinical scanning proton beam into beamlets. Using a multilayer ionization chamber (MLIC) as the imager, the proposed pCT system and the reconstruction approach can minimize extra ambient neutron dose and simplify the beamline design by eliminating an additional collimator to confine the proton beam. Monte Carlo simulation was applied to digitally simulate the IDDs of the exiting proton beams detected by the MLIC. A forward model was developed to model each IDD into a weighted sum of percentage depth doses of the constituent beamlets separated laterally by 1mm. The water equivalent path lengths (WEPLs) of the beamlets were determined by iteratively minimizing the squared L2-norm between the forward projected and simulated IDDs. The final WEPL values were reconstructed to pCT images, that is, proton stopping power ratio (SPR) maps, through simultaneous algebraic reconstruction technique with total variation regularization. The reconstruction process was tested with a digital cylindrical water-based phantom and an ICRP adult reference computational phantom. The mean of SPR within regions of interest (ROIs) and the WEPL along a 4mm-wide beam ( ) were compared with the reference values. The spatial resolution was analyzed at the edge of a cortical insert of the cylindrical phantom. The percentage deviations from reference SPR were within±1% in all selected ROIs. The mean absolute error of the reconstructed SPR was 0.33%, 0.19%, and 0.27% for the cylindrical phantom, the adult phantom at the head and lung region, respectively. The corresponding percentage deviations from reference were 0.48±0.64%, 0.28±0.48%, and 0.22±0.49%. The full width at half maximum of the line spread function (LSF) derived from the radial edge spread function (ESF) of a cortical insert was 0.13cm. The frequency at 10% of the modulation transfer function (MTF) was 6.38cm-1 . The mean signal-to-noise ratio (SNR) of all the inserts was 2.45. The mean imaging dose was 0.29 and 0.25cGy at the head and lung region of the adult phantom, respectively. A new pCT reconstruction approach was developed by modeling the IDDs of the uncollimated scanning proton beams in the pencil beam geometry. SPR accuracy within±1%, spatial resolution of better than 2mm at 10% MTF, and imaging dose at the magnitude of mGy were achieved. Potential side effects caused by neutron dose were eliminated by removing the extra beam collimator.

Optimizing calibration settings for accurate water equivalent path length assessment using flat panel proton radiography

Objective: Proton range uncertainties can compromise the effectiveness of proton therapy treatments. Water equivalent path length (WEPL) assessment by flat panel detector proton radiography (FP-PR) can provide means of range uncertainty detection. Since WEPL accuracy intrinsically relies on the FP-PR calibration parameters, the purpose of this study is to establish an optimal calibration procedure that ensures high accuracy of WEPL measurements. To that end, several calibration settings were investigated. Approach: FP-PR calibration datasets were obtained simulating PR fields with different proton energies, directed towards water-equivalent material slabs of increasing thickness. The parameters investigated were the spacing between energy layers (ΔE) and the increment in thickness of the water-equivalent material slabs (ΔX) used for calibration. 30 calibrations were simulated, as a result of combining ΔE = 9, 7, 5, 3, 1 MeV and ΔX = 10, 8, 5, 3, 2, 1 mm. FP-PRs through a CIRS electron density phantom were simulated, and WEPL images corresponding to each calibration were obtained. Ground truth WEPL values were provided by range probing multi-layer ionization chamber simulations on each insert of the phantom. Relative WEPL errors between FP-PR simulations and ground truth were calculated for each insert. Mean relative WEPL errors and standard deviations across all inserts were computed for WEPL images obtained with each calibration. Main results: Large mean and standard deviations were found in WEPL images obtained with large ΔE values (ΔE = 9 or 7 MeV), for any ΔX. WEPL images obtained with ΔE ≤ 5 MeV and ΔX ≤ 5 mm resulted in a WEPL accuracy with mean values within ±0.5% and standard deviations around 1%. Significance: An optimal FP calibration in the framework of this study was established, characterized by 3 MeV ≤ ΔE ≤ 5 MeV and 2 mm ≤ ΔX ≤ 5 mm. Within these boundaries, highly accurate WEPL acquisitions using FP-PR are feasible and practical, holding the potential to assist future online range verification quality control procedures.

Investigating beam range uncertainty in proton prostate treatment using pelvic-like biological phantoms

This study aims to develop a method for verifying site-specific and/or beam path specific proton beam range, which could reduce range uncertainty margins and the associated treatment complications. It investigates the range uncertainties from both CT HU to relative stopping power conversion and patient positioning errors for prostate treatment using pelvic-like biological phantoms. Three 25 × 14 × 12 cm3 phantoms, made of fresh animal tissues mimicking the pelvic anatomies of prostate patients, were scanned with a general electric CT simulator. A 22 cm circular passive scattering beam with 29 cm range and 8 cm modulation width was used to measure the water equivalent path lengths (WEPL) through the phantoms at multiple points using the dose extinction method with a MatriXXPT detector. The measured WEPLs were compared to those predicted by TOPAS simulations and ray-tracing WEPL calculations. For the three phantoms, the WEPL differences between measured and theoretical prediction (WDMT) are below 1.8% for TOPAS, and 2.5% for ray-tracing. WDMT varies with phantom anatomies by about 0.5% for both TOPAS and ray-tracing. WDMT also correlates with the tissue types of a specific treated region. For the regions where the proton beam path is parallel to sharp bone edges, the WDMTs of TOPAS and ray-tracing respectively reach up to 1.8% and 2.5%. For the region where proton beams pass through just soft tissues, the WDMT is mostly less than 1% for both TOPAS and ray-tracing. For prostate treatments, range uncertainty depends on the tissue types within a specific treated region, patient anatomies and the range calculation methods in the planning algorithms. Our study indicates range uncertainty is less than 2.5% for the whole treated region with both ray-tracing and TOPAS, which suggests the potential to reduce the current 3.5% range uncertainty margin used in the clinics by at least 1% even for single-energy CT data.

In-vivo proton range verification for reducing the risk of permanent alopecia in medulloblastoma treatment

ObjectiveThe purpose of this work is the clinical commissioning of a recently developed in-vivo range verification system for the head treatment of pediatric medulloblastoma patients. Inaccurate beam range for such treatment could lead to either inadequate dose coverage of the target volume or excessive dose to the head skin resulting permanent alopecia. MethodsThe in-vivo range verification system is designed to perform pre-treatment range verification and adjustment. An array of Si-diode detectors is to be placed on the patient immobilization mask in the exit direction of a whole-brain field; signal is analyzed, and the extracted water equivalent path length (WEPL) is compared to the expected one, revealing if a range correction is needed. The method was tested in solid water and anthropomorphic head phantom, with validation based on independent WEPL measurements. The measured WEPL were compared to those computed by the treatment planning system (TPS). ResultsThe accuracy for the WEPL measurements by the diode system in both solid water and anthropomorphic head phantom were on average within a millimeter from more accurate measurement by the dose-extinction technique, with the error range for the two phantoms as (0–1 ​mm) and (0–1.3 ​mm), respectively. When compared to the WEPL calculated by the treatment planning system, the measured values were on average within 1 ​% (range 0–3%) of the beam range. The accuracy of dose measurements by the diodes in the fall-off part of the depth dose profile was validated against the reference Markus chamber. No need for further correction (due to different beam parameters and detector dose ageing effects) was found. ConclusionsThe range verification workflow was successfully tested in the anthropomorphic head phantom. The performance of the in-vivo range verification system and related workflow meet the clinical requirements in terms of the needed WEPL accuracy for pretreatment range verification.

Investigating particle track topology for range telescopes in particle radiography using convolutional neural networks

Background Proton computed tomography (pCT) and radiography (pRad) are proposed modalities for improved treatment plan accuracy and in situ treatment validation in proton therapy. The pCT system of the Bergen pCT collaboration is able to handle very high particle intensities by means of track reconstruction. However, incorrectly reconstructed and secondary tracks degrade the image quality. We have investigated whether a convolutional neural network (CNN)-based filter is able to improve the image quality. Material and methods The CNN was trained by simulation and reconstruction of tens of millions of proton and helium tracks. The CNN filter was then compared to simple energy loss threshold methods using the Area Under the Receiver Operating Characteristics curve (AUROC), and by comparing the image quality and Water Equivalent Path Length (WEPL) error of proton and helium radiographs filtered with the same methods. Results The CNN method led to a considerable improvement of the AUROC, from 74.3% to 97.5% with protons and from 94.2% to 99.5% with helium. The CNN filtering reduced the WEPL error in the helium radiograph from 1.03 mm to 0.93 mm while no improvement was seen in the CNN filtered pRads. Conclusion The CNN improved the filtering of proton and helium tracks. Only in the helium radiograph did this lead to improved image quality.

Technical Note: First report on an in vivo range probing quality control procedure for scanned proton beam therapy in head and neck cancer patients.

The capability of proton therapy to provide highly conformal dose distributions is impaired by range uncertainties. The aim of this work is to apply range probing (RP), a form of a proton radiography-based quality control (QC) procedure for range accuracy assessment in head and neck cancer (HNC) patients in a clinical setting. This study included seven HNC patients. RP acquisition was performed using a multi-layer ionization chamber (MLIC). Per patient, two RP frames were acquired within the first two weeks of treatment, on days when a repeated CT scan was obtained. Per RP frame, integral depth dose (IDD) curves of 81 spots around the treatment isocenter were acquired. Range errors are determined as a discrepancy between calculated IDDs in the treatment planning system and measured residual ranges by the MLIC. Range errors are presented relative to the water equivalent path length of individual proton spots. In addition to reporting results for complete measurement frames, an analysis, excluding range error contributions due to anatomical changes, is presented. Discrepancies between measured and calculated ranges are smaller when performing RP calculations on the day-specific patient anatomy rather than the planning CT. The patient-specific range evaluation shows an agreement between calculated and measured ranges for spots in anatomically consistent areas within 3% (1.5 standard deviation). The results of an RP-based QC procedure implemented in the clinical practice for HNC patients have been demonstrated. The agreement of measured and simulated proton ranges confirms the 3% uncertainty margin for robust optimization. Anatomical variations show a predominant effect on range accuracy, motivating efforts towards the implementation of adaptive radiotherapy.

Helium radiography with a digital tracking calorimeter—a Monte Carlo study for secondary track rejection

Radiation therapy using protons and heavier ions is a fast-growing therapeutic option for cancer patients. A clinical system for particle imaging in particle therapy would enable online patient position verification, estimation of the dose deposition through range monitoring and a reduction of uncertainties in the calculation of the relative stopping power of the patient.Several prototype imaging modalities offer radiography and computed tomography using protons and heavy ions. A Digital Tracking Calorimeter (DTC), currently under development, has been proposed as one such detector. In the DTC 43 longitudinal layers of laterally stacked ALPIDE CMOS monolithic active pixel sensor chips are able to reconstruct a large number of simultaneously recorded proton tracks.In this study, we explored the capability of the DTC for helium imaging which offers favorable spatial resolution over proton imaging. Helium ions exhibit a larger cross section for inelastic nuclear interactions, increasing the number of produced secondaries in the imaged object and in the detector itself. To that end, a filtering process able to remove a large fraction of the secondaries was identified, and the track reconstruction process was adapted for helium ions.By filtering on the energy loss along the tracks, on the incoming angle and on the particle ranges, 97.5% of the secondaries were removed. After passing through 16 cm water, 50.0% of the primary helium ions survived; after the proposed filtering 42.4% of the primaries remained; finally after subsequent image reconstruction 31% of the primaries remained. Helium track reconstruction leads to more track matching errors compared to protons due to the increased available focus strength of the helium beam. In a head phantom radiograph, the Water Equivalent Path Length error envelope was 1.0 mm for helium and 1.1 mm for protons. This accuracy is expected to be sufficient for helium imaging for pre-treatment verification purposes.

Accuracy Evaluation of EPL and ETAR Algorithms in the Treatment Planning Systems using CIRS Thorax Phantom.

Background: It is recommended for each set of radiation data and algorithm that subtle deliberation is done regarding dose calculation accuracy. Knowing the errors in dose calculation for each treatment plan will result in an accurate estimate of the actual dose achieved by the tumor.Objective: This study aims to evaluate the equivalent path length (EPL) and equivalent tissue air ratio (ETAR) algorithms in radiation dose calculation.Material and Methods: In this experimental study, the TEC-DOC 1583 guideline was used. Measurements and calculations were obtained for each algorithm at specific points in thorax CIRS phantom for 6 and 18 MVs and results were compared.Results: In the EPL, calculations were in agreement with measurements for 27 points and differences between them ranged from 0.1% to 10.4% at 6 MV. The calculations were in agreement with measurements for 21 points and differences between them ranged from 0.4% to 13% at 18 MV. In ETAR, calculations were also in consistent with measurements for 21 points, and differences between them ranged from 0.1% to 9% at 6 MV. Moreover, for 18 MV, the calculations were in agreement with measurements for 17 points and differences between them ranged from 0% to 11%. Conclusion: For the EPL algorithm, more dose points were in consistent with acceptance criteria. The errors in the ETAR were 1% to 2% less than the EPL. The greatest calculation error occurs in low-density lung tissue with inhomogeneities or in high-density bone. Errors were larger in shallow depths. The error in higher energy was more than low energy beam.

OpenPR - A computational tool for CT conversion assessment with proton radiography.

One of the main sources of uncertainty in proton therapy is the conversion of the Hounsfield Units of the planning CT to (relative) proton stopping powers. Proton radiography provides range error maps but these can be affected by other sources of errors as well as the CT conversion (e.g., residual misalignment). To better understand and quantify range uncertainty, it is desirable to measure the individual contributions and particularly those associated to the CT conversion. A workflow is proposed to carry out an assessment of the CT conversion solely on the basis of proton radiographs of real tissues measured with a multilayer ionization chamber (MLIC). The workflow consists of a series of four stages: (a) CT and proton radiography acquisitions, (b) CT and proton radiography registration in postprocessing, (c) sample-specific validation of the semi-empirical model both used in the registration and to estimate the water equivalent path length (WEPL), and (d) WEPL error estimation. The workflow was applied to a pig head as part of the validation of the CT calibration of the proton therapy center PARTICLE at UZ Leuven, Belgium. The CT conversion-related uncertainty computed based on the well-established safety margin rule of 1.2mm+2.4% were overestimated by 71% on the pig head. However, the range uncertainty was very much underestimated where cavities were encountered by the protons. Excluding areas with cavities, the overestimation of the uncertainty was 500%. A correlation was found between these localized errors and HUs between -1000 and -950, suggesting that the underestimation was not a consequence of an inaccurate conversion but was probably rather due to the resolution of the CT leading to material mixing at interfaces. To reduce these errors, the CT calibration curve was adapted by increasing the HU interval corresponding to the air up to -950. The application of the workflow as part of the validation of the CT conversion to RSPs showed an overall overestimation of the expected uncertainty. Moreover, the largest WEPL errors were found to be related to the presence of cavities which nevertheless are associated with low WEPL values. This suggests that the use of this workflow on patients or in a generalized study on different types of animal tissues could shed sufficient light on how the contributions to the CT conversion-related uncertainty add up to potentially reduce up to several millimeters the uncertainty estimations taken into account in treatment planning. All the algorithms required to perform the workflow were implemented in the computational tool named openPR which is part of openREGGUI, an open-source image processing platform for adaptive proton therapy.

PO-1461: Benchmarking proton therapy water equivalent path length calculations against TPS algorithms

PO-1461: Benchmarking proton therapy water equivalent path length calculations against TPS algorithms

A High-Granularity Digital Tracking Calorimeter Optimized for Proton CT

A typical proton CT (pCT) detector comprises a tracking system, used to measure the proton position before and after the imaged object, and an energy/range detector to measure the residual proton range after crossing the object. The Bergen pCT collaboration was established to design and build a prototype pCT scanner with a high granularity digital tracking calorimeter (DTC) used as both tracking and energy/range detector. In this work the conceptual design and the layout of the mechanical and electronics implementation, along with a Monte Carlo (MC) simulation of the new pCT system are reported. The DTC is a multilayer structure with a lateral aperture of 27 cm × 15 cm, made of 41 detector/absorber sandwich layers (calorimeter), with aluminum (3.5 mm) used both as absorber and carrier, and 2 additional layers used as tracking system (rear trackers) positioned downstream of the imaged object; no tracking upstream the object is included. The rear tracker’s structure only differs from the calorimeter layers for the carrier made of ~200 μm carbon fleece and carbon paper (carbon-Epoxy sandwich), to minimize scattering. Each sensitive layer consists of 108 ALPIDE chip sensors (developed for ALICE, CERN) bonded on a polyimide flex and subsequently bonded to a larger flexible printed circuit board. Beam tests tailored to the pCT operation have been performed using high-energetic (50–220 MeV/u) proton and ion beams at the Heidelberg Ion-Beam Therapy Center (HIT) in Germany. These tests proved the ALPIDE response independent of occupancy and proportional to the particle energy deposition, making the distinction of different ion tracks possible. The read-out electronics is able to handle enough data to acquire a single 2D image in few seconds making the system fast enough to be used in a clinical environment. For the reconstructed images in the modeled MC simulation, the water equivalent path length error is lower than 2 mm, and the relative stopping power accuracy is better than 0.4%. Thanks to its ability to detect different types of radiation and its specific design, the pCT scanner can be employed for additional online applications during the treatment, such as in-situ proton range verification.

Evaluation of an a priori scatter correction algorithm for cone-beam computed tomography based range and dose calculations in proton therapy.

Background and purposeScatter correction of cone-beam computed tomography (CBCT) projections may enable accurate online dose-delivery estimations in photon and proton-based radiotherapy. This study aimed to evaluate the impact of scatter correction in CBCT-based proton range/dose calculations, in scans acquired in both proton and photon gantries.Material and methodsCBCT projections of a Catphan and an Alderson phantom were acquired on both a proton and a photon gantry. The scatter corrected CBCTs (corrCBCTs) and the clinical reconstructions (stdCBCTs) were compared against CTs rigidly registered to the CBCTs (rigidCTs). The CBCTs of the Catphan phantom were segmented by materials for CT number analysis. Water equivalent path length (WEPL) maps were calculated through the Alderson phantom while proton plans optimized on the rigidCT and recalculated on all CBCTs were compared in a gamma analysis.ResultsIn medium and high-density materials, the corrCBCT CT numbers were much closer to those of the rigidCT than the stdCBCTs. E.g. in the 50% bone segmentations the differences were reduced from above 300 HU (with stdCBCT) to around 60–70 HU (with corrCBCT). Differences in WEPL from the rigidCT were typically well below 5 mm for the corrCBCTs, compared to well above 10 mm for the stdCBCTs with the largest deviations in the head and thorax regions. Gamma pass rates (2%/2mm) when comparing CBCT-based dose re-calculations to rigidCT calculations were improved from around 80% (with stdCBCT) to mostly above 90% (with corrCBCT).ConclusionScatter correction leads to substantial artefact reductions, improving accuracy of CBCT-based proton range/dose calculations.

Feasibility of energy-resolved dose imaging technique in pencil beam scanning mode

Purpose: Proton energy-resolved dose imaging (pERDI) is a recently proposed technique to generate water equivalent path length (WEPL) images using a single detector. Owing to its simplicity in instrumentation, analysis and the possibility of using the in-room x-ray flat panels as detectors, this technique offers a promising avenue towards a clinically usable imaging system for proton therapy using scanned beams. The purpose of this study is to estimate the achievable accuracy in WEPL and Relative Stopping Power (RSP) using the pERDI technique and to assess the minimum dose required to achieve such accuracy. The novelty of this study is the first demonstration of the feasibility of pERDI technique in the pencil beam scanning (PBS) mode. Methods: A solid water wedge was placed in front of a 2D detector (Lynx). A library of energy-resolved dose functions (ERDF) was generated from the dose deposited in the detector by 50 PBS layers of energy varying from 100 MeV to 225 MeV. This set-up is further used to image the following configurations using the pERDI technique: stair-case shaped solid water phantom (configuration 1), electron density phantom (configuration 2) and head phantom (configuration 3). The result from configuration 1 was used to determine the achievable WEPL accuracy. The result from configuration 2 was used to estimate the relative uncertainty in RSP. Configuration 3 was used to evaluate the effect of range mixing on the WEPL. In all three cases, the variation of the accuracy with respect to dose, by varying the number of scanning layers, was also studied. Results: An accuracy of 1 mm in WEPL was achieved using the Lynx detector with an imaging field of 10 PBS layers or more, which is equivalent to a total dose of 5 cGy. The RSP is measured with a precision better than 2% for all homogeneous inserts of tissue surrogates. The pERDI technique failed for tissues surrogates with total WEPL outside the calibration window (WEPL < 70 mm) like in the case of lung exhale and lung inhale. The imaging of an anthropomorphic head phantom, in the same condition, produced a WEPL radiograph and compared to the WEPL derived from CT using gamma index analysis. The gamma index failed in the heterogeneous areas due to range mixing. Conclusions: The pERDI technique is a promising clinically usable imaging modality for reducing range uncertainties and set-up errors in proton therapy. The first results have demonstrated that WEPL and RSP can be estimated with clinically acceptable accuracy using the Lynx detector. Similar accuracy is also expected with in-room flat-panel detectors but at significantly reduced imaging dose. Though the issue of range mixing is still to be addressed, we expect that a statistical moment analysis of the ERDFs can be implemented to filter out the regions with high gradient of range mixing.

Evolutions of microstructure and magnetic properties of heatproof domain-refined silicon steel during annealing and its application

Heatproof domain-refined grain-oriented (DRGO) silicon steel is an ideal material for the manufacture of S15-type ultra-high efficiency transformer with three-dimensional wound-core (TDWC). In this paper, the microstructure, grain orientation, and magnetic performance of 23ZDMH80 heatproof silicon steel during stress-relief annealing were investigated by scanning electron microscopy (SEM) and electron backscattered diffraction (EBSD) techniques. Based on heatproof DRGO steel, the Epstein-SST (single sheet tester) conversion factor and equivalent magnetic path length were calculated, and then an ultra-high energy efficiency transformer was developed. Results show that, the average size of micrograins in notched region increased from 42.3 μm to 68.2 μm after annealing at 850 ℃ for 0–8 h. The orientations of micrograins grown on large-size Goss grain included not only the non-〈0 0 1〉 orientations, such as {210} 〈−241〉, {215} < 1–20>, and {1 1 0} < 1–12}, which were harmful to the magnetic properties, but also the favorable {1 0 0} 〈0 0 1〉 and Goss. Correspondingly, the core loss of 5 pairs of Epstein specimens first decreased then increased slowly (the growth rates below 1.2%) during annealing. The Epstein-SST conversion factor δP for 23ZDMH80 heatproof steel was 8.6% (higher than 5.0% the value recommended by IEC standard), and the equivalent magnetic path length was 0.489 m. Finally, the developed 10 kV/400 kVA transformer with TDWC presented ultra low no-load loss of 289 W, load loss of 3072 W, and noise of 35.4 dB, it means significant advantage of energy conservation and environmental protection.

Assessment of range uncertainty in lung-like tissue using a porcine lung phantom and proton radiography

Thoracic tumours are increasingly considered indications for pencil beam scanned proton therapy (PBS-PT) treatments. Conservative robustness settings have been suggested due to potential range straggling effects caused by the lung micro-structure. Using proton radiography (PR) and a 4D porcine lung phantom, we experimentally assess range errors to be considered in robust treatment planning for thoracic indications. A human-chest-size 4D phantom hosting inflatable porcine lungs and corresponding 4D computed tomography (4DCT) were used. Five PR frames were planned to intersect the phantom at various positions. Integral depth-dose curves (IDDs) per proton spot were measured using a multi-layer ionisation chamber (MLIC). Each PR frame consisted of 81 spots with an assigned energy of 210 MeV (full width at half maximum (FWHM) 8.2 mm). Each frame was delivered five times while simultaneously acquiring the breathing signal of the 4D phantom, using an ANZAI load cell. The synchronised ANZAI and delivery log file information was used to retrospectively sort spots into their corresponding breathing phase. Based on this information, IDDs were simulated by the treatment planning system (TPS) Monte Carlo dose engine on a dose grid of 1 mm. In addition to the time-resolved TPS calculations on the 4DCT phases, IDDs were calculated on the average CT. Measured IDDs were compared with simulated ones, calculating the range error for each individual spot. In total, 2025 proton spots were individually measured and analysed. The range error of a specific spot is reported relative to its water equivalent path length (WEPL). The mean relative range error was 1.2% (1.5 SD 2.3 %) for the comparison with the time-resolved TPS calculations, and 1.0% (1.5 SD 2.2 %) when comparing to TPS calculations on the average CT. The determined mean relative range errors justify the use of 3% range uncertainty for robust treatment planning in a clinical setting for thoracic indications.

Similar-cases-based planning approaches with beam angle optimizations using water equivalent path length for lung stereotactic body radiation therapy.

This study aimed to propose automated treatment planning approaches based on similar cases with beam angle optimizations using water equivalent path length (WEPL) to avoid lung and rib doses for lung stereotactic body radiation therapy (SBRT). Similar cases to an objective case were defined as cases, which were close to the objective case with respect to the Euclidean distances based on geometrical features. Initial similar-case-based (ISC) plans were generated by applying lung SBRT plans of similar cases to objective cases. Similar cases were selected using the Euclidean distances based on lung shape and geometrical features from a radiation treatment planning database with 174 cases. Beam angles of the ISC plans were optimized using a greedy algorithm based on a cost function to include absorbed doses in the lung and ribs in the WEPL. The 12 dose evaluation indices for the planning target volume, lung, spinal cord, and ribs were evaluated in the original plans, ISC plans, and optimized similar-case-based (OSC) plans with and without WEPL for 20 test cases to investigate its dosimetric impact. These findings revealed that V10 and the mean dose for the lung and V20, V30, and V40 for the ribs in the OSC plan with WEPL improved more significantly than those in the original and ISC plans. This study indicates a potential of similar cases, whose beam angle configurations were optimized with WEPL to avoid lung and rib doses in lung SBRT plans.