Engineering Papers⌕ Search

Engineering topics

Liu, C

Publications and source records attributed to Liu, C.

At least 19 records

Linear magnetoconductivity as a probe of time-reversal symmetry breaking

Several optical experiments have shown that in magnetic materials, the principal axes of response tensors can rotate as an odd function of an applied magnetic field. Here we offer a microscopic explanation of this effect, and we propose a closely related dc transport phenomenon—an off-diagonal symmetric conductivity, linear and odd in a magnetic field, which we refer to as linear magnetoconductivity (LMC). Although LMC has the same functional dependence on a magnetic field as the Hall effect, its origin is fundamentally different: LMC requires time-reversal symmetry to be broken even before a magnetic field is applied, and is therefore a sensitive probe of magnetism. We demonstrate LMC in three different ways: via a tight-binding toy model, a density functional theory calculation on MnPSe3, and a semiclassical treatment. The third approach identifies two distinct mechanisms yielding LMC: momentum-dependent band magnetization and Berry curvature. Finally, we propose an experimental geometry suitable for detecting LMC, and we demonstrate its applicability using Landauer-Büttiker simulations. Our results emphasize the importance of measuring the full conductivity tensor in magnetic materials, and they introduce LMC as a new transport probe of symmetry.

Sunko, V↗

DIRECT RF SAMPLING BASED LLRF CONTROL SYSTEM FOR C-BAND LINEAR ACCELERATOR

Low Level RF (LLRF) control systems of linear accel- erators (LINACs) are typically implemented with hetero- dyne based architectures, which have complex analog RF mixers for up and down conversion. The Gen 3 Radio Fre- quency System-on-Chip (RFSoC) device from AMD Xilinx integrates data converters with maximum RF frequency of 6 GHz. This enables direct RF sampling of C-band LLRF signal typically operated at 5.712 GHz without any analogue mixers, which can significantly simplify the system architec- ture. The data converters sample RF signals in higher order Nyquist zones and then up or down convert digitally by the integrated data path in RFSoC. The closed-loop feedback control firmware implemented in FPGA integrated in RF- SoC can process the base-band signal from the ADC data path and calculate the updated phase and amplitude to be up- mixed by the DAC data path. We have developed a C-band LLRF control RFSoC platform with direct RF sampling, which targets Cool Copper Collider (𝐶3) and other C or S band LINAC research and development projects. In this paper, the architecture of the platform will be described. We have optimized the configuration of the data converter and characterized performance of them with RF pulses. The test results for some of the key performance parameters for the LLRF platform with our custom solid-state amplifier, such as phase and amplitude stability, will be discussed in this paper.

Liu, C↗

SU-C-BRC-04: Efficient Dose Calculation Algorithm for FFF IMRT with a Simplified Bivariate Gaussian Source Model

Purpose: To develop an efficient and accurate independent dose calculation algorithm with a simplified analytical source model for the quality assurance and safe delivery of Flattening Filter Free (FFF)-IMRT on an Elekta Versa HD. Methods: The source model consisted of a point source and a 2D bivariate Gaussian source, respectively modeling the primary photons and the combined effect of head scatter, monitor chamber backscatter and collimator exchange effect. The in-air fluence was firstly calculated by back-projecting the edges of beam defining devices onto the source plane and integrating the visible source distribution. The effect of the rounded MLC leaf end, tongue-and-groove and interleaf transmission was taken into account in the back-projection. The in-air fluence was then modified with a fourth degree polynomial modeling the cone-shaped dose distribution of FFF beams. Planar dose distribution was obtained by convolving the in-air fluence with a dose deposition kernel (DDK) consisting of the sum of three 2D Gaussian functions. The parameters of the source model and the DDK were commissioned using measured in-air output factors (Sc) and cross beam profiles, respectively. A novel method was used to eliminate the volume averaging effect of ion chambers in determining the DDK. Planar dose distributions of five head-and-neck FFF-IMRT plans were calculated and compared against measurements performed with a 2D diode array (MapCHECK™) to validate the accuracy of the algorithm. Results: The proposed source model predicted Sc for both 6MV and 10MV with an accuracy better than 0.1%. With a stringent gamma criterion (2%/2mm/local difference), the passing rate of the FFF-IMRT dose calculation was 97.2±2.6%. Conclusion: The removal of the flattening filter represents a simplification of the head structure which allows the use of a simpler source model for very accurate dose calculation. The proposed algorithm offers an effective way to ensure the safe delivery of FFF-IMRT.

60 APPLIED LIFE SCIENCES↗

SU-C-202-05: Pilot Study of Online Treatment Evaluation and Adaptive Re-Planning for Laryngeal SBRT

Purpose: We have instigated a phase I trial of 5-fraction stereotactic body radiotherapy (SBRT) for advanced-stage laryngeal cancer. We conducted this pilot dosimetric study to confirm the potential utility of online adaptive re-planning to preserve treatment quality. Methods: Ten cases of larynx cancer were evaluated. Baseline and daily SBRT treatment plans were generated per trial protocol. Daily volumetric images were acquired prior to every fraction of treatment. Reference simulation CT images were deformably registered to daily volumetric images using Eclipse. Planning contours were then deformably propagated to daily images. Reference SBRT plans were directly copied to calculate delivered dose distributions on deformed reference CT images. In-house software platform has been developed to calculate cumulative dose over a course of treatment in four steps: 1) deforming delivered dose grid to reference CT images using deformation information exported from Eclipse; 2) generating tetrahedrons using deformed dose grid as vertices; 3) resampling dose to a high resolution within every tetrahedron; 4) calculating dose-volume histograms. Our inhouse software was benchmarked with a commercial software, Mirada. Results: In all ten cases including 49 fractions of treatments, delivered daily doses were completely evaluated and treatment could be re-planned within 10 minutes. Prescription dose coverage of PTV was less than intended in 53% of fractions of treatment (mean: 94%, range: 84%–98%) while minimum coverage of CTV and GTV was 94% and 97%, respectively. Maximum bystander point dose limits to arytenoids, parotids, and spinal cord remained respected in all cases, although variances in carotid artery doses were observed in a minority of cases. Conclusion: Although GTV and CTV coverage is preserved by in-room 3D image guidance of larynx SBRT, PTV coverage can vary significantly from intended plans. Online adaptive treatment evaluation and re-planning is potentially necessary and our procedure is clinically applicable to fully preserve treatment quality. This project is supported by CPRIT Individual Investigator Research Award RP150386.

60 APPLIED LIFE SCIENCES↗

SU-D-204-07: Retrospective Correlation of Dose Accuracy with Regions of Local Failure for Early Stage Lung Cancer Patients Treated with Stereotactic Body Radiotherapy

Purpose: To correlate dose distributions computed using six algorithms for recurrent early stage non-small cell lung cancer (NSCLC) patients treated with stereotactic body radiotherapy (SBRT), with outcome (local failure). Methods: Of 270 NSCLC patients treated with 12Gyx4, 20 were found to have local recurrence prior to the 2-year time point. These patients were originally planned with 1-D pencil beam (1-D PB) algorithm. 4D imaging was performed to manage tumor motion. Regions of local failures were determined from follow-up PET-CT scans. Follow-up CT images were rigidly fused to the planning CT (pCT), and recurrent tumor volumes (Vrecur) were mapped to the pCT. Dose was recomputed, retrospectively, using five algorithms: 3-D PB, collapsed cone convolution (CCC), anisotropic analytical algorithm (AAA), AcurosXB, and Monte Carlo (MC). Tumor control probability (TCP) was computed using the Marsden model (1,2). Patterns of failure were classified as central, in-field, marginal, and distant for Vrecur ≥95% of prescribed dose, 95–80%, 80–20%, and ≤20%, respectively (3). Results: Average PTV D95 (dose covering 95% of the PTV) for 3-D PB, CCC, AAA, AcurosXB, and MC relative to 1-D PB were 95.3±2.1%, 84.1±7.5%, 84.9±5.7%, 86.3±6.0%, and 85.1±7.0%, respectively. TCP values for 1-D PB, 3-D PB, CCC, AAA, AcurosXB, and MC were 98.5±1.2%, 95.7±3.0, 79.6±16.1%, 79.7±16.5%, 81.1±17.5%, and 78.1±20%, respectively. Patterns of local failures were similar for 1-D and 3D PB plans, which predicted that the majority of failures occur in centraldistal regions, with only ∼15% occurring distantly. However, with convolution/superposition and MC type algorithms, the majority of failures (65%) were predicted to be distant, consistent with the literature. Conclusion: Based on MC and convolution/superposition type algorithms, average PTV D95 and TCP were ∼15% lower than the planned 1-D PB dose calculation. Patterns of failure results suggest that MC and convolution/superposition type algorithms predict different outcomes for patterns of failure relative to PB algorithms. Work supported in part by Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗

SU-F-R-38: Impact of Smoothing and Noise On Robustness of CBCT Textural Features for Prediction of Response to Radiotherapy Treatment of Head and Neck Cancers

Purpose: To examine the impact of image smoothing and noise on the robustness of textural information extracted from CBCT images for prediction of radiotherapy response for patients with head/neck (H/N) cancers. Methods: CBCT image datasets for 14 patients with H/N cancer treated with radiation (70 Gy in 35 fractions) were investigated. A deformable registration algorithm was used to fuse planning CT’s to CBCT’s. Tumor volume was automatically segmented on each CBCT image dataset. Local control at 1-year was used to classify 8 patients as responders (R), and 6 as non-responders (NR). A smoothing filter [2D Adaptive Weiner (2DAW) with 3 different windows (ψ=3, 5, and 7)], and two noise models (Poisson and Gaussian, SNR=25) were implemented, and independently applied to CBCT images. Twenty-two textural features, describing the spatial arrangement of voxel intensities calculated from gray-level co-occurrence matrices, were extracted for all tumor volumes. Results: Relative to CBCT images without smoothing, none of 22 textural features extracted showed any significant differences when smoothing was applied (using the 2DAW with filtering parameters of ψ=3 and 5), in the responder and non-responder groups. When smoothing, 2DAW with ψ=7 was applied, one textural feature, Information Measure of Correlation, was significantly different relative to no smoothing. Only 4 features (Energy, Entropy, Homogeneity, and Maximum-Probability) were found to be statistically different between the R and NR groups (Table 1). These features remained statistically significant discriminators for R and NR groups in presence of noise and smoothing. Conclusion: This preliminary work suggests that textural classifiers for response prediction, extracted from H&N CBCT images, are robust to low-power noise and low-pass filtering. While other types of filters will alter the spatial frequencies differently, these results are promising. The current study is subject to Type II errors. A much larger cohort of patients is needed to confirm these results. This work was supported in part by a grant from Varian Medical Systems (Palo Alto, CA)

60 APPLIED LIFE SCIENCES↗

SU-F-R-41: Regularized PCA Can Model Treatment-Related Changes in Head and Neck Patients Using Daily CBCTs

Purpose: To use daily cone beam CTs (CBCTs) to develop regularized principal component analysis (PCA) models of anatomical changes in head and neck (H&N) patients, to guide replanning decisions in adaptive radiation therapy (ART). Methods: Known deformations were applied to planning CT (pCT) images of 10 H&N patients to model several different systematic anatomical changes. A Pinnacle plugin was used to interpolate systematic changes over 35 fractions, generating a set of 35 synthetic CTs for each patient. Deformation vector fields (DVFs) were acquired between the pCT and synthetic CTs and random fraction-to-fraction changes were superimposed on the DVFs. Standard non-regularized and regularized patient-specific PCA models were built using the DVFs. The ability of PCA to extract the known deformations was quantified. PCA models were also generated from clinical CBCTs, for which the deformations and DVFs were not known. It was hypothesized that resulting eigenvectors/eigenfunctions with largest eigenvalues represent the major anatomical deformations during the course of treatment. Results: As demonstrated with quantitative results in the supporting document regularized PCA is more successful than standard PCA at capturing systematic changes early in the treatment. Regularized PCA is able to detect smaller systematic changes against the background of random fraction-to-fraction changes. To be successful at guiding ART, regularized PCA should be coupled with models of when anatomical changes occur: early, late or throughout the treatment course. Conclusion: The leading eigenvector/eigenfunction from the both PCA approaches can tentatively be identified as a major systematic change during radiotherapy course when systematic changes are large enough with respect to random fraction-to-fraction changes. In all cases the regularized PCA approach appears to be more reliable at capturing systematic changes, enabling dosimetric consequences to be projected once trends are established early in the treatment course. This work is supported in part by a grant from Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗

SU-F-J-11: Radiobiologically Optimized Patient Localization During Prostate External Beam Localization

Purpose: To present a novel positioning strategy which optimizes radiation delivery with radiobiological response knowledge, and to evaluate its application during prostate external beam radiotherapy. Methods: Ten patients with low or intermediate risk prostate cancer were evaluated retrospectively in this IRB-approved study. For each patient, a VMAT plan was generated on the planning CT (PCT) to deliver 78 Gy in 39 fractions with PTV = prostate + 7 mm margin, except for 5mm in the posterior direction. Five representative pretreatment CBCT images were selected for each patient, and prostate, rectum, and bladder were delineated on all CBCT images. Each CBCT was auto-registered to the corresponding PCT. Starting from this auto-matched position (AM-position), a search for optimal treatment position was performed utilizing a score function based on radiobiological and dosimetric indices (D98-DTV, NTCP-rectum, and NTCP-bladder) for the daily target volume (DTV), rectum, and bladder. DTV was defined as prostate + 4 mm margin to account for intra-fraction motion as well as contouring variability on CBCT. We termed the optimal treatment position the radiobiologically optimized couch shift position (ROCS-position). Results: The indices, averaged over the 10 patients’ treatment plans, were (mean±SD): 77.7±0.2 Gy (D98-PTV), 12.3±2.7% (NTCP-rectum), and 53.2±11.2% (NTCP-bladder). The corresponding values calculated on all 50 CBCT images at the AM-positions were 72.9±11.3 Gy (D98-DTV), 15.8±6.4% (NTCP-rectum), and 53.0±21.1% (NTCP-bladder), respectively. In comparison, calculated on CBCT at the ROCS-positions, the indices were 77.0±2.1 Gy (D98-DTV), 12.1±5.7% (NTCP-rectum), and 60.7±16.4% (NTCP-bladder). Compared to autoregistration, ROCS-optimization recovered dose coverage to target volume and lowered the risk to rectum. Moreover, NTCPrectum for one patient remained high after ROCS-optimization and therefore could potentially benefit from adaptive planning. Conclusion: These encouraging results illustrate the potential utility of applying radiobiologically optimized correction for online image-guided radiotherapy of prostate patients.

60 APPLIED LIFE SCIENCES↗

SU-F-J-38: Dose Rates and Preliminary Evaluation of Contouring Similarity Metrics Using 4D Cone Beam CT

Purpose: 4D imaging modalities require detailed characterization for clinical optimization. The On-Board Imager mounted on the linear accelerator was used to investigate dose rates in a tissue mimicking phantom using 4D-CBCT and assess variability of contouring similarity metrics between 4D-CT and 4D-CBCT retrospective reconstructions. Methods: A 125 kVp thoracic protocol was used. A phantom placed on a motion platform simulated a patient’s breathing cycle. An ion chamber was affixed inside the phantom’s tissue mimicking cavities (i.e. bone, lung, and soft tissue). A sinusoidal motion waveform was executed with a five second period and superior-inferior motion. Dose rates were measured at six ion chamber positions. A preliminary workflow for contouring similarity between 4D-CT and 4D-CBCT was established using a single lung SBRT patient’s historical data. Average intensity projection (Ave-IP) and maximum intensity projection (MIP) reconstructions generated offline were compared between the 4D modalities. Similarity metrics included Dice similarity coefficient (DSC), Hausdorff distance, and center of mass (COM) deviation. Two isolated lesions were evaluated in the patient’s scans: one located in the right lower lobe (ITVRLL) and one located in the left lower lobe (ITVLLL). Results: Dose rates ranged from 2.30 (lung) to 5.18 (bone) E-3 cGy/mAs. For fixed acquisition parameters, cumulative dose is inversely proportional to gantry speed. For ITVRLL, DSC were 0.70 and 0.68, Hausdorff distances were 6.11 and 5.69 mm, and COM deviations were 1.24 and 4.77 mm, for Ave-IP and MIP respectively. For ITVLLL, DSC were 0.64 and 0.75, Hausdorff distances were 10.74 and 8.00 mm, and COM deviations were 7.55 and 4.3 mm, for Ave-IP and MIP respectively. Conclusion: While the dosimetric output of 4D-CBCT is low, characterization is necessary to assure clinical optimization. A basic workflow for comparison of simulation and treatment 4D image-based contours was established. This work was partially supported by a Research Scholar Grant (RSG-15-137-01-CCE) from the American Cancer Society.

60 APPLIED LIFE SCIENCES↗

SU-F-J-68: Deformable Dose Accumulation for Voxel-Based Dose Tracking of PTV Cold Spots for Adaptive Radiotherapy of the Head and Neck

Purpose: To utilize deformable dose accumulation (DDA) to determine how cold spots within the PTV change over the course of fractionated head and neck (H&N) radiotherapy. Methods: Voxel-based dose was tracked using a DDA platform. The DDA process consisted of B-spline-based deformable image registration (DIR) and dose accumulation between planning CT’s and daily cone-beam CT’s for 10 H&N cancer patients. Cold spots within the PTV (regions receiving less than the prescription, 70 Gy) were contoured on the cumulative dose distribution. These cold spots were mapped to each fraction, starting from the first fraction to determine how they changed. Spatial correlation between cold spot regions over each fraction, relative to the last fraction, was computed using the Jaccard index Jk (Mk,N), where N is the cold spot within the PTV at the end of the treatment, and Mk the same region for fraction k. Results: Figure 1 shows good spatial correlation between cold spots, and highlights expansion of the cold spot region over the course of treatment, as a result of setup uncertainties, and anatomical changes. Figure 2 shows a plot of Jk versus fraction number k averaged over 10 patients. This confirms the good spatial correlation between cold spots over the course of treatment. On average, Jk reaches ∼90% at fraction 22, suggesting that possible intervention (e.g. reoptimization) may mitigate the cold spot region. The cold spot, D99, averaged over 10 patients corresponded to a dose of ∼65 Gy, relative to the prescription dose of 70 Gy. Conclusion: DDA-based tracking provides spatial dose information, which can be used to monitor dose in different regions of the treatment plan, thereby enabling appropriate mid-treatment interventions. This work is supported in part by Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗

SU-F-T-273: Using a Diode Array to Explore the Weakness of TPS DoseCalculation Algorithm for VMAT and Sliding Window Techniques

Purpose: To identify the weakness of dose calculation algorithm in a treatment planning system for volumetric modulated arc therapy (VMAT) and sliding window (SW) techniques using a two-dimensional diode array. Methods: The VMAT quality assurance(QA) was implemented with a diode array using multiple partial arcs that divided from a VMAT plan; each partial arc has the same segments and the original monitor units. Arc angles were less than ± 30°. Multiple arcs delivered through consecutive and repetitive gantry operating clockwise and counterclockwise. The source-toaxis distance setup with the effective depths of 10 and 20 cm were used for a diode array. To figure out dose errors caused in delivery of VMAT fields, the numerous fields having the same segments with the VMAT field irradiated using different delivery techniques of static and step-and-shoot. The dose distributions of the SW technique were evaluated by creating split fields having fine moving steps of multi-leaf collimator leaves. Calculated doses using the adaptive convolution algorithm were analyzed with measured ones with distance-to-agreement and dose difference of 3 mm and 3%.. Results: While the beam delivery through static and step-and-shoot techniques showed the passing rate of 97 ± 2%, partial arc delivery of the VMAT fields brought out passing rate of 85%. However, when leaf motion was restricted less than 4.6 mm/°, passing rate was improved up to 95 ± 2%. Similar passing rate were obtained for both 10 and 20 cm effective depth setup. The calculated doses using the SW technique showed the dose difference over 7% at the final arrival point of moving leaves. Conclusion: Error components in dynamic delivery of modulated beams were distinguished by using the suggested QA method. This partial arc method can be used for routine VMAT QA. Improved SW calculation algorithm is required to provide accurate estimated doses.

60 APPLIED LIFE SCIENCES↗

SU-F-T-511: Feasibility Study of Using Flattening-Filter-Free Photon Beams to Deliver Conventional Flat Beam

Purpose: Various dosimetric benefits such as increased dose rate, and reduced leakage and out of field dose has led to the growth of FFF technology in the clinic. In this study, we concentrate on investigating the feasibility of using flattening-filter-free (FFF) beams to deliver conventional flat beams (CFB), since completely getting rid of the flattening-filter module from the gantry head can not only simplify the gantry design but also decrease the workload on machine maintain and quality assurance. Methods: The sliding window based IMRT technique was utilized to generate the CFB from the FFF beam for various beam configurations on the Elekta Versa HD. The flat beam reproducibility and MU efficiency were compared for each beam configuration among FFF planning, delivery and CFB planning. Results: Compared to the CFB plan, the 3%3mm passing rates of the FFF beams from both measurement and plan are 100% and 95%(or better) for 15×15 cm{sup 2} or smaller field size and for any field size greater than 15×15 cm{sup 2}at 10 cm depth, respectively. The largest discrepancy is about 5% and typically appears at the field shoulder area. The MU increase average was 80% for FFF compared to CFB, however has a minimal impact on treatment delivery time. Conclusion: The ability to deliver conventional flat treatments is not absent when operating in FFF mode. With proper TPS manipulation and beam modulation, FFF mode can achieve reasonable flat profiles and comparable dose coverage as CFB does for various conventional treatment techniques, such as four field box, or long spine treatment techniques. The ability to deliver most clinical treatments from the same treatment unit, will allow for less quality assurance as well as maintenance, and completely eliminate the need for the flattening filter on modern linacs.

60 APPLIED LIFE SCIENCES↗

SU-F-T-518: Development and Characterization of a Gated Treatment System Implemented with An In-House Optical Tracking System and the Elekta Response Interface

Purpose: To report the development and characterization of the first in-house gating system implemented with an optical tracking system (OTS) and the Elekta Response™ interface. Methods: The Response™ connects a patient tracking system with a linac, enabling the tracking system to control radiation delivery. The developed system uses an in-house OTS to monitor patient breathing. The OTS consists of two infrared-based cameras, tracking markers affixed on patient. It achieves gated or breath-held (BH) treatment by calling beam ON/OFF functions in the Response™ dynamic-link library (DLL). A 4D motion phantom was used to evaluate its dosimetric and time delay characteristics. Two FF- and two FFF-IMRT beams were delivered in non-gated, BH and gated mode. The sinusoidal gating signal had a 6 sec period and 15 mm amplitude. The duty cycle included 10%, 20%, 30% and 50%. The BH signal was adapted from the sinusoidal wave by inserting 15 sec BHs. Each delivery was measured with a 2D diode array (MapCHECK™) and compared with the non-gated delivery using gamma analysis (3%). The beam ON/OFF time was captured using the service graphing utility of the linac. Results: The gated treatments were successfully delivered except the 10% duty cycle. The BH delivery had perfect agreement (100%) with non-gated delivery; the agreement of gated delivery decreased from 99% to 88% as duty cycle reduced from 50% to 20%. The beam on/off delay was on average 0.25/0.06 sec. The delivery time for the 50%, 30% and 20% duty cycle increased by 29%, 71% and 139%, respectively. No dosimetric or time delay difference was noticed between FF- and FFF-IMRT beams. Conclusion: The in-house gating system was successfully developed with dosimetric and time delay characteristics in line with published results for commercial systems. It will be an important platform for further research and clinical development of gated treatment.

43 PARTICLE ACCELERATORS↗