Search NASA⌕ Search

SEARCH · Search NASA

Results for “PMB”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Scattering Observables from One- and Two-body Densities: Formalism and Application to $\pmb \gamma $ ${}^3\hbox {He}$ Scattering

We introduce the transition-density formalism, an efficient and general method for calculating the interaction of external probes with light nuclei. One- and two-body transition densities that encode the nuclear structure of the target are evaluated once and stored. They are then convoluted with an interaction kernel to produce amplitudes, and hence observables. Here, by choosing different kernels, the same densities can be used for any reaction in which a probe interacts perturbatively with the target. The method therefore exploits the factorisation between nuclear structure and interaction kernel that occurs in such processes. We study in detail the convergence in the number of partial waves for matrix elements relevant in elastic Compton scattering on 3 He. The results are fully consistent with our previous calculations in Chiral Effective Field Theory. But the new approach is markedly more computationally efficient, which facilitates the inclusion of more partial-wave channels in the calculation. We also discuss the usefulness of the transition-density method for other nuclei and reactions. Calculations of elastic Compton scattering on heavier targets like 4 He are straightforward extensions of this study, since the same interaction kernels are used. And the generality of the formalism means that our 3 He densities can be used to evaluate any 3 He elastic-scattering observable with contributions from one- and two-body operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Charged-particle pseudorapidity density at mid-rapidity in p–Pb collisions at $\pmb {\sqrt{s_{\scriptscriptstyle {\mathrm{NN}}}}}$ = 8.16 TeV

The pseudorapidity density of charged particles, d$N$ ch /d$η$, in p–Pb collisions has been measured at a centre-of-mass energy per nucleon–nucleon pair of $\sqrt{s_{NN}}$ = 8.16 TeV at mid-pseudorapidity for non-single-diffractive events. The results cover 3.6 units of pseudorapidity, |$η$| < 1.8. The d$N$ ch /d$η$ value is 19.1 ± 0.7 at |$η$| < 0.5. This quantity divided by (N part )/2 is 4.73 ± 0.20 , where (N part ) is the average number of participating nucleons, is 9.5% higher than the corresponding value for p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Measurements are compared with models based on different mechanisms for particle production. All models agree within uncertainties with data in the Pb-going side, while HIJING overestimates, showing a symmetric behaviour, and EPOS underestimates the p-going side of the d$N$ ch /d$η$ distribution. Saturation-based models reproduce the distributions well for < -1.3. The d$N$ ch /d$η$ is also measured for different centrality estimators, based both on the charged-particle multiplicity and on the energy deposited in the Zero-Degree Calorimeters. A study of the implications of the large multiplicity fluctuations due to the small number of participants for systems like p–Pb in the centrality calculation for multiplicity-based estimators is discussed, demonstrating the advantages of determining the centrality with energy deposited near beam rapidity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Azimuthal correlations of heavy-flavor hadron decay electrons with charged particles in pp and p–Pb collisions at $\pmb {\sqrt{s_{\mathrm{{NN}}}}}$ = 5.02 TeV

The azimuthal (Δφ ) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p–Pb collisions at $\sqrt{s_{NN}}$=5.02 TeV. Results are reported for electrons with transverse momentum 4 < p T < 16 GeV/c and pseudorapidity |η| <0.6. The associated charged particles are selected with transverse momentum 1< p T < 7 GeV/c, and relative pseudorapidity separation with the leading electron |Δη| <1. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p–Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The Δφ distribution and the peak observables in pp and p–Pb collisions are compared with calculations from various Monte Carlo event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Measurements of $\mathbf{t\bar t} \pmb{\gamma}$ , $\mathbf{t\bar t + Z}$, and $\mathbf{t\bar t + W}$ production

Measurements of the production of top quark pairs in association withelectroweak vector bosons from the CMS, ATLAS and CDF collaborationsare presented. The CDF collaboration measures the cross section oftop pair plus photon production and the ratio of this cross sectionthe to inclusive top pair cross section using 6 fb$^{-1}$ of dataat $\sqrt{s} = 1.96$ TeV with the results $\sigma_{t\bar{t}\gamma}= 0.18 \pm 0.07 (stat.) \pm 0.04 (sys.) \pm 0.01 (lumi.)$ pb and $\sigma_{t\bar{t}\gamma}/\sigma_{t\bar{t}= 0.024 \pm 0.009}$. The ATLAS collaboration measures the top pairplus photon cross section at $\sqrt{s} = 7$ TeV using 1.04 fb $^{-1}$of data with the result $\sigma_{t\bar{t}\gamma} = 2.0 \pm 0.5 (stat.)\pm 0.7 (sys.) \pm 0.08 (lumi.)$. The ATLAS collaboration also setsan upper limit on $\sigma_{t\bar{t}Z}$ of 0.71 pb using 4.0 fb $^{-1}$of data at $\sqrt{s} = 7$ TeV. The CMS collaboration measures the$\sigma_{t\bar{t}Z}$ and $\sigma_{t\bar{t}V (V= W, Z)}$ cross sectionsusing 5.0 fb $^{-1}$ of data at $\sqrt{s} = 7$ TeV with the results$\sigma_{t\bar{t}V (V= W, Z)} = 0.43^{+0.17}_{-0.15} (stat.) ^{+0.09}_{-0.07}(sys.)$ pb and $\sigma_{t\bar{t}Z} = 0.28^{+0.14}_{-0.11} (stat.)^{+0.06}_{-0.03} (sys.)$ pb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Desert dust induces TLR signaling to trigger Th2-dominant lung allergic inflammation via a MyD88-dependent signaling pathway

Asian sand dust (ASD) is known to exacerbate asthma, although its mechanism is not yet well understood. In this study, when the effects on inflammatory response by LPS present in ASD was investigated by measuring the gene expression of cytokines and chemokines in RAW264.7 cells treated with ASD and/or polymyxin B (PMB), the ASD effects were attenuated by PMB, but not completely. When an in vitro study was performed using bone marrow-derived macrophages (BMDMs) from WT, TLR2{sup −/−}, TLR4{sup −/−}, and MyD88{sup −/−} BALB/c mice and BMDMs from WT, TLR2{sup −/−}, TLR4{sup −/−}, TLR2/4{sup −/−}, TLR7/9{sup −/−}, and MyD88{sup −/−} C57BL/6J mice, cytokine (IL-6, IL-12) production in BMDMs was higher in ASD-stimulated TLR2{sup −/−} cells than in TLR4{sup −/−} cells, whereas it was lower or undetectable in TLR2/4{sup −/−} and MyD88{sup −/−} cells. These results suggest that ASD causes cytokine production predominantly in a TLR4/MyD88-dependent pathway. When WT and TLRs 2{sup −/−}, 4{sup −/−}, and MyD88{sup −/−} BALB/c mice were intratracheally challenged with OVA and/or ASD, ASD caused exacerbation of lung eosinophilia along with Th2 cytokine and eosinophil-relevant chemokine production. Serum OVA-specific IgE and IgG1 similar to WT was observed in TLRs 2{sup −/−}, 4{sup −/−} mice, but not in MyD88{sup −/−} mice. The Th2 responses in TLR2{sup −/−} mice were attenuated remarkably by PMB. These results indicate that ASD exacerbates lung eosinophilia in a MyD88-dependent pathway. TLRs 2 and 4 signaling may be important in the increase in lung eosinophilia. Also, the TLR4 ligand LPS and TLR2 ligand like β-glucan may be strong candidates for exacerbation of lung eosinophilia. - Highlights: • ASD enhanced Th2 response in TLR2{sup −/−}, TLR4{sup −/−} and WT mice, but not in MyD88{sup −/−}. • Th2 responses in TLR2{sup −/−} mice were attenuated by LPS inhibitor polymyxin B. • TLR2 and TLR4 signaling is important in allergic lung disease aggravation by ASD. • MyD88 is the key adapter molecule in the signaling pathway for Th2 activation. • ASD-bound LPS and β-glucan may be strong candidates for the aggravating substances.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Weight-adapted ultra-low-dose pancreatic perfusion CT: radiation dose, image quality, and perfusion parameters

Purpose: We evaluate the reliability and feasibility of weight-adapted ultra-low-dose pancreatic perfusion CT. Methods: A total of 100 (47 men, 53 women) patients were enrolled prospectively and were assigned to five groups (A, B, C, D, and E) with different combination of tube voltage and tube current according to their body weight. Radiation dose parameters including volume CT dose index (CTDI) and dose-length product (DLP) were recorded. Image quality was evaluated both subjectively and objectively (noise, signal-to-noise ratio, contrast-to-noise ratio). Perfusion parameters including blood flow (BF), blood volume (BV), and permeability (PMB) were measured. The dose, image quality measurements, and perfusion parameters were compared between the five groups using one-way analysis of variance (ANOVA). Results: Radiation dose reached 8.7 mSv in patients under 50 kg and was 18.9 mSv in patients above 80 kg. The mean subjective image quality score was above 4.45 on a 5-point scale with good agreement between two radiologists. Groups A–D had equivalent performance on objective image quality (P > 0.05), while Group E performed even better (P < 0.05). No significant differences emerged in comparison with perfusion parameters (BF, BV, PMB) of normal pancreas parenchyma between the five groups. Conclusion: Weight-adapted ultra-low-dose pancreatic perfusion CT can effectively reduce radiation dose without prejudice to image quality, and the perfusion parameters of normal parenchyma are accurate and reliable.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Peri-Net-Pro: the neural processes with quantified uncertainty for crack patterns

Abstract This paper develops a deep learning tool based on neural processes (NPs) called the Peri-Net-Pro, to predict the crack patterns in a moving disk and classifies them according to the classification modes with quantified uncertainties. In particular, image classification and regression studies are conducted by means of convolutional neural networks (CNNs) and NPs. First, the amount and quality of the data are enhanced by using peridynamics to theoretically compensate for the problems of the finite element method (FEM) in generating crack pattern images. Second, case studies are conducted with the prototype microelastic brittle (PMB), linear peridynamic solid (LPS), and viscoelastic solid (VES) models obtained by using the peridynamic theory. The case studies are performed to classify the images by using CNNs and determine the suitability of the PMB, LBS, and VES models. Finally, a regression analysis is performed on the crack pattern images with NPs to predict the crack patterns. The regression analysis results confirm that the variance decreases when the number of epochs increases by using the NPs. The training results gradually improve, and the variance ranges decrease to less than 0.035. The main finding of this study is that the NPs enable accurate predictions, even with missing or insufficient training data. The results demonstrate that if the context points are set to the 10th, 100th, 300th, and 784th, the training information is deliberately omitted for the context points of the 10th, 100th, and 300th, and the predictions are different when the context points are significantly lower. However, the comparison of the results of the 100th and 784th context points shows that the predicted results are similar because of the Gaussian processes in the NPs. Therefore, if the NPs are employed for training, the missing information of the training data can be supplemented to predict the results.

Mathematics↗
DOE OSTI2019

Simulating the Fracture of Notched Mortar Beams through Extended Finite-Element Method and Peridynamics

This paper simulates fracture in notched mortar beams under three-point bending using extended finite element method (XFEM) and peridynamics. A three-phase microstructure (i.e., cement paste, aggregates, and paste-aggregate interface) is used for constitutive modeling of the mortar to obtain the elastic properties for simulation. In the XFEM approach, the simulated homogenized elastic modulus is used along with the total fracture energy of the cement mortar in a damage model to predict the fracture response of the mortar including crack propagation and its fracture parameters (Mode I stress intensity factor, K IC and critical crack tip opening displacement, CTOD C ). The damage model incorporates a maximum principal stress-based damage initiation criteria and a traction-separation law for damage evolution. In the peridynamics approach, a bond-based model involving a prototype microelastic brittle (PMB) material model is used. The elastic properties and fracture energy release rates are used as inputs in the PMB model, along with the choice of peridynamic horizon size. Comparison with experimental fracture properties ( K IC , CTOD C ) as well as crack propagation paths from digital image correlation show that both the approaches yield satisfactory results, particularly for K IC and crack extension. Furthermore, both these methods can be adopted for fracture simulation of cement-based materials.

42 ENGINEERING↗
DOE OSTI2016

Effective of the q-deformed pseudoscalar magnetic field on the charge carriers in graphene

In this paper, we have obtained exact analytical solutions of the time-independent Dirac-Weyl equation for the charge carriers with q-deformed pseudoscalar magnetic barrier (PMB) in graphene by using the ansatz method. We have also found a solution that describes the left propagating wave function to calculation of the reflection and transmission coefficients using the Riemann’s equation. This allows us to conclude about the Dirac-Weyl equation with PMB and to understand quantum behavior of the Dirac fermions. Finally, some of the numerical results are shown, too.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Experimental and numerical investigation of flame stabilization and pollutant formation in matrix stabilized ammonia-hydrogen combustion

Ammonia (NH 3 ) is a carbon-free fuel that offers an attractive alternative for reducing greenhouse gas emissions. However, the slow flame speed, low heating value, and emissions of nitrogen-containing pollutants present significant issues for practical combustion applications. To address these issues, we investigate the use of matrix stabilized combustion. In this type of burner, combustion is performed within an inert porous ceramic foam, heat is recirculated by solid conduction and radiation, which enhances flame speed and combustion stabilization, thereby permitting combustion over a wide range of equivalence ratio conditions. We present a new porous media burner (PMB) capable of stabilizing NH 3 /air flames at ambient conditions. An extensive experimental characterization of the stability of this burner is conducted with up to 30% by volume of hydrogen (H 2 ) in the fuel stream. A 15:1 turndown ratio is demonstrated, with a high thermal power density of 62 MW m -3 . Concentrations of NO, unburnt NH 3 , and H 2 in the exhaust stream are measured. Two regimes are identified for low NO operation: rich and very lean. For rich conditions, NO emissions decrease with increasing equivalence ratio and decreasing H 2 blending. Unburnt NH 3 emissions follow opposite trends. These measurements are complemented by simulations in which the burner is represented by a coupled solid-gas reactor network. In conclusion, this model captures the burner’s pollutant emissions to good accuracy and is used to analyze the mechanisms of pollutant formation.

42 ENGINEERING↗
DOE OSTI2017

Automated glycan assembly of galactosylated xyloglucan oligosaccharides and their recognition by plant cell wall glycan-directed antibodies

We report the automated glycan assembly of oligosaccharides related to the plant cell wall hemicellulosic polysaccharide xyloglucan. The synthesis of galactosylated xyloglucan oligosaccharides was enabled by introducing p-methoxybenzyl (PMB) as a temporary protecting group for automated glycan assembly. The generated oligosaccharides were printed as microarrays, and the binding of a collection of xyloglucan-directed monoclonal antibodies (mAbs) to the oligosaccharides was assessed. We also demonstrated that the printed glycans can be further enzymatically modified while appended to the microarray surface by Arabidopsis thaliana xyloglucan xylosyltransferase 2 (AtXXT2).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

SU-E-J-159: Intra-Patient Deformable Image Registration Uncertainties Quantified Using the Distance Discordance Metric

Purpose: The quantitative evaluation of deformable image registration (DIR) is currently challenging due to lack of a ground truth. In this study we test a new method proposed for quantifying multiple-image based DIRrelated uncertainties, for DIR of pelvic images. Methods: 19 patients were analyzed, each with 6 CT scans, who previously had radiotherapy for prostate cancer. Manually delineated structures for rectum and bladder, which served as ground truth structures, were delineated on the planning CT and each subsequent scan. For each patient, voxel-by-voxel DIR-related uncertainties were evaluated, following B-spline based DIR, by applying a previously developed metric, the distance discordance metric (DDM; Saleh et al., PMB (2014) 59:733). The DDM map was superimposed on the first acquired CT scan and DDM statistics were assessed, also relative to two metrics estimating the agreement between the propagated and the manually delineated structures. Results: The highest DDM values which correspond to greatest spatial uncertainties were observed near the body surface and in the bowel due to the presence of gas. The mean rectal and bladder DDM values ranged from 1.1–11.1 mm and 1.5–12.7 mm, respectively. There was a strong correlation in the DDMs between the rectum and bladder (Pearson R = 0.68 for the max DDM). For both structures, DDM was correlated with the ratio between the DIR-propagated and manually delineated volumes (R = 0.74 for the max rectal DDM). The maximum rectal DDM was negatively correlated with the Dice Similarity Coefficient between the propagated and the manually delineated volumes (R= −0.52). Conclusion: The multipleimage based DDM map quantified considerable DIR variability across different structures and among patients. Besides using the DDM for quantifying DIR-related uncertainties it could potentially be used to adjust for uncertainties in DIR-based accumulated dose distributions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-D-BRE-04: Modeling Optimal Concurrent Chemotherapy Schedules

Purpose: Concurrent chemo-radiation therapy (CCRT) has become a more common cancer treatment option with a better tumor control rate for several tumor sites, including head and neck and lung cancer. In this work, possible optimal chemotherapy schedules were investigated by implementing chemotherapy cell-kill into a tumor response model of RT. Methods: The chemotherapy effect has been added into a published model (Jeong et al., PMB (2013) 58:4897), in which the tumor response to RT can be simulated with the effects of hypoxia and proliferation. Based on the two-compartment pharmacokinetic model, the temporal concentration of chemotherapy agent was estimated. Log cell-kill was assumed and the cell-kill constant was estimated from the observed increase in local control due to concurrent chemotherapy. For a simplified two cycle CCRT regime, several different starting times and intervals were simulated with conventional RT regime (2Gy/fx, 5fx/wk). The effectiveness of CCRT was evaluated in terms of reduction in radiation dose required for 50% of control to find the optimal chemotherapy schedule. Results: Assuming the typical slope of dose response curve (γ50=2), the observed 10% increase in local control rate was evaluated to be equivalent to an extra RT dose of about 4 Gy, from which the cell-kill rate of chemotherapy was derived to be about 0.35. Best response was obtained when chemotherapy was started at about 3 weeks after RT began. As the interval between two cycles decreases, the efficacy of chemotherapy increases with broader range of optimal starting times. Conclusion: The effect of chemotherapy has been implemented into the resource-conservation tumor response model to investigate CCRT. The results suggest that the concurrent chemotherapy might be more effective when delayed for about 3 weeks, due to lower tumor burden and a larger fraction of proliferating cells after reoxygenation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-04: Integration of Interstitial High Intensity Therapeutic Ultrasound Applicators On a Clinical MRI-Guided High Intensity Focused Ultrasound Treatment Planning Software Platform

Purpose: Interstitial high intensity therapeutic ultrasound (HITU) applicators can be used to ablate tissue percutaneously, allowing for minimally-invasive treatment without ionizing radiation [1,2]. The purpose of this study was to evaluate the feasibility and usability of combining multielement interstitial HITU applicators with a clinical magnetic resonance imaging (MRI)-guided focused ultrasound software platform. Methods: The Sonalleve software platform (Philips Healthcare, Vantaa, Finland) combines anatomical MRI for target selection and multi-planar MRI thermometry to provide real-time temperature information. The MRI-compatible interstitial US applicators (Acoustic MedSystems, Savoy, IL, USA) had 1–4 cylindrical US elements, each 1 cm long with either 180° or 360° of active surface. Each applicator (4 Fr diameter, enclosed within a 13 Fr flexible catheter) was inserted into a tissue-mimicking agar-silica phantom. Degassed water was circulated around the transducers for cooling and coupling. Based on the location of the applicator, a virtual transducer overlay was added to the software to assist targeting and to allow automatic thermometry slice placement. The phantom was sonicated at 7 MHz for 5 minutes with 6–8 W of acoustic power for each element. MR thermometry data were collected during and after sonication. Results: Preliminary testing indicated that the applicator location could be identified in the planning images and the transducer locations predicted within 1 mm accuracy using the overlay. Ablation zones (thermal dose ≥ 240 CEM43) for 2 active, adjacent US elements ranged from 18 mm × 24 mm (width × length) to 25 mm × 25 mm for the 6 W and 8 W sonications, respectively. Conclusion: The combination of interstitial HITU applicators and this software platform holds promise for novel approaches in minimally-invasive MRI-guided therapy, especially when bony structures or air-filled cavities may preclude extracorporeal HIFU.[1] Diederich et al. IEEE UFFFC 46.5 (1999): 1218.[2] Chopra et al. PMB 50.21 (2005): 4957. Funding support was provided by Philips Healthcare and in-kind support from Acoustic MedSystems Inc. Ari Partanen is a paid employee of Philips Healthcare. Goutam Ghoshal and Everette Clif Burdette are paid employees of Acoustic MedSystems Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-R-48: Early Prediction of Pathological Response of Locally Advanced Rectal Cancer Using Perfusion CT:A Prospective Clinical Study

Purpose: To prospectively evaluate the tumor vascularity assessed by perfusion CT for prediction of chemo-radiation treatment (CRT) response in locally advanced rectal cancer (LARC). Methods: Eighteen consecutive patients (61.9±8.8 years, from March–June 2015) diagnosed with LARC who underwent 6–8 weeks CRT followed by surgery were included. The pre-treatment perfusion CT was acquired after a 5s delay of contrast agent injection for 45s with 1s interval. A total of 7-cm craniocaudal range covered the tumor region with 3-mm slice thickness. The effective radiation dose is around 15mSv, which is about 1.5 the conventional abdomen/pelvis CT dose. The parametric map of blood flow (BF), blood volume (BV), mean transit time (MTT), permeability (PMB), and maximum intensity map (MIP) were obtained from commercial software (Syngo-CT 2011A, Siemens). An experienced radiation oncologist outlined the tumor based on the pre-operative MR and pathologic residual region, but was blinded with regards to pathological tumor stage. The perfusion parameters were compared to histopathological response quantified by tumor regression grade as good-responder (GR, TRG 0-1) vs. non-good responder (non-GR). Furthermore, the predictive value for pathological complete response (pCR) was also investigated. Results: Both BV (p=0.02) and MTT (P=0.02) was significantly higher and permeambility was lower (p=0.004) in the good responders. The BF was higher in GR group but not statistically significant. Regarding the discrimination of pCR vs non-pCR, the BF was higher in the pCR group (p=0.08) but none of those parameters showed statistically significant differences. Conclusion: BV and MTT can discriminate patients with a favorable response from those that fail to respond well, potentially selecting high-risk patients with resistant tumors that may benefit from an aggressive preoperative treatment approach. However, future studies with more patient data are needed to verify the prognostic value of perfusion CT especially for pCR prediction. This work is supported by the National High-tech R&D program for Young Scientists by the Ministry of Science and Technology of China (Grant No. 2015AA020917), Natural Science Foundation of China (NSFC Grant No. 81201091).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-FG-CAMPUS-TeP3-02: Benchmarks of a Proton Relative Biological Effectiveness (RBE) Model for DNA Double Strand Break (DSB) Induction in the FLUKA, MCNP, TOPAS, and RayStation™ Treatment Planning System

Purpose: Validate implementation of a published RBE model for DSB induction (RBEDSB) in several general purpose Monte Carlo (MC) code systems and the RayStation™ treatment planning system (TPS). For protons and other light ions, DSB induction is a critical initiating molecular event that correlates well with the RBE for cell survival. Methods: An efficient algorithm to incorporate information on proton and light ion RBEDSB from the independently tested Monte Carlo Damage Simulation (MCDS) has now been integrated into MCNP (Stewart et al. PMB 60, 8249–8274, 2015), FLUKA, TOPAS and a research build of the RayStation™ TPS. To cross-validate the RBEDSB model implementation LET distributions, depth-dose and lateral (dose and RBEDSB) profiles for monodirectional monoenergetic (100 to 200 MeV) protons incident on a water phantom are compared. The effects of recoil and secondary ion production ({sub 2}H{sub +}, {sub 3}H{sub +}, {sub 3}He{sub 2+}, {sub 4}He{sub 2+}), spot size (3 and 10 mm), and transport physics on beam profiles and RBEDSB are examined. Results: Depth-dose and RBEDSB profiles among all of the MC models are in excellent agreement using a 1 mm distance criterion (width of a voxel). For a 100 MeV proton beam (10 mm spot), RBEDSB = 1.2 ± 0.03 (− 2–3%) at the tip of the Bragg peak and increases to 1.59 ± 0.3 two mm distal to the Bragg peak. RBEDSB tends to decrease as the kinetic energy of the incident proton increases. Conclusion: The model for proton RBEDSB has been accurately implemented into FLUKA, MCNP, TOPAS and the RayStation™TPS. The transport of secondary light ions (Z > 1) has a significant impact on RBEDSB, especially distal to the Bragg peak, although light ions have a small effect on (dosexRBEDSB) profiles. The ability to incorporate spatial variations in proton RBE within a TPS creates new opportunities to individualize treatment plans and increase the therapeutic ratio. Dr. Erik Traneus is employed full-time as a Research Scientist at RaySearch Laboratories. The research build of the RayStation used in the study was made available to the University of Washington free of charge. RaySearch Laboratories did not provide any monetary support for the reported studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-FG-BRB-03: Basis Vector Model Based Method for Proton Stopping Power Estimation From Experimental Dual Energy CT Data

Purpose: This work aims at reducing the uncertainty in proton stopping power (SP) estimation by a novel combination of a linear, separable basis vector model (BVM) for stopping power calculation (Med Phys 43:600) and a statistical, model-based dual-energy CT (DECT) image reconstruction algorithm (TMI 35:685). The method was applied to experimental data. Methods: BVM assumes the photon attenuation coefficients, electron densities, and mean excitation energies (I-values) of unknown materials can be approximated by a combination of the corresponding quantities of two reference materials. The DECT projection data for a phantom with 5 different known materials was collected on a Philips Brilliance scanner using two scans at 90 kVp and 140 kVp. The line integral alternating minimization (LIAM) algorithm was used to recover the two BVM coefficient images using the measured source spectra. The proton stopping powers are then estimated from the Bethe-Bloch equation using electron densities and I-values derived from the BVM coefficients. The proton stopping powers and proton ranges for the phantom materials estimated via our BVM based DECT method are compared to ICRU reference values and a post-processing DECT analysis (Yang PMB 55:1343) applied to vendorreconstructed images using the Torikoshi parametric fit model (tPFM). Results: For the phantom materials, the average stopping power estimations for 175 MeV protons derived from our method are within 1% of the ICRU reference values (except for Teflon with a 1.48% error), with an average standard deviation of 0.46% over pixels. The resultant proton ranges agree with the reference values within 2 mm. Conclusion: Our principled DECT iterative reconstruction algorithm, incorporating optimal beam hardening and scatter corrections, in conjunction with a simple linear BVM model, achieves more accurate and robust proton stopping power maps than the post-processing, nonlinear tPFM based DECT analysis applied to conventional reconstructions of low and high energy scans. Funding Support: NIH R01CA 75371; NCI grant R01 CA 149305.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-202-11: Radiobiological Modeling of Tumor Response During Radiotherapy Based On Pre-Treatment Dynamic PET Imaging Data

Purpose: To evaluate the ability of a multiscale radiobiological model of tumor response to predict mid-treatment hypoxia images, based on pretreatment imaging of perfusion and hypoxia with [18-F]FMISO dynamic PET and glucose metabolism with [18-F]FDG PET. Methods: A mechanistic tumor control probability (TCP) radiobiological model describing the interplay between tumor cell proliferation and hypoxia (Jeong et al., PMB 2013) was extended to account for intra-tumor nutrient heterogeneity, dynamic cell migration due to nutrient gradients, and stromal cells. This extended model was tested on 10 head and neck cancer patients treated with chemoradiotherapy, randomly drawn from a larger MSKCC protocol involving baseline and mid-therapy dynamic PET scans. For each voxel, initial fractions of proliferative and hypoxic tumor cells were obtained by finding an approximate solution to a system of linear equations relating cell fractions to voxel-level FDG uptake, perfusion (FMISO K{sub 1}) and hypoxia (FMISO k{sub 3}). The TCP model then predicted their evolution over time up until the mid treatment scan. Finally, the linear model was reapplied to predict each lesion’s median hypoxia level (k{sub 3}[med,sim]) which in turn was compared to the FMISO k{sub 3}[med] measured at mid-therapy. Results: The average k3[med] of the tumors in pre-treatment scans was 0.0035 min{sup −1}, with an inter-tumor standard deviation of σ[pre]=0.0034 min{sup −1}. The initial simulated k{sub 3}[med,sim] of each tumor agreed with the corresponding measurements within 0.1σ[pre]. In 7 out of 10 lesions, the mid-treatment k{sub 3}[med,sim] prediction agreed with the data within 0.3σ[pre]. The remaining cases corresponded to the most extreme relative changes in k{sub 3}[med]. Conclusion: This work presents a method to personalize the prediction of a TCP model using pre-treatment kinetic imaging data, and validates the modeling of radiotherapy response by predicting changes in median hypoxia values during treatment. Variations from predicted response may be a useful biomarker, which should be further explored. Partially supported by NIH grant #1 R01 CA157770-01A1 and a grant from Varian Corporation.

60 APPLIED LIFE SCIENCES↗