Functionally weighted airway sparing (FWAS): a functional avoidance method for preserving post-treatment ventilation in lung radiotherapy.

Recent changes to the guidelines for screening and early diagnosis of lung cancer have increased the interest in preserving post-radiotherapy lung function. Current investigational approaches are based on spatially mapping functional regions and generating regional avoidance plans that preferentially spare highly ventilated/perfused lung. A potentially critical, yet overlooked, aspect of functional avoidance is radiation injury to peripheral airways, which serve as gas conduits to and from functional lung regions. Dose redistribution based solely on regional function may cause irreparable damage to the 'supply chain'. To address this deficiency, we propose the functionally weighted airway sparing (FWAS) method. FWAS (i) maps the bronchial pathways to each functional sub-lobar lung volume; (ii) assigns a weighting factor to each airway based on the relative contribution of the sub-volume to overall lung function; and (iii) creates a treatment plan that aims to preserve these functional pathways. To evaluate it, we used four cases from a retrospective cohort of SAbR patients treated for lung cancer. Each patient's airways were auto-segmented from a diagnostic-quality breath-hold CT using a research virtual bronchoscopy software. A ventilation map was generated from the planning 4DCT to map regional lung function. For each terminal airway, as resolved by the segmentation software, the total ventilation within the sub-lobar volume supported by that airway was estimated and used as a function-based weighting factor. Upstream airways were weighted based on the cumulative volumetric ventilation supported by corresponding downstream airways. Using a previously developed model for airway radiosensitivity, dose constraints were determined for each airway corresponding to a <5% probability of airway collapse. Airway dose constraints, ventilation scores, and clinical dose constraints were input to a swarm optimization-based inverse planning engine to create a 3D conformal SAbR plan (CRT). The FWAS plans were compared to the patients' prescribed CRT clinical plans and the inverse-optimized clinical plans. Depending on the size and location of the tumour, the FWAS plan showed superior preservation of ventilation due to airflow preservation through open pathways (i.e. cumulative ventilation score from the sub-lobar volumes of open pathways). Improvements ranged between 3% and 23%, when comparing to the prescribed clinical plans, and between 3% and 35%, when comparing to the inverse-optimized clinical plans. The three plans satisfied clinical requirements for PTV coverage and OAR dose constraints. These initial results suggest that by sparing pathways to high-functioning lung subregions it is possible to reduce post-SAbR loss of respiratory function.

Quantifying the impact of respiratory-gated 4D CT acquisition on thoracic image quality: a digital phantom study.

PURPOSE: Prospective respiratory-gated 4D CT has been shown to reduce tumor image artifacts by up to 50% compared to conventional 4D CT. However, to date no studies have quantified the impact of gated 4D CT on normal lung tissue imaging, which is important in performing dose calculations based on accurate estimates of lung volume and structure. To determine the impact of gated 4D CT on thoracic image quality, the authors developed a novel simulation framework incorporating a realistic deformable digital phantom driven by patient tumor motion patterns. Based on this framework, the authors test the hypothesis that respiratory-gated 4D CT can significantly reduce lung imaging artifacts. METHODS: Our simulation framework synchronizes the 4D extended cardiac torso (XCAT) phantom with tumor motion data in a quasi real-time fashion, allowing simulation of three 4D CT acquisition modes featuring different levels of respiratory feedback: (i) "conventional" 4D CT that uses a constant imaging and couch-shift frequency, (ii) "beam paused" 4D CT that interrupts imaging to avoid oversampling at a given couch position and respiratory phase, and (iii) "respiratory-gated" 4D CT that triggers acquisition only when the respiratory motion fulfills phase-specific displacement gating windows based on prescan breathing data. Our framework generates a set of ground truth comparators, representing the average XCAT anatomy during beam-on for each of ten respiratory phase bins. Based on this framework, the authors simulated conventional, beam-paused, and respiratory-gated 4D CT images using tumor motion patterns from seven lung cancer patients across 13 treatment fractions, with a simulated 5.5 cm(3) spherical lesion. Normal lung tissue image quality was quantified by comparing simulated and ground truth images in terms of overall mean square error (MSE) intensity difference, threshold-based lung volume error, and fractional false positive/false negative rates. RESULTS: Averaged across all simulations and phase bins, respiratory-gating reduced overall thoracic MSE by 46% compared to conventional 4D CT (p ∼ 10(-19)). Gating leads to small but significant (p < 0.02) reductions in lung volume errors (1.8%-1.4%), false positives (4.0%-2.6%), and false negatives (2.7%-1.3%). These percentage reductions correspond to gating reducing image artifacts by 24-90 cm(3) of lung tissue. Similar to earlier studies, gating reduced patient image dose by up to 22%, but with scan time increased by up to 135%. Beam paused 4D CT did not significantly impact normal lung tissue image quality, but did yield similar dose reductions as for respiratory-gating, without the added cost in scanning time. CONCLUSIONS: For a typical 6 L lung, respiratory-gated 4D CT can reduce image artifacts affecting up to 90 cm(3) of normal lung tissue compared to conventional acquisition. This image improvement could have important implications for dose calculations based on 4D CT. Where image quality is less critical, beam paused 4D CT is a simple strategy to reduce imaging dose without sacrificing acquisition time.

Application of Doppler spectroscopy in H2 to the prediction of experimental D(d,n)3He reaction rates in an inertial electrostatic confinement device.

In our previous paper, we developed an optical measurement for absolute densities of fast (tens of keV) H ions in an abnormal hollow cathode discharge of hydrogen in the units to tens of mTorr pressure range. We apply this method to a cylindrically symmetric inertial electrostatic confinement (IEC) discharge of hydrogen using the Doppler spectrum of Halpha. We predict neutron production rates for an equivalent discharge of deuterium and compare these with experimental values under the separate conditions of constant deuterium gas pressure (approximately 6 mTorr) and voltage (30 kV). Our predictions capture the variation of production rates with dc current (10-50 mA) and agree with experiment to within an order of magnitude. The applicability of this diagnostic to the cylindrical IEC discharge is thus demonstrated, with our results supporting the theory that the discharge is dominated by energetic neutrals emerging from the cathode apertures.

Absolute densities of energetic hydrogen ion species in an abnormal hollow cathode discharge.

We develop an optical measurement for densities of fast (units to tens of keV) hydrogen ions in an abnormal hollow cathode discharge in units to tens of mTorr pressure range. This method combines results from previous collisional-radiative models, comparing the intensity of Balmer H_{alpha} due to dissociative excitation of H2 by fast electrons to Doppler-shifted emission arising from charge exchange of energetic ions. The method requires only two inputs: the current density due to fast electrons and a single H_{alpha} spectrum of the characteristic emission channel at the anode. We model in particular the cylindrical interelectrode discharge of an inertial electrostatic confinement device. Experimentally, we find that the density of fast ions emerging from the cathode (bias -5 kV at 20 mTorr ) is in the order 10;{14}m;{-3} , increasing approximately linearly with current in the 10-30 mA range. Calculated densities agree with values obtained in similar apparatus using Langmuir probes and analysis of dust particle motion.