Search NASA⌕ Search

SEARCH · Search NASA

Results for “Open Source Open Science”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Simulations of the WFIRST Supernova Survey and Forecasts of Cosmological Constraints

The Wide Field InfraRed Survey Telescope (WFIRST) was the highest-ranked large space-based mission of the 2010 New Worlds, New Horizons decadal survey. It is now a NASA mission in formulation with a planned launch in the mid 2020s. A primary mission objective is to precisely constrain the nature of dark energy through multiple probes, including Type Ia supernovae (SN Ia). Here, we present the first realistic simulations of the WFIRST SN survey based on current hardware specifications and using open-source tools. We simulate SN light curves and spectra as viewed by the WFIRST wide-field channel (WFC) imager and integral field channel (IFC) spectrometer, respectively. We examine 11 survey strategies with different time allocations between the WFC and IFC, two of which are based upon the strategy described by the WFIRST Science Definition Team, which measures SN distances exclusively from IFC data. We propagate statistical and, crucially, systematic uncertainties to predict the Dark Energy Task Force figure of merit (FoM) for each strategy. Of the strategies investigated, we find the most successful to be WFC focused. However, further work in constraining systematics is required to fully optimize the use of the IFC. Even without improvements to other cosmological probes, the WFIRST SN survey has the potential to increase the FoM by more than an order of magnitude from the current values. Although the survey strategies presented here have not been fully optimized, these initial investigations are an important step in the development of the final hardware design and implementation of the WFIRST mission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

BMX: Biological modelling and interface exchange

Abstract High performance computing has a great potential to provide a range of significant benefits for investigating biological systems. These systems often present large modelling problems with many coupled subsystems, such as when studying colonies of bacteria cells. The aim to understand cell colonies has generated substantial interest as they can have strong economic and societal impacts through their roles in in industrial bioreactors and complex community structures, called biofilms, found in clinical settings. Investigating these communities through realistic models can rapidly exceed the capabilities of current serial software. Here, we introduce BMX, a software system developed for the high performance modelling of large cell communities by utilising GPU acceleration. BMX builds upon the AMRex adaptive mesh refinement package to efficiently model cell colony formation under realistic laboratory conditions. Using simple test scenarios with varying nutrient availability, we show that BMX is capable of correctly reproducing observed behavior of bacterial colonies on realistic time scales demonstrating a potential application of high performance computing to colony modelling. The open source software is available from the zenodo repository https://doi.org/10.5281/zenodo.8084270 under the BSD-2-Clause licence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

International Neuroscience Initiatives through the Lens of High-Performance Computing

Neuroscience initiatives aim to develop new technologies and tools to measure and manipulate neuronal circuits. In this paper, to deal with the massive amounts of data generated by these tools, the authors envision the co-location of open data repositories in standardized formats together with high-performance computing hardware utilizing open source optimized analysis codes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Engineering-scale Modeling of High-Temperature Creep and Creep Crack Growth in Alloy 316H

This document demonstrates completion of the goals described in the technical narrative of the Department of Energy’s Industry Funding Opportunity Announcement (iFOA) project entitled “Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing,” which relates to the development and demonstration of capabilities for conducting engineering-scale simulations of Alloy 316H components under high-temperature conditions. This work encompassed two major aspects: integrating and testing constitutive models for the creep response of 316H at high temperatures, and developing tools for modeling creep crack growth (CCG) in 316H. Two classes of constitutive models were used in this effort: phenomenological models based on behavior observed at the engineering scale, and reduced-order models (ROMs) that represent the nonlinear response of mesoscale models that capture the sensitivity to material microstructure and processing. Likewise, the approaches employed for CCG modeling considered both simplified engineering approaches and detailed simulations of creep and damage ahead of the crack tip. These developments, which were performed by utilizing the Grizzly code as well as the open-source libraries it depends on, strengthen Grizzly’s ability to support licensing and safety analyses of high-temperature reactor components.

316H↗
DOE OSTI2025

Source Attribution of Methane Enhancements at the Southern Great Plains Atmospheric Radiation Measurement Site Field Campaign Report

Tower data from 2017 to 2020 revealed that the U.S. Department of Energy’s Southern Great Plains (SGP) Atmospheric Radiation Measurements (ARM) observatory in Oklahoma exhibited significant CH 4 enhancements compared to other U.S. tower sites. On average, near-surface CH 4 at the ARM site was higher, had larger diurnal and seasonal variations, and had sharper near-surface increases at night compared to other tower sites in the U.S. A field campaign was conducted in June 2024 to identify and quantify the surrounding CH 4 emission sources by two researchers from the University of Oklahoma. During the campaign, a mobile measurement was conducted using a LI-COR 7810 tracer gas analyzer mounted on a vehicle during June 23-26, 2024 near the SGP site (within 7 km distance). We sampled CH 4 along predefined routes surrounding a suspected source, i.e., a large animal feeding operation (AFO) farm 5.8 km northwest of the SGP ARM site. Measurements were collected between 04:30 and 09:00 local time when the planetary boundary layer (PBL) was shallow, increasing sensitivity to surface emissions. Driving speed was kept below 5 mph to minimize pressure-related artifacts. CH 4 concentrations downwind of the AFO exceeded 6000 ppb on the early morning of June 26, corresponding to an emission rate of ~95 kg⋅hr -1 estimated using the mass balance method. Additional sources identified during the campaign include nearby open-range cattle and leaks from a nearby natural gas pipeline.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Multi-scale Molecular Systems Biology: Reconstruction and Model Optimization (Final Report)

Within the Biological and Environmental Research program, the goal of the Genomic Science program is to achieve a predictive, systems-level understanding of plants, microbes, and biological communities to enable biobased solutions to DOE mission challenges in energy and environment. To achieve this goal, an ultimate objective is to develop the knowledgebase, computational infrastructure, and modeling capabilities to advance predictive understanding and manipulation of biological systems. Our project has increased the modeling capabilities and computational infrastructure required to advance predictive understanding and manipulation of biological systems. In particular, we have developed novel algorithms, for modelling biochemical reaction kinetics at genome-scale. These new algorithms have been publicly disseminated as open source software with extensive documentation and tutorials that enable widespread application of this advanced software to practical modelling scenarios. The technical effectiveness of all novel software has been conclusively demonstrated and the quality of this software is maintained by a dedicated continuous integration system. Ultimately, the results of this project will benefit the public because it will enable biologists to more accurately computationally model biological systems that need to be understood to address DOE challenges in energy and environment.

09 BIOMASS FUELS↗
DOE OSTI2018

Hybrid Polymer/Metal–Organic Framework Films for Colorimetric Water Sensing over a Wide Concentration Range [Hybrid Polymer/MOF Films for Colorimetric Water Sensing over a Wide Concentration Range]

Because of their extraordinary surface areas and tailorable porosity, metal–organic frameworks (MOFs) have the potential to be excellent sensors of gas-phase analytes. MOFs with open metal sites are particularly attractive for detecting Lewis basic atmospheric analytes, such as water. Here, we demonstrate that thin films of the MOF HKUST-1 can be used to quantitatively determine the relative humidity (RH) of air using a colorimetric approach. HKUST-1 thin films are spin-coated onto rigid or flexible substrates and are shown to quantitatively determine the RH within the range of 0.1–5% RH by either visual observation or a straightforward optical reflectivity measurement. At high humidity (>10% RH), a polymer/MOF bilayer is used to slow the transport of H 2 O to the MOF film, enabling quantitative determination of RH using time as the distinguishing metric. Lastly, the sensor is combined with an inexpensive light-emitting diode light source and Si photodiode detector to demonstrate a quantitative humidity detector for low humidity environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LandScan Mosaic Rapid Population Update: Jamaica After Hurricane Melissa (V1)

During a natural disaster such as Hurricane Melissa, understanding where people are located is critical for situational awareness, operational planning and humanitarian support and consequence assessment. Traditional population datasets focus on mapping populations based on residential, or "business-as-usual" scenarios. However, natural disasters can create disruptions in daily routines of population in addition to the magnitude of the population displacement, depending on the type, duration, context, and location of the event. The Geospatial Science and Human Security Division at Oak Ridge National Laboratory (ORNL) produced this latest LandScan Mosaic Rapid Population Update for Jamaica following Hurricane Melissa, a category 5 hurricane that made landfall on Jamaica on October 28 2025. This Rapid Population Update captures the immediate population displacement following the hurricane using a combination of open-source building damage assessment data from Microsoft, flood exposure data from the Global Flood Monitoring service, reported population displacement information, and humanitarian shelter locations from the Jamaican Office of Disaster Preparedness and Emergency Management and the underlying LandScan Mosaic Jamaica as a base population.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Genome-Scale 13 C Fluxomics Modeling for Metabolic Engineering of Saccharomyces cerevisiae .

Synthetic biology is a rapidly developing field that pursues the application of engineering principles and development approaches to biological engineering. Synthetic biology is poised to change the way biology is practiced, and has important practical applications: for example, building genetically engineered organisms to produce biofuels, medicines, and other chemicals. Traditionally, synthetic biology has focused on manipulating a few genes (e.g., in a single pathway or genetic circuit), but its combination with systems biology holds the promise of creating new cellular architectures and constructing complex biological systems from the ground up. Enabling this merge of synthetic and systems biology will require greater predictive capability for modeling the behavior of cellular systems, and more comprehensive data sets for building and calibrating these models. The so-called "-omics" data sets can now be generated via high throughput techniques in the form of genomic, proteomic, transcriptomic, and metabolomic information on the engineered biological system. Of particular interest with respect to the engineering of microbes capable of producing biofuels and other chemicals economically and at scale are metabolomic datasets, and their insights into intracellular metabolic fluxes. Metabolic fluxes provide a rapid and easy to understand picture of how carbon and energy flow throughout the cell. Furthermore, we present a detailed guide to performing metabolic flux analysis and modeling using the open source JBEI Quantitative Metabolic Modeling (jQMM) library. This library allows the user to transform metabolomics data in the form of isotope labeling data from a 13 C labeling experiment into a determination of cellular fluxes that can be used to develop genetic engineering strategies for metabolic engineering.The jQMM library presents a complete toolbox for performing a range of different tasks of interest in metabolic engineering. Various different types of flux analysis and modeling can be performed such as flux balance analysis, 13 C metabolic flux analysis, and two-scale 13 C metabolic flux analysis (2S- 13 C MFA). 2S- 13 C MFA is a novel method that determines genome-scale fluxes without the need of every single carbon transition in the metabolic network. In addition to several other capabilities, the jQMM library can make model based predictions for how various genetic engineering strategies can be incorporated toward bioengineering goals: it can predict the effects of reaction knockouts on metabolism using both the MoMA and ROOM methodologies. In this chapter, we will illustrate the use of the jQMM library through a step-by-step demonstration of flux determination and knockout prediction in a complex eukaryotic model organism: Saccharomyces cerevisiae ( S. cerevisiae ). Included with this chapter is a digital Jupyter Notebook file that provides a computable appendix showing a self-contained example of jQMM usage, which can be changed to fit the user's specific needs. As an open source software project, users can modify and extend the code base to make improvements at will, allowing them to share their development work and contribute back to the jQMM modeling community.

-Omics data↗
DOE OSTI2023

Entangled Photon Correlations Allow a Continuous-Wave Laser Diode to Measure Single-Photon, Time-Resolved Fluorescence

Fluorescence lifetime experiments are a standard approach for measuring excited state dynamics and local environment effects. Here, we show that entangled photon pairs produced from a continuous-wave (CW) laser diode can replicate pulsed laser experiments without phase modulation. As a proof of principle, picosecond fluorescence lifetimes of indocyanine green are measured in multiple environments. Here, the use of entangled photons has three unique advantages. First, low power CW laser diodes and entangled photon source design lead to straightforward on-chip integration for a direct path to distributable fluorescence lifetime measurements. Second, the entangled pair wavelength is easily tuned by temperature or electric field, allowing a single source to cover octave bandwidths. Third, femtosecond temporal resolutions can be reached without requiring major advances in source technology or external phase modulation. Entangled photons could therefore provide increased accessibility to time-resolved fluorescence while also opening new scientific avenues in photosensitive and inherently quantum systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The QICK (Quantum Instrumentation Control Kit): Readout and control for qubits and detectors

We introduce a Xilinx RF System-on-Chip (RFSoC)-based qubit controller (called the Quantum Instrumentation Control Kit, or QICK for short), which supports the direct synthesis of control pulses with carrier frequencies of up to 6 GHz. The QICK can control multiple qubits or other quantum devices. The QICK consists of a digital board hosting an RFSoC field-programmable gate array, custom firmware, and software and an optional companion custom-designed analog front-end board. We characterize the analog performance of the system as well as its digital latency, important for quantum error correction and feedback protocols. We benchmark the controller by performing standard characterizations of a transmon qubit. We achieve an average gate fidelity of ℱ avg =99.93%. All of the schematics, firmware, and software are open-source.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

OPEN leaf : an open‐source cloud‐based phenotyping system for tracking dynamic changes at leaf‐specific resolution in Arabidopsis

SUMMARY The first draft of the Arabidopsis genome was released more than 20 years ago and despite intensive molecular research, more than 30% of Arabidopsis genes remained uncharacterized or without an assigned function. This is in part due to gene redundancy within gene families or the essential nature of genes, where their deletion results in lethality (i.e., thedark genome). High‐throughput plant phenotyping (HTPP) offers an automated and unbiased approach to characterize subtle or transient phenotypes resulting from gene redundancy or inducible gene silencing; however, access to commercial HTPP platforms remains limited. Here we describe the design and implementation ofOPEN leaf, an open‐source phenotyping system with cloud connectivity and remote bilateral communication to facilitate data collection, sharing and processing.OPEN leaf, coupled with our SMART imaging processing pipeline was able to consistently document and quantify dynamic changes at the whole rosette level and leaf‐specific resolution when plants experienced changes in nutrient availability. Our data also demonstrate that VIS sensors remain underutilized and can be used in high‐throughput screens to identify and characterize previously unidentified phenotypes in a leaf‐specific time‐dependent manner. Moreover, the modular and open‐source design ofOPEN leafallows seamless integration of additional sensors based on users and experimental needs.

Plant Sciences↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.

ATS↗
DOE OSTI2024

Contactless Position Measurement System for Remote Alignment of Highly Reflective Objects

The Contactless Position Measurement System (CPMS) permits automatic alignment of highly-reflecting components that cannot or should not be touched by human hands, such as those used in SRF cavities. CPMS also has application during maintenance of fusion reactors, when the components in the reactor are radioactive. The length of time for which humans can handle the components in order to carry out maintenance or repairs is limited. We would like to measure the relative positions of such objects so we can bring them into alignment with motorized stages, thus eliminating the need for human contact. The CPMS is a computer vision system that measures the location of highly reflective components to within a fraction of a millimeter by finding the edges of silhouette images. This is done by placing a uniform, infrared backlight behind the components, and viewing each component with two low-aberration, infrared cameras. Each of these stereoscopic cameras are located within a coordinate system we set up with reference platforms, light sources, and survey cameras distributed around the perimeter of the string assembly room, or fusion reactor. We obtain stereoscopic silhouette images of each component, and we use these to determine the position of the component within our string assembly coordinate system. The CPMS combines the silhouette images with knowledge of the dimensions of the objects to obtain the relative positions of components in its field of view. These position measurements then allow us to mechanically maneuver the components into contact. Once we know where they are, we can move the components to where they are supposed to be with motorized stages or other mechanical means, check they are in the right place, and bolt them together. In Phase I, we built a prototype infrared backlights and cameras, tested it, wrote an analysis program to fit the images. At a range of 50 cm, we are able to measure the relative positions of two stainless steel flanges with an accuracy of 150 μm rms. The next steps are to test the analysis program with motorized stages, build a full-size prototype, and implement it in an cleanroom with SRF cavity assembly. This can be done in a cost effective way given the materials with which the system is designed. As Open Source Instruments achieved more in Phase I than anticipated, despite not receiving a Phase II SBIR grant, we hope to move forward with product development and sales with a testing partner.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Upper-air soundings collected during the CROCUS Urban Canyons 2024 campaign in Chicago, Illinois USA

Funded by the Department of Energy’s Office of Science, Biological and Environmental Research program, Community Research on Climate and Urban Science (CROCUS) studies urban climate change and the impact it has on communities, with particular focus on disinvested, under-resourced communities. This information leads to new insights on urban climate challenges and informs future actions for mitigating and adapting to climate change at the street, neighborhood and regional levels.As part of the CROCUS effort, the Urban Canyons 2024 project was undertaken to study conditions at unprecedented detail over various neighborhoods in Chicago, Illinois. This dataset consists of upper air soundings that were collected as part of this effort. Soundings were launched during two intensive observing periods, IOP1 occurred on 22-23 July 2024, while IOP2 occurred on 27-28 July 2024. For IOP1, soundings were launched at coordinated times from three sites, Shedd Aquarium in Downtown Chicago, Abizu Campus High School in Humboldt Park, and Gary Comer Youth Center in West Woodlawn. For IOP2, the Gary Comer site was replaced by a neighborhood site in West Woodlawn, Chicago. The Abizu Campos site was operated by Valparaiso University and used iMET-4 rawinsondes, the other sites were operated by the University of Illinois Urbana-Champaign and used GRAW DFM-19 sondes.This dataset contains netCDF files containing quality-controlled temperature, dewpoint, geopotential height, pressure, and vector wind measurements at 1 second intervals following launch. These files are readable by the open-source netCDF software libraries available in many software packages (i.e., python, R, fortran, C++, etc.). The dataset also contains quicklook plots of each launch on a skew-T log-p thermodynamic diagram. These are in png format viewable by most web browsers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Tuscaloosa Marine Shale Laboratory

The Tuscaloosa Marine Shale (TMS) in Louisiana and Mississippi is an Upper Cretaceous source rock formation sandwiched between the sands of the upper and lower Tuscaloosa sections. The TMS is believed to be the source rock for underlying prolific Tuscaloosa sand formation. The TMS has an unproven estimate of 7,000,000,000 bbls of recoverable oil while its current total average production is about 3,000 bbls of oil per day in 2017. In 2013 and 2014, more than 80 wells were drilled horizontally into the TMS that were fractured using multi-stage fracturing technology. The results from this have been mixed, but recent production for several wells show an appealing initial oil production rate of more than 1000 bbl/day. The preliminary core analysis by industry partners and a few literature studies shows that the TMS is one of the most clay-rich and sensitive shales to water. Due to these and other technical problems, there is high risk for the economic development of TMS compared to other shale plays. The experiences of major industrial players in the TMS show the necessity of open and collaborative efforts to better understand the critical gaps in the development of this challenging and potentially highly economic shale play to enable more cost-efficient and environmentally-sound recovery from this unconventional liquid-rich shale play. The overall objective of this project is to form a consortium of science and industry partners to address the following six major objectives using scientific and technical approaches: 1. To improve wellbore integrity by better understanding the sources of the wellbore instability issues, proposing innovative mud and cement design for the TMS. 2. To improve formation evaluation using laboratory techniques for the evaluation of mineralogical composition, organic content, and produced-water chemistry as well as well log and geophysical analysis. 3. To determine the role of geologic discontinuities on fracture growth and shale creep behavior using digital image correlation technique. 4. To investigate the application of stable CO 2 foam and super-hydrophobic proppants for improved reservoir stimulation. 5. To better understand the nature of water/hydrocarbon/CO 2 flow in clay and organic-rich formation and the role of water/fluid interaction on recovery. 6. To prepare better socio-economic environment for TMS development by community engagement. Subsequently, the TMS virtual laboratory conducted testing and analysis of various properties of rock and formation fluids from the TMS, including but not limited to the following: Analyzing reports and logs to better understand the source of wellbore instability in TMS wells; Experiments to design a customized cement based on TMS requirements; Experiments to obtain the mineralogical and geochemical composition of TMS samples; Seismic analysis of TMS geophysical data to better predict total organic carbon (TOC) content and brittleness in TMS; Well log analysis to better estimate the TOC and geo-mechanical properties of TMS; Experiments on formation water to understand the chemistry of produced water; Experiments to determine the role of lamination and natural fractures on fracture propagation or rock deformation using digital image correlation technique in in-direct tensile tests, semi-circular bend test and creep tests Experiments to determine the stability and rheological properties of nanoparticle-stabilized CO 2 foam in TMS rock samples; Experiments to determine fluid dynamics in un-propped TMS fractures and the role of nano-coating of proppants on fluid dynamics in fractures with proppants; Micro-fluidics experiments to enhance the understanding of fluid dynamics in tight liquidrich pores with high clay content; Socio-economic studies to better engage communities in TMS development.

58 GEOSCIENCES↗
DOE OSTI2020

On-the-fly active learning of interpretable Bayesian force fields for atomistic rare events

Abstract Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important rare events, such as chemical reactions and diffusion. We present an adaptive Bayesian inference method for automating the training of interpretable, low-dimensional, and multi-element interatomic force fields using structures drawn on the fly from molecular dynamics simulations. Within an active learning framework, the internal uncertainty of a Gaussian process regression model is used to decide whether to accept the model prediction or to perform a first principles calculation to augment the training set of the model. The method is applied to a range of single- and multi-element systems and shown to achieve a favorable balance of accuracy and computational efficiency, while requiring a minimal amount of ab initio training data. We provide a fully open-source implementation of our method, as well as a procedure to map trained models to computationally efficient tabulated force fields.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient generation of grids and traversal graphs in compositional spaces towards exploration and path planning

Abstract Diverse disciplines across science and engineering deal with problems related to compositions, which exist in non-Euclidean simplex spaces, rendering many standard tools inaccurate or inefficient. This work explores such spaces conceptually in the context of materials discovery, quantifies their computational feasibility, and implements several essential methods specific to simplex spaces through a new high-performance open-source library . Most significantly, we derive and implement an algorithm for constructing a novel n-dimensional simplex graph data structure, containing all discretized compositions and possible neighbor-to-neighbor transitions. Critically, no distance or neighborhood calculations are performed, instead leveraging pure combinatorics and order in procedurally generated simplex grids, keeping the algorithm $${\mathcal{O}}(N)$$ O ( N ) , with minimal memory, enabling rapid construction of graphs with billions of transitions in seconds. Additionally, we demonstrate how such graph representations can be combined to homogeneously express complex path-planning problems, while facilitating efficient deployment of existing high-performance gradient descent, graph traversal, and other optimization algorithms.

Krajewski, Adam M. (ORCID:0000000222660099)↗