Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed data structures”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Detecting hidden layers from spreading dynamics on complex networks

When dealing with spreading processes on networks it can be of the utmost importance to test the reliability of data and identify potential unobserved spreading paths. In this paper we address these problems and propose methods for hidden layer identification and reconstruction. We also explore the interplay between difficulty of the task and the structure of the multilayer network describing the whole system where the spreading process occurs. Our methods stem from an exact expression for the likelihood of a cascade in the susceptible-infected model on an arbitrary graph. We then show that by imploring statistical properties of unimodal distributions and simple heuristics describing joint likelihood of a series of cascades one can obtain an estimate of both existence of a hidden layer and its content with success rates far exceeding those of a null model. Furthermore, we conduct our analyses on both synthetic and real-world networks providing evidence for the viability of the approach presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Correction of confidence intervals in excess relative risk models using Monte Carlo dosimetry systems with shared errors

Introduction: Measurement errors are ubiquitous in epidemiological studies, especially for doses arising from environmental and occupational exposures, which are difficult to assess. Correction methods for independent measurement errors, have been comprehensively described elsewhere, including regression calibration, simulation extrapolation (SIMEX), etc. However, when exposure estimates are constructed using complicated physical and biological models, uncertainties in the dosimetry system can become very complex. In such systems, uncertainties in some parameters may affect a large group of study participants simultaneously, and hence are “shared” as opposed to independent dosimetric uncertainties. Representing the structure of the uncertainty in the dosimetry by providing multiple realizations of possible dose has been described for several studies. The goal of this paper is to implement, and explore the properties of, a generalized estimating equation-based approach to shared uncertainties based on principles described in an earlier paper. Methods: The effects of shared and unshared dosimetry error on parameter estimates and confidence intervals in a generalized excess relative risk model are examined. A post-hoc correction method is implemented that constructs corrected confidence intervals based on an underlying variance structure described in an earlier paper. In this correction method it is assumed that the estimated slope parameter in a linear excess relative risk model is distributed as a mixture of normal and lognormal components. A simulation study of the behavior of this correction was developed based on multiple realizations provided by the Monte Carlo dosimetry system used for the plutonium exposures in the Mayak Worker Cohort and a risk model with close similarity to that used in analyses of lung cancer data in that cohort. Results: In the simulation study we found that shared dose errors severely degraded the coverage probabilities of confidence limits for the slope parameter in the linear excess relative risk model relating uncertain exposure to lung cancer risk. Using our correction method to adjust the confidence intervals provided much improved coverage probabilities for this parameter. A smaller amount of degradation in confidence interval coverage for other parameters (e.g. those in the baseline and dose modification terms) was noted, and the proposed correction method also improved coverage for these parameters as well. Discussion: Compared to other recently described methods our proposed correction method relies on relatively familiar software to fit dose response models. The basic modeling philosophy described here first uses existing Poisson regression methods to estimate parameters in the excess relative risk model; in this first stage mean dose (computed from the multiple realizations) is used as the independent variable in the regression. Once the model is fit we then correct (widen) the standard errors for all parameters using the proposed correction method. The simulation studies indicate that the methods described here are an important contribution to the available tools for dealing with shared measurement errors in complex risk analysis.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Ducted Assembly Steady-State Heat Transfer Software (DASSH)

The Ducted Assembly Steady-State Heat transfer software (DASSH) is being developed at Argonne National Laboratory to perform a steady-state thermal fluids calculation to determine the coolant flow and temperature distribution for a hexagonal reactor core configuration with ducted assemblies. DASSH is intended to be used in the early stages of the reactor design process when assembly components can be undefined or are undergoing considerable design changes. DASSH provides a rapid assessment of the temperature and flow distribution to allow quick characterization of the system and identification of problem areas. This document is a guide for DASSH users. It summarizes the core capabilities of DASSH and highlights differences between DASSH and SE2-ANL. Instructions for obtaining and installing DASSH are provided along with some guidelines and recommendations for newer Python users. It covers the structure of the input file, provides directions for running DASSH, describes the various output files produced by the code, and provides instructions for visualizing data. Examples demonstrating various DASSH options are included. DASSH is under active development, so it is anticipated that this guide will grow and evolve as the code is updated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Realized ecological forecast through an interactive Ecological Platform for Assimilating Data (EcoPAD, v1.0) into models

Abstract. Predicting future changes in ecosystem services is not only highly desirable but is also becoming feasible as several forces (e.g., available big data, developed data assimilation (DA) techniques, and advanced cyber-infrastructure) are converging to transform ecological research into quantitative forecasting. To realize ecological forecasting, we have developed an Ecological Platform for Assimilating Data (EcoPAD, v1.0) into models. EcoPAD (v1.0) is a web-based software system that automates data transfer and processing from sensor networks to ecological forecasting through data management, model simulation, data assimilation, forecasting, and visualization. It facilitates interactive data–model integration from which the model is recursively improved through updated data while data are systematically refined under the guidance of model. EcoPAD (v1.0) relies on data from observations, process-oriented models, DA techniques, and the web-based workflow. We applied EcoPAD (v1.0) to the Spruce and Peatland Responses Under Climatic and Environmental change (SPRUCE) experiment in northern Minnesota. The EcoPAD-SPRUCE realizes fully automated data transfer, feeds meteorological data to drive model simulations, assimilates both manually measured and automated sensor data into the Terrestrial ECOsystem (TECO) model, and recursively forecasts the responses of various biophysical and biogeochemical processes to five temperature and two CO 2 treatments in near-real time (weekly). Forecasting with EcoPAD-SPRUCE has revealed that mismatches in forecasting carbon pool dynamics are more related to model (e.g., model structure, parameter, and initial value) than forcing variables, opposite to forecasting flux variables. EcoPAD-SPRUCE quantified acclimations of methane production in response to warming treatments through shifted posterior distributions of the CH 4 :CO 2 ratio and the temperature sensitivity ( Q 10 ) of methane production towards lower values. Different case studies indicated that realistic forecasting of carbon dynamics relies on appropriate model structure, correct parameterization, and accurate external forcing. Moreover, EcoPAD-SPRUCE stimulated active feedbacks between experimenters and modelers to identify model components to be improved and additional measurements to be taken. It has become an interactive model–experiment (ModEx) system and opens a novel avenue for interactive dialogue between modelers and experimenters. Altogether, EcoPAD (v1.0) acts to integrate multiple sources of information and knowledge to best inform ecological forecasting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Neutral Pion Electroproduction and development of a Neutral Particle Spectrometer

Protons and neutrons, i.e. nucleons, are the basic building blocks of the matter in the visible universe. The strong force binds the nucleons to form nuclei. The electromagnetic force forms the atoms by binding the electrons with the nuclei. The electromagnetic interaction is well understood by Quantum Electrodynamics (QED), which shows the most precise predictability amongst all the theories in physics. In QED, charges interact with each other by exchanging photons. The Quantum Chromodynamics (QCD) describes the strong interaction. Its degrees of freedom are quarks and gluons, the fundamental constituents of the nucleons. The quarks interact with each other by exchanging gluons. However, unlike QED, the gluons interact amongst themselves. This feature of self bindings of the gluons confines the quarks and gluons in the nucleons/hadrons, never to be seen as free. In order to study some of the features of QCD, such as confinement or the structure of the nucleon, one usually needs to rely on experiments. Electromagnetic probes, governed by the well-understood QED, are excellent tools to probe the nucleon. In general, different scales, e.g. electron beam energies, probe different regions of the nucleon. At low energy, of the order of a few GeV, the electron probes the nucleon in the valence quark region. As its energy increases, the electron probes the sea quark and gluon regions. The study of the nucleon structure in all these regions is needed to fully understand QCD. Form factors and parton distribution functions measured from elastic scattering and deep inelastic scattering of leptons off nucleons have provided a partial view of the internal structure of the nucleon. In the mid-1990s, Generalized Parton Distributions (GPDs) were developed. These new objects are a generalization of the form factors and parton distribution functions, but contain richer information on the nucleon internal structure. GPDs are accessible experimentally by deep exclusive reactions. Deeply virtual Compton scattering (DVCS) and deeply virtual meson production (DVMP) are some examples. The first dedicated DVCS/DVMP experiment took place in 2004 in Hall A at The Thomas Jefferson National Accelerator Facility, i.e. Jefferson Lab, in Virginia, U.S.A. A new DVCS/DVMP experiment, after the beam energy upgrade of Jefferson Lab, was carried out in Hall A in a wider kinematic range. Its data were taken from 2014 to 2016. In Hall C at Jefferson Lab, the next DVCS/DVMP experiment will take place. The Hall C experiment will further exploit the kinematic range with higher precision. A Neutral Particle Spectrometer (NPS) is in development to measure DVCS/DVMP events under high background conditions. Jefferson Lab will provide the highest precision data in the valence quark region for various exclusive reactions. The Electron-Ion Collider (EIC) is a future experimental facility currently planned to start operations around 2030 in the U.S.A. Its high energy and high luminosity will probe the sea quark and gluon regions providing answers to the outstanding questions of QCD, in particular in the region where matter is dominated by gluons. First of all, this document describes the data analysis and results of the Hall A neutral pion electroproduction off the proton, from the data taken in 2014-2016. Later, some of the developments towards the construction of the electromagnetic calorimeter of the NPS for the upcoming DVCS/DVMP experiment in Hall C are presented. Finally, one of the candidate materials for the EIC calorimeter, a glass scintillator, will be briefly introduced. I have participated to all these projects, in collaboration with many colleagues. I present in this thesis my contributions to each of these projects. My contributions to the neutral pion data analysis were focused on background subtractions on the calorimeter, acceptance calculations, and the estimation of the systematic uncertainty associated to the event selection cuts. Some necessary information on calibrations of the detectors and data analysis methods are also described. In the NPS project, I performed background dose calculations and energy and position resolution studies of the calorimeter, all using Monte Carlo simulations, with realistic geometries of the experimental apparatus. Characterization of the crystals of the calorimeter was also done. Additionally, I measured the radiation hardness of some glass scintillator in its early stage of development. In order to have a future reference when the glass calorimeter prototype will be tested, I simulated the energy resolution of the prototype.

Ko, Ho-San↗
DOE OSTI2025

Modeling Cyber Supply Chain Incidents with Multilayered Graph Motifs

As noted within the literature, supply chain includes people and organizations---manufacturers, integrators, and third-party vendors---that are involved in one or more stages of a product lifecycle. Since supply chains, by definition, include organizations and people, supply chain risk management activities must consider dependencies between an organization's business processes and third-party resources. Just as adversarial tactics can be implemented via techniques implemented via networked computer systems, so can such tactics be expressed via legal business relationships. A cyber incident may have an exponential impact downstream, for example, by leveraging a product's distribution channel (e.g. malicious updates in SolarWinds, buggy updates in CrowdStrike). Similarly, legitimate and legal business relationships also affect the attack surface exposure of systems, enabling long-term persistence and/or unknown impacts to product quality that are hard to detect. This paper catalogs several recent digital supply chain incidents and applies a multilayered network formalism to develop structural indicators (graph motifs) that reflect potentially-adversarial behavior. Finally, we compare and contrast the characteristics of adversarial tactics (e.g. Loss of Availability, Data Collection) that leverage cyber-physical dependencies to those that leverage legal organizational relationships.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Statistical analysis for the neutrinoless double-β-decay matrix element of 48 Ca

Neutrinoless double-β-decay (0⁢vββ) nuclear matrix elements (NME) are the object of many theoretical calculation methods, and are very important for analysis and guidance of a large number of experimental efforts. However, there are large discrepancies between the NME values provided by different methods. Here, in this paper, we propose a statistical analysis of the 48 Ca 0v⁢ββ NME using the interacting shell model, emphasizing the range of the NME probable values and their correlations with observables that can be obtained from the existing nuclear data. Based on this statistical analysis with three independent effective Hamiltonians, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (0.45–0.95) at 90% confidence level, and a mean value of 0.68.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Morphogenesis of respiratory syncytial virus in human primary nasal ciliated epithelial cells occurs at surface membrane microdomains that are distinct from cilia

The distribution of cilia and the respiratory syncytial virus (RSV) nucleocapsid (N) protein, fusion (F) protein, attachment (G) protein, and M2-1 protein in human ciliated nasal epithelial cells was examined at between 1 and 5 days post-infection (dpi). All virus structural proteins were localized at cell surface projections that were distinct from cilia. The F protein was also trafficked into the cilia, and while its presence increased as the infection proceeded, the N protein was not detected in the cilia at any time of infection. The presence of the F protein in the cilia correlated with cellular changes in the cilia and reduced cilia function. At 5 dpi extensive cilia loss and further reduced cilia function was noted. These data suggested that although RSV morphogenesis occurs at non-cilia locations on ciliated nasal epithelial cells, RSV infection induces changes in the cilia body that leads to extensive cilia loss. - Highlights: • Respiratory syncytial virus (RSV) infects nasal ciliated epithelial cells. • Virus morphogenesis occurs within filamentous projections distinct from cilia. • The RSV N protein was not detected in the cilia at any time during infection. • Trafficking of the F protein into the cilia occurred early in infection. • Presence of the F protein in cilia correlated with impaired cilia function.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Decoding the proton’s gluonic density with lattice QCD-informed machine learning

We present a first machine learning-based decoding of the gluonic structure of the proton from lattice QCD using a variational autoencoder inverse mapper (VAIM). Harnessing the power of generative AI, we predict the parton distribution function (PDF) of the gluon given information on the reduced pseudo-Ioffe-time distributions (RpITDs) as calculated from an ensemble with lattice spacing a ≈ 0.09 fm and a pion mass of M π ≈ 310 MeV. The resulting gluon PDF is consistent with phenomenological global fits within uncertainties, particularly in the intermediate-to-high-x region where lattice data are most constraining. A subsequent correlation analysis confirms that the VAIM learns a meaningful latent representation, highlighting the potential of generative AI to bridge lattice QCD and phenomenological extractions within a unified analysis framework.

Gluon parton distribution function↗
DOE OSTI2018

Structure, Chemistry, and Charge Transfer Resistance of the Interface between Li 7 La 3 Zr 2 O 12 Electrolyte and LiCoO 2 Cathode

All-solid-state batteries promise significant safety and energy density advantages over liquid-electrolyte batteries. The interface between the cathode and the solid electrolyte is an important contributor to charge transfer resistance. Strong bonding of solid oxide electrolytes and cathodes requires sintering at elevated temperatures. Knowledge of the temperature dependence of the composition and charge transfer properties of this interface is important for determining the ideal sintering conditions. To understand the interfacial decomposition processes and their onset temperatures, model systems of LiCoO 2 (LCO) thin films deposited on cubic Al-doped Li 7 La 3 Zr 2 O 12 (LLZO) pellets were studied as a function of temperature using interface-sensitive techniques. X-ray photoelectron spectroscopy, secondary ion mass spectroscopy, and energy-dispersive X-ray spectroscopy data indicated significant cation interdiffusion and structural changes starting at temperatures as low as 300 °C. La 2 Zr 2 O 7 and Li 2 CO 3 were identified as decomposition products after annealing at 500 °C by synchrotron X-ray diffraction. X-ray absorption spectroscopy results indicate the presence of also LaCoO 3 in addition to La 2 Zr 2 O 7 and Li 2 CO 3 . On the basis of electrochemical impedance spectroscopy and depth profiling of the Li distribution upon potentiostatic hold experiments on symmetric LCO|LLZO|LCO cells, the interfaces exhibited significantly increased impedance, up to 8 times that of the as-deposited samples after annealing at 500 °C. Here, our results indicate that lower-temperature processing conditions, shorter annealing time scales, and CO 2 -free environments are desirable for obtaining ceramic cathode|electrolyte interfaces that enable fast Li transfer and high capacity.

25 ENERGY STORAGE↗
DOE OSTI2024

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G 0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G 0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Distributed Acoustic Sensing Using Dark Fiber for Near-Surface Characterization and Broadband Seismic Event Detection

We present one of the first case studies demonstrating the use of distributed acoustic sensing deployed on regional unlit fiber-optic telecommunication infrastructure (dark fiber) for broadband seismic monitoring of both near-surface soil properties and earthquake seismology. We recorded 7 months of passive seismic data on a 27 km section of dark fiber stretching from West Sacramento, CA to Woodland, CA, densely sampled at 2 m spacing. This dataset was processed to extract surface wave velocity information using ambient noise interferometry techniques; the resulting VS profiles were used to map both shallow structural profiles and groundwater depth, thus demonstrating that basin-scale variations in hydrological state could be resolved using this technique. The same array was utilized for detection of regional and teleseismic earthquakes and evaluated for long period response using records from the M8.1 Chiapas, Mexico 2017, Sep 8th event. The combination of these two sets of observations conclusively demonstrates that regionally extensive fiber-optic networks can effectively be utilized for a host of geoscience observation tasks at a combination of scale and resolution previously inaccessible.

58 GEOSCIENCES↗
DOE OSTI2021

Galaxy zoo builder: Morphological dependence of spiral galaxy pitch angle

ABSTRACT Spiral structure is ubiquitous in the Universe, and the pitch angle of arms in spiral galaxies provide an important observable in efforts to discriminate between different mechanisms of spiral arm formation and evolution. In this paper, we present a hierarchical Bayesian approach to galaxy pitch angle determination, using spiral arm data obtained through the Galaxy Builder citizen science project. We present a new approach to deal with the large variations in pitch angle between different arms in a single galaxy, which obtains full posterior distributions on parameters. We make use of our pitch angles to examine previously reported links between bulge and bar strength and pitch angle, finding no correlation in our data (with a caveat that we use observational proxies for both bulge size and bar strength which differ from other work). We test a recent model for spiral arm winding, which predicts uniformity of the cotangent of pitch angle between some unknown upper and lower limits, finding our observations are consistent with this model of transient and recurrent spiral pitch angle as long as the pitch angle at which most winding spirals dissipate or disappear is larger than 10°.

Lingard, Timothy↗
DOE OSTI2023

Robust clustering of the local Milky Way stellar kinematic substructures with Gaia eDR3

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm HDBSCAN on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the HDBSCAN member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Comparative Genomics and Functional Analysis of Carbohydrate Utilization Networks in Unicyanobacterial Consortia Derived from Hypersaline Lake Microbial Mats

Hot Lake is a meromictic, hypersaline lake located in Washington state, USA. It harbors a phototrophic microbial mat community stable in dramatic seasonal environmental changes. The complexity of natural community hinders to analyze metabolic interactions between its members. To manage this obstacle, unicyanobacterial consortia (UCC) were sampled from the community and sequenced. UCC is a microbial community in which one cyanobacterium serves as the sole primary producer of carbon sources to several heterotrophic bacteria. Recently available metagenomes for two UCC revealed a shared set of 19 microbial operational taxonomic units (OTU), whose individual genomes were reconstructed by co‐assembly and abundance profile binning. We used this genomic data and a comparative genomics approach to infer carbohydrate utilization pathways and their transcription regulation. We scanned all UCC proteomes against the KEGG Orthology, Pfam, GO and SEED databases to reveal genes potentially involved in sugar utilization. Then we analyzed the genomic and functional contexts of gene loci encoding the obtained proteins and reconstructed the respective metabolic pathways and transcriptional regulons. As result, we observed catabolic pathways for 32 various carbohydrates and their derivatives. Distribution of the inferred catabolic pathways has highly mosaic structure across the UCC genomes. The largest number of pathways was observed in two representatives of the Halomonas genus (14–15 pathways) and one Rhodobacteriaceae sp. (15 pathways). In other UCC genomes, we found up to 10 carbohydrate utilization pathways per genome. However in 7 genomes we didn't observe any pathway. The most widely distributed pathways in UCC were the alpha‐glucosides and the DeLey‐Doudoroff galactose utilization pathways. Both were found in 7 genomes. Majority of pathways are distributed among 1–3 genomes. The reconstructed catabolic pathways were partially validated by phenotypic analyses of UCC community members that were isolated in a pure culture. In the Rhodobacteriaceae sp., we found a novel catabolic pathway that was proposed to be involved in mannoheptulose utilization. To our knowledge, the proposed mannoheptulose catabolic pathway is the first case of a pathway for heptose utilization and will require further experimental validation. The identified pathways are potentially regulated by local transcription factors (TFs) that are often co‐localized with the respective sugar catabolic gene loci. For 63 out of the 80 identified TFs we identified their cognate DNA binding motifs and reconstructed regulons using the bioinformatics approach. Support or Funding Information This research was supported by the Genomic Science Program (GSP), Office of Biological and Environmental Research (OBER), and U.S. Department of Energy (DOE) and is a contribution of the Pacific Northwest National Laboratory (PNNL) Foundational Scientific Focus Area.

Leyn, Semen A.↗
DOE OSTI2015

SU-F-BRD-14: Dose Weighted Linear Energy Transfer Analysis of Critical Structures in Proton Therapy of Pediatric Brain Tumor Patients

Purpose: Dose weighted Linear Energy Transfer (LETd) analysis of critical structures may be useful in understanding the side effects of the proton therapy. The objective is to analyze the differences between LETd and dose distributions in brain tumor patients receiving double scattering proton therapy, to quantify LETd variation in critical organs, and to identify beam arrangements contributing to high LETd in critical organs. Methods: Monte Carlo simulations of 9 pediatric brain tumor patients were performed. The treatment plans were reconstructed with the TOPAS Monte Carlo code to calculate LETd and dose. The beam data were reconstructed proximal to the aperture of the double scattering nozzle. The dose and LETd to target and critical organs including brain stem, optic chiasm, lens, optic nerve, pituitary gland, and hypothalamus were computed for each beam. Results: Greater variability in LETd compared to dose was observed in the brainstem for patients with a variety of tumor types including 5 patients with tumors located in the posterior fossa. Approximately 20%–44% brainstem volume received LETd of 5kev/µm or greater from beams within gantry angles 180°±30° for 5 patients treated with a 3 beam arrangement. Critical organs received higher LETd when located in the vicinity of the beam distal edge. Conclusion: This study presents a novel strategy in the evaluation of the proton treatment impact on critical organs. While the dose to critical organs is confined below the required limits, the LETd may have significant variation. Critical organs in the vicinity of beam distal edge receive higher LETd and depended on beam arrangement, e.g. in posterior fossa tumor treatment, brainstem receive higher LETd from posterior-anterior beams. This study shows importance of the LETd analysis of the radiation impact on the critical organs in proton therapy and may be used to explain clinical imaging observations after therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Machine learning applied to classifying neutron resonances

The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonance structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. In this project, we apply a machine learning technique, namely decision trees, to automate the quantum number assignments. The tree is trained from simulated data generated to mimic the errors found in real data. We explore the use of several physics-motivated features for training our trees, including the nearest neighbor spacing distribution, cumulative level distribution, and channel width distributions. Initial results using random matrix theory motivated fits which demonstrated that we can determine resonance spin groups somewhat reliably. If we use these fits as features in our trees, we can train them to spot outliers corresponding to misassigned resonances. We found that with the large number of features used in this project that the decision tree tended to over t training data resulting in poor performance with respect to the test data. By reducing the number of features, we can achieve nearly perfect assignment of quantum numbers with our training data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗