Search

Purpose: Historically, spot scanning proton therapy (SSPT) treatment planning utilizes dose volume constraints and linear-energy-transfer (LET) volume constraints separately to balance tumor control and organs-at-risk (OARs) protection. We propose a novel dose-LET volume constraint (DLVC)-based robust optimization (DLVCRO) method for SSPT in treating prostate cancer to obtain a desirable joint dose and LET distribution to minimize adverse events (AEs). Methods: DLVCRO treats DLVC as soft constraints controlling the joint distribution of dose and LET. Ten prostate cancer patients were included with rectum and bladder as OARs. DLVCRO was compared with the conventional robust optimization (RO) method using the worst-case analysis method. Besides the dose-volume histogram (DVH) indices, the analogous LETVH and extra-biological-dose (xBD)-volume histogram indices were also used. The Wilcoxon signed rank test was used to measure statistical significance. Results: In nominal scenario, DLVCRO significantly improved dose, LET and xBD distributions to protect OARs (rectum: V70Gy: 3.07\% vs. 2.90\%, p = .0063, RO vs. DLVCRO; $\text{LET}_{\max}$ (keV/um): 11.53 vs. 9.44, p = .0101; $\text{xBD}_{\max}$ (Gy$\cdot$keV/um): 420.55 vs. 398.79, p = .0086; bladder: V65Gy: 4.82\% vs. 4.61\%, p = .0032; $\text{LET}_{\max}$ 8.97 vs. 7.51, p = .0047; $\text{xBD}_{\max}$ 490.11 vs. 476.71, p = .0641). The physical dose distributions in targets are comparable (D2%: 98.57\% vs. 98.39\%; p = .0805; CTV D2% - D98%: 7.10\% vs. 7.75\%, p = .4624). In the worst-case scenario, DLVCRO robustly enhanced OAR while maintaining the similar plan robustness in target dose coverage and homogeneity. Conclusion: DLVCRO upgrades 2D DVH-based to 3D DLVH-based treatment planning to adjust dose/LET distributions simultaneously and robustly. DLVCRO is potentially a powerful tool to improve patient outcomes in SSPT.

A survey was designed to inquire about the practice of proton SBRT treatment for prostate cancer. The survey was distributed to all 30 proton therapy centers in the United States that participate in the National Clinical Trial Network in Feb. 2023. The survey focused on usage, patient selection criteria, prescriptions, target contours, dose constraints, treatment plan optimization and evaluation methods, patient-specific QA, and IGRT methods. Results: We received responses from 25 centers (83% participation). Only 8 respondent proton centers (32%) reported performing SBRT of the prostate. The remaining 17 centers cited three primary reasons for not offering this treatment: no clinical need, lack of volumetric imaging, and/or lack of clinical evidence. Only 1 center cited the reduction in overall reimbursement as a concern for not offering prostate SBRT. Several common practices among the 8 centers offering SBRT for the prostate were noted, such as using Hydrogel spacers, fiducial markers, and MRI for target delineation. Most proton centers (87.5%) utilized pencil beam scanning (PBS) delivery and completed Imaging and Radiation Oncology Core (IROC) phantom credentialing. Treatment planning typically used parallel opposed lateral beams, and consistent parameters for setup and range uncertainties were used for plan optimization and robustness evaluation. Measurements-based patient-specific QA, beam delivery every other day, fiducial contours for IGRT, and total doses of 35-40 GyRBE were consistent across all centers. However, there was no consensus on the risk levels for patient selection. Conclusion: Prostate SBRT is used in about 1/3 of proton centers in the US. There was a significant consistency in practices among proton centers treating with proton SBRT. It is possible that the adoption of proton SBRT may become more common if proton SBRT is more commonly offered in clinical trials.

Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally-fractionated PBSPT due to concerns of amplified uncertainties at the larger dose per fraction. NRG Oncology and Particle Therapy Cooperative Group (PTCOG) Thoracic Subcommittee surveyed US proton centers to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Amongst other points, the recommendations highlight the need for volumetric image guidance and multiple CT-based robust optimization and robustness tools to minimize further the impact of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques., Comment: 36 pages, 4 figures, 4 tables

We present the Radiation Oncology NLP Database (ROND), the first dedicated Natural Language Processing (NLP) dataset for radiation oncology, an important medical specialty that has received limited attention from the NLP community in the past. With the advent of Artificial General Intelligence (AGI), there is an increasing need for specialized datasets and benchmarks to facilitate research and development. ROND is specifically designed to address this gap in the domain of radiation oncology, a field that offers many opportunities for NLP exploration. It encompasses various NLP tasks including Logic Reasoning, Text Classification, Named Entity Recognition (NER), Question Answering (QA), Text Summarization, and Patient-Clinician Conversations, each with a distinct focus on radiation oncology concepts and application cases. In addition, we have developed an instruction-tuning dataset consisting of over 20k instruction pairs (based on ROND) and trained a large language model, CancerChat. This serves to demonstrate the potential of instruction-tuning large language models within a highly-specialized medical domain. The evaluation results in this study could serve as baseline results for future research. ROND aims to stimulate advancements in radiation oncology and clinical NLP by offering a platform for testing and improving algorithms and models in a domain-specific context. The ROND dataset is a joint effort of multiple U.S. health institutions. The data is available at https://github.com/zl-liu/Radiation-Oncology-NLP-Database., Comment: 10 pages, 7 figures, 6 tables

Purpose: Prior AI-based dose prediction studies in photon and proton therapy often neglect underlying physics, limiting their generalizability to handle outlier clinical cases, especially for pencil beam scanning proton therapy (PBSPT). Our aim is to design a physics-aware and generalizable AI-based PBSPT dose prediction method that has the underlying physics considered to achieve high generalizability to properly handle the outlier clinical cases. Methods and Materials: This study analyzed PBSPT plans of 103 prostate and 78 lung cancer patients from our institution,with each case comprising CT images, structure sets, and plan doses from our Monte-Carlo dose engine (serving as the ground truth). Three methods were evaluated in the ablation study: the ROI-based method, the beam mask and sliding window method, and the noisy probing dose method. Twelve cases with uncommon beam angles or prescription doses tested the methods' generalizability to rare treatment planning scenarios. Performance evaluation used DVH indices, 3D Gamma passing rates (3%/2mm/10%), and dice coefficients for dose agreement. Results: The noisy probing dose method showed improved agreement of DVH indices, 3D Gamma passing rates, and dice coefficients compared to the conventional methods for the testing cases. The noisy probing dose method showed better generalizability in the 6 outlier cases than the ROI-based and beam mask-based methods with 3D Gamma passing rates (for prostate cancer, targets: 89.32%$\pm$1.45% vs. 93.48%$\pm$1.51% vs. 96.79%$\pm$0.83%, OARs: 85.87%$\pm$1.73% vs. 91.15%$\pm$1.13% vs. 94.29%$\pm$1.01%). The dose predictions were completed within 0.3 seconds. Conclusions: We've devised a novel noisy probing dose method for PBSPT dose prediction in prostate and lung cancer patients. With more physics included, it enhances the generalizability of dose prediction in handling outlier clinical cases.

We propose an oAPT workflow that incorporates all these functionalities and validate its clinical implementation feasibility with prostate patients. AI-based auto-segmentation tool AccuContourTM (Manteia, Xiamen, China) was seamlessly integrated into oAPT. Initial spot arrangement tool on the vCT for re-optimization was implemented using raytracing. An LET-based biological effect evaluation tool was developed to assess the overlap region of high dose and high LET in selected OARs. Eleven prostate cancer patients were retrospectively selected to verify the efficacy and efficiency of the proposed oAPT workflow. The time cost of each component in the workflow was recorded for analysis. The verification plan showed significant degradation of the CTV coverage and rectum and bladder sparing due to the interfractional anatomical changes. Re-optimization on the vCT resulted in great improvement of the plan quality. No overlap regions of high dose and high LET distributions were observed in bladder or rectum in re-plans. 3D Gamma analyses in PSQA confirmed the accuracy of the re-plan doses before delivery (Gamma passing rate = 99.57%), and after delivery (98.59%). The robustness of the re-plans passed all clinical requirements. The average time for the complete execution of the workflow was 9.12minutes, excluding manual intervention time. The AI-facilitated oAPT workflow was demonstrated to be both efficient and effective by generating a re-plan that significantly improved the plan quality in prostate cancer treated with PBSPT.

Purpose: To study the effect of proton linear energy transfer (LET) on rib fracture in breast cancer patients treated with pencil-beam scanning proton therapy (PBS) using a novel tool of dose-LET volume histogram (DLVH). Methods: From a prospective registry of patients treated with post-mastectomy proton therapy to the chest wall and regional lymph nodes for breast cancer between 2015 and 2020, we retrospectively identified rib fracture cases detected after completing treatment. Contemporaneously treated control patients that did not develop rib fracture were matched to patients 2:1 considering prescription dose, boost location, reconstruction status, laterality, chest wall thickness, and treatment year. The DLVH index, V(d, l), defined as volume(V) of the structure with at least dose(d) and LET(l), was calculated. DLVH plots between the fracture and control group were compared. Conditional logistic regression (CLR) model was used to establish the relation of V(d, l) and the observed fracture at each combination of d and l. The p-value derived from CLR model shows the statistical difference between fracture patients and the matched control group. Using the 2D p-value map, the DLVH features associated with the patient outcomes were extracted. Results: Seven rib fracture patients were identified, and fourteen matched patients were selected for the control group. The median time from the completion of proton therapy to rib fracture diagnosis was 12 months (range 5 to 14 months). Two patients had grade 2 symptomatic rib fracture while the remaining 5 were grade 1 incidentally detected on imaging. The derived p-value map demonstrated larger V(0-36 Gy[RBE], 4.0-5.0 keV/um) in patients experiencing fracture (p<0.1). Conclusions: In breast cancer patients receiving PBS, a larger volume of chest wall receiving moderate dose and high LET may result in increased risk of rib fracture., Comment: 1 Table and 3 Figures

Purpose: To enhance an in-house graphic-processing-unit (GPU) accelerated virtual particle (VP)-based Monte Carlo (MC) proton dose engine (VPMC) to model aperture blocks in both dose calculation and optimization for pencil beam scanning proton therapy (PBSPT)-based stereotactic radiosurgery (SRS). Methods and Materials: A block aperture module was integrated into VPMC. VPMC was validated by an opensource code, MCsquare, in eight water phantom simulations with 3cm thick brass apertures: four were with aperture openings of 1, 2, 3, and 4cm without a range shifter, while the other four were with same aperture opening configurations with a range shifter of 45mm water equivalent thickness. VPMC was benchmarked with MCsquare and RayStation MC for 10 patients with small targets (average volume 8.4 cc). Finally, 3 patients were selected for robust optimization with aperture blocks using VPMC. Results: In the water phantoms, 3D gamma passing rate (2%/2mm/10%) between VPMC and MCsquare were 99.71$\pm$0.23%. In the patient geometries, 3D gamma passing rates (3%/2mm/10%) between VPMC/MCsquare and RayStation MC were 97.79$\pm$2.21%/97.78$\pm$1.97%, respectively. The calculation time was greatly decreased from 112.45$\pm$114.08 seconds (MCsquare) to 8.20$\pm$6.42 seconds (VPMC), both having statistical uncertainties of about 0.5%. The robustly optimized plans met all the dose-volume-constraints (DVCs) for the targets and OARs per our institutional protocols. The mean calculation time for 13 influence matrices in robust optimization by VPMC was 41.6 seconds. Conclusion: VPMC has been successfully enhanced to model aperture blocks in dose calculation and optimization for the PBSPT-based SRS., Comment: 3 tables, 3 figures

Purpose: To develop a DL-based PBSPT dose prediction workflow with high accuracy and balanced complexity to support on-line adaptive proton therapy clinical decision and subsequent replanning. Methods: PBSPT plans of 103 prostate cancer patients and 83 lung cancer patients previously treated at our institution were included in the study, each with CTs, structure sets, and plan doses calculated by the in-house developed Monte-Carlo dose engine. For the ablation study, we designed three experiments corresponding to the following three methods: 1) Experiment 1, the conventional region of interest (ROI) method. 2) Experiment 2, the beam mask (generated by raytracing of proton beams) method to improve proton dose prediction. 3) Experiment 3, the sliding window method for the model to focus on local details to further improve proton dose prediction. A fully connected 3D-Unet was adopted as the backbone. Dose volume histogram (DVH) indices, 3D Gamma passing rates, and dice coefficients for the structures enclosed by the iso-dose lines between the predicted and the ground truth doses were used as the evaluation metrics. The calculation time for each proton dose prediction was recorded to evaluate the method's efficiency. Results: Compared to the conventional ROI method, the beam mask method improved the agreement of DVH indices for both targets and OARs and the sliding window method further improved the agreement of the DVH indices. For the 3D Gamma passing rates in the target, OARs, and BODY (outside target and OARs), the beam mask method can improve the passing rates in these regions and the sliding window method further improved them. A similar trend was also observed for the dice coefficients. In fact, this trend was especially remarkable for relatively low prescription isodose lines. The dose predictions for all the testing cases were completed within 0.25s.

We present the first study to investigate Large Language Models (LLMs) in answering radiation oncology physics questions. Because popular exams like AP Physics, LSAT, and GRE have large test-taker populations and ample test preparation resources in circulation, they may not allow for accurately assessing the true potential of LLMs. This paper proposes evaluating LLMs on a highly-specialized topic, radiation oncology physics, which may be more pertinent to scientific and medical communities in addition to being a valuable benchmark of LLMs. We developed an exam consisting of 100 radiation oncology physics questions based on our expertise at Mayo Clinic. Four LLMs, ChatGPT (GPT-3.5), ChatGPT (GPT-4), Bard (LaMDA), and BLOOMZ, were evaluated against medical physicists and non-experts. ChatGPT (GPT-4) outperformed all other LLMs as well as medical physicists, on average. The performance of ChatGPT (GPT-4) was further improved when prompted to explain first, then answer. ChatGPT (GPT-3.5 and GPT-4) showed a high level of consistency in its answer choices across a number of trials, whether correct or incorrect, a characteristic that was not observed in the human test groups. In evaluating ChatGPTs (GPT-4) deductive reasoning ability using a novel approach (substituting the correct answer with "None of the above choices is the correct answer."), ChatGPT (GPT-4) demonstrated surprising accuracy, suggesting the potential presence of an emergent ability. Finally, although ChatGPT (GPT-4) performed well overall, its intrinsic properties did not allow for further improvement when scoring based on a majority vote across trials. In contrast, a team of medical physicists were able to greatly outperform ChatGPT (GPT-4) using a majority vote. This study suggests a great potential for LLMs to work alongside radiation oncology experts as highly knowledgeable assistants.

Charge accumulation was measured under calibration conditions in the spread-out Bragg peak (SOBP) using the calibration bias as well as a range of voltages from 10V to 500V and a Farmer-style ion chamber. Collection efficiency was determined by extrapolating to infinite voltage. Similar measurements were taken in an identical dose distribution with a much shorter spot duration. The impact of each of the three models on calibration was then quantified using the TRS-398 protocol. The collection efficiency for the standard calibration was determined to agree well with the prediction of a continuous beam recombination correction. The standard calibration field was found to persistently agree with a continuous beam recombination correction for much lower operating biases. The collection efficiency result for the short spot duration field did not agree with either the continuous or pulsed-beam correction. Using the incorrect recombination model under the standard calibration conditions resulted in a 0.5% calibration difference. We have determined that our spot scanning system would be most appropriately calibrated using a recombination correction with continuous beam model. Physicists responsible for the calibration of such systems are advised to take measurements described here to correctly identify the applicable recombination model for their clinics., Comment: Submitted December 16, 2015 to Medical Physics

We discuss an analytic approach for modeling structure formation in sheets, filaments and knots. This is accomplished by combining models of triaxial collapse with the excursion set approach: sheets are defined as objects which have collapsed along only one axis, filaments have collapsed along two axes, and halos are objects in which triaxial collapse is complete. In the simplest version of this approach, which we develop here, large scale structure shows a clear hierarchy of morphologies: the mass in large-scale sheets is partitioned up among lower mass filaments, which themselves are made-up of still lower mass halos. Our approach provides analytic estimates of the mass fraction in sheets, filaments and halos, and its evolution, for any background cosmological model and any initial fluctuation spectrum. In the currently popular $\Lambda$CDM model, our analysis suggests that more than 99% of the cosmic mass is in sheets, and 72% in filaments, with mass larger than $10^{10} M_{\odot}$ at the present time. For halos, this number is only 46%. Our approach also provides analytic estimates of how halo abundances at any given time correlate with the morphology of the surrounding large-scale structure, and how halo evolution correlates with the morphology of large scale structure., Comment: 22 pages, 7 figures, Accepted for publication in ApJ

Using an eigenspectrum decomposition technique, we separate the host galaxy from the broad line active galactic nucleus (AGN) in a set of 4666 spectra from the Sloan Digital Sky Survey (SDSS), from redshifts near zero up to about 0.75. The decomposition technique uses separate sets of galaxy and quasar eigenspectra to efficiently and reliably separate the AGN and host spectroscopic components. The technique accurately reproduces the host galaxy spectrum, its contributing fraction, and its classification. We show how the accuracy of the decomposition depends upon S/N, host galaxy fraction, and the galaxy class. Based on the eigencoefficients, the sample of SDSS broad-line AGN host galaxies spans a wide range of spectral types, but the distribution differs significantly from inactive galaxies. In particular, post-starburst activity appears to be much more common among AGN host galaxies. The luminosities of the hosts are much higher than expected for normal early-type galaxies, and their colors become increasingly bluer than early-type galaxies with increasing host luminosity. Most of the AGNs with detected hosts are emitting at between 1% and 10% of their estimated Eddington luminosities, but the sensitivity of the technique usually does not extend to the Eddington limit. There are mild correlations among the AGN and host galaxy eigencoefficients, possibly indicating a link between recent star formation and the onset of AGN activity. The catalog of spectral reconstruction parameters is available as an electronic table., Comment: 18 pages; accepted for publication in AJ