Search NASA⌕ Search

SEARCH · Search NASA

Results for “Applied Probability”

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 361 records · Page 20

Centralized Cryptographic Key Management and Critical Risk Assessment - (CRADA Final Report)

The Department of Energy Office of Electricity Delivery and Energy Reliability (DOE-OE) Cyber Security for Energy Delivery Systems (CSEDS) industry led program (DE-FOA-0000359) entitled "Innovation for Increasing Cyber Security for Energy Delivery Systems (12CSEDS)," awarded a contract to Sypris Electronics LLC to develop a Cryptographic Key Management System for the smart grid (Scalable Key Management Solutions for Critical Infrastructure Protection). Oak Ridge National Laboratory (ORNL) and Sypris Electronics, LLC as a result of that award entered into a CRADA (NFE-11-03562) between ORNL and Sypris Electronics, LLC. ORNL provided its Cyber Security Econometrics System (CSES) as a tool to be modified and used as a metric to address risks and vulnerabilities in the management of cryptographic keys within the Advanced Metering Infrastructure (AMI) domain of the electric sector. ORNL concentrated our analysis on the AMI domain of which the National Electric Sector Cyber security Organization Resource (NESCOR) Working Group 1 (WG1) has documented 29 failure scenarios. The computational infrastructure of this metric involves system stakeholders, security requirements, system components and security threats. To compute this metric, we estimated the stakes that each stakeholder associates with each security requirement, as well as stochastic matrices that represent the probability of a threat to cause a component failure and the probability of a component failure to cause a security requirement violation. We applied this model to estimate the security of the AMI, by leveraging the recently established National Institute of Standards and Technology Interagency Report (NISTIR) 7628 guidelines for smart grid security and the International Electrotechnical Commission (IEC) 63351, Part 9 to identify the life cycle for cryptographic key management, resulting in a vector that assigned to each stakeholder an estimate of their average loss in terms of dollars per day of system operation. To further address probabilities of threats, information security analysis can be performed using game theory implemented in dynamic Agent Based Game Theoretic (ABGT) simulations. Such simulations can be verified with the results from game theory analysis and further used to explore larger scale, real world scenarios involving multiple attackers, defenders, and information assets. The strategy for the game was developed by analyzing five electric sector representative failure scenarios contained in the AMI functional domain from NESCOR WG1. From these five selected scenarios, we characterized them into three specific threat categories affecting confidentiality, integrity and availability (CIA). The analysis using our ABGT simulation demonstrated how to model the AMI functional domain using a set of rationalized game theoretic rules decomposed from the failure scenarios in terms of how those scenarios might impact the AMI network with respect to CIA.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A COMPLETE SPECTROSCOPIC SURVEY OF THE MILKY WAY SATELLITE SEGUE 1: DARK MATTER CONTENT, STELLAR MEMBERSHIP, AND BINARY PROPERTIES FROM A BAYESIAN ANALYSIS

We introduce a comprehensive analysis of multi-epoch stellar line-of-sight velocities to determine the intrinsic velocity dispersion of the ultrafaint satellites of the Milky Way. Our method includes a simultaneous Bayesian analysis of both membership probabilities and the contribution of binary orbital motion to the observed velocity dispersion within a 14-parameter likelihood. We apply our method to the Segue 1 dwarf galaxy and conclude that Segue 1 is a dark-matter-dominated galaxy at high probability with an intrinsic velocity dispersion of 3.7{sup +1.4}{sub -1.1} km s{sup -1}. The dark matter halo required to produce this dispersion must have an average density of {rho}-bar{sub 1/2} = 2.5{sub -1.9}{sup +4.1}, M{sub sun}, pc{sup -3} within a sphere that encloses half the galaxy's stellar luminosity. This is the highest measured density of dark matter in the Local Group. Our results show that a significant fraction of the stars in Segue 1 may be binaries with the most probable mean period close to 10 years, but also consistent with the 180 year mean period seen in the solar vicinity at about 1{sigma}. Despite this binary population, the possibility that Segue 1 is a bound star cluster with the observed velocity dispersion arising from the orbital motion of binary stars is disfavored by the multi-epoch stellar velocity data at greater than 99% C.L. Finally, our treatment yields a projected (two-dimensional) half-light radius for the stellar profile of Segue 1 of R{sub 1/2} = 28{sup +5}{sub -4} pc, in excellent agreement with photometric measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

A Semi-Quantitative, Synteny-Based Method to Improve Functional Predictions for Hypothetical and Poorly Annotated Bacterial and Archaeal Genes

During microbial evolution, genome rearrangement increases with increasing sequence divergence. If the relationship between synteny and sequence divergence can be modeled, gene clusters in genomes of distantly related organisms exhibiting anomalous synteny can be identified and used to infer functional conservation. We applied the phylogenetic pairwise comparison method to establish and model a strong correlation between synteny and sequence divergence in all 634 available Archaeal and Bacterial genomes from the NCBI database and four newly assembled genomes of uncultivated Archaea from an acid mine drainage (AMD) community. In parallel, we established and modeled the trend between synteny and functional relatedness in the 118 genomes available in the STRING database. By combining these models, we developed a gene functional annotation method that weights evolutionary distance to estimate the probability of functional associations of syntenous proteins between genome pairs. The method was applied to the hypothetical proteins and poorly annotated genes in newly assembled acid mine drainage Archaeal genomes to add or improve gene annotations. This is the first method to assign possible functions to poorly annotated genes through quantification of the probability of gene functional relationships based on synteny at a significant evolutionary distance, and has the potential for broad application.

59 BASIC BIOLOGICAL SCIENCES↗

Exploring the capabilities of optical pump X-ray probe NEXAFS spectroscopy to track photo-induced dynamics mediated by conical intersections

X-ray spectroscopy is gaining a growing interest in the scientific community, as it represents a versatile and powerful experimental toolbox for probing the dynamics of both core and valence electronic excitations, nuclear motions and material structure, with element and site specificity. Among the various X-ray based techniques, near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, which investigates the energy and probability of resonant core-to-valence transitions, has started to be applied to organic molecules: a recent UV-pump X-ray probe time-resolved NEXAFS experiment [Wolf et al., Nat. Commun., 2017, 8, 1] has shown the capability of the technique to provide information about the ultrafast internal conversion between the bright ππ* and the dark nπ* electronic states of the nucleobase thymine. In the present contribution we introduce an accurate theoretical approach for the simulation of NEXAFS spectra of organic molecules, employing azobenzene as a test case. The electronic structure calculations, which provide both energy levels and transition probabilities of core-to-valence excitations, were here performed with a high level multiconfigurational method, the restricted active space self consistent field (RASSCF/RASPT2). GS- and nπ*-NEXAFS spectra were obtained on the top of key molecular geometries (as the optimized cis, trans and conical intersection(s) structures) as well as along the fundamental isomerization coordinates (namely, symmetric and asymmetric bendings of the phenyl rings, and torsion around the central dihedral). We eventually characterize and explain the origin of the simulated signals, highlighting the specific signatures that make it possible to follow the excited state evolution from the nπ* Franck–Condon point, towards the conical intersection(s).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗

Use of Several Thermal Analysis Techniques on a Hypalon Paint Coating for the Solid Rocket Booster (SRB) of the Space Shuttle

White Hypalon paint is brush-applied as a moisture barrier coating over cork surfaces on each of the two Space Shuttle SRBS. Fine cracks have been observed in the Hypalon coating three times historically on laboratory witness panels, but never on flight hardware. Recent samples of the cracked and standard ("good") Hypalon were removed from cork surfaces and were tested by Thermal Gravimetric Analysis (TGA), Thermomechanical (TMA) and Differential Scanning Calorimetry (DSC) thermal analysis techniques. The TGA data showed that at 700 C, where only paint pigment solids remain, the cracked material had about 9 weight percent more material remaining than the standard material, probably indicating incomplete mixing of the paint before it was brush-applied to produce the cracked material. Use of the TMA film tension method showed that the average static modulus vs. temperature was about 3 to 6 times higher for the cracked material than for the standard material, indicating a much higher stiffness for the cracked Hypalon. The TMA data also showed than an increased coating thickness for the cracked Hypalon was not a factor in the anomaly.

Wingard, Charles D.↗

Use of Several Thermal Analysis Techniques on a Hypalon Paint Coating for the Solid Rocket Booster (SRB) of the Space Shuttle

White Hypalon paint is brush-applied as a moisture barrier coating over cork surfaces on each of the two Space Shuttle SRBs. Fine cracks have been observed in the Hypalon coating three times historically on laboratory witness panels, but never on flight hardware. Samples of the cracked and standard ("good") Hypalon were removed from witness panel cork surfaces, and were tested in 1998 by Thermogravimetric Analysis (TGA), TMA and Differential Scanning Calorimetry (DSC) thermal analysis techniques. The TGA data showed that at 700C, where only paint pigment solids remain, the cracked material had about 9 weight percent more material remaining than the standard material, probably indicating incomplete mixing of the paint before it was brush-applied to produce the cracked material. Use of the TMA film/fiber technique showed that the average modulus (stiffness) vs. temperature was about 3 to 6 times higher for the cracked material than for the standard material. The TMA data also showed that an increase in coating thickness for the cracked Hypalon was not a factor in the anomaly.

Wingard, Charles D.↗

Multiscale analysis of the CMB temperature derivatives

We study the Planck CMB temperature at different scales through its derivatives up to second order, which allows one to characterize the local shape and isotropy of the field. The problem of having an incomplete sky in the calculation and statistical characterization of the derivatives is addressed in the paper. The analysis confirms the existence of a low variance in the CMB at large scales, which is also noticeable in the derivatives. Moreover, deviations from the standard model in the gradient, curvature and the eccentricity tensor are studied in terms of extreme values on the data. As it is expected, the Cold Spot is detected as one of the most prominent peaks in terms of curvature, but additionally, when the information of the temperature and its Laplacian are combined, another feature with similar probability at the scale of 10{sup o} is also observed. However, the p -value of these two deviations increase above the 6% when they are referred to the variance calculated from the theoretical fiducial model, indicating that these deviations can be associated to the low variance anomaly. Finally, an estimator of the directional anisotropy for spinorial quantities is introduced, which is applied to the spinors derived from the field derivatives. An anisotropic direction whose probability is <1% is detected in the eccentricity tensor.

79 ASTRONOMY AND ASTROPHYSICS↗

Probabilistic Risk Assessment of Spent Fuel Pool Cooling System Failure Considering Safety and Security Initiating Events

In order to quantify the probable risk associated with a safety system, nuclear safety risk evaluation process is applied which is generally performed in three steps. The first step is to calculate the damage frequency of a critical system. This step is followed by the radioactive source term estimation due to the potential system damage. The final step is the assessment of off-site consequences. The safety risk formula is shown by Eq. (1). P{sub F} is the probability of safety system failure (the probability used describes a random, independent event), and C is the consequence value. R{sub SAF} = P{sub F} * C (1) The security risk associated with an attack on critical infrastructure and/or special nuclear materials theft is usually evaluated using three- term risk equation (probability of attack, probability of adversary success, and consequence). A security threat and the related initiating events are not considered in a pure safety risk analysis. The vulnerability of the physical protection system is performed independently even though the consequence resulting from a security initiating event can be similar to a safety related initiating event and thus gets double counted. The security risk can be calculated using Eq. (2) R{sub SEC} = P{sub A} * P{sub S} * C (2) PA is the probability of attack, P{sub S} is the probability of adversary success (P{sub S} = 1 - P{sub I} P{sub N}). P{sub I} being the probability of interruption which accounts for the probabilities of detection of all timely detection elements along an adversary path. PN being the probability of neutralization: the result of engagement of the response force with the adversaries after the interruption has occurred. C is the consequence magnitude from a successful attack. The two security probabilities do not describe random, independent events. As a result, because of the differences in the randomness and independence of the respective safety and security probabilities, the safety and the security risks are currently calculated separately. The consequence term in nuclear safety and security risk evaluation have the same definition. An effort is made here to study the effect on the overall risk by considering security initiating events along with the subsequent failure of the intermediate safety element. The target selected is a spent fuel storage pool for analyzing the public risk arising from the interplay between nuclear security and safety initiating events. EASI (Estimate of Adversary Sequence Interruption) model is used to calculate the probability of interruption related to a security system that leads to the calculation of the security initiating event frequency. Further progression of this initiating security event to a top event is analyzed using a Probabilistic Safety Assessment (PSA) tool, SAPHIRE (version 8) to obtain a joint pool damage frequency. For the progression of events leading to a top event, the failure of the intermediate safety elements is assumed to be external random events such as heavy load drop and fire events. The objectives of this study is to (1) determine the double accounting of frequency of failure, risk and allocated resources in the scenario of a non-random security initiating event leading to a safety random event chain (2) compare the damage frequency from pure security with the frequency of combined security-safety failure event and pure safety failure events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SU-E-T-531: Large--Scale DVH Quality Study: Correlated Aims Lead Relaxations

Purpose: Intensity modulated radiation therapy plans are designed to optimally target a tumor while sparing surrounding tissue. Desired dose distributions are iteratively approached via inverse planning. This leads to tradeoffs between clinical objectives for the planning target volume (PTV), organs at risk, and normal tissues. Dose volume histogram (DVHs) related aims are followed that are either institutional or internationally recommended. We analyze common goals and identify potential reasons that often lead to tradeoffs. Methods: 524 IMRT plans for various tumor sites were analyzed based on the main institutional DVH goal for PTV (D95) and the recommendations by ICRU-83 (D2, D50, and D98). Robust statistical tools are developed and applied to ensure that the results are immune to data uncertainties. The probability of violation was measured for each of the DVH goals based on the frequency of not meeting recommended doses. Conditional probabilities for satisfying and/or violating DVH aims were computed to test the hypothesized pair-wise relations between DVH aims. For example, for plans that satisfied D50, the probability of violating D98 was computed via P(D98 < 95% | 98% ≤ D50 ≤ 102%). The equality constraint D50 = 100% was relaxed to encompass the range [98,102]%. Results: A large majority of cases (88%) satisfied the institutional goal for PTV of D95 ≤ 95%. Similar consensus existed for D98. 51% of cases satisfied D2 ≥ 107%. However, only 18% of cases satisfied D50. The conditional probability showed correlations amongst the studied DVH goals. In fact, a negative correlation was revealed between D50 and D95 (and D98), suggesting that these competing goals cannot be satisfied concurrently. Conclusion: The majority of plans followed the institutional guidelines. The reason for their deviation from international recommendations seems to be that the latter goals are competing and cannot be satisfied concurrently in clinical practice.

60 APPLIED LIFE SCIENCES↗

Determining the Probability of Violating Upper-Level Wind Constraints for the Launch of Minuteman III Ballistic Missiles at Vandenberg Air Force Base

The 30th Operational Support Squadron Weather Flight (30 OSSWF) provides comprehensive weather services to the space program at Vandenberg Air Force Base (VAFB) in California. One of their responsibilities is to monitor upper-level winds to ensure safe launch operations of the Minuteman Ill ballistic missile. The 30 OSSWF tasked the Applied Meteorology Unit (AMU) to analyze VAFB sounding data with the goal of determining the probability of violating (PoV) their upper-level thresholds for wind speed and shear constraints specific to this launch vehicle, and to develop a tool that will calculate the PoV of each constraint on the day of launch. In order to calculate the probability of exceeding each constraint, the AMU collected and analyzed historical data from VAFB. The historical sounding data were retrieved from the National Oceanic and Atmospheric Administration Earth System Research Laboratory archive for the years 1994-2011 and then stratified into four sub-seasons: January-March, April-June, July-September, and October-December. The AMU determined the theoretical distributions that best fit the maximum wind speed and maximum wind shear datasets and applied this information when calculating the averages and standard deviations needed for the historical and real-time PoV calculations. In addition, the AMU included forecast sounding data from the Rapid Refresh model. This information provides further insight for the launch weather officers (LWOs) when determining if a wind constraint violation will occur over the next few hours on the day of launch. The AMU developed an interactive graphical user interface (GUI) in Microsoft Excel using Visual Basic for Applications. The GUI displays the critical sounding data easily and quickly for LWOs on day of launch. This tool will replace the existing one used by the 30 OSSWF, assist the LWOs in determining the probability of exceeding specific wind threshold values, and help to improve the overall upper winds forecast for the launch customer. This presentation will describe how the AMU calculated the historical and real-time PoV values for the specific upper-level wind launch constraints and outline the development of the interactive GUI display.

Shafer, Jaclyn A.↗

Superconducting surface impedance under radiofrequency field

Based on BCS theory with moving Cooper pairs, the electron states distribution at 0K and the probability of electron occupation with finite temperature have been derived and applied to anomalous skin effect theory to obtain the surface impedance of a superconductor under radiofrequency (RF) field. We present the numerical results for Nb and compare these with representative RF field-dependent effective surface resistance measurements from a 1.5 GHz resonant structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗

IMPROVED Ni I log(gf) VALUES AND ABUNDANCE DETERMINATIONS IN THE PHOTOSPHERES OF THE SUN AND METAL-POOR STAR HD 84937

Atomic transition probability measurements for 371 Ni I lines in the UV through near-IR are reported. Branching fractions from data recorded using a Fourier transform spectrometer and a new echelle spectrograph are combined with published radiative lifetimes to determine these transition probabilities. Generally good agreement is found in comparisons to previously reported Ni I transition probability measurements. Use of the new echelle spectrograph, independent radiometric calibration methods, and independent data analysis routines enable a reduction of systematic errors and overall improvement in transition probability uncertainty over previous measurements. The new Ni I data are applied to high-resolution visible and UV spectra of the Sun and metal-poor star HD 84937 to derive new, more accurate Ni abundances. Lines covering a wide range of wavelength and excitation potential are used to search for non-LTE effects.

79 ASTRONOMY AND ASTROPHYSICS↗

The detective quantum efficiency of photon-counting x-ray detectors using cascaded-systems analyses

Purpose: Single-photon counting (SPC) x-ray imaging has the potential to improve image quality and enable new advanced energy-dependent methods. The purpose of this study is to extend cascaded-systems analyses (CSA) to the description of image quality and the detective quantum efficiency (DQE) of SPC systems. Methods: Point-process theory is used to develop a method of propagating the mean signal and Wiener noise-power spectrum through a thresholding stage (required to identify x-ray interaction events). The new transfer relationships are used to describe the zero-frequency DQE of a hypothetical SPC detector including the effects of stochastic conversion of incident photons to secondary quanta, secondary quantum sinks, additive noise, and threshold level. Theoretical results are compared with Monte Carlo calculations assuming the same detector model. Results: Under certain conditions, the CSA approach can be applied to SPC systems with the additional requirement of propagating the probability density function describing the total number of image-forming quanta through each stage of a cascaded model. Theoretical results including DQE show excellent agreement with Monte Carlo calculations under all conditions considered. Conclusions: Application of the CSA method shows that false counts due to additive electronic noise results in both a nonlinear image signal and increased image noise. There is a window of allowable threshold values to achieve a high DQE that depends on conversion gain, secondary quantum sinks, and additive noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-F-T-129: Impact of Radial Fluctuations in RBE for Therapeutic Proton Beams

Purpose: To evaluate the off axis relative biological effectiveness (RBE) for actively scanned proton beams and determine if a constant radial RBE can be assumed. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary protons were simulated for a 0.6cm diameter pencil beam. Beam energies corresponding to Bragg Peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially for annuli of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm outer radius. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model, for human submandibular gland (HSG) cells, to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. Results: RBE was generally seen to increase as distance from the central axis (CAX) increased. However, this increase was only seen in low dose regions and its overall effects on the transverse biological dose remains low. In the entrance region of the phantom (10mm depth), minimum and maximum calculated RBEs varied between 15.22 and 18.88% for different energies. At the Bragg peak, this difference ranged from 3.15 to 26.77%. Despite these rather large variations the dose-weighted RBE and the CAX RBE varied by less than 0.14% at 10mm depth and less than 0.16% at the Bragg peak. Similarly small variations were found at all depths proximal of the Bragg peak. Conclusion: Although proton RBE does vary radially, its overall effect on biological dose is minimal and the use of a radially constant RBE in treatment planning for scanned proton beams would not produce large errors.

60 APPLIED LIFE SCIENCES↗

WE-H-BRA-09: Application of a Modified Microdosimetric-Kinetic Model to Analyze Relative Biological Effectiveness of Ions Relevant to Light Ion Therapy Using the Particle Heavy Ion Transport System

Purpose: To evaluate the dose and biological effectiveness of various ions that could potentially be used for actively scanned particle therapy. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary particles for 1H beams and ten million particles for 4He, 7Li, 10B, 12C, 14N, 16O, and 20Ne were simulated for 0.6cm diameter pencil beams. Beam energies corresponding to Bragg peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially in annuli with outer radius of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model for five different cell types to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. The product of the calculated RBEs and the simulated physical dose was taken to create biological dose and comparisons were then made between the various ions. Results: Transversely, the 10B beam was seen to minimize relative biological dose in both the constant and accelerated dose change regions, proximal to the Bragg Peak, for all beams traveling greater than 50mm. For the 50mm beam, 7Li was seen to provide the most optimal biological dose profile. Radially small fluctuations (<4.2%) were seen in RBE while physical dose was greater than 1% for all beams. Conclusion: Even with the growing usage of 12C, it may not be the most optimal ion in all clinical situations. Boron was calculated to have slightly enhanced RBE characteristics, leading to lower relative biological doses.

60 APPLIED LIFE SCIENCES↗

Anharmonic Bloch oscillations of electrons in electrically biased superlattices

The phenomenon of anharmonic Bloch oscillations (i.e., oscillations with a frequency which is a multiple of the Bloch frequency) is considered. The energy-band structure of silicon-carbide polytypes where these oscillations are observed is calculated ab initio. A one-dimensional model potential making it possible to calculate the Stark wave functions in a biased superlattice is constructed for these polytypes. The transfer matrix method is used to find the set of energy levels (the so-called Stark ladder) and calculate the electron wave functions. It is shown that both radiative transitions between neighboring Stark levels (at the Bloch frequency) and transitions in the form of jumps via several levels, accompanied by emission at frequencies that are multiples of the Bloch frequency, may occur. The calculated probabilities of transitions between Stark-ladder levels increase with increasing applied field strength.

36 MATERIALS SCIENCE↗