Publications

Abstract not provided.

Concentrating solar mirrors need to have correct geometry to accurately focus sunlight onto a receiver. SOFAST is a tool created at Sandia National Laboratories which measures a mirror’s surface normal at a dense sampling of points, enabling evaluation of a given mirror’s expected performance. SOFAST supports general mirrors, including concentrating solar heliostats, troughs, and dishes. The examples presented in this report predominantly focus on heliostats.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The objective of the project is to demonstrate that critical control functions can be implemented using little resources in modern microelectronics. A finite state machine (Figure 1) is implemented onto a field programmable gate array (FPGA). The functionality of the system is demonstrated by sending binary instructions to the controller. The controller transmits patterns through an LED, controls an electromechanical device, and uses pulse-width modulation (PWM) for radar functions.

This report summarizes the international collaborations conducted by Sandia funded by the US Department of Energy Office (DOE) of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-24SN010303063. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss granular salt reconsolidation (KOMPASS/MEASURES), engineered barriers (RANGERS), numerical model comparison (DECOVALEX), an NEA Salt Club working group on the development of scenarios as part of the performance assessment development process, and progress on seal percolation experiments with German colleagues at Gesellschaft für Anlagen- und Reaktorsicherheit (GRS). Finally, we summarize events related to the US/German Workshop on Repository Research, Design and Operations.

This study conducts a comparative analysis, using non-equilibrium Green’s functions (NEGF), of two state-of-the-art two-well (TW) Terahertz Quantum Cascade Lasers (THz QCLs) supporting clean 3-level systems. The devices have nearly identical parameters and the NEGF calculations with an abrupt-interface roughness height of 0.12 nm predict a maximum operating temperature (T_max) of ~ 250 K for both devices. However, experimentally, one device reaches a T_max of ~ 250 K and the other a T_max of only ~ 134 K. Both devices were fabricated and measured under identical conditions in the same laboratory, with high quality processes as verified by reference devices. The main difference between the two devices is that they were grown in different MBE reactors. Our NEGF-based analysis considered all parameters related to MBE growth, including the maximum estimated variation in aluminum content, growth rate, doping density, background doping, and abrupt-interface roughness height. From our NEGF calculations it is evident that the sole parameter to which a drastic drop in T_max could be attributed is the abrupt-interface roughness height. We can also learn from the simulations that both devices exhibit high-quality interfaces, with one having an abrupt-interface roughness height of approximately an atomic layer and the other approximately a monolayer. However, these small differences in interface sharpness are the cause of the large performance discrepancy. This underscores the sensitivity of device performance to interface roughness and emphasizes its strategic role in achieving higher operating temperatures for THz QCLs. We suggest Atom Probe Tomography (APT) as a path to analyze and measure the (graded)-interfaces roughness (IFR) parameters for THz QCLs, and subsequently as a design tool for higher performance THz QCLs, as was done for mid-IR QCLs. Our study not only addresses challenges faced by other groups in reproducing the record T_max of ~ 250 K and ~ 261 K but also proposes a systematic pathway for further improving the temperature performance of THz QCLs beyond the state-of-the-art.

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

Understanding and accurately characterizing energy dissipation mechanisms in civil structures during earthquakes is an important element of seismic assessment and design. The most commonly used model is attributed to Rayleigh. This paper proposes a systematic approach to quantify the uncertainty associated with Rayleigh's damping model. Bayesian calibration with embedded model error is employed to treat the coefficients of the Rayleigh model as random variables using modal damping ratios. Through a numerical example, we illustrate how this approach works and how the calibrated model can address modeling uncertainty associated with the Rayleigh damping model.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Topological insulator–magnetic insulator (TI–MI) heterostructures hold significant promise in the field of spintronics, offering the potential for manipulating magnetization through topological surface state–enabled spin-orbit torque. However, many TI–MI interfaces are plagued by issues such as contamination within the magnetic insulator layer and the presence of a low-density transitional region of the topological insulator. These interfacial challenges often obscure the intrinsic behavior of the TI–MI system. In this study, we addressed these challenges by depositing sputtered Bi₂⁢Te₃ on liquid phase epitaxy grown Y₃⁢Fe₅⁢O₁₂/Gd₃⁢Ga₅⁢O₁₂. The liquid phase epitaxy grown Y₃⁢Fe₅⁢O₁₂ has been previously shown to have exceptional interface quality, without an extended transient layer derived from interdiffusion processes of the substrate or impurity ions, thereby eliminating rare-earth impurity-related losses in the MI at low temperatures. At the TI–MI interface, high-resolution depth-sensitive polarized neutron reflectometry confirmed the absence of a low-density transitional growth region of the TI. By overcoming these undesirable interfacial effects, we isolate and probe the intrinsic low-temperature magnetization dynamics and transport properties of the TI–MI interface. In conclusion, our findings revealed strong spin pumping at low temperatures, accompanied by an additional in-plane anisotropy. The enhanced spin pumping at low temperatures is correlated with the observed suppression of bulk conduction and the weak antilocalization in the TI film, highlighting the interplay between the transport and spin pumping behavior in the TI–MI system.

We present a comprehensive benchmarking framework for evaluating machine-learning approaches applied to phase-field problems. This framework focuses on four key analysis areas crucial for assessing the performance of such approaches in a systematic and structured way. Firstly, interpolation tasks are examined to identify trends in prediction accuracy and accumulation of error over simulation time. Secondly, extrapolation tasks are also evaluated according to the same metrics. Thirdly, the relationship between model performance and data requirements is investigated to understand the impact on predictions and robustness of these approaches. Finally, systematic errors are analyzed to identify specific events or inadvertent rare events triggering high errors. Quantitative metrics evaluating the local and global description of the microstructure evolution, along with other scalar metrics representative of phase-field problems, are used across these four analysis areas. This benchmarking framework provides a path to evaluate the effectiveness and limitations of machine-learning strategies applied to phase-field problems, ultimately facilitating their practical application.

Solvent expulsion away from an intervening region between two approaching particles plays important roles in particle aggregation yet remains poorly understood. In this work, we use metadynamics molecular simulations to study the free energy landscape of removing water molecules from gibbsite and pyrophyllite slit pores representing the confined spaces between two approaching particles. For gibbsite, removing water from the intervening region is both entropically and enthalpically unfavorable. The closer the particles approach each other, the harder it is to expel water molecules. For pyrophyllite, water expulsion is spontaneous, which is different from the gibbsite system. A smaller pore makes the water removal more favorable. When water is being drained from the intervening region, single chains of water molecules are observed in gibbsite pore, while in pyrophyllite pore water cluster is usually observed. Water-gibbsite hydrogen bonds help stabilize water chains, while water forms clusters in pyrophyllite pore to maximize the number of hydrogen bonds among themselves. This work provides the first assessment into the energetics and structure of water being drained from the intervening region between two approaching particles during oriented attachment and aggregation.

This study investigates high performance electrochromic windows used on a passive house and residential dwelling to IECC 2021 (i.e., IECC dwelling). In the lab, the electrochromic film switches transmitted solar heat gain coefficient (SHGC) from 0.09 to 0.7 and visible transmittance from 0.15 to 0.82 with power consumption of 1.23 W/m2 during switching times less than 3 minutes. We extrapolate these results to a window assembly. Building energy models of the houses were evaluated in Santa Fe, New Mexico. A Monte Carlo analysis for 2020, 2040, 2060, and 2080 was conducted for Shared Socioeconomic Pathways 2-4.5, 3-7.0, and 5-8.5. Cases with and without the electrochromic windows and with and without electricity were used to determine energy use intensity and hours beyond thermal safety thresholds. The passive house showed 1.3-3.1% mean energy savings and the IECC dwelling 4.4-5.1% with electrochromic efficiency benefits growing into the future for both cases. Even so, overall savings decrease into the future for the passive house, due to growth in cooling load being dominant, conversely overall energy savings increase into the future for the IECC dwelling due to heating loads being dominant. For thermal resilience, the passive house exhibited a mean percent decrease of 0.02-0.31% hours in the extreme caution (i.e., > 32.2 ∘C, ≤ 39.4 ∘C) range while the IECC dwelling exhibited 0.38-4.38%. The study therefore shows that electrochromic windows will have smaller benefits for the passive house in comparison to the IECC dwelling. The relationship between electrochromic windows is shown to have a complex relationship between house efficiency and climate change by these results.

Here, it is demonstrated that single-crystalline and highly doped GaAs membranes are excellent candidates for realizing infrared-transparent shields of electromagnetic interference at millimeter frequencies. Measured optical transmittance spectra for the semiconductor membranes show resonant features between 750 and 2500 nm, with a 100% maximum transmittance. The shielding effectiveness of the membranes is extracted from measured scattering parameters between 65 and 85 GHz. Selected GaAs membranes and membranes/polyamide films exhibit shielding effectiveness ranging from 22 to 40 dB, which are suitable values to ensure the safe operation of infrared devices for commercial applications. Theoretical calculations based on a plane wave model show that the interplay of primary reflection and multiple internal reflections of the radio-frequency waves results in broadband shielding capabilities of the membrane between 10 and 300 GHz.

This report provides documentation for the Sandia Toolkit (STK) modules. STK modules are intended to provide infrastructure that assists the development of computational engineering software such as finite-element analysis applications. STK includes modules for unstructured-mesh data structures, reading/writing mesh files, geometric proximity search, transfers, MPMD coupling support, and various other utilities. This document contains a chapter for each module, and each chapter contains overview descriptions and usage examples. Usage examples are primarily code listings which are generated from working test programs that are included in the STK code-base. A goal of this approach is to ensure that the usage examples will not fall out of date.

The laser system is the most complex component of a light-pulse atom interferometer (LPAI), controlling frequencies and intensities of multiple laser beams to configure quantum gravity and inertial sensors. Its main functions include cold-atom generation, state preparation, state-selective detection, and generating a coherent two-photon process for the light-pulse sequence. To achieve substantial miniaturization and ruggedization, we integrate key laser system functions onto a photonic integrated circuit. Our study focuses on a high-performance silicon photonic suppressed-carrier single-sideband (SC-SSB) modulator at 1560 nanometers, capable of dynamic frequency shifting within the LPAI. By independently controlling radio frequency (RF) channels, we achieve 30-decibel carrier suppression and unprecedented 47.8-decibel sideband suppression at peak conversion efficiency of -6.846 decibels (20.7%). We investigate imbalances in both amplitudes and phases between the RF signals. Using this modulator, we demonstrate cold-atom generation, state-selective detection, and atom interferometer fringes to estimate gravitational acceleration, g ≈ 9.77 ± 0.01 meters per second squared, in a rubidium (87Rb) atom system.

Helium-4-based scintillation detector technology is emerging as a strong alternative to pulse-shape discrimination-capable organic scintillators for fast neutron detection and spectroscopy, particularly in extreme gamma-ray environments. The ⁴He detector is intrinsically insensitive to gamma radiation, as it has a relatively low cross-section for gamma-ray interactions, and the stopping power of electrons in the ⁴He medium is low compared to that of ⁴He recoil nuclei. Consequently, gamma rays can be discriminated by simple energy deposition thresholding instead of the more complex pulse shape analysis. The energy resolution of ⁴He scintillation detectors has not yet been well-characterized over a broad range of energy depositions, which limits the ability to deconvolve the source spectra. In this work, an experiment was performed to characterize the response of an Arktis S670 ⁴He detector to nuclear recoils up to 9 MeV. The ⁴He detector was positioned in the center of a semicircular array of organic scintillation detectors operated in coincidence. Deuterium–deuterium and deuterium–tritium neutron generators provided monoenergetic neutrons, yielding geometrically constrained nuclear recoils ranging from 0.0925 to 8.87 MeV. The detector response provides evidence for scintillation linearity beyond the previously reported energy range. Finally, the measured response was used to develop an energy resolution function applicable to this energy range for use in high-fidelity detector simulations needed by future applications.

In global sensitivity analysis (GSA) of a model, a proper convergence analysis of metrics is essential for ensuring a level of confidence or trustworthiness in sensitivity results obtained, yet is somewhat deficient in practice. The level of confidence in sensitivity measures, particularly in relation to their influence and support for decisions from scientific, social and policy perspectives, is heavily reliant on the convergence of GSA. We review the literature and summarize the available methods for monitoring and assessing convergence of sensitivity measures based on application purposes. The aim is to expose the various choices for convergence assessment and encourage further testing of available methods to clarify their level of robustness. Furthermore, the review identifies a pressing need for comparative studies on convergence assessment methods to establish a clear hierarchy of effectiveness and encourages the adoption of systematic approaches for enhanced robustness in sensitivity analysis.

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln3+) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln3+ coordinated to ligands in water and embedded in metal-organic frameworks (MOFs), and ΔΔE for Ln3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln3+ in the periodic table. Static DFT calculations with a full Ln3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. We also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

Boron nitride nanotubes (BNNTs) are high-strength, high-modulus nanotubes with high thermal and oxidative stabilities. Two hybrid composites were prepared with satin weave carbon fiber (CF) and resole-type phenolic resin: one with surface layers of BNNTs and one with alternating interlayers of BNNTs. The samples were subjected to hot jet tests that simulate realistic high-pressure-temperature conditions to understand the behavior of BNNTs under high-pressure erosion. Adding BNNTs to CF/phenolic laminates enhanced the ablation resistance by reinforcing the char material and mitigated localized thermal damage. Hybrid laminates exhibited up to 14% lower weight loss, 55% increase in flexural modulus, higher thermal diffusivity, and improved char yield and microstructure compared to CF/phenolic samples. The surface layer hybrid had many surviving nanotubes reinforcing the char and crystalline oxide structures that could mitigate further oxygen diffusion. Further, various characterization methods were used to deduce possible mechanisms and their products, indicating that BNNTs could serve as growth templates for direct crystalline boron oxide formation. Overall, hybrid BNNT/CF/phenolic laminates displayed better ablation resistance and favorable microstructure evolution under high-pressure conditions.

Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

Analytical and semi–analytical models for stream depletion with transient stream stage drawdown induced by groundwater pumping are developed to address a deficiency in existing models, namely, the use of a fixed stream stage condition at the stream–aquifer interface. Here field data are presented to demonstrate that stream stage drawdown does indeed occur in response to groundwater pumping near aquifer–connected streams. A model that predicts stream depletion with transient stream drawdown is developed based on stream channel mass conservation and finite stream channel storage. The resulting models are shown to reduce to existing fixed–stage models in the limit as stream channel storage becomes infinitely large, and to the confined aquifer flow with a no–flow boundary at the streambed in the limit as stream storage becomes vanishingly small. The model is applied to field measurements of aquifer and stream drawdown, giving estimates of aquifer hydraulic parameters, streambed conductance, and a measure of stream channel storage. The results of the modeling and data analysis presented herein have implications for sustainable groundwater management.

This dataset is comprised of a library of atomistic structure files and corresponding X-ray diffraction (XRD) profiles and vibrational density of states (VDoS) profiles for bulk single crystal silicon (Si), gold (Au), magnesium (Mg), and iron (Fe) with and without disorder introduced into the atomic structure and with and without mechanical loading. Included with the atomistic structure files are descriptor files that measure the stress state, phase fractions, and dislocation content of the microstructures. All data was generated via molecular dynamics or molecular statics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. This dataset can inform the understanding of how local or global changes to a materials microstructure can alter their spectroscopic and diffraction behavior across a variety of initial structure types (cubic diamond, face-centered cubic (FCC), hexagonal close-packed (HCP), and body-centered cubic (BCC) for Si, Au, Mg, and Fe, respectively) and overlapping changes to the microstructure (i.e., both disorder insertion and mechanical loading).

Efficient carbon capture requires engineered porous systems that selectively capture CO2 and have low energy regeneration pathways. Porous liquids (PLs), solvent-based systems containing permanent porosity through the incorporation of a porous host, increase the CO2 adsorption capacity. A proposed mechanism of PL regeneration is the application of isostatic pressure in which the dissolved nanoporous host is compressed to alter the stability of gases in the internal pore. This regeneration mechanism relies on the flexibility of the porous host, which can be evaluated through molecular simulations. Here, the flexibility of porous organic cages (POCs) as representative porous hosts was evaluated, during which pore windows decreased by 10-40% at 6 GPa. POCs with sterically smaller functional groups, such as the 1,2-ethane in the CC1 POC resulted in greater imine cage flexibility relative to those with sterically larger functional groups, such as the cyclohexane in the CC3 POC that protected the imine cage from the application of pressure. Structural changes in the POC also caused CO2 adsorption to be thermodynamically unfavorable beginning at ∼2.2 GPa in the CC1 POC, ∼1.1 GPa in the CC3 POC, and ∼1.0 GPa in the CC13 POC, indicating that the CO2 would be expelled from the POC at or above these pressures. Energy barriers for CO2 desorption from inside the POC varied based on the geometry of the pore window and all the POCs had at least one pore window with a sufficiently low energy barrier to allow for CO2 desorption under ambient temperatures. The results identified that flexibility of the CC1, CC3, or CC13 POCs under compression can result in the expulsion of captured gas molecules.

Abstract not provided.

Abstract not provided.

High-throughput image segmentation of atomic resolution electron microscopy data poses an ongoing challenge for materials characterization. In this paper, we investigate the application of the polyhedral template matching (PTM) method, a technique widely employed for visualizing three-dimensional (3D) atomistic simulations, to the analysis of two-dimensional (2D) atomic resolution electron microscopy images. This technique is complementary with other atomic resolution data reduction techniques, such as the centrosymmetry parameter, that use the measured atomic peak positions as the starting input. Furthermore, since the template matching process also gives a measure of the local rotation, the method can be used to segment images based on local orientation. We begin by presenting a 2D implementation of the PTM method, suitable for atomic resolution images. We then demonstrate the technique's application to atomic resolution scanning transmission electron microscopy images from close-packed metals, providing examples of the analysis of twins and other grain boundaries in FCC gold and martensite phases in 304 L austenitic stainless steel. Finally, we discuss factors, such as positional errors in the image peak locations, that can affect the accuracy and sensitivity of the structural determinations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Presentation for the GeoMIP 2024 Meeting

Abstract not provided.

Abstract not provided.