Search NASA⌕ Search

SEARCH · Search NASA

Results for “data intensive computing”

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,171 records · Page 65

SU-F-P-09: A Global Medical Physics Collaboration for Implementation of Modern Radiotherapy in Botswana

Purpose: The global burden of cancer is considerable, particularly in low and middle-income countries. Massachusetts General Hospital (MGH) and Botswana-Harvard AIDS Institute have partnered with the oncology community and government of Botswana to form BOTSOGO (BOTSwana Oncology Global Outreach) to address the rising burden of cancer in Botswana. Currently, radiation therapy (RT) is only available at a single linear accelerator (LINAC) in Gaborone Private Hospital (GPH). BOTSOGO worked to limit the absence of RT during a LINAC upgrade and ensure a safe transition to modern radiotherapy techniques. Methods: The existing Elekta Precise LINAC was decommissioned in November 2015 and replaced with a new Elekta VERSA-HD with IMRT/VMAT/CBCT capability. Upgraded treatment planning and record-and-verify systems were also installed. Physicists from GPH and MGH collaborated during an intensive on-site visit in Botswana during the commissioning process. Measurements were performed using newly purchased Sun Nuclear equipment. Photon beams were matched with an existing model to minimize the time needed for beam modeling and machine down time. Additional remote peer review was also employed. Independent dosimetry was performed by irradiating OSLDs, which were subsequently analyzed at MGH. Results: Photon beam quality agreed with reference data within 0.2%. Electron beam data agreed with example clinical data within 3%. Absolute dose calibration was performed using both IAEA and AAPM protocols. Absolute dose measurements with OSLDs agreed within 5%. Quentry cloud-based software was installed to facilitate remote review of treatment plans. Patient treatments resumed in February 2016. The time without RT was reduced, therefore likely resulting in reduced patient morbidity/mortality. Conclusion: A global physics collaboration was utilized to commission a modern LINAC in a resource-constrained setting. This can be a useful model in other areas with limited resources. Further use of technology and on-site exchanges will facilitate the introduction of more advanced techniques in Botswana. We acknowledge funding support from the AAPM International Educational Activities Committee and the NCI Federal Share Proton Beam Program Income Grant.

60 APPLIED LIFE SCIENCES↗

Data center growth in the United States: Decoupling the demand for services from electricity use

Data centers are energy intensive buildings that have grown in size and number to meet the increasing demands of a digital economy. This paper presents a bottom-up model to estimate data center electricity demand in the United States over a 20 year period and examines observed and projected electricity use trends in the context of changing data center operations. Results indicate a rapidly increasing electricity demand at the turn of the century that has significantly subsided to a nearly steady annual electricity use of about 70 billion kWh in recent years. While data center workloads continue to grow exponentially, comparable increases in electricity demand have been avoided through the adoption of key energy efficiency measures and a shift towards large cloud-based service providers. Alternative projections from the model illustrate the wide range in potential electricity that could be consumed to support data centers, with the US data center workload demand estimated for 2020 requiring a total electricity use that varies by about 135 billion kWh, depending on the adoption rate of efficiency measures during this decade. While recent improvements in data center energy efficiency have been a success, the growth of data center electricity use beyond 2020 is uncertain, as modeled trends indicate that the efficiency measures of the past may not be enough for the data center workloads of the future. The results show that successful stabilization of data center electricity will require new innovations in data center efficiency to further decouple electricity demand from the ever-growing demand for data center services.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

SU-E-T-372: Evaluation and Comparison of Second-Check Monitor Unit Calculation Software with Pinnacle Treatment Planning System

Purpose: The purpose of this study was to evaluate and compare the accuracy of dose calculation algorithms in the second check software programs Radcalc, Diamond, IMSure, and MUcheck, against the Pinnacle3 treatment planning system (TPS). Methods: Baseline accuracy of the second check software was established by comparison against Pinnacle TPS data using open square fields of 5, 10, 20, 30 and 40cm in a SAD setup. 18 previously treated patients’ files were exported from the Pinnacle3 TPS to each of the four second check softwares, consisting of 146 step and shoot intensity modulated radiotherapy (IMRT) beams and 60 Smart Arcs. Monitor unit (MU) calculated in each of the software were compared with the TPS and the values were represented as a percent difference. Open fields were calculated as a baseline for each software’s accuracy using 5×5, 10×10, 20×20, 30×30, and 40×40 fields. Box plots, Pearson correlation, and Bland-Altman analysis were used for comparison of the results. Results: The baseline accuracy was established to within 0.6%, −1.4%, −0.2%, and −1.0% for Diamond, IMSure,MUcheck, and Radcalc, respectively. In the clinical data, the dose difference represented as mean ± 1 standard deviation were 0.7%±0.1%, −0.3%±0.1%, −1.5%±0.1%, and 0.4%±0.0% for Diamond, IMSure, MUcheck, and Radcalc, respectively Conclusion: The implementation of Clarkson algorithm for the dose calculation between each of the software in question can vary considerably. The currently used second check software, Radcalc has shown the best agreement on average, variance, and smallest percent range from Pinnacle3 TPS values. The closest in average percent difference from the TPS data was the IMSure software, but has significantly larger variance and percent range. The mean percent differences in Diamond and MUcheck were significantly larger than Radcalc and IMSure.

60 APPLIED LIFE SCIENCES↗

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗

Highly resolved Eulerian simulations of fuel spray transients in single and multi-hole injectors: Nozzle flow and near-exit dynamics

In high pressure fuel injectors, needle opening and closing transients cause complex off-design fluid dynamics behaviors that profoundly impact the spray and mixture formation processes. These dynamics are completely different from what is known to occur in steady state conditions. In this study, diesel spray transients have been investigated in single-hole and 3-hole nozzles, encompassing internal and external nozzle flow and including needle motion, performing highly resolved (2.5 μm) computational fluid dynamics (CFD) simulations. We focused on end-of-injection (EOI) and start-of-injection (SOI) processes, in order to provide insights in to the physics. The liquid fuel, vapor and gas species are modeled with a single-fluid multiphase mixture approach, with diffuse interface, and with large eddy simulations (LES) of the turbulence. Occurrence of phase change due to cavitation is accounted for, and the spray dispersion is described with a turbulent dispersion model. Detailed needle motion data and orifice internal surface are available from xray synchrotron source measurements carried out at Argonne National Laboratory, and shared through the Engine Combustion Network (ECN) community. Simulations are compared against x-ray phase contrast imaging and radiography of the internal and near-exit flow, in addition to optical microscopy data of the near-exit sprays. Simulation results are found to agree well with available experimental data, and are able to realistically capture local and global features. The simulations allow to gain insight into the physics of gas ingestion and dribbles at EOI, for different hole diameters, operating conditions and number of holes. At SOI, timing of liquid appearance out of the injector and spray tip penetration are adequately predicted, by using the EOI flow field as in-nozzle initialization, and by prescribing the measured tip needle displacement with an informed effective valve opening point inferred from the x-ray observations. Lastly, the variation of spreading angle over time is also discussed in detail for the multi-hole case, including hole-to-hole variations. Due to real geometry features and asymmetric needle motion with eccentric components, it is found that the three holes exhibit swirling flows of increasing intensity as the lift decreases, causing the near cone angle to open and spread, in a quasi-hollow cone structure. These features are not observed in axial single-hole injectors because of their relative simplicity and intrinsic symmetry.

42 ENGINEERING↗

Control-Oriented Modeling of Lithium-Ion Batteries

Battery management systems (BMSs), which monitor and optimize performance while ensuring safety, require control-oriented models, i.e., models tailored to the design and implementation of estimation and control algorithms. Physics-based electrochemical models describe detailed battery phenomena, but are too computationally intensive for use in estimation and control applications. Single particle models (SPMs), which retain some of the physics of electrochemical models, are often used for control-oriented battery modeling since they are computationally efficient; however, they are only valid over very low frequency ranges and C-rates. Empirical equivalent circuit models (ECMs) are also used for control-oriented battery modeling since they are computationally efficient and can describe battery behavior over wide frequency ranges; however, they provide no physical understanding of the battery and, therefore, have limited applicability. Further, fractional order terms (e.g., Warburg impedances) are often employed, making the models unwieldy for use in the time domain. This work provides a control-oriented battery model that combines the benefits of SPM and ECM models, while overcoming their limitations. The proposed model incorporates some of the battery physics found in electrochemical models, can easily be used in both the time and frequency domains, and describes battery behavior over its entire frequency range. A linearized single particle model, which incorporates key electrochemical parameters, is used for modeling battery physics at very low frequencies. For low frequencies, integer-order linear systems are used to approximate diffusion physics described by Warburg impedances, and high frequency behavior is modeled by the double layer capacitance effect. The proposed battery model is more computationally efficient than full electrochemical models since it does not require the solution of PDEs, is accurate for a wider frequency range than the SPM considered in this paper, and does not suffer from the unwieldiness and limited applicability of empirical ECMs. Finally, the model is validated in the time and frequency domains via a comparison to pseudo-two-dimensional (P2D) model simulations and experimental data.

25 ENERGY STORAGE↗

End-to-End Pipeline for Trigger Detection on Hit and Track Graphs

There has been a surge of interest in applying deep learning in particle and nuclear physics to replace labor-intensive offline data analysis with automated online machine learning tasks. This paper details a novel AI-enabled triggering solution for physics experiments in Relativistic Heavy Ion Collider and future Electron-Ion Collider. The triggering system consists of a comprehensive end-to-end pipeline based on Graph Neural Networks that classifies trigger events versus background events, makes online decisions to retain signal data, and enables efficient data acquisition. Here, the triggering system first starts with the coordinates of pixel hits lit up by passing particles in the detector, applies three stages of event processing (hits clustering, track reconstruction, and trigger detection), and labels all processed events with the binary tag of trigger versus background events. By switching among different objective functions, we train the Graph Neural Networks in the pipeline to solve multiple tasks: the edge-level track reconstruction problem, the edge-level track adjacency matrix prediction, and the graph-level trigger detection problem. We propose a novel method to treat the events as track-graphs instead of hit-graphs. This method focuses on intertrack relations and is driven by underlying physics processing. As a result, it attains a solid performance (around 72% accuracy) for trigger detection and outperforms the baseline method using hit-graphs by 2% higher accuracy.

97 MATHEMATICS AND COMPUTING↗

PARLO: PArallel Run-Time Layout Optimization for Scientific Data Explorations with Heterogeneous Access Patterns

The size and scope of cutting-edge scientific simulations are growing much faster than the I/O and storage capabilities of their run-time environments. The growing gap is exacerbated by exploratory, data-intensive analytics, such as querying simulation data with multivariate, spatio-temporal constraints, which induces heterogeneous access patterns that stress the performance of the underlying storage system. Previous work addresses data layout and indexing techniques to improve query performance for a single access pattern, which is not sufficient for complex analytics jobs. We present PARLO a parallel run-time layout optimization framework, to achieve multi-level data layout optimization for scientific applications at run-time before data is written to storage. The layout schemes optimize for heterogeneous access patterns with user-specified priorities. PARLO is integrated with ADIOS, a high-performance parallel I/O middleware for large-scale HPC applications, to achieve user-transparent, light-weight layout optimization for scientific datasets. It offers simple XML-based configuration for users to achieve flexible layout optimization without the need to modify or recompile application codes. Experiments show that PARLO improves performance by 2 to 26 times for queries with heterogeneous access patterns compared to state-of-the-art scientific database management systems. Compared to traditional post-processing approaches, its underlying run-time layout optimization achieves a 56% savings in processing time and a reduction in storage overhead of up to 50%. PARLO also exhibits a low run-time resource requirement, while also limiting the performance impact on running applications to a reasonable level.

Computer Science↗

De-noising drift chambers in CLAS12 using convolutional auto encoders

Modern Nuclear Physics experimental setups run experiments with higher beam intensity resulting in increased noise in detector components used for particle track reconstruction. Increased uncorrelated signals (noise) result in decreased particle reconstruction efficiency. In this paper, we investigate the usage of Machine Learning, specifically Convolutional Neural Network Auto-Encoders (CAE), for de-noising raw hits from drift chambers in the CLAS12 detector. To the best of our knowledge, this is the first time CAE is employed to perform such an operation in this field. During the de-noising phase, it is important to remove as much noise as possible while retaining the valid hits to avoid losing crucial information about the experiment. Here, we show that using CAE, it is possible to remove noise hits while retaining up to 94% of valid tracks for a beam current of 110nA while for lower beam currents (45-55nA), we get up to 98% efficiency. Studies on experimental conditions with increasing noise show that CAE performs better than conventional tracking algorithms in isolating hits belonging to tracks. Specifically, the de-noising algorithm results in tracking efficiency improvements greater than 15%, in real data production procedures with nominal conditions, and up to two times better efficiency in synthetically generated data with high luminosity conditions (90-110nA), indicating that machine learning can lead to significantly shorter times for conducting physics experiments.

97 MATHEMATICS AND COMPUTING↗

Computational Imaging for Intelligence in Highly Scattering Aerosols (Final Report)

Natural and man-made degraded visual environments pose major threats to national security. The random scattering and absorption of light by tiny particles suspended in the air reduces situational awareness and causes unacceptable down-time for critical systems and operations. To improve the situation, we have developed several approaches to interpret the information contained within scattered light to enhance sensing and imaging in scattering media. These approaches were tested at the Sandia National Laboratory Fog Chamber facility and with tabletop fog chambers. Computationally efficient light transport models were developed and leveraged for computational sensing. The models are based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes. After the new model was experimentally validated, statistical approaches for detection, localization, and imaging of objects hidden in fog were developed and demonstrated. A binary hypothesis test and the Neyman-Pearson lemma provided the highest theoretically possible probability of detection for a specified false alarm rate and signal-to-noise ratio. Maximum likelihood estimation allowed estimation of the fog optical properties as well as the position, size, and reflection coefficient of an object in fog. A computational dehazing approach was implemented to reduce the effects of scatter on images, making object features more readily discernible. We have developed, characterized, and deployed a new Tabletop Fog Chamber capable of repeatably generating multiple unique fog-analogues for optical testing in degraded visual environments. We characterized this chamber using both optical and microphysical techniques. In doing so we have explored the ability of droplet nucleation theory to describe the aerosols generated within the chamber, as well as Mie scattering theory to describe the attenuation of light by said aerosols, and correlated the aerosol microphysics to optical properties such as transmission and meteorological optical range (MOR). This chamber has proved highly valuable and has supported multiple efforts inclusive to and exclusive of this LDRD project to test optics in degraded visual environments. Circularly polarized light has been found to maintain its polarization state better than linearly polarized light when propagating through fog. This was demonstrated experimentally in both the visible and short-wave infrared (SWIR) by imaging targets made of different commercially available retroreflective films. It was found that active circularly polarized imaging can increase contrast and range compared to linearly polarized imaging. We have completed an initial investigation of the capability for machine learning methods to reduce the effects of light scattering when imaging through fog. Previously acquired experimental long-wave images were used to train an autoencoder denoising architecture. Overfitting was found to be a problem because of lack of variability in the object type in this data set. The lessons learned were used to collect a well labeled dataset with much more variability using the Tabletop Fog Chamber that will be available for future studies. We have developed several new sensing methods using speckle intensity correlations. First, the ability to image moving objects in fog was shown, establishing that our unique speckle imaging method can be implemented in dynamic scattering media. Second, the speckle decorrelation over time was found to be sensitive to fog composition, implying extensions to fog characterization. Third, the ability to distinguish macroscopically identical objects on a far-subwavelength scale was demonstrated, suggesting numerous applications ranging from nanoscale defect detection to security. Fourth, we have shown the capability to simultaneously image and localize hidden objects, allowing the speckle imaging method to be effective without prior object positional information. Finally, an interferometric effect was presented that illustrates a new approach for analyzing speckle intensity correlations that may lead to more effective ways to localize and image moving objects. All of these results represent significant developments that challenge the limits of the application of speckle imaging and open important application spaces. A theory was developed and simulations were performed to assess the potential transverse resolution benefit of relative motion in structured illumination for radar systems. Results for a simplified radar system model indicate that significant resolution benefits are possible using data from scanning a structured beam over the target, with the use of appropriate signal processing.

58 GEOSCIENCES↗

Comparison of Gas Phase Fragmentation Behaviors of Nuclear Fuel Cycle Ligands in Lanthanide and Americium Metal Ligand Nitrate Clusters

Introduction (120 words): Transport of metal ions across the aqueous-organic phase boundary is an essential step in a hydrometallurgical nuclear fuel reprocessing strategy. The study of transport agents for nuclear fuel elements is imperative to guide the design of ligands that boost the separation efficiency of the recovery process from fission products. However, limited studies have been made on the chemistry of these transport agents when complexing with transuranic elements in gas-phase where all surrounding factors are essentially excluded. This work investigates the reagent ligand complexations to transuranic and other metals and their dissociations in the gas phase. Comparisons are made between 4f and 5f elements and between ligands. Methods (120 words): (N,N-diisobutylcarbamoylmethyl)phenyloctylphosphine oxide (CMPO) and N,N,N',N'-tetraoctyldiglycolamide (TODGA) have been selected to complex with metal nitrates. The actinide americium and lanthanides neodymium, samarium, and europium were investigated as part of this work. The lanthanides were selected to act as size and electron configuration analogues of the minor actinides. Metal complexes with two ligands and two nitrates ([M(NO3)2(CMPO)2]+, for example) are studied in Bruker micrOTOF-Q II mass spectrometer equipped with collision-induced dissociation capability. The comparisons of the mass spectra are made in groups of homogenous ligands and mixed TODGA-CMPO ligands clusters. Comparisons are also made based on the complexed metals (Am and lanthanides). Preliminary data (300 words): Collision-induced dissociation mass spectrometry data are collected on two ligands complexed with metal nitrates where the two ligands are homogenous, with (CMPO)2 or (TODGA)2, or heterogeneous, with (TODGA)(CMPO). Several fragmentation patterns are observed among complexes with the CMPO ligand whereas the TODGA ligand commonly dissociates intact from the complex. Most of the metal complexes exhibit similar fragmentation patterns, but there are a few notable deviations in fragmentation patterns between the Am and Ln-bearing complexes. For the [M(CMPO)2(NO3)2]+ complexes, the initial loss of nitrate in the form of nitric acid is observed in all four complexes. However, [Am(CMPO)2(NO3)2]+ exhibits an additional fragmentation not found in the lanthanide complexes. Also, a significantly different ratio of the second nitric acid loss is found in the Am complex. These deviations may indicate the different interaction behaviors between actinides and lanthanides. The [M(TODGA)2(NO3)2]+ complexes exhibit the fragmentation as the loss of one TODGA ligand as an intact form and the loss of nitrate as nitric acid. The Am complex exhibits an additional fragmentation after losing the TOGDA ligand, which is not observed among the Ln complexes. The heterogeneous [M(TODGA)(CMPO)(NO3)2]+ complexes exhibit both similarities and differences between the Am and Ln complexes. For example, the heterogenous Am complex does not exhibit the loss of an intact TODGA ligand while all three Ln complexes do. This indicates that TODGA may bind more strongly to Am than Ln. Additionally, the intensity of the loss of CMPO ligand (as partially or whole) is found to be significantly larger than that of the loss of TODGA (as partially of whole) indicating that TODGA is bound to the metal significantly stronger than CMPO. Planned computational analysis will help understand the deviation in fragmentation behaviors between americium and lanthanide metal centers, or between TODGA and CMPO ligands. Novel aspect (20 words): Gas-phase actinide and lanthanide complex formation and fragmentation provide insight into the coordination environment differences of f-element metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of Gas Phase Fragmentation Behaviors of Nuclear Fuel Cycle Ligands in Lanthanide and Americium Metal Ligand Nitrate Clusters

Title (20 word): Comparison of Gas Phase Fragmentation Behaviors of Nuclear Fuel Cycle Ligands in Lanthanide and Americium Metal Ligand Nitrate Clusters Introduction (120 words): Transport of metal ions across the aqueous-organic phase boundary is an essential step in a hydrometallurgical nuclear fuel reprocessing strategy. The study of transport agents for nuclear fuel elements is imperative to guide the design of ligands that boost the separation efficiency of the recovery process from fission products. However, limited studies have been made on the chemistry of these transport agents when complexing with transuranic elements in gas-phase where all surrounding factors are essentially excluded. This work investigates the reagent ligand complexations to transuranic and other metals and their dissociations in the gas phase. Comparisons are made between 4f and 5f elements and between ligands. Methods (120 words): (N,N-diisobutylcarbamoylmethyl)phenyloctylphosphine oxide (CMPO) and N,N,N',N'-tetraoctyldiglycolamide (TODGA) have been selected to complex with metal nitrates. The actinide americium and lanthanides neodymium, samarium, and europium were investigated as part of this work. The lanthanides were selected to act as size and electron configuration analogues of the minor actinides. Metal complexes with two ligands and two nitrates ([M(NO3)2(CMPO)2]+, for example) are studied in Bruker micrOTOF-Q II mass spectrometer equipped with collision-induced dissociation capability. The comparisons of the mass spectra are made in groups of homogenous ligands and mixed TODGA-CMPO ligands clusters. Comparisons are also made based on the complexed metals (Am and lanthanides). Preliminary data (300 words): Collision-induced dissociation mass spectrometry data are collected on two ligands complexed with metal nitrates where the two ligands are homogenous, with (CMPO)2 or (TODGA)2, or heterogeneous, with (TODGA)(CMPO). Several fragmentation patterns are observed among complexes with the CMPO ligand whereas the TODGA ligand commonly dissociates intact from the complex. Most of the metal complexes exhibit similar fragmentation patterns, but there are a few notable deviations in fragmentation patterns between the Am and Ln-bearing complexes. For the [M(CMPO)2(NO3)2]+ complexes, the initial loss of nitrate in the form of nitric acid is observed in all four complexes. However, [Am(CMPO)2(NO3)2]+ exhibits an additional fragmentation not found in the lanthanide complexes. Also, a significantly different ratio of the second nitric acid loss is found in the Am complex. These deviations may indicate the different interaction behaviors between actinides and lanthanides. The [M(TODGA)2(NO3)2]+ complexes exhibit the fragmentation as the loss of one TODGA ligand as an intact form and the loss of nitrate as nitric acid. The Am complex exhibits an additional fragmentation after losing the TOGDA ligand, which is not observed among the Ln complexes. The heterogeneous [M(TODGA)(CMPO)(NO3)2]+ complexes exhibit both similarities and differences between the Am and Ln complexes. For example, the heterogenous Am complex does not exhibit the loss of an intact TODGA ligand while all three Ln complexes do. This indicates that TODGA may bind more strongly to Am than Ln. Additionally, the intensity of the loss of CMPO ligand (as partially or whole) is found to be significantly larger than that of the loss of TODGA (as partially of whole) indicating that TODGA is bound to the metal significantly stronger than CMPO. Planned computational analysis will help understand the deviation in fragmentation behaviors between americium and lanthanide metal centers, or between TODGA and CMPO ligands. Novel aspect (20 words): Gas-phase actinide and lanthanide complex formation and fragmentation provide insight into the coordination environment differences of f-element metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Verification of the ENDF/B-VII.1 Based MC 2 -3 Library Rev.1

The MC 2 -3 code, developed by Argonne National Laboratory under the DOE-NE NEAMS program, is a multigroup cross section generation code for fast reactor applications. Last year, the ENDF/B-VII.0 (E70) MC 2 -3 library, which has been extensively used, verified, and validated over a long period, was intensively reverified and updated to support the commercial grade dedication (CGD) requirement of the TerraPower Natrium project. This year, the ENDF/B-VII.1 (E71) MC 2 -3 library, the preliminary version of which was generated several years ago, was regenerated and rigorously verified to support the Natrium project as well as the completion of verification of the E71 library. The E71 library was verified using the process developed during the verification of the E70 library, including comparisons of cross sections with the NJOY-generated cross sections, comparisons of the resolved resonance cross sections with those using the PEDNF library, and comparison of total cross sections with the sum of partial cross sections. Additional verifications were conducted to ensure that the benchmark problem solutions with the E71 library are reasonable compared to the corresponding Monte Carlo solutions. Furthermore, the E71 gamma library was generated, which includes data for prompt gamma, delayed gamma, and delayed beta as well as neutron and gamma heating. The gamma library was verified at the level of individual isotopes. The EBR-II core solutions from MC 2 -3/ DIF3D and MCNP were compared, demonstrating that those solutions in terms of k-effective and assembly powers were in good agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Direct Observation of Halide Migration and its Effect on the Photoluminescence of Methylammonium Lead Bromide Perovskite Single Crystals

Optoelectronic devices based on hybrid perovskites have demonstrated outstanding performance within a few years of intense study. However, commercialization of these devices requires barriers to their development to be overcome, such as their chemical instability under operating conditions. In order to investigate this instability and its consequences, the electric field applied to single crystals of methylammonium lead bromide (CH 3 NH 3 PbBr 3 ) is varied, and changes are mapped in both their elemental composition and photo-luminescence. Synchrotron-based nanoprobe X-ray fluorescence (nano-XRF) with 250 nm resolution reveals quasi-reversible field-assisted halide migration, with corresponding changes in photoluminescence. It is observed that higher local bromide concentration is correlated to superior optoelectronic performance in CH 3 NH 3 PbBr 3 . A lower limit on the electromigration rate is calculated from these experiments and the motion is interpreted as vacancy-mediated migration based on nudged elastic band density functional theory (DFT) simulations. The XRF mapping data provide direct evidence of field-assisted ionic migration in a model hybrid-perovskite thin single crystal, while the link with photoluminescence proves that the halide stoichiometry plays a key role in the optoelectronic properties of the perovskite.

36 MATERIALS SCIENCE↗

Change-point Detection and Image Segmentation for Time Series of Astrophysical Images

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multiband image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all of the change points (in time domain), as well as all of the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real data sets show that the proposed method enjoys very promising empirical properties. Applications to two real data sets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiation-induced refraction artifacts in the optical CT readout of polymer gel dosimeters

Purpose: The objective of this work is to demonstrate imaging artifacts that can occur during the optical computed tomography (CT) scanning of polymer gel dosimeters due to radiation-induced refractive index (RI) changes in polyacrylamide gels. Methods: A 1 L cylindrical polyacrylamide gel dosimeter was irradiated with 3 × 3 cm{sup 2} square beams of 6 MV photons. A prototype fan-beam optical CT scanner was used to image the dosimeter. Investigative optical CT scans were performed to examine two types of rayline bending: (i) bending within the plane of the fan-beam and (ii) bending out the plane of the fan-beam. To address structured errors, an iterative Savitzky–Golay (ISG) filtering routine was designed to filter 2D projections in sinogram space. For comparison, 2D projections were alternatively filtered using an adaptive-mean (AM) filter. Results: In-plane rayline bending was most notably observed in optical CT projections where rays of the fan-beam confronted a sustained dose gradient that was perpendicular to their trajectory but within the fan-beam plane. These errors caused distinct streaking artifacts in image reconstructions due to the refraction of higher intensity rays toward more opaque regions of the dosimeter. Out-of-plane rayline bending was observed in slices of the dosimeter that featured dose gradients perpendicular to the plane of the fan-beam. These errors caused widespread, severe overestimations of dose in image reconstructions due to the higher-than-actual opacity that is perceived by the scanner when light is bent off of the detector array. The ISG filtering routine outperformed AM filtering for both in-plane and out-of-plane rayline errors caused by radiation-induced RI changes. For in-plane rayline errors, streaks in an irradiated region (>7 Gy) were as high as 49% for unfiltered data, 14% for AM, and 6% for ISG. For out-of-plane rayline errors, overestimations of dose in a low-dose region (∼50 cGy) were as high as 13 Gy for unfiltered data, 10 Gy for AM, and 3.1 Gy for ISG. The ISG routine also addressed unrelated artifacts that previously needed to be manually removed in sinogram space. However, the ISG routine blurred reconstructions, causing losses in spatial resolution of ∼5 mm in the plane of the fan-beam and ∼8 mm perpendicular to the fan-beam. Conclusions: This paper reveals a new category of imaging artifacts that can affect the optical CT readout of polyacrylamide gel dosimeters. Investigative scans show that radiation-induced RI changes can cause significant rayline errors when rays confront a prolonged dose gradient that runs perpendicular to their trajectory. In fan-beam optical CT, these errors manifested in two ways: (1) distinct streaking artifacts caused by in-plane rayline bending and (2) severe overestimations of opacity caused by rays bending out of the fan-beam plane and missing the detector array. Although the ISG filtering routine mitigated these errors better than an adaptive-mean filtering routine, it caused unacceptable losses in spatial resolution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Global Carbon Budget 2017

Here an accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the global carbon budget – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO 2 emissions from fossil fuels and industry ( E FF ) are based on energy statistics and cement production data, respectively, while emissions from land-use change ( E LUC ), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO 2 concentration is measured directly and its rate of growth ( G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink ( S OCEAN ) and terrestrial CO 2 sink ( S LAND ) are estimated with global process models constrained by observations. The resulting carbon budget imbalance ( B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1 σ . For the last decade available (2007–2016), E FF was 9.4 ± 0.5 GtC yr –1 , E LUC 1.3 ± 0.7 GtC yr –1 , G ATM 4.7 ± 0.1 GtC yr –1 , S OCEAN 2.4 ± 0.5 GtC yr –1 , and S LAND 3.0 ± 0.8 GtC yr –1 , with a budget imbalance B IM of 0.6 GtC yr –1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in E FF was approximately zero and emissions remained at 9.9 ±\ 0.5 GtC yr –1 . Also for 2016, E LUC was 1.3 ± .7 GtC yr –1 , G ATM was 6.1 ± 0.2 GtC yr –1 , S OCEAN was 2.6 ± 0.5 GtC yr –1 , and S LAND was 2.7 ± 1.0 GtC yr –1 , with a small B IM of –0.3 GtC. G ATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the high fossil emissions and the small S LAND consistent with El Niño conditions. The global atmospheric CO 2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6–9 months indicate a renewed growth in E FF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set.

54 ENVIRONMENTAL SCIENCES↗

Refinement for single-nanoparticle structure determination from low-quality single-shot coherent diffraction data

With the emergence of X-ray free-electron lasers, it is possible to investigate the structure of nanoscale samples by employing coherent diffractive imaging in the X-ray spectral regime. In this work, we developed a refinement method for structure reconstruction applicable to low-quality coherent diffraction data. The method is based on the gradient search method and considers the missing region of a diffraction pattern and the small number of detected photons. We introduced an initial estimate of the structure in the method to improve the convergence. The present method is applied to an experimental diffraction pattern of an Xe cluster obtained in an X-ray scattering experiment at the SPring-8 Angstrom Compact free-electron LAser (SACLA) facility. It is found that the electron density is successfully reconstructed from the diffraction pattern with a large missing region, with a good initial estimate of the structure. The diffraction pattern calculated from the reconstructed electron density reproduced the observed diffraction pattern well, including the characteristic intensity modulation in each ring. Our refinement method enables structure reconstruction from diffraction patterns under difficulties such as missing areas and low diffraction intensity, and it is potentially applicable to the structure determination of samples that have low scattering power.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗