Publications

Abstract not provided.

This project sought to develop a fundamental understanding of the mechanisms underlying a newly observed enhanced germanium (Ge) diffusion process in silicon germanium (SiGe) semiconductor nanostructures during thermal oxidation. Using a combination of oxidationdiffusion experiments, high resolution imaging, and theoretical modeling, a model for the enhanced Ge diffusion mechanism was proposed. Additionally, a nanofabrication approach utilizing this enhanced Ge diffusion mechanism was shown to be applicable to arbitrary 3D shapes, leading to the fabrication of stacked silicon quantum dots embedded in SiGe nanopillars. A new wet etch-based method for preparing 3D nanostructures for highresolution imaging free of obscuring material or damage was also developed. These results enable a new method for the controlled and scalable fabrication of on-chip silicon nanostructures with sub-10 nm dimensions needed for next generation microelectronics, including low energy electronics, quantum computing, sensors, and integrated photonics.

Abstract not provided.

Abstract not provided.

Deep neural networks (NNs) typically outperform traditional machine learning (ML) approaches for complicated, non-linear tasks. It is expected that deep learning (DL) should offer superior performance for the important non-proliferation task of predicting explosive device configuration based upon observed optical signature, a task which human experts struggle with. However, supervised machine learning is difficult to apply in this mission space because most recorded signatures are not associated with the corresponding device description, or “truth labels.” This is challenging for NNs, which traditionally require many samples for strong performance. Semi-supervised learning (SSL), low-shot learning (LSL), and uncertainty quantification (UQ) for NNs are emerging approaches that could bridge the mission gaps of few labels and rare samples of importance. NN explainability techniques are important in gaining insight into the inferential feature importance of such a complex model. In this work, SSL, LSL, and UQ are merged into a single framework, a significant technical hurdle not previously demonstrated. Exponential Average Adversarial Training (EAAT) and Pairwise Neural Networks (PNNs) are chosen as the SSL and LSL methods of choice. Permutation feature importance (PFI) for functional data is used to provide explainability via the Variable importance Explainable Elastic Shape Analysis (VEESA) pipeline. A variety of uncertainty quantification approaches are explored: Bayesian Neural Networks (BNNs), ensemble methods, concrete dropout, and evidential deep learning. Two final approaches, one utilizing ensemble methods and one utilizing evidential learning, are constructed and compared using a well-quantified synthetic 2D dataset along with the DIRSIG Megascene.

Abstract not provided.

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia's Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia's ElectroMagnetic Plasma In Realistic Environments ( EMPIRE ) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia's HPC users.

Abstract not provided.

Microstructures and corresponding properties of compacted powders ultimately depend on the mechanical response of individual particles. In principle, computational simulations can predict the results of powder compaction processes, but the selection of appropriate models for both particle–particle interactions and particle deformations across all relevant length scales remain nontrivial tasks, especially in material systems lacking detailed mechanical property information. The work presented here addresses these issues by conducting uniaxial compressions in situ inside of a scanning electron microscope to characterize the mechanical response of individual micron-sized particles of a molecular crystal, hexanitrohexaazaisowurtzitane (CL-20). This experimental approach enabled the collection of quantitative force and displacement data alongside simultaneous imaging to capture morphology changes. The results reveal information about elastic deformation, yield, plastic deformation, creep, and fracture phenomena. Accordingly, this work demonstrates a generalizable approach for assessing the mechanical response of individual micron-sized molecular crystal particles and utilizing those responses in particle-level models. Graphic abstract: [Figure not available: see fulltext.].

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report provides a demonstration of MELCOR and MELCOR Accident Consequence Code System (MACCS) capabilities to perform a dose assessment for a Molten Salt Reactor (MSR) off-gas system. A primary containment system salt spill is used as the off-normal scenario, along with a normal operation dose assessment for comparison. This report discusses the tools, methods, and information used in this assessment so that it may be utilized as a starting point for future advanced reactor consequence analyses. This report also highlights several gaps, to include the need for reactor inventory information specific to advanced reactors, and the need for specific atmospheric transport models that take into account the unique deposition behaviors of tritium and carbon-14, and makes recommendations for closing these gaps. This report satisfies the DOE NE Milestone M4RD-21SN0601062.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This document describes the Cybersecurity Research Development and Demonstration (RD&D) Program, established by the Department of Energy Office of Nuclear Energy (NE) to provide sciencebased methods and technologies necessary for cost-effective, cyber-secure digital instrumentation, control and communication in collaboration with nuclear energy stakeholders. It provides an overview of program goals, objectives, linkages to organizational strategies, management structure, and stakeholder and cross-program interfaces.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The DOE-NE NWM Cloud was designed to be a generic set of tools and applications for any nuclear waste management program. As policymakers continue to consider approaches that emphasize consolidated interim storage and transportation of spent nuclear fuel, a gap analysis of the tools and applications provided for spent nuclear fuel and high-level radioactive waste disposal in comparison those needed for siting, licensing, and developing a consolidated interim storage facility and/or for a transportation campaign will help prepare DOE for implementing such potential policy direction. This report evaluates the points of alignment and potential gaps between the applications on the NWM Cloud that supported SNF disposal project, and the applications needed to address QA requirements and for other project support needs of an SNF storage project.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

While computer systems, software applications, and operational technology (OT)/Industrial Control System (ICS) devices are regularly updated through automated and manual processes, there are several unique challenges associated with distributed energy resource (DER) patching. Millions of DER devices from dozens of vendors have been deployed in home, corporate, and utility network environments that may or may not be internet-connected. These devices make up a growing portion of the electric power critical infrastructure system and are expected to operate for decades. During that operational period, it is anticipated that critical and noncritical firmware patches will be regularly created to improve DER functional capabilities or repair security deficiencies in the equipment. The SunSpec/Sandia DER Cybersecurity Workgroup created a Patching Subgroup to investigate appropriate recommendations for the DER patching, holding fortnightly meetings for more than nine months. The group focused on DER equipment, but the observations and recommendations contained in this report also apply to DERMS tools and other OT equipment used in the end-to-end DER communication environment. The group found there were many standards and guides that discuss firmware lifecycles, patch and asset management, and code-signing implementations, but did not singularly cover the needs of the DER industry. This report collates best practices from these standards organizations and establishes a set of best practices that may be used as a basis for future national or international patching guides or standards.

Seven generation III+ and generation IV nuclear reactor types, based on twelve reactor concepts surveyed, are examined using functional decomposition to extract relevant operational technology (OT) architecture information. This information is compared to existing nuclear power plants (NPPs) OT architectures to highlight novel and emergent cyber risks associated with next generation NPPs. These insights can help inform operational technology architecture requirements that will be unique to a given reactor type. Next generation NPPs have streamlined OT architectures relative to the current generation II commercial NPP fleet. Overall, without compensatory measures that provide sufficient and efficient cybersecurity controls, next generation NPPs will have increased cyber risk. Verification and validation of cyber-physical testbeds and cyber risk assessment methodologies may be an important next step to reduce cyber risk in the OT architecture design and testing phase. Coordination with safety requirements can result in OT architecture design being an iterative process.

This report presents the results of the “Foundations of Rigorous Cyber Experimentation” (FORCE) Laboratory Directed Research and Development (LDRD) project. This project is a companion project to the “Science and Engineering of Cyber security through Uncertainty quantification and Rigorous Experimentation” (SECURE) Grand Challenge LDRD project. This project leverages the offline, controlled nature of cyber experimentation technologies in general, and emulation testbeds in particular, to assess how uncertainties in network conditions affect uncertainties in key metrics. We conduct extensive experimentation using a Firewheel emulation-based cyber testbed model of Invisible Internet Project (I2P) networks to understand a de-anonymization attack formerly presented in the literature. Our goals in this analysis are to see if we can leverage emulation testbeds to produce reliably repeatable experimental networks at scale, identify significant parameters influencing experimental results, replicate the previous results, quantify uncertainty associated with the predictions, and apply multi-fidelity techniques to forecast results to real-world network scales. The I2P networks we study are up to three orders of magnitude larger than the networks studied in SECURE and presented additional challenges to identify significant parameters. The key contributions of this project are the application of SECURE techniques such as UQ to a scenario of interest and scaling the SECURE techniques to larger network sizes. This report describes the experimental methods and results of these studies in more detail. In addition, the process of constructing these large-scale experiments tested the limits of the Firewheel emulation-based technologies. Therefore, another contribution of this work is that it informed the Firewheel developers of scaling limitations, which were subsequently corrected.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Physical protection of public buildings has long been a concern of police and security services where a balance of facility security and personnel safety is vital. Due to the nature of public spaces, the use of permanently installed and deploy-on-demand physical barrier systems must be safe for the legitimate occupants and visitors of that space. Such systems must seek to mitigate the personal and organizational consequences of unintentionally seriously injuring or killing an innocent bystander by slamming a heavy, rigid, and quick-deploying barrier into place. Consideration and implementation of less-than-lethal technologies is necessary to reduce risk to visitors and building personnel. One potential barrier solution is a fast-acting, high-strength, composite airbag barrier system for doorways and hallways to quickly deploy a less-than-lethal barrier at entry points as well as isolate intruders who have already gained access. This system is envisioned to be stored within an architecturally attractive selectively frangible shell that could be permanently installed at a facility or installed in remote or temporary locations as dictated by risk. The system would be designed to be activated remotely (hardwired or wireless) from a Central Alarm Station (CAS) or other secure location.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We present our research findings on the novel NDN protocol. In this work, we defined key attack scenarios for possible exploitation and detail software security testing procedures to evaluate the security of the NDN software. This work was done in the context of distributed energy resources (DER). The software security testing included an execution of unit tests and static code analyses to better understand the software rigor and the security that has been implemented. The results from the penetration testing are presented. Recommendations are discussed to provide additional defense for secure end-to-end NDN communications.

Abstract not provided.

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.2 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Cybersecurity for industrial control systems is an important consideration that advance reactor designers will need to consider. How cyber risk is managed is the subject of on-going research and debate in the nuclear industry. This report seeks to identify potential cyber risks for advance reactors. Identified risks are divided into absorbed risk and licensee managed risk to clearly show how cyber risks for advance reactors can potentially be transferred. Absorbed risks are risks that originate external to the licensee but may unknowingly propagate into the plant. Insights include (1) the need for unification of safety, physical security, and cybersecurity risk assessment frameworks to ensure optimal coordination of risk, (2) a quantitative risk assessment methodology in conjunction with qualitative assessments may be useful in efficiently and sufficiently managing cyber risks, and (3) cyber risk management techniques should align with a risked informed regulatory framework for advance reactors.

Highlights: Battery energy storage may improve energy efficiency and reliability of hybrid energy systems composed by diesel and solar photovoltaic power generators serving isolated communities.In projects aiming update of power plants serving electrically isolated communities with redundant diesel generation, battery energy storage can improve overall economic performance of power supply system by reducing fuel usage, decreasing capital costs by replacing redundant diesel generation units, and increasing generator system life by shortening yearly runtime.Fast-acting battery energy storage systems with grid-forming inverters might have potential for improving drastically the reliability indices of isolated communities currently supplied by diesel generation. Abstract: This paper will highlight unique challenges and opportunities with regard to energy storage utilization in remote, self-sustaining communities. The energy management of such areas has unique concerns. Diesel generation is often the go-to power source in these scenarios, but these systems are not devoid of issues. Without dedicated maintenance crews as in large, interconnected network areas, minor interruptions can be frequent and invasive not only for those who lose power, but also for those in the community that must then correct any faults. Although the immediate financial benefits are perhaps not readily apparent, energy storage could be used to address concerns related to reliability, automation, fuel supply concerns, generator degradation, solar utilization, and, yes, fuel costs to name a few. These ideas are shown through a case study of the Levelock Village of Alaska. Currently, the community is faced with high diesel prices and a difficult supply chain, which makes temporary loss of power very common and reductions in fuel consumption very impactful. This study will investigate the benefits that an energy storage system could bring to the overall system life, fuel costs, and reliability of the power supply. The variable efficiency of the generators, impact of startup/shutdown process, and low-load operation concerns are considered. The technological benefits of the combined system will be explored for various scenarios of future diesel prices and technology maintenance/replacement costs as well as for the avoidance of power interruptions that are so common in the community currently. Graphic abstract: [Figure not available: see fulltext.] Discussion: In several cases, energy storage can provide a means to promote energy equity by improving remote communities’ power supply reliability to levels closer to what the average urban consumer experiences at a reduced cost compared to transmission buildout. Furthermore, energy equity represents a hard-to-quantify benefit achieved by the integration of energy storage to isolated power systems of under-served communities, which suggests that the financial aspects of such projects should be questioned as the main performance criterion. To improve battery energy storage system valuation for diesel-based power systems, integration analysis must be holistic and go beyond fuel savings to capture every value stream possible.

Abstract not provided.

This paper presents a practical methodology for propagating and processing uncertainties associated with random measurement and estimation errors (that vary from test-to-test) and systematic measurement and estimation errors (uncertain but similar from test-to-test) in inputs and outputs of replicate tests to characterize response variability of stochastically varying test units. Also treated are test condition control variability from test-to-test and sampling uncertainty due to limited numbers of replicate tests. These aleatory variabilities and epistemic uncertainties result in uncertainty on computed statistics of output response quantities. The methodology was developed in the context of processing experimental data for “real-space” (RS) model validation comparisons against model-predicted statistics and uncertainty thereof. The methodology is flexible and sufficient for many types of experimental and data uncertainty, offering the most extensive data uncertainty quantification (UQ) treatment of any model validation method the authors are aware of. It handles both interval and probabilistic uncertainty descriptions and can be performed with relatively little computational cost through use of simple and effective dimension- and order-adaptive polynomial response surfaces in a Monte Carlo (MC) uncertainty propagation approach. A key feature of the progressively upgraded response surfaces is that they enable estimation of propagation error contributed by the surrogate model. Sensitivity analysis of the relative contributions of the various uncertainty sources to the total uncertainty of statistical estimates is also presented. The methodologies are demonstrated on real experimental validation data involving all the mentioned sources and types of error and uncertainty in five replicate tests of pressure vessels heated and pressurized to failure. Simple spreadsheet procedures are used for all processing operations.

Abstract not provided.

State-of-the-art engineering and science codes have grown in complexity dramatically over the last two decades. Application teams have adopted more sophisticated development strategies, leveraging third party libraries, deploying comprehensive testing, and using advanced debugging and profiling tools. In today's environment of diverse hardware platforms, these applications also desire performance portability-avoiding the need to duplicate work for various platforms. The Kokkos EcoSystem provides that portable software stack. Based on the Kokkos Core Programming Model, the EcoSystem provides math libraries, interoperability capabilities with Python and Fortran, and Tools for analyzing, debugging, and optimizing applications. In this article, we overview the components, discuss some specific use cases, and highlight how codesigning these components enables a more developer friendly experience.

Abstract not provided.

To date, disinformation research has focused largely on the production of false information ignoring the suppression of select information. We term this alternative form of disinformation information suppression. Information suppression occurs when facts are withheld with the intent to mislead. In order to detect information suppression, we focus on understanding the actors who withhold information. In this research, we use knowledge of human behavior to find signatures of different gatekeeping behaviors found in text. Specifically, we build a model to classify the different types of edits on Wikipedia using the added text alone and compare a human-informed feature engineering approach to a featureless algorithm. Being able to computationally distinguish gatekeeping behaviors is a first step towards identifying when information suppression is occurring.

Abstract not provided.

Graph partitioning has been an important tool to partition the work among several processors to minimize the communication cost and balance the workload. While accelerator-based supercomputers are emerging to be the standard, the use of graph partitioning becomes even more important as applications are rapidly moving to these architectures. However, there is no distributed-memory-parallel, multi-GPU graph partitioner available for applications. We developed a spectral graph partitioner, Sphynx, using the portable, accelerator-friendly stack of the Trilinos framework. In Sphynx, we allow using different preconditioners and exploit their unique advantages. We use Sphynx to systematically evaluate the various algorithmic choices in spectral partitioning with a focus on the GPU performance. We perform those evaluations on two distinct classes of graphs: regular (such as meshes, matrices from finite element methods) and irregular (such as social networks and web graphs), and show that different settings and preconditioners are needed for these graph classes. The experimental results on the Summit supercomputer show that Sphynx is the fastest alternative on irregular graphs in an application-friendly setting and obtains a partitioning quality close to ParMETIS on regular graphs. When compared to nvGRAPH on a single GPU, Sphynx is faster and obtains better balance and better quality partitions. Sphynx provides a good and robust partitioning method across a wide range of graphs for applications looking for a GPU-based partitioner.

Zero-emissions hydrogen fuel cell electrical vehicles (FCEVs) have become more popular in recent years. However, the limited availability of hydrogen fueling stations is considered a critical barrier to sustainable adoption of hydrogen FCEV. To enable the widespread deployment and commercialization of hydrogen FCEV, the availability of hydrogen refueling stations needs to improve. One of the consequences of the lack of hydrogen refueling infrastructure is that consumers can suffer from “range anxiety”, meaning consumers would get anxious of running out of fuel during long-distance trip [4]. A practical solution is to provide a compact emergency hydrogen refueler that can be used if the consumer runs out of hydrogen before reaching the nearest hydrogen refueling station. A safe, compact, and user-friendly hydrogen refueler would give consumers the flexibility they need to feel comfortable using their hydrogen FCEV when planning a long-distance trip. Offering this product would alleviate range anxiety, and it would make Hydrogen FCEV a more attractive alternative to gasoline vehicles. The emergency hydrogen refueler consists of a lithium hydride bed that reacts with liquid water to produce hydrogen gas and lithium hydroxide.

Abstract not provided.

The predictive understanding of catalytic surface reactions requires accurate microkinetic models, and while decades of work has been devoted to the elucidation of the reaction steps in these models, many open questions remain. One key issue is a lack of approaches enabling the local spatially resolved assessment of catalytic activity over a surface. In this report, we detail efforts to develop a new diagnostic approach to solve this problem. The approach is based upon laser resonance enhanced multiphoton ionization of reaction products emitted into the gas phase followed by spatially resolved imaging of the resultant ions or electrons. Ion imaging is pursued with a velocity-selected spatially resolved ion imaging microscope, while electron imaging was attempted in a low energy electron microscope. Successful demonstration of the ion imaging microscope coupled with the development of transport simulations shows promise for a revolutionary new tool to assess local catalytic activity

Laser powder bed fusion (LPBF) Additive manufacturing (AM) has attracted interest as an agile method of building production metal parts to reduce design-build-test cycle times for systems. However, predicting part performance is difficult due to inherent process variabilities. This makes qualification challenging. Computational process models have attempted to address some of these challenges, including mesoscale, full physics models and reduced fidelity conduction models. The goal of this work is credible multi-fidelity modeling of the LPBF process by investigating methods for estimating the error between models of two different fidelities. Two methods of error estimation are investigated, adjoint-based error estimation and Bayesian calibration. Adjoint-based error estimation is found to effectively bounding the error between the two models, but with very conservative bounds, making predictions highly uncertain. Bayesian parameter calibration applied to conduction model heat source parameters is found to effectively bound the observed error between the models for melt pool morphology quantities of interest. However, the calibrations do not effectively bound the error in heat distribution.

Abstract not provided.

Abstract not provided.

Arithmetic Coding (AC) using Prediction by Partial Matching (PPM) is a compression algorithm that can be used as a machine learning algorithm. This paper describes a new algorithm, NGram PPM. NGram PPM has all the predictive power of AC/PPM, but at a fraction of the computational cost. Unlike compression-based analytics, it is also amenable to a vector space interpretation, which creates the ability for integration with other traditional machine learning algorithms. AC/PPM is reviewed, including its application to machine learning. Then NGram PPM is described and test results are presented, comparing them to AC/PPM.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report is a preliminary test plan of the seismic shake table test. The final report will be developed when all decisions regarding the test hardware, instrumentation, and shake table inputs are made. A new revision of this report will be issued in spring of 2022. The preliminary test plan documents the free-field ground motions that will be used as inputs to the shake table, the test hardware, and instrumentation. It also describes the facility at which the test will take place in late summer of 2022.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Metallic lattice structures are being considered for shock mitigation applications due to their superior mechanical properties, energy absorption capability and lightweight characteristics inherent of the additive manufacturing process. In this study, shock compression experiments coupled to x-ray phase contrast imaging (PCI) were conducted on 316L stainless steel lattices. Meso-scale simulations incorporating the as-built lattice structure characterized by computed tomography were used to simulate PCI radiographs in CTH for direct comparison to experimental data. The methodology presented here offers robust validation for constitutive properties to further our understanding of lattice compaction at application-relevant strain rates.

Abstract not provided.

This report summarizes the findings and outcomes of the LDRD-express project with title “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”. The primary motivation of this project was the computational/mathematical exploration of the ideas advanced aiming to improve the state-of-the-art on numerical methods for the one-fluid Euler-Poisson models and gain some understanding on the Euler-Maxwell model. Euler-Poisson and Euler-Maxwell, by themselves are not the most technically relevant PDE plasma-models. However, both of them are elementary building blocks of PDE-models used in actual technical applications and include most (if not all) of their mathematical difficulties. Outside the classical ideal MHD models, rigorous mathematical and numerical understanding of one-fluid models is still a quite undeveloped research area, and the treatment/understanding of boundary conditions is minimal (borderline non-existent) at this point in time. This report focuses primarily on bulk-behaviour of Euler-Poisson’s model, touching boundary conditions only tangentially.

Supplementing an existing high-quality seismic monitoring network with openly available station data could improve coverage and decrease magnitudes of completeness; however, this can present challenges when varying levels of data quality exist. Without discerning the quality of openly available data, using it poses significant data management, analysis, and interpretation issues. Incorporating additional stations without properly identifying and mitigating data quality problems can degrade overall monitoring capability. If openly available stations are to be used routinely, a robust, automated data quality assessment for a wide range of quality control (QC) issues is essential. To meet this need, we developed Pycheron, a Python-based library for QC of seismic waveform data. Pycheron was initially based on the Incorporated Research Institutions for Seismology's Modular Utility for STAtistical kNowledge Gathering but has been expanded to include more functionality. Pycheron can be implemented at the beginning of a data processing pipeline or can process stand-alone data sets. Its objectives are to (1) identify specific QC issues; (2) automatically assess data quality and instrumentation health; (3) serve as a basic service that all data processing builds on by alerting downstream processing algorithms to any quality degradation; and (4) improve our ability to process orders of magnitudes more data through performance optimizations. This article provides an overview of Pycheron, its features, basic workflow, and an example application using a synthetic QC data set.

Abstract not provided.

Abstract not provided.

A series of nanopillar compression tests were performed on tungsten as a function of temperature using in situ transmission electron microscopy with localized laser heating. Surface oxidation was observed to form on the pillars and grow in thickness with increasing temperature. Deformation between 850◦C and 1120◦C is facilitated by long-range diffusional transport from the tungsten pillar onto adjacent regions of the Y2O3-stabilized ZrO2 indenter. The constraint imposed by the surface oxidation is hypothesized to underly this mechanism for localized plasticity, which is generally the so-called whisker growth mechanism. The results are discussed in context of the tungsten fuzz growth mechanism in He plasma-facing environments. The two processes exhibit similar morphological features and the conditions under which fuzz evolves appear to satisfy the conditions necessary to induce whisker growth.

With the rapid proliferation of additive manufacturing and 3D printing technologies, architected cellular solids including truss-like 3D lattice topologies offer the opportunity to program the effective material response through topological design at the mesoscale. The present report summarizes several of the key findings from a 3-year Laboratory Directed Research and Development Program. The program set out to explore novel lattice topologies that can be designed to control, redirect, or dissipate energy from one or multiple insult environments relevant to Sandia missions, including crush, shock/impact, vibration, thermal, etc. In the first 4 sections, we document four novel lattice topologies stemming from this study: coulombic lattices, multi-morphology lattices, interpenetrating lattices, and pore-modified gyroid cellular solids, each with unique properties that had not been achieved by existing cellular/lattice metamaterials. The fifth section explores how unintentional lattice imperfections stemming from the manufacturing process, primarily sur face roughness in the case of laser powder bed fusion, serve to cause stochastic response but that in some cases such as elastic response the stochastic behavior is homogenized through the adoption of lattices. In the sixth section we explore a novel neural network screening process that allows such stocastic variability to be predicted. In the last three sections, we explore considerations of computational design of lattices. Specifically, in section 7 using a novel generative optimization scheme to design novel pareto-optimal lattices for multi-objective environments. In section 8, we use computational design to optimize a metallic lattice structure to absorb impact energy for a 1000 ft/s impact. And in section 9, we develop a modified micromorphic continuum model to solve wave propagation problems in lattices efficiently.

The properties of materials can change dramatically at the nanoscale new and useful properties can emerge. An example is found in the paramagnetism in iron oxide magnetic nanoparticles. Using magnetically sensitive nitrogen-vacancy centers in diamond, we developed a platform to study electron spin resonance of nanoscale materials. To implement the platform, diamond substrates were prepared with nitrogen vacancy centers near the surface. Nanoparticles were placed on the surface using a drop casting technique. Using optical and microwave pulsing techniques, we demonstrated T1 relaxometry and double electron-electron resonance techniques for measuring the local electron spin resonance. The diamond NV platform developed in this project provides a combination of good magnetic field sensitivity and high spatial resolution and will be used for future investigations of nanomaterials and quantum materials.

Abstract not provided.

Abstract not provided.

In 2006, the United States Congress mandated for the Department of Homeland Security (DHS) to screen all cargo containers to protect against terrorist acts and ensure the safety and security of the nation. Containers are screened by gamma and neutron detectors to ensure threat material is not smuggled into the country. However, because commerce is radioactive, detecting the presence of radioactive materials is not sufficient in ensuring the safety and security of the nation. Radioactive materials must also be identified in real time, thus distinguishing threat sources (Uranium-235 and Plutonium-230) from non-threat sources (kitty litter, pot ash, medical isotopes). Screening cargo containers can be considered a two-step process (1) alerting to the presence of radioactive material when gamma counts exceed a threshold setting, and (2) once alerted, identifying the type of radioactive material, which is done by collecting a gamma spectra and analyzing it with an analysis tool/algorithm. For this reason, it is important to evaluate not only emerging technology in neutron and gamma detection, but also investigate new advances in algorithm development for radioisotope identification (RIID). New candidates in detection and on-board algorithm analysis might offer opportunities to make the scanning, detection, and identification process more efficient while still ensuring the health and safety of the public. This research will investigate emerging technology in radiation detection focused on gamma spectroscopy capabilities and RIID algorithms for DHS applications.

Abstract not provided.

This project developed prototype germanium telluride switches, which can be used in RF applications to improve SWAP (size, weight, and power) and signal quality in RF systems. These switches can allow for highly reconfigurable systems, including antennas, communications, optical systems, phased arrays, and synthetic aperture radar, which all have high impact on current National Security goals for improved communication systems and communication technology supremacy. The final result of the project was the demonstration of germanium telluride RF switches, which could act as critical elements necessary for a single chip RF communication system that will demonstrate low SWAP and high reconfigurability

Abstract not provided.

Abstract not provided.

Abstract not provided.

All energy production systems need efficient energy conversion systems. Current Rankine cycles use water to generate steam at temperatures where efficiency is limited to around 40%. As existing fossil and nuclear power plants are decommissioned due to end of effective life and/or societies’ desire for cleaner generation options, more efficient energy conversion is needed to keep up with increasing electricity demands. Modern energy generation technologies, such as advanced nuclear reactors and concentrated solar, coupled to high efficiency sCO₂ conversion systems provide a solution to efficient, clean energy systems. Leading R&D communities worldwide agree that the successful development of sCO₂ Brayton power cycle technology will eventually bring about large-scale changes to existing multi-billion-dollar global markets and enable power applications not currently possible or economically justifiable. However, all new technologies face challenges in the path to commercialization and the electricity sector is distinctively risk averse. The Sandia sCO₂ Brayton team needs to better understand what the electricity sector needs in terms of new technology risk mitigation, generation efficiency, reliability improvements above current technology, and cost requirements which would make new technology adoption worthwhile. Relying on the R&D community consensus that a sCO₂ power cycle will increase the revenue of the electrical industry, without addressing the electrical industry’s concerns, significantly decreases the potential for adoption at commercial scale. With a clear understanding of the market perspectives on technology adoption, including military, private sector, and utilities customers, the Sandia Brayton Team can resolve industry concerns for smoother development and faster transition to commercialization. An extensive customer discovery process, similar to that executed through the NSF’s I-Corp program, is necessary in order to understand the pain points of the market and articulate the value proposition of Brayton systems in terms that engage decision makers and facilitate commercialization of the technology.

Abstract not provided.

Abstract not provided.

Identification and characterization of underground events from surface or remote data requires a thorough understanding of the rock material properties. However, material properties usually come from borehole data, which is expensive and not always available. A potential alternative is to use topographic characteristics to approximate the strength, but this has never been done before quantitatively. Here we present the results from the first steps towards this goal. We have found that there are strong correlations between compressive and tensile strengths and slopes, but these correlations vary depending on data analysis details. Rugosity may be better correlated to strength than slope values. More comprehensive analyses are needed to fully understand the best method of predicting strength from topography for this area. We also found that misalignment of multiple GIS datasets can have a large influence on the ability to make interpretations. Lastly, these results will require further study in a variety of climatic conditions before being applicable to other sites.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This project focused on providing a fundamental physico-chemical understanding of the coupling mechanisms of corrosion- and radiation-induced degradation at material-salt interfaces in Ni-based alloys operating in emulated Molten Salt Reactor(MSR) environments through the use of a unique suite of aging experiments, in-situ nanoscale characterization experiments on these materials, and multi-physics computational models. The technical basis and capabilities described in this report bring us a step closer to accelerate the deployment of MSRs by closing knowledge gaps related to materials degradation in harsh environments.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report details the results of a three-fold investigation of sensitivity analysis (SA) for machine learning (ML) explainability (MLE): (1) the mathematical assessment of the fidelity of an explanation with respect to a learned ML model, (2) quantifying the trustworthiness of a prediction, and (3) the impact of MLE on the efficiency of end-users through multiple users studies. We focused on the cybersecurity domain as the data is inherently non-intuitive. As ML is being using in an increasing number of domains, including domains where being wrong can elicit high consequences, MLE has been proposed as a means of generating trust in a learned ML models by end users. However, little analysis has been performed to determine if the explanations accurately represent the target model and they themselves should be trusted beyond subjective inspection. Current state-of-the-art MLE techniques only provide a list of important features based on heuristic measures and/or make certain assumptions about the data and the model which are not representative of the real-world data and models. Further, most are designed without considering the usefulness by an end-user in a broader context. To address these issues, we present a notion of explanation fidelity based on Shapley values from cooperative game theory. We find that all of the investigated MLE explainability methods produce explanations that are incongruent with the ML model that is being explained. This is because they make critical assumptions about feature independence and linear feature interactions for computational reasons. We also find that in deployed, explanations are rarely used due to a variety of reason including that there are several other tools which are trusted more than the explanations and there is little incentive to use the explanations. In the cases when the explanations are used, we found that there is the danger that explanations persuade the end users to wrongly accept false positives and false negatives. However, ML model developers and maintainers find the explanations more useful to help ensure that the ML model does not have obvious biases. In light of these findings, we suggest a number of future directions including developing MLE methods that directly model non-linear model interactions and including design principles that take into account the usefulness of explanations to the end user. We also augment explanations with a set of trustworthiness measures that measure geometric aspects of the data to determine if the model output should be trusted.

Recent work has shown that artificial opsonins stimulate the targeted destruction of bacteria by phagocyte immune cells. Artificial opsonization has the potential to direct the innate immune system to target novel antigens, potentially even viral pathogens. Furthermore, the engagement of innate immunity presents a potential solution for the spread of pandemics in a scenario when a vaccine is unavailable or ineffective. Funded by the LDRD late start bioscience pandemic response program, we tested whether artificial opsonins can be developed to target viral pathogens using phage MS2 and a SARS-CoV-2 surrogate. To direct opsonization against these viruses we purified antibody derived viral targeting motifs and attempted the same chemical conjugation strategies that produced bacterial targeting artificial opsonins. However, the viral targeting motifs proved challenging to conjugate using these methods, frequently resulting in precipitation and loss of product. Future studies may be successful with this approach if a smaller and more soluble viral-targeting peptide could be used.

Abstract not provided.

We present an experimental and numerical study of a terahertz metamaterial with a nonlinear response that is controllable via the relative structural arrangement of two stacked split ring resonator arrays. The first array is fabricated on an n-doped GaAs substrate, and the second array is fabricated vertically above the first using a polyimide spacer layer. Due to GaAs carrier dynamics, the on-resonance terahertz transmission at 0.4 THz varies in a nonlinear manner with incident terahertz power. The second resonator layer dampens this nonlinear response. In samples where the two layers are aligned, the resonance disappears, and the total nonlinear modulation of the on-resonance transmission decreases. The nonlinear modulation is restored in samples where an alignment offset is imposed between the two resonator arrays. Structurally tunable metamaterials and metasurfaces can therefore act as a design template for tunable nonlinear THz devices by controlling the coupling of confined electric fields to nonlinear phenomena in a complex material substrate or inclusion.

This SAND report documents Exploratory Express LDRD Project 223790, "Structure-preserving numerical discretizations for domains with boundaries", which developed a method to incorporate consistent treatment of domain boundaries and arbitrary boundary conditions in discrete exterior calculus (DEC) for arbitrary polygonal (2D) and tensor-product structure prism (3D) grids. The new DEC required the development of novel discrete exterior derivatives, boundary operators, wedge products and Hodge stars. This was accomplished through the use of boundary extension and the blending of known 2D operators on the interior with 1D operators on the boundary. The Hodge star was based on the Voronoi Hodge star, and retained the limitation of a triangular circumcentric primal or dual grid along with low-order accuracy. In addition to the new DEC, two related software packages were written: one for the study of DEC operators on arbitrary polygonal and polyhedral grids using both symbolic and numerical approaches and one for a (thermal) shallow water testbed using TRiSK-type numerics. Immediately relevant (already funded, through CANGA) followup work is the development of a high-order, geometrically flexible Hodge star and structure-preserving, high-order, oscillation-limiting transport operators (using WENO) for n-forms on arbitrary 2D and 3D grids. This will provide all of the machinery required for a high-order version of TRiSK with boundaries on arbitrary 2D and tensor-product 3D grids, which is applicable to both the atmospheric (CRM in E3SM-MMF) and oceanic (MPAS-O) components of E3SM.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic ℎ-adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.

Abstract not provided.

Denoising contaminated seismic signals for later processing is a fundamental problem in seismic signals analysis. The most straightforward denoising approach, using spectral filtering, is not effective when noise and seismic signal occupy the same frequency range. Neural network approaches have shown success denoising local signal when trained on short-time Fourier transform spectrograms (Zhu et al 2018; Tibi et al 2021). Scalograms, a wavelet-based transform, achieved ~15% better reconstruction as measured by dynamic time warping on a seismic waveform test set than spectrograms, suggesting their use as an alternative for denoising. We train a deep neural network on a scalogram dataset derived from waveforms recorded by the University of Utah Seismograph Stations network. We find that initial results are no better than a spectrogram approach, with additional overhead imposed by the significantly larger size of scalograms. A robust exploration of neural network hyperparameters and network architecture was not performed, which could be done in follow on work.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

There is a dearth in the literature on how to capture the uncertainty generated by material surface evolution in thermal modeling. This leads to inadequate or highly variable uncertainty representations for material properties, specifically emissivity when minimal information is available. Inaccurate understandings of prediction uncertainties may lead decision makers to incorrect conclusions, so best engineering practices should be developed for this domain. In order to mitigate the aforementioned issues, this study explores different strategies to better capture the thermal uncertainty response of engineered systems exposed to fire environments via defensible emissivity uncertainty characterizations that can be easily adapted to a variety of use cases. Two unique formulations (one physics-informed and one mathematically based) are presented. The formulations and methodologies presented herein are not exhaustive but more so are a starting point and give the reader a basis for how to customize their uncertainty definitions for differing fire scenarios and materials. Finally, the impact of using this approach versus other commonly used strategies and the usefulness of adding rigor to material surface evolution uncertainty is demonstrated.

Abstract not provided.

Abstract not provided.

Understanding the fundamental mechanisms underpinning shock initiation is critical to predicting energetic material (EM) safety and performance. Currently, the timescales and pathways by which shock-excited lattice modes transfer energy into specific chemical bonds remains an open question. Towards understanding these mechanisms, our group has previously measured the vibrational energy transfer (VET) pathways in several energetic thin films using broadband, femtosecond transient absorption spectroscopy. However, new technologies are needed to move beyond these thin film surrogates and measure broadband VET pathways in realistic EM morphologies. Herein, we describe a new broadband, femtosecond, attenuated total reflectance spectroscopy apparatus. Performance of the system is benchmarked against published data and the first VET results from a pressed EM pellet are presented. This technology enables fundamental studies of VET dynamics across sample configurations and environments (pressure, temperature, etc .) and supports the potential use of VET studies in the non-destructive surveillance of EM components.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report summarizes the 2021 fiscal year (FY21) status of ongoing borehole heater tests in salt funded by the disposal research and development (R&D) program of the Office of Spent Fuel & Waste Science and Technology (SFWST) of the US Department of Energy’s Office of Nuclear Energy’s (DOE-NE) Office of Spent Fuel and Waste Disposition (SFWD). This report satisfies SFWST milestone M2SF- 21SN010303052 by summarizing test activities and data collected during FY21. The Brine Availability Test in Salt (BATS) is fielded in a pair of similar arrays of horizontal boreholes in an experimental area at the Waste Isolation Pilot Plant (WIPP). One array is heated, the other unheated. Each array consists of 14 boreholes, including a central borehole with gas circulation to measure water production, a cement seal exposure test, thermocouples to measure temperature, electrodes to infer resistivity, a packer-isolated borehole to add tracers, fiber optics to measure temperature and strain, and piezoelectric transducers to measure acoustic emissions. The key new data collected during FY21 include a series of gas tracer tests (BATS phase 1b), a pair of liquid tracer tests (BATS phase 1c), and data collected under ambient conditions (including a period with limited access due to the ongoing pandemic) since BATS phase 1a in 2020. A comparison of heated and unheated gas tracer test results clearly shows a decrease in permeability of the salt upon heating (i.e., thermal expansion closes fractures, which reduces permeability).

University partnerships play an essential role in sustaining Sandia’s vitality as a national laboratory. The SAA is an element of Sandia’s broader University Partnerships program, which facilitates recruiting and research collaborations with dozens of universities annually. The SAA program has two three-year goals. SAA aims to realize a step increase in hiring results, by growing the total annual inexperienced hires from each out-of-state SAA university. SAA also strives to establish and sustain strategic research partnerships by establishing several federally sponsored collaborations and multi-institutional consortiums in science & technology (S&T) priorities such as autonomy, advanced computing, hypersonics, quantum information science, and data science. The SAA program facilitates access to talent, ideas, and Research & Development facilities through strong university partnerships. Earlier this year, the SAA program and campus executives hosted John Myers, Sandia’s former Senior Director of Human Resources (HR) and Communications, and senior-level staff at Georgia Tech, U of Illinois, Purdue, UNM, and UT Austin. These campus visits provided an opportunity to share the history of the partnerships from the university leadership, tours of research facilities, and discussions of ongoing technical work and potential recruiting opportunities. These visits also provided valuable feedback to HR management that will help Sandia realize a step increase in hiring from SAA schools. The 2020-2021 Collaboration Report is a compilation of accomplishments in 2020 and 2021 from SAA and Sandia’s valued SAA university partners.

Develop, verify, and document model capabilities sufficient for comparing field wake measurements from SWiFT with synthetic lidar wake measurements from Nalu-Wind (hereafter referred to as `Nalu').

Layered van der Waals heterostructures provide extraordinary opportunities for applications such as thermoelectrics and allow for tunability of optical and electronic properties. The performance of devices made from these heterostructures will depend on their properties, which are sensitive to the nanoarchitecture (constituent layer thicknesses, layer sequence, etc.). However, performance will also be impacted by defects, which will vary in concentration and identity with the nanoarchitecture and preparation conditions. Here, we identify several types of defects and propose mechanisms for their formation, focusing on compounds in the ([SnSe]1+δ)m(TiSe2)n system prepared using the modulated elemental reactants method. The defects were observed by atomic resolution high-angle annular dark-field scanning transmission electron microscopy and can be broadly categorized into those that form domain boundaries as a result of rotational disorder from the self-assembly process and those that are layer-thickness-related and result from local or global deviations in the amount of material deposited. Defect type and density were found to depend on the nanoarchitecture of the heterostructure. Categorizing the defects provides insights into defect formation in these van der Waals layered heterostructures and suggests strategies for controlling their concentrations. Strategies for controlling defect type and concentration are proposed, which would have implications for transport properties for applications in thermoelectrics.

This report summarizes the FY21 Activities for EBS International Collaborations Work Package. The international collaborations work packages aim to leverage knowledge, expertise, and tools from the international nuclear waste community, as deemed relevant according to SFWST “roadmap” priorities. This report describes research and development (R&D) activities conducted during fiscal year 2021(FY21) specifically related to the Engineered Barrier System (EBS) R&D Work Package in the Spent Fuel and Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Energy (DOE). It fulfills the SFWST Campaign deliverable M4SF- 21SN010308062. The R&D activities described in this report focus on understanding EBS component evolution and interactions within the EBS, as well as interactions between the host media and the EBS. A primary goal is to advance the development of process models that can be implemented directly within the Generic Disposal System Analysis (GDSA) platform or that can contribute to the safety case in some manner such as building confidence, providing further insight into the processes being modeled, establishing better constraints on barrier performance, etc. Sandia National Laboratories is participating in THM modeling in the international projects EBS Task Force and DECOVALEX 2023. EBS Task Force, Task 11 is on modeling of laboratory-scale High Temperature Column Test conducted at Lawrence Berkeley National Laboratory. DECOVALEX 2023, Task C is on THM modeling of the full-scale emplacement experiment (FE experiment) at the Mont Terri Underground Rock Laboratory, Switzerland. This report summarizes Sandia’s progress in the modeling studies of DECOVALEX 2023, Task C. Modeling studies related to the High Temperature Column Test will be documented in future reports.

The final review for the FY21 Advanced Simulation and Computing (ASC) Computational Systems and Software Environments (CSSE) L2 Milestone #7840 was conducted on August 25th, 2021 at Sandia National Laboratories in Albuquerque, New Mexico. The review committee/panel unanimously agreed that the milestone has been successfully completed, exceeding expectations on several of the key deliverables.

Steam cracking of ethane, a non-catalytic thermochemical process, remains the dominant means of ethylene production. The severe reaction conditions and energy expenditure involved in this process incentivize the search for alternative reaction pathways and reactor designs which maximize ethylene yield while minimizing cost and energy input. Herein, we report a comparison of catalytic and non-catalytic non-oxidative dehydrogenation of ethane. We achieve ethylene yields as high as 67 % with an open tube quartz reactor without the use of a catalyst at residence times ∼4 s. The open tube reactor design promotes simplicity, low cost, and negligible coke formation. Pristine quartz tubes were most effective, since coke formation was detected when defects were introduced by scratching the surface of the quartz. Surprisingly, the addition of solids to the quartz tube, such as quartz sand, alumina powder, or even Pt-based intermetallic catalysts, led to lower ethylene yield. Pt alloy catalysts are effective at lower temperatures, such as at 575 °C, but conversion is limited due to thermodynamic constraints. When operated at industrially relevant temperatures, such as 700 °C and above, these catalysts were not stable in our tests, causing ethylene yield to drop below that of the open tube. These results suggest that future research on non-oxidative dehydrogenation should be directed at optimizing reactor designs to improve the conversion of ethane to ethylene, since this approach shows promise for decentralized production of ethylene from natural gas deposits.

The synthesis and electrical properties of a new misfit compound containing BiSe, Bi2Se3, and MoSe2 constituent layers are reported. The reaction pathway involves competition between the formation of (BiSe)1+x(Bi2Se3)1+y(BiSe)1+x(MoSe2) and [(Bi2Se3)1+y]2(MoSe2). Excess Bi and Se are required in the precursor to synthesize (BiSe)1+x(Bi2Se3)1+y(BiSe)1+x(MoSe2). High-angle annular dark field-scanning transmission electron microscopy (HAADF-STEM) confirm the stacking sequence of the heterostructure. Small grains of both 2H-and 1T-MoSe2 are observed in the MoSe2 layers. X-ray photoelectron spectroscopy (XPS) indicates that there is a significantly higher percentage of 1T-MoSe2 in (BiSe)1+x(Bi2Se3)1+y(BiSe)1+x(MoSe2) than in (BiSe)0.97(MoSe2), suggesting that more charge transfer to MoSe2 occurs due to the additional BiSe layer. The additional charge transfer results in (BiSe)1+x(Bi2Se3)1+y(BiSe)1+x(MoSe2) having a low resistivity (14-19 μω m) with metallic temperature dependence. The heterogeneous mix of MoSe2 polytypes observed in the XPS complicates the interpretation of the Hall data as two bands contribute to the electrical continuity.

Atomic precision advanced manufacturing (APAM) leverages the highly reactive nature of Si dangling bonds relative to H- or Cl-passivated Si to selectively adsorb precursor molecules into lithographically defined areas with sub-nanometer resolution. Due to the high reactivity of dangling bonds, this process is confined to ultra-high vacuum (UHV) environments, which currently limits its commercialization and broad-based appeal. In this work, we explore the use of halogen adatoms to preserve APAM-derived lithographic patterns outside of UHV to enable facile transfer into real-world commercial processes. Specifically, we examine the stability of H-, Cl-, Br-, and I-passivated Si(100) in inert N2 and ambient environments. Characterization with scanning tunneling microscopy and x-ray photoelectron spectroscopy (XPS) confirmed that each of the fully passivated surfaces were resistant to oxidation in 1 atm of N2 for up to 44 h. Varying levels of surface degradation and contamination were observed upon exposure to the laboratory ambient environment. Characterization by ex situ XPS after ambient exposures ranging from 15 min to 8 h indicated the Br– and I–passivated Si surfaces were highly resistant to degradation, while Cl–passivated Si showed signs of oxidation within minutes of ambient exposure. As a proof-of-principle demonstration of pattern preservation, a H–passivated Si sample patterned and passivated with independent Cl, Br, I, and bare Si regions was shown to maintain its integrity in all but the bare Si region post-exposure to an N2 environment. The successful demonstration of the preservation of APAM patterns outside of UHV environments opens new possibilities for transporting atomically-precise devices outside of UHV for integrating with non-UHV processes, such as other chemistries and commercial semiconductor device processes.

Lithium metal is considered the "holy grail"material to replace typical Li-ion anodes due to the absence of a host structure coupled with a high theoretical capacity. The absence of a host structure results in large volumetric changes when lithium is electrodeposited/dissolved, making the lithium prone to stranding and parasitic reactions with the electrolyte. Lithium research is focused on enabling highly reversible lithium electrodeposition/dissolution, which is important to achieving long cycle life. Understanding the various mechanisms of self-discharge is also critical for realizing practical lithium metal batteries but is often overlooked. In contrast to previous work, it is shown here that self-discharge via galvanic corrosion is negligible, particularly when lithium is cycled to relevant capacities. Rather, the continued electrochemical cycling of lithium metal results in self-discharge when periodic rest is applied during cycling. The extent of self-discharge can be controlled by increasing the capacity of plated lithium, tuning electrolyte chemistry, incorporating regular rest, or introducing lithiophilic materials. The Coulombic losses that occur during periodic rest are largely reversible, suggesting that the dominant self-discharge mechanism in this work is not an irreversible chemical process but rather a morphological process.

Porous nanoscale carbonaceous materials are widely employed for catalysis, separations, and electrochemical devices where device performance often relies upon specific and well-defined regular feature sizes. The use of block polymers as templates has enabled affordable and scalable production of diverse porous carbons. However, popular carbon preparations use equilibrating micelles which can change dimensions in response to the processing environment. Thus, polymer methods have not yet demonstrated carbon nanomaterials with constant average template diameter and tailored wall thickness. In contrast, persistent micelle templates (PMTs) use kinetic control to preserve constant micelle template diameters, and thus PMT has enabled constant pore diameter metrics. With PMT, the wall thickness is independently adjustable via the amount of material precursor added to the micelle templates. Previous PMT demonstrations relied upon thermodynamic barriers to inhibit chain exchange while in solution, followed by rapid evaporation and cross-linking of material precursors to mitigate micelle reorganization once the solvent evaporated. It is shown here that this approach, however, fails to deliver kinetic micelle control when used with slowly cross-linking material precursors such as those for porous carbons. A new modality for kinetic control over micelle templates, glassy-PMTs, is shown using an immobilized glassy micelle core composed of polystyrene (PS). Although PS based polymers have been used to template carbon materials before, all prior reports included plasticizers that prevented kinetic micelle control. Here the key synthetic conditions for carbon materials with glassy-PMT control are enumerated, including dependencies upon polymer block selection, block molecular mass, solvent selection, and micelle processing timeline. The use of glassy-PMTs also enables the direct observation of micelle cores by TEM which are shown to be commensurate with template dimensions. Glassy-PMTs are thus robust and insensitive to material processing kinetics, broadly enabling tailored nanomaterials with diverse chemistries.

Over the last year, the ECP xSDK-multiprecision effort has made tremendous progress in developing and deploying new mixed precision technology and customizing the algorithms for the hardware deployed in the ECP flagship supercomputers. The effort also has succeeded in creating a cross-laboratory community of scientists interested in mixed precision technology and now working together in deploying this technology for ECP applications. In this report, we highlight some of the most promising and impactful achievements of the last year. Among the highlights we present are: Mixed precision IR using a dense LU factorization and achieving a 1.8× speedup on Spock; results and strategies for mixed precision IR using a sparse LU factorization; a mixed precision eigenvalue solver; Mixed Precision GMRES-IR being deployed in Trilinos, and achieving a speedup of 1.4× over standard GMRES; compressed Basis (CB) GMRES being deployed in Ginkgo and achieving an average 1.4× speedup over standard GMRES; preparing hypre for mixed precision execution; mixed precision sparse approximate inverse preconditioners achieving an average speedup of 1.2×; and detailed description of the memory accessor separating the arithmetic precision from the memory precision, and enabling memory-bound low precision BLAS 1/2 operations to increase the accuracy by using high precision in the computations without degrading the performance. We emphasize that many of the highlights presented here have also been submitted to peer-reviewed journals or established conferences, and are under peer-review or have already been published.

In this study, model derivations are carried out of a dynamical system under base excitations with a piezoelectric energy harvesting absorber as the tuned-mass-damper. Additionally, amplitude stoppers are included to the absorber in order to create a broadband resonant response, increasing the window of operational use for energy harvesting and system's control. This study is unique in the fact that the energy harvester is coupled to the source of its excitation. A nonlinear reduced-order model is developed using Euler–Lagrange principle and the Galerkin method to accurately estimate the energy harvesting absorber's displacement, harvested power, and the oscillating response of the primary structure. The nonlinear interaction of the energy harvesting absorber and the amplitude stoppers are the focus of this study, where an in-depth investigation of bifurcation points of the primary structure and energy harvesting absorber responses is performed. Due to a transfer of energy between the primary structure and the absorber, it is shown that a soft stopper with stiffness $5 \times {10^3}\,{\text{N}}\;{{\text{m}}^{ - 1}}\,$ has great control of the primary structure with 60% of the uncontrolled amplitude being reduced, as well as an increase of the harvested energy. Medium stoppers with small initial gaps size and hard stoppers do not control the primary structure and show a decrease in the energy harvesting capabilities due to the activation of the nonlinear contact-impact interactions. Finally, these stoppers also generate aperiodic regions due to the possible presence of grazing bifurcations.

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heatgenerating nuclear waste disposition in deep-seated clay/shale/argillaceous rock. International collaboration activities such as heater tests, continuous field data monitoring, and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum and discrete) to tackle issues related to flow and transport at various scales of the host-rock, its interactions with barrier materials, and EBS design concept.

As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. We conclude with a discussion of expectations for future developments.

The following trade study was done to answer the following task from the Sandia JPL Collaboration for Europa Lander Statement of Work: Survey facility infrastructure SNL may have for performing aseptic assembly and integration of S/C and assess its suitability for PP applications.

The goal of this work is to develop a Bayesian framework to characterize the uncertainty of material response when using a nonlocal, homogenized model to describe wave propagation through heterogeneous, disordered materials. Our approach is based on an operator regression technique combined with Bayesian optimization, through which the nonlocal kernel for a specific disordered microstructure is investigated.

This report is a functional review of the radionuclide containment strategies of fluoride-salt-cooled high temperature reactor (FHR), molten salt reactor (MSR) and high temperature gas reactor (HTGR) systems. This analysis serves as a starting point for further, more in-depth analyses geared towards identifying phenomenological gaps that still exist, hindering the creation of a mechanistic source term for these reactor types. As background information to this review, an overview of how a mechanistic source term is created and used for consequence assessment necessary for licensing is provided. How a mechanistic source term is used within the Licensing Modernization Project (LMP) is also provided. Lastly, the characteristics of non-LWR mechanistic source terms are examined. This report does not assess the viability of any software system for use with advanced reactor designs, but instead covers system function requirements. Future work within the Nuclear Energy Advanced Modeling and Simulations (NEAMS) program will address such gaps. This document is an update of SAND 2020-6730. An additional chapter is included as well as edits to original content.

Information from 2015 annual report highlighting several tasks, including: Task 7: Research of microspectrophotometry for inspection and validation of laser color markings. Task 8: Investigate new laser fabrication techniques that produce color markings with improved corrosion resistance. Task 9: Research new methods for laser marking curved surfaces (and large areas). Task 10: Complete model simulations of laser-induced ripple formation-involves an ElectroMagnetic field solver.

This report summarizes initial results from a series of gun experiments which were conducted at the DICE facility. The target of these experiments was a modified metal slug composed of a tantalum/tungsten alloy (Ta-10W). The general geometry of the slug was a right circular cylinder with a through-hole cut normal to the cylinder's axis. In all experiments, hardened steel impactors were used, the desired impact velocity was 200 m/s, the slug was preheated to a target temperature of 175° C, photon doppler velocimetry (PDV) was used to measure the projectile velocity before and after impact, and the impact event was recorded with high-speed video. In two of the impacts the slug was oriented perpendicular to the projectile, while in the remaining two it was tilted 8° from normal. Initial high-speed speed video results showed slug failure in the tilted impact case, while the slug survived normal impacts. Recovery fixtures were used to preserve impacted slugs for future postmortem analysis. Discussions are included regarding improvements to potential future experiments involving these slugs.

Sandia will provide technical assistance to New Mexico Department of Health to provide analysis of SafeGraph mobility data (for which Sandia already has the data and a Data Use Agreement in place with the data provider). Sandia will produce analysis to determine the contribution of travel to SARS-CoV-2 spread within New Mexico.

Sandia will provide technical assistance to Helpful Engineering to develop and test the Universal Citizen Protection Device (UCPD) which is a UV-based filter-less PPE concept that aims to keep the Sars-CoV-2 virus out of eyes, nose and mouth with a 99%+ reliability. The heart of the device would be a concealed UV Chamber that decontaminates all air going in and out of the PPE. Helpful Engineering’s goal is to build this device to be reusable and cost less than $100 to construct and can be worn for 8 hours. The UCPD is an open source project, and once developed, prototyped, tested and approved, it will be shared with interested manufacturers globally.

MELCOR is a fully integrated, engineering-level computer code for modeling the progression of severe accidents in light water reactors (LWR) at nuclear power plants and nuclear fuel cycle facilities. Originally developed to assess severe accidents following Three Mile Island, MELCOR’s flexible modeling framework has enabled it to be applied to safety assessments of a much broader range of nuclear power reactor designs and other types of nuclear facilities processing radioactive material. Further, MELCOR can model a broad spectrum of severe accident phenomena such as thermal-hydraulic response in a reactor coolant system; core heat-up, degradation, and relocation; and transport behavior in both boiling water and pressurized water reactors.

Historically, nuclear component manufacturing vendors, from small businesses through large conglomerates, have felt compelled to obtain an American Society for Mechanical Engineers (ASME) Nuclear Certification, known colloquially as an "N-stamp", to assure supply chain quality standards that will be acceptable to regulators and safety concerns. Since the N-stamp quality standard is a U.S.-origin code, combined with the apparent decline in the U.S. nuclear industry alongside the growth of the Asian nuclear industry, there is the question of whether the rest of the world, including new entrants to the nuclear industry, also regard N-stamp as a needed certification. This study addresses this question through analysis of the entire N-stamp database of holders, and former holders, of N-stamp certificates of all types and for all regions worldwide from 1989-2020 (the dates available in the database). From this 30 years of data, we find that actually U.S.-based vendors still consistently obtain the largest number of N-stamps worldwide over all time periods, but also find that the countries participating in the N-stamp certification process has broadly expanded beyond just North America, Japan, S. Korea and Western Europe (the primary N-stamp recipients before the mid-2000's). We produced global heats maps and bar charts to illustrate our findings, as well as further investigation into why the data shows changes over time and region. We note that nuclear entities involved with Soviet-type reactors do not participate in the N-stamp process, but instead pursue the Russian version PNAE, which is substantially similar to the ASME code. We conclude that at least from the N-stamp database, the United States nuclear component manufacturing industry is alive and well, although there have been some consolidations, and that the ASME N-stamp appears to still be a valued certificate worldwide, including in China which now ranks second only to the United States in obtaining N-stamp certificates in recent years. We further note that the vendors of new reactor types, in particular High Temperature Gas-Cooled Reactors (HTGRs) and Small Modular Reactors (SMRs), are actively engaged with ASME (and other U.S.-based nuclear standards bodies such as the American Nuclear Society and Nuclear Energy Institute) to coordinate updates to the ASME N-Stamp criteria to ensure applicability of the code for these new designs. Implications of these findings include the following: The global use of the U.S.-origin N-stamp certification supports the view that, despite the decline of the U.S. nuclear industry, the United States remains an esteemed global leader in the area of nuclear safety. As the U.S. Government works to revitalize the U.S. nuclear industry, especially in the area of exports, it may be beneficial to leverage the global standing of the N-stamp certification. The findings indicate that the N-stamp database would be a useful tool for the U.S. Government to use to track the growth of the civil nuclear industry in foreign countries, under certain circumstances. The Excel-format N-stamp database produced as part of this study may be a useful tool for this purpose. N-stamp data may be a useful tool for foreign governments to use to identify nuclear manufacturers within their own country, especially to identify "targets" for outreach on nuclear export control compliance. The U.S. Government could carry this message to foreign partners during bilateral engagements or as part of Nuclear Suppliers Group (NSG) discussions on industry outreach.

The advanced materials team investigated the use of additively manufactured metallic lattice structures for mitigating impact response in a Davis gun earth penetrator impact experiment. High-fidelity finite element models were developed and validated with quasistatic experiments. These models were then used to simulate the response of such lattices when subjected to the acceleration loads expected in the Davis gun experiment. Results reveal how the impact mitigation performance of lattices can change drastically at a certain relative density. Based on these observations, an experiment deck was designed to probe the response of lattices with different relative densities during the Davis gun phase 2 shots. The expected performance of these lattices is predicted before testing based on simulation results. The results of the Davis gun phase 2 shots are expected to provide data which will be used to assess the predictive capability of the finite element simulations in such a complex impact environment.

We present an approach for constructing a surrogate from ensembles of information sources of varying cost and accuracy. The multifidelity surrogate encodes connections between information sources as a directed acyclic graph, and is trained via gradient-based minimization of a nonlinear least squares objective. While the vast majority of state-of-the-art assumes hierarchical connections between information sources, our approach works with flexibly structured information sources that may not admit a strict hierarchy. The formulation has two advantages: (1) increased data efficiency due to parsimonious multifidelity networks that can be tailored to the application; and (2) no constraints on the training data—we can combine noisy, non-nested evaluations of the information sources. Finally, numerical examples ranging from synthetic to physics-based computational mechanics simulations indicate the error in our approach can be orders-of-magnitude smaller, particularly in the low-data regime, than single-fidelity and hierarchical multifidelity approaches.

Currently, traditional methods such as short-term average/long-term average (STA/LTA) are used to detect arrivals in three-component seismic waveform data. Accurately establishing the identity and arrival of these waves is helpful in detecting and locating seismic events. Convolutional Neural Networks (CNNs) have been shown to significantly improve performance at local distances. This work will expand the use of CNNs to more remote distances and lower magnitudes. Sandia National Labs (SNL) will explore the advantages and limits of a particular approach and investigate requirements for expanding this technique to different types, distances, and magnitudes of events in the future. The team will describe detailed performance results of this method tuned on a curated dataset from Utah with its expert-defined arrival picks.

This Storm Water Pollution Prevention Plan (SWPPP) has been prepared for the Sandia National Laboratories Water Line Project – Northern Portion, in Livermore, CA. The project, located at 7011 East Avenue, and will entail the portion of the site north of the Arroyo Section. The project is comprised of 19,584 linear feet of water line improvements totaling approximately 9.0 acres. The property is owned by the U.S. Department of Energy, and managed and operated by National Technology & Engineering Solutions of Sandia, LLC with this project being developed by NTESS for Sandia National Laboratories.

Automatic mechanism generation is used to determine mechanisms for the CO₂ hydrogenation on Ni(111) in a two-stage process while considering the correlated uncertainty in DFT-based energetic parameters systematically. In a coarse stage, all the possible chemistry is explored with gas-phase products down to the ppb level, while a refined stage discovers the core methanation submechanism. Five thousand unique mechanisms were generated, which contain minor perturbations in all parameters. Global uncertainty assessment, global sensitivity analysis, and degree of rate control analysis are performed to study the effect of this parametric uncertainty on the microkinetic model predictions. Comparison of the model predictions with experimental data on a Ni/SiO₂ catalyst find a feasible set of microkinetic mechanisms within the correlated uncertainty space that are in quantitative agreement with the measured data, without relying on explicit parameter optimization. Global uncertainty and sensitivity analyses provide tools to determine the pathways and key factors that control the methanation activity within the parameter space. Together, these methods reveal that the degree of rate control approach can be misleading if parametric uncertainty is not considered. The procedure of considering uncertainties in the automated mechanism generation is not unique to CO₂ methanation and can be easily extended to other challenging heterogeneously catalyzed reactions.

As the demand for higher-performance batteries has increased, so has the body of research on theoretical high-capacity anode materials. However, the research has been hindered because the high-capacity anode material properties and interactions are not well understood, largely due to the difficulty of observing cycling in situ. Using electrochemical scanning transmission electron microscopy (ec-STEM), we report the real-time observation and electrochemical analysis of pristine tin (Sn) and titanium dioxide-coated Sn (TiO₂@Sn) electrodes during lithiation/delithiation. As expected, we observed a volume expansion of the pristine Sn electrodes during lithiation, but we further observed that the expansion was followed by Sn detachment from the current collector. Remarkably, although the TiO₂@Sn electrodes also exhibited similar volume expansion during lithiation, they showed no evidence of Sn detachment. We found that the TiO₂ surface layer acted as an electrochemically activated artificial solid-electrolyte interphase that serves to conduct Li ions. As a physical coating, it mechanically prevented Sn detachment following volume changes during cycling, providing significant degradation resistance and 80% Coulombic efficiency for a complete lithiation/delithiation cycle. Interestingly, upon delithiation, TiO₂@Sn electrode displayed a self-healing mechanism of small pore formation in the Sn particle followed by agglomeration into several larger pores as delithiation continued.

We describe a method to automatically generate an ion implantation recipe, a set of energies and fluences, to produce a desired defect density profile in a solid using the fewest required energies. We simulate defect density profiles for a range of ion energies, fit them with an appropriate function, and interpolate to yield defect density profiles at arbitrary ion energies. Given N energies, we then optimize a set of N energy-fluence pairs to match a given target defect density profile. Finally, we find the minimum N such that the error between the target defect density profile and the defect density profile generated by the N energy-fluence pairs is less than a given threshold. Inspired by quantum sensing applications with nitrogen-vacancy centers in diamond, we apply our technique to calculate optimal ion implantation recipes to create uniform-density 1 μm surface layers of 15N or vacancies (using 4He).

Nanostructures with a high density of interfaces, such as in nanoporous materials and nanowires, resist radiation damage by promoting the annihilation and migration of defects. This study details the size effect and origins of the radiation damage mechanisms in nanowires and nanoporous structures in model face-centered (gold) and body-centered (niobium) cubic nanostructures using accelerated multi-cascade atomistic simulations and in-situ ion irradiation experiments. Our results reveal three different size-dependent mechanisms of damage accumulation in irradiated nanowires and nanoporous structures: sputtering for very small nanowires and ligaments, the formation and accumulation of point defects and dislocation loops in larger nanowires, and a face-centered-cubic to hexagonal-close-packed phase transformation for a narrow range of wire diameters in the case of gold nanowires. Smaller nanowires and ligaments have a net effect of lowering the radiation damage as compared to larger wires that can be traced back to the fact that smaller nanowires transition from a rapid accumulation of defects to a saturation and annihilation mechanism at a lower dose than larger nanowires. These irradiation damage mechanisms are accompanied with radiation-induced surface roughening resulting from defect-surface interactions. Comparisons between nanowires and nanoporous structures show that the various mechanisms seen in nanowires provide adequate bounds for the defect accumulation mechanisms in nanoporous structures with the difference attributed to the role of nodes connecting ligaments in nanoporous structures. Taken together, our results shed light on the compounded, size-dependent mechanisms leading to the radiation resistance of nanowires and nanoporous structures.

The trapping of paraffins is beneficial compared to selective olefin adsorption for adsorptive olefin purification from a process engineering point of view. Here we demonstrate the use of a series of Zn2(X-bdc)2(dabco) (where X-bdc2−is bdc2−= 1,4-benzenedicarboxylate with substituting groups X, DM-bdc2−= 2,5-dimethyl-1,4-benzenedicarboxylate or TM-bdc2−= 2,3,5,6-tetramethyl-1,4-benzenedicarboxylate and dabco = diazabicyclo[2.2.2.]octane) metal-organic frameworks (MOFs) for the adsorptive removal of ethane from ethylene streams. The best performing material from this series is Zn2(TM-bdc)2(dabco) (DMOF-TM), which shows a high ethane uptake of 5.31 mmol g−1at 110 kPa, with a good IAST selectivity of 1.88 towards ethane over ethylene. Through breakthrough measurements a high productivity of 13.1 L kg−1per breakthrough is revealed with good reproducibility over five consecutive cycles. Molecular simulations show that the methyl groups of DMOF-TM are forming a van der Waals trap with the methylene groups from dabco, snuggly fitting the ethane. Further, rarely used high pressure coadsorption measurements, in pressure regimes that most scientific studies on hydrocarbon separation on MOFs ignore, reveal an increase in ethane capacity and selectivity for binary mixtures with increased pressures. The coadsorption measurements reveal good selectivity of 1.96 at 1000 kPa, which is verified also through IAST calculations up to 3000 kPa. This study overall showcases the opportunities that pore engineering by alkyl group incorporation and pressure increase offer to improve hydrocarbon separation in reticular materials.

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Fuel Cycle Technology (FCT) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). Two high priorities for SFWST disposal R&D are design concept development and disposal system modeling. These priorities are directly addressed in the SFWST Geologic Disposal Safety Assessment (GDSA) control account, which is charged with developing a geologic repository system modeling and analysis capability, and the associated software, GDSA Framework, for evaluating disposal system performance for nuclear waste in geologic media. GDSA Framework is supported by SFWST Campaign and its predecessor the Used Fuel Disposition (UFD) campaign. This report fulfills the GDSA Uncertainty and Sensitivity Analysis Methods work package (SF-21SN01030404) level 3 milestone, Uncertainty and Sensitivity Analysis Methods and Applications in GDSA Framework (FY2021) (M3SF-21SN010304042). It presents high level objectives and strategy for development of uncertainty and sensitivity analysis tools, demonstrates uncertainty quantification (UQ) and sensitivity analysis (SA) tools in GDSA Framework in FY21, and describes additional UQ/SA tools whose future implementation would enhance the UQ/SA capability of GDSA Framework. This work was closely coordinated with the other Sandia National Laboratory GDSA work packages: the GDSA Framework Development work package (SF-21SN01030405), the GDSA Repository Systems Analysis work package (SF-21SN01030406), and the GDSA PFLOTRAN Development work package (SF-21SN01030407). This report builds on developments reported in previous GDSA Framework milestones, particularly M3SF 20SN010304032.

Organic materials are an attractive choice for structural components due to their light weight and versatility. However, because they decompose at low temperatures relative to tradiational materials they pose a safety risk due to fire and loss of structural integrity. To quantify this risk, analysts use chemical kinetics models to describe the material pyrolysis and oxidation using thermogravimetric analysis. This process requires the calibration of many model parameters to closely match experimental data. Previous efforts in this field have largely been limited to finding a single best-fit set of parameters even though the experimental data may be very noisy. Furthermore the chemical kinetics models are often simplified representations of the true de- composition process. The simplification induces model-form errors that the fitting process cannot capture. In this work we propose a methodology for calibrating decomposition models to thermogravimetric analysis data that accounts for uncertainty in the model-form and experimental data simultaneously. The methodology is applied to the decomposition of a carbon fiber epoxy composite with a three-stage reaction network and Arrhenius kinetics. The results show a good overlap between the model predictions and thermogravimetric analysis data. Uncertainty bounds capture devia- tions of the model from the data. The calibrated parameter distributions are also presented. In conclusion, the distributions may be used in forward propagation of uncertainty in models that leverage this material.

Designing polymers with controlled nanoscale morphologies and scalable synthesis is of great interest in the development of fluorine-free materials for proton-exchange membranes in fuel cells. This study focuses on a precision polyethylene with phenylsulfonic acid branches at every fifth carbon, p5PhSA, with a high ion-exchange capacity (4.2 mmol/g). The polymers self-assemble into hydrophilic and hydrophobic co-continuous nanoscale domains. In the hydrated state, the hydrophilic domain, composed of polar sulfonic acid moieties and water, serves as a pathway for efficient mesoscopic proton conductivity. The morphology and proton transport of p5PhSA are evaluated under hydrated conditions using in situ X-ray scattering and electrochemical impedance spectroscopy techniques. At 40 °C and 95% relative humidity, the proton conductivity of p5PhSA is 0.28 S/cm, which is four times greater than Nafion 117 under the same conditions. Atomistic molecular dynamics (MD) simulations are also used to elucidate the interplay between the structure and the water dynamics. The MD simulations show strong nanophase separation between the percolated hydrophilic and hydrophobic domains over a wide range of water contents. The percolated hydrophilic nanoscale domain facilitates the rapid proton transport in p5PhSA and demonstrates the potential of precise hydrocarbon-based polymers as processible and effective proton-exchange membranes.

Continued creation of harmful emissions such as NO_x and soot from compression-ignition engines utilizing mixing-controlled combustion systems (i.e., diesel engines) remains a problem and is the subject of on-going research. The inherently high efficiency, relatively low cost, and numerous other desirable attributes of such engines, coupled with a widely supported infrastructure, motivates their continued advancement. Recently, a scientifically distinct and mechanically simple technology called ducted fuel injection (DFI) has shown a robust ability to allow such engines to operate with simultaneously low engine-out soot and NO_x emissions when it is employed with simulated exhaust-gas recirculation. To better understand the property ranges of sustainable, oxygenated-fuel blending stocks that will most improve engine performance, two oxygenated blendstocks were separately blended with a commercial diesel base fuel and tested within a heavy-duty diesel optical engine equipped with a four-duct DFI configuration. Conventional and crank-angle-resolved optical diagnostics were used to elucidate the effects of fuel ignition quality, oxygenate molecular structure, and overall oxygen content on engine performance.

The representation of material heterogeneity (also referred to as "spatial variation") plays a key role in the material failure simulation method used in ALEGRA. ALEGRA is an arbitrary Lagrangian-Eulerian shock and multiphysics code developed at Sandia National Laboratories and contains several methods for incorporating spatial variation into simulations. A desirable property of a spatial variation method is that it should produce consistent stochastic behavior regardless of the mesh used (a property referred to as "mesh independence"). However, mesh dependence has been reported using the Weibull distribution with ALEGRA's spatial variation method. This report describes efforts towards providing additional insight into both the theory and numerical experiments investigating such mesh dependence. In particular, we have implemented a discrete minimum order statistic model with properties that are theoretically mesh independent.

This paper describe our team's experience using minimega, a network emulation system using node and network virtualization, to support evaluation of a set of networked and distributed systems for topology discovery, traffic classification and engineering in the DARPA Searchlight program [18]. We present the methodology we developed to encode network and traffic definitions into an experiment description model, and how our tools compile this model onto the underlying minimega API. We then present three cases studies which demonstrate the ability of our EDM to support experiments with diverse network topologies, diverse traffic mixes, and networks with specialized layer-2 connectivity requirements. We conclude with the overall takeaways from using minimega to support our evaluation process.

Heteroepitaxy of β-phase gallium oxide (β-Ga₂O₃) thin films on foreign substrates shows promise for the development of next-generation deep ultraviolet solar blind photodetectors and power electronic devices. In this work, the influences of the film thickness and crystallinity on the thermal conductivity of ($\bar{2}01$)-oriented β-Ga₂O₃ heteroepitaxial thin films were investigated. Unintentionally doped β-Ga₂O₃ thin films were grown on c-plane sapphire substrates with off-axis angles of 0° and 6° toward $\langle$$11\bar{2}0$$\rangle$ via metal–organic vapor phase epitaxy (MOVPE) and low-pressure chemical vapor deposition. The surface morphology and crystal quality of the β-Ga₂O₃ thin films were characterized using scanning electron microscopy, X-ray diffraction, and Raman spectroscopy. The thermal conductivities of the β-Ga₂O₃ films were measured via time-domain thermoreflectance. The interface quality was studied using scanning transmission electron microscopy. The measured thermal conductivities of the submicron-thick β-Ga₂O₃ thin films were relatively low as compared to the intrinsic bulk value. The measured thin film thermal conductivities were compared with the Debye–Callaway model incorporating phononic parameters derived from first-principles calculations. The comparison suggests that the reduction in the thin film thermal conductivity can be partially attributed to the enhanced phonon-boundary scattering when the film thickness decreases. They were found to be a strong function of not only the layer thickness but also the film quality, resulting from growth on substrates with different offcut angles. Growth of β-Ga₂O₃ films on 6° offcut sapphire substrates was found to result in higher crystallinity and thermal conductivity than films grown on on-axis c-plane sapphire. However, the β-Ga₂O₃ films grown on 6° offcut sapphire exhibit a lower thermal boundary conductance at the β-Ga₂O₃/sapphire heterointerface. In addition, the thermal conductivity of MOVPE-grown ($\bar{2}01$)-oriented β-(Al_xGa_1–x)₂O₃ thin films with Al compositions ranging from 2% to 43% was characterized. Because of phonon-alloy disorder scattering, the β-(Al_xGa_1–x)₂O₃ films exhibit lower thermal conductivities (2.8–4.7 W/m∙K) than the β-Ga₂O₃ thin films. The dominance of the alloy disorder scattering in β-(Al_xGa_1–x)₂O₃ is further evidenced by the weak temperature dependence of the thermal conductivity. This work provides fundamental insight into the physical interactions that govern phonon transport within heteroepitaxially grown β-phase Ga₂O₃ and (Al_xGa_1–x)₂O₃ thin films and lays the groundwork for the thermal modeling and design of β-Ga₂O₃ electronic and optoelectronic devices.

Metallic enclosures are commonly used to protect electronic circuits against unwanted electromagnetic (EM) interactions. However, these enclosures may be sealed with imperfect mechanical seams or joints. These joints form narrow slots that allow external EM energy to couple into the cavity and then to the internal circuits. This coupled EM energy can severely affect circuit operations, particularly at the cavity resonance frequencies when the cavity has a high Q factor. To model these slots and the corresponding EM coupling, a thin-slot sub-cell model [1] , developed for slots in infinite ground plane and extended to numerical modeling of cavity-backed apertures, was successfully implemented in Sandia's electromagnetic code EIGER [2] and its next-generation counterpart Gemma [3]. However, this thin-slot model only considers resonances along the length of the slot. At sufficiently high frequencies, the resonances due to the slot depth must also be considered. Currently, slots must be explicitly meshed to capture these depth resonances, which can lead to low-frequency instability (due to electrically small mesh elements). Therefore, a slot sub-cell model that considers resonances in both length and depth is needed to efficiently and accurately capture the slot coupling.

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (Tg) and chain length and thus have the potential for phase separation. One crosslinker has a high Tg corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low Tg. A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

We study the deformation of tantalum under extreme loading conditions. Experimental velocity data are drawn from both ramp loading experiments on Sandia's Z-machine and gas gun compression experiments. The drive conditions enable the study of materials under pressures greater than 100 GPa. We provide a detailed forward model of the experiments including a model of the magnetic drive for the Z-machine. Utilizing these experiments, we simultaneously infer several different types of physically motivated parameters describing equation of state, plasticity, and anelasticity via the computational device of Bayesian model calibration. Characteristics of the resulting calculated posterior distributions illustrate relationships among the parameters of interest via the degree of cross correlation. The calibrated velocity traces display good agreement with the experiments up to experimental uncertainty as well as improvement over previous calibrations. Examining the Z-shots and gun-shots together and separately reveals a trade-off between accuracy and transferability across different experimental conditions. Implications for model calibration, limitations from model form, and suggestions for improvements are discussed.

Electric vehicles (EVs) represent an important socio-economic development opportunity for islands and remote locations because they can lead to reduced fuel imports, electricity storage, grid services, and environmental and health benefits. This paper presents an overview of opportunities, challenges, and examples of EVs in islands and remote power systems, and is meant to provide background to researchers, utilities, energy offices, and other stakeholders interested in the impacts of electrification of transportation. The impact of uncontrolled EV charging on the electric grid operation is discussed, as well as several mitigation strategies. Of particular importance in many islands and remote systems is taking advantage of local resources by combining renewable energy and EV charging. Policy and economic issues are presented, with emphasis on the need for an overarching energy policy to guide the strategies for EVs growth. The key conclusion of this paper is that an orderly transition to EVs, one that maximizes benefits while addressing the challenges, requires careful analysis and comprehensive planning.

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

Superionic phases of bulk anhydrous salts based on large cluster-like polyhedral (carba)borate anions are generally stable only well above room temperature, rendering them unsuitable as solid-state electrolytes in energy-storage devices that typically operate at close to room temperature. To unlock their technological potential, strategies are needed to stabilize these superionic properties down to subambient temperatures. One such strategy involves altering the bulk properties by confinement within nanoporous insulators. In the current study, the unique structural and ion dynamical properties of an exemplary salt, NaCB11H12, nanodispersed within porous, high-surface-area silica via salt-solution infiltration were studied by differential scanning calorimetry, X-ray powder diffraction, neutron vibrational spectroscopy, nuclear magnetic resonance, quasielastic neutron scattering, and impedance spectroscopy. Combined results hint at the formation of a nanoconfined phase that is reminiscent of the high-temperature superionic phase of bulk NaCB11H12, with dynamically disordered CB11H12- anions exhibiting liquid-like reorientational mobilities. However, in contrast to this high-temperature bulk phase, the nanoconfined NaCB11H12 phase with rotationally fluid anions persists down to cryogenic temperatures. Moreover, the high anion mobilities promoted fast-cation diffusion, yielding Na+ superionic conductivities of ∼0.3 mS/cm at room temperature, with higher values likely attainable via future optimization. It is expected that this successful strategy for conductivity enhancement could be applied as well to other related polyhedral (carba)borate-based salts. Thus, these results present a new route to effectively utilize these types of superionic salts as solid-state electrolytes in future battery applications.

Poly(carbon monofluoride), or (CF)n, is a layered fluorinated graphite material consisting of nanosized platelets. Here, we present experimental multidimensional solid-state NMR spectra of (CF)n, supported by density functional theory (DFT) calculations of NMR parameters, which overhauls our understanding of structure and bonding in the material by elucidating many ways in which disorder manifests. We observe strong 19F NMR signals conventionally assigned to elongated or "semi-ionic"C-F bonds and find that these signals are in fact due to domains where the framework locally adopts boat-like cyclohexane conformations. We calculate that C-F bonds are weakened but are not elongated by this conformational disorder. Exchange NMR suggests that conformational disorder avoids platelet edges. We also use a new J-resolved NMR method for disordered solids, which provides molecular-level resolution of highly fluorinated edge states. The strings of consecutive difluoromethylene groups at edges are relatively mobile. Topologically distinct edge features, including zigzag edges, crenellated edges, and coves, are resolved in our samples by solid-state NMR. Disorder should be controllable in a manner dependent on synthesis, affording new opportunities for tuning the properties of graphite fluorides.

Finding dense regions of graphs is fundamental in graph mining. We focus on the computation of dense hierarchies and regions with graph nuclei - -a generalization of k-cores and trusses. Static computation of nuclei, namely through variants of 'peeling', are easy to understand and implement. However, many practically important graphs undergo continuous change. Dynamic algorithms, maintaining nucleus computations on dynamic graph streams, are nuanced and require significant effort to port between nuclei, e.g., from k-cores to trusses. We propose a unifying framework to maintain nuclei in dynamic graph streams. First, we show no dynamic algorithm can asymptotically beat re-computation, highlighting the need to experimentally understand variability. Next, we prove equivalence between k-cores on a special hypergraph and nuclei. Our algorithm splits the problem into maintaining the special hypergraph and maintaining k-cores on it. We implement our algorithm and experimentally demonstrate improvements up to 108 x over re-computation. We show algorithmic improvements on k-cores apply to trusses and outperform truss-specific implementations.

Recently, lithium nitride (Li3N) has been proposed as a chemical warfare agent (CWA) neutralization reagent for its ability to produce nucleophilic ammonia molecules and hydroxide ions in aqueous solution. Quantum chemical calculations can provide insight into the Li3N neutralization process that has been studied experimentally. Here, we calculate reaction-free energies associated with the Li3N-based neutralization of the CWA VX using quantum chemical density functional theory and ab initio methods. We find that alkaline hydrolysis is more favorable to either ammonolysis or neutral hydrolysis for initial P-S and P-O bond cleavages. Reaction-free energies of subsequent reactions are calculated to determine the full reaction pathway. Notably, products predicted from favorable reactions have been identified in previous experiments.

Costs to permit Marine Energy projects are poorly understood. In this paper we examine environmental compliance and permitting costs for 19 projects in the U.S., covering the last 2 decades. Guided discussions were conducted with developers over a 3-year period to obtain historical and ongoing project cost data relative to environmental studies (e.g., baseline or pre-project site characterization as well as post-installation effects monitoring), stakeholder outreach, and mitigation, as well as qualitative experience of the permitting process. Data are organized in categories of technology type, permitted capacity, pre-and post-installation, geographic location, and funding types. We also compare our findings with earlier logic models created for the Department of Energy (i.e., Reference Models). Environmental studies most commonly performed were for Fish and Fisheries, Noise, Marine Habitat/Benthic Studies and Marine Mammals. Studies for tidal projects were more expensive than those performed for wave projects and the range of reported project costs tended to be wider than ranges predicted by logic models. For eight projects reporting full project costs, from project start to FERC or USACE permit, the average amount for environmental permitting compliance was 14.6%.

PCalc is a software tool that computes travel-time predictions, ray path geometry and model queries. This software has a rich set of features, including the ability to use custom 3D velocity models to compute predictions using a variety of geometries. The PCalc software is especially useful for research related to seismic monitoring applications.

Automated vehicles (AV) hold great promise for improving safety, as well as reducing congestion and emissions. In order to make automated vehicles commercially viable, a reliable and highperformance vehicle-based computing platform that meets ever-increasing computational demands will be key. Given the state of existing digital computing technology, designers will face significant challenges in meeting the needs of highly automated vehicles without exceeding thermal constraints or consuming a large portion of the energy available on vehicles, thus reducing range between charges or refills. The accompanying increases in energy for AV use will place increased demand on energy production and distribution infrastructure, which also motivates increasing computational energy efficiency.

The generalized singular value decomposition (GSVD) is a valuable tool that has many applications in computational science. However, computing the GSVD for large-scale problems is challenging. Motivated by applications in hyper-differential sensitivity analysis (HDSA), we propose new randomized algorithms for computing the GSVD which use randomized subspace iteration and weighted QR factorization. Detailed error analysis is given which provides insight into the accuracy of the algorithms and the choice of the algorithmic parameters. We demonstrate the performance of our algorithms on test matrices and a large-scale model problem where HDSA is used to study subsurface flow.

Underground explosions nonlinearly deform the surrounding earth material and can interact with the free surface to produce spall. However, at typical seismological observation distances the seismic wavefield can be accurately modeled using linear approximations. Although nonlinear algorithms can accurately simulate very near field ground motions, they are computationally expensive and potentially unnecessary for far field wave simulations. Conversely, linearized seismic wave propagation codes are orders of magnitude faster computationally and can accurately simulate the wavefield out to typical observational distances. Thus, devising a means of approximating a nonlinear source in terms of a linear equivalent source would be advantageous both for scenario modeling and for interpretation of seismic source models that are based on linear, far-field approximations. This allows fast linear seismic modeling that still incorporates many features of the nonlinear source mechanics built into the simulation results so that one can have many of the advantages of both types of simulations without the computational cost of the nonlinear computation. In this report we first show the computational advantage of using linear equivalent models, and then discuss how the near-source (within the nonlinear wavefield regime) environment affects linear source equivalents and how well we can fit seismic wavefields derived from nonlinear sources.

In situ TEM/SEM microscopy to investigate the structural evolution that occurs due to various extreme environments

Based on the latest DOE (Department of Energy) milestones, Sandia needs to convert to IPv6 (Internet Protocol version 6)-only networks over the next 5 years. Our original IPv6 migration plan did not include migrating to IPv6-only networks at any point within the next 10 years, so it must necessarily change. To be successful in this endeavor, we need to evaluate technologies that will enable us to deploy IPv6-only networks early without creating system stability or security issues. We have set up a test environment using technology representative of our production network where we configured and evaluated industry standard translation technologies and techniques. Based on our results, bidirectional translation between IPv4 (Internet Protocol version 4) and IPv6 is achievable with our current equipment, but due to the complexity of the configuration, may not scale well to our production environment.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

LocOO3D is a software tool that computes geographical locations for seismic events at regional to global scales. This software has a rich set of features, including the ability to use custom 3D velocity models, correlated observations and master event locations. The LocOO3D software is especially useful for research related to seismic monitoring applications, since it allows users to easily explore a variety of location methods and scenarios and is compatible with the CSS3.0 data format used in monitoring applications. The LocOO3D software, User's Manual, and Examples are available on the web at: https://github.com/sandialabs/LocOO3D For additional information on GeoTess, SALSA3D, RSTT, and other related software, please see: https://github.com/sandialabs/GeoTessJava, www.sandia.gov/geotess, www.sandia.gov/salsa3d, and www.sandia.gov/rstt