To detect a specific radio-frequency (rf) magnetic field, rf optically pumped magnetometers (OPMs) require a static magnetic field to set the Larmor frequency of the atoms equal to the frequency of interest. However, unshielded and variable magnetic field environments (e.g., an rf OPM on a moving platform) pose a problem for rf OPM operation. Here, we demonstrate the use of a natural-abundance rubidium vapor to make a comagnetometer to address this challenge. Our implementation builds upon the simultaneous application of several OPM techniques within the same vapor cell. First, we use a modified implementation of an OPM variometer based on 87Rb to detect and actively cancel unwanted external fields at frequencies ≲60 Hz using active feedback to a set of field control coils. In this experiment, we exploit this stabilized field environment to implement a high-sensitivity rf magnetometer using 85Rb. Using this approach, we demonstrate the ability to measure rf fields with a sensitivity of approximately 9 fT Hz-1/2 inside a magnetic shield in the presence of an applied field of approximately 20 μT along three mutually orthogonal directions. This demonstration opens up a path toward completely unshielded operation of a high-sensitivity rf OPM.
The long-standing problem of predicting the electronic structure of matter on ultra-large scales (beyond 100,000 atoms) is solved with machine learning.
Physics-constrained machine learning is emerging as an important topic in the field of machine learning for physics. One of the most significant advantages of incorporating physics constraints into machine learning methods is that the resulting model requires significantly less data to train. By incorporating physical rules into the machine learning formulation itself, the predictions are expected to be physically plausible. Gaussian process (GP) is perhaps one of the most common methods in machine learning for small datasets. In this paper, we investigate the possibility of constraining a GP formulation with monotonicity on three different material datasets, where one experimental and two computational datasets are used. The monotonic GP is compared against the regular GP, where a significant reduction in the posterior variance is observed. The monotonic GP is strictly monotonic in the interpolation regime, but in the extrapolation regime, the monotonic effect starts fading away as one goes beyond the training dataset. Imposing monotonicity on the GP comes at a small accuracy cost, compared to the regular GP. The monotonic GP is perhaps most useful in applications where data are scarce and noisy, and monotonicity is supported by strong physical evidence.
The recently-developed ability to control phosphorous-doping of silicon at an atomic level using scanning tunneling microscopy, a technique known as atomic precision advanced manufacturing (APAM), has allowed us to tailor electronic devices with atomic precision, and thus has emerged as a way to explore new possibilities in Si electronics. In these applications, critical questions include where current flow is actually occurring in or near APAM structures as well as whether leakage currents are present. In general, detection and mapping of current flow in APAM structures are valuable diagnostic tools to obtain reliable devices in digital-enhanced applications. In this report, we used nitrogen-vacancy (NV) centers in diamond for wide-field magnetic imaging (with a few-mm field of view and micron-scale resolution) of magnetic fields from surface currents flowing in an APAM test device made of a P delta-doped layer on a Si substrate, a standard APAM witness material. We integrated a diamond having a surface NV ensemble with the device (patterned in two parallel mm-sized ribbons), then mapped the magnetic field from the DC current injected in the APAM device in a home-built NV wide-field microscope. The 2D magnetic field maps were used to reconstruct the surface current densities, allowing us to obtain information on current paths, device failures such as choke points where current flow is impeded, and current leakages outside the APAM-defined P-doped regions. Analysis on the current density reconstructed map showed a projected sensitivity of ~0.03 A m-1, corresponding to a smallest-detectable current in the 200 μm wide APAM ribbon of ~6 μA. These results demonstrate the failure analysis capability of NV wide-field magnetometry for APAM materials, opening the possibility to investigate other cutting-edge microelectronic devices.
Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.
Cryogenic plasma focused ion beam (PFIB) electron microscopy analysis is applied to visualizing ex situ (surface industrial) and in situ (subsurface geologic) carbonation products, to advance understanding of carbonation kinetics. Ex situ carbonation is investigated using NIST fly ash standard #2689 exposed to aqueous sodium bicarbonate solutions for brief periods of time. In situ carbonation pathways are investigated using volcanic flood basalt samples from Schaef et al. (2010) exposed to aqueous CO2 solutions by them. The fly ash reaction products at room temperature show small amounts of incipient carbonation, with calcite apparently forming via surface nucleation. Reaction products at 75° C show beginning stages of an iron carbonate phase, e.g., siderite or ankerite, common phases in subsurface carbon sequestration environments. This may suggest an alternative to calcite in carbonation low calcium-bearing fly ashes. Flood basalt carbonation reactions show distinct zonation with high calcium and calcium-magnesium bearing zones alternating with high iron-bearing zones. The calcium-magnesium zones are notable with occurrence of localized pore space. Oscillatory zoning in carbonate minerals is distinctly associated with far-from-equilibrium conditions where local chemical environments fluctuate via a coupling of reaction with transport. The high porosity zones may reflect a precursor phase (e.g., aragonite) with higher molar volume that then “ripens” to the high-Mg calcite phase-plus-porosity. These observations reveal that carbonation can proceed with evolving local chemical environments, formation and disappearance of metastable phases, and evolving reactive surface areas. Together this work shows that future application of cryo-PFIB in carbonation studies would provide advanced understanding of kinetic mechanisms for optimizing industrial-scale and commercial-scale applications.
Laser-induced photoemission of electrons offers opportunities to trigger and control plasmas and discharges. However, the underlying mechanisms are not sufficiently characterized to be fully utilized. Photoemission is highly nonlinear, achieved through multiphoton absorption, above threshold ionization, photo-assisted tunneling, etc., where the dominant process depends on the work function of the material, photon energy and associated fields, surface heating, background fields, etc. To characterize the effects of photoemission on breakdown, breakdown experiments were performed and interpreted using a 0D plasma discharge circuit model and quantum model of photoemission.
Single particle aerosol mass spectrometry (SPAMS), an analytical technique for measuring the size and composition of individual micron-scale particles, is capable of analyzing atmospheric pollutants and bioaerosols much more efficiently and with more detail than conventional methods which require the collection of particles onto filters for analysis in the laboratory. Despite SPAMS’ demonstrated capabilities, the primary mechanisms of ionization are not fully understood, which creates challenges in optimizing and interpreting SPAMS signals. In this paper, we present a well-stirred reactor model for the reactions involved with the laser-induced vaporization and ionization of an individual particle. The SPAMS conditions modeled in this paper include a 248 nm laser which is pulsed for 8 ns to vaporize and ionize each particle in vacuum. The ionization of 1 μm, spherical Al particles was studied by approximating them with a 0-dimensional plasma chemistry model. The primary mechanism of absorption of the 248 nm photons was pressure-broadened direct photoexcitation to Al(y2D). Atoms in this highly excited state then undergo superelastic collisions with electrons, heating the electrons and populating the lower energy excited states. We found that the primary ionization mechanism is electron impact ionization of various excited state Al atoms, especially Al(y2D). Because the gas expands rapidly into vacuum, its temperature decreases rapidly. The rate of three-body recombination (e- + e- + Al+ → Al + e-) increases at low temperature, and most of the electrons and ions produced recombine within several μs of the laser pulse. The importance of the direct photoexcitation indicates that the relative peak heights of different elements in SPAMS mass spectra may be sensitive to the available photoexcitation transitions. We also discuss the effects of laser intensity, particle diameter, and expansion dynamics.
Monitoring of cooling tower performance in a nuclear reactor facility is necessary to ensure safe operation; however, instrumentation for measuring performance characteristics can be difficult to install and may malfunction or break down over long duration experiments. This paper describes employing a thermodynamic approach to quantify cooling tower performance, the Merkel model, which requires only five parameters, namely, inlet water temperature, outlet water temperature, liquid mass flowrate, gas mass flowrate, and wet bulb temperature. Using this model, a general method to determine cooling tower operation for a nuclear reactor was developed in situations when neither the outlet water temperature nor gas mass flowrate are available, the former being a critical piece of information to bound the Merkel integral. Furthermore, when multiple cooling tower cells are used in parallel (as would be in the case of large-scale cooling operations), only the average outlet temperature of the cooling system is used as feedback for fan speed control, increasing the difficulty of obtaining the outlet water temperature for each cell. To address these shortcomings, this paper describes a method to obtain individual cell outlet water temperatures for mechanical forced-air cooling towers via parametric analysis and optimization. In this method, the outlet water temperature for an individual cooling tower cell is acquired as a function of the liquid-to-gas ratio (L/G). Leveraging the tight tolerance on the average outlet water temperature, an error function is generated to describe the deviation of the parameterized L/G to the highly controlled average outlet temperature. The method was able to determine the gas flowrate at rated conditions to be within 3.9% from that obtained from the manufacturer’s specification, while the average error for the four individual cooling cell outlet water temperatures were 1.6 °C, -0.5 °C, -1.0 °C, and 0.3 °C.
Stangebye, Sandra S.; Lei, Changhui L.; Kinghorn, Aubri K.; robertson, ian m.; Kacher, Josh K.; Hattar, Khalid M.
We report the dynamics of the gold–silicon eutectic reaction in limited dimensions were studied using in situ transmission electron microscopy and scanning transmission electron microscopy heating experiments. The phase transformation, viewed in both plan-view and cross-section of the film, occurs through a complex combination of dislocation and grain boundary motion and diffusion of silicon along gold grain boundaries, which results in a dramatic change in the microstructure of the film. The conversion observed in cross-section shows that the eutectic mixture forms at the Au–Si interface and proceeds into the Au film at a discontinuous growth rate. This complex process can lead to a variety of microstructures depending on sample geometry, heating temperature, and the ratio of gold to silicon which was found to have the largest impact on the eutectic microstructure. The eutectic morphology varied from dendrites to hollow rectangular structures to Au–Si eutectic agglomerates with increasing silicon to gold ratio.
Yarritu, Kevin A.; Hongkyu, Yoon H.; Roesler, Erika R.
Reliable climate predictions are important for making robust decisions in response to the changing climate. This project aims to reduce mis-modeling uncertainties arising from the representation of the land-atmosphere coupling in the Energy Exascale Earth System Model (E3SM) by using a machine learning approach. This approach will use an encoder-decoder architecture to represent the information that is developed in the land model and given to the atmosphere model. The simulated data will be taken from the E3SM simulation. However, the incorporation of observed data into the simulated dataset reduces mis-modeling uncertainties.
Many distributed applications implement complex data flows and need a flexible mechanism for routing data between producers and consumers. Recent advances in programmable network interface cards, or SmartNICs, represent an opportunity to offload data-flow tasks into the network fabric, thereby freeing the hosts to perform other work. System architects in this space face multiple questions about the best way to leverage SmartNICs as processing elements in data flows. In this paper, we advocate the use of Apache Arrow as a foundation for implementing data-flow tasks on SmartNICs. We report on our experiences adapting a partitioning algorithm for particle data to Apache Arrow and measure the on-card processing performance for the BlueField-2 SmartNIC. Our experiments confirm that the BlueField-2’s (de)compression hardware can have a significant impact on in-transit workflows where data must be unpacked, processed, and repacked.
Sanz-Matias, Ana S.; Roychoudhury, Subhayan x.; Feng, Xuefei F.; Yang, Feipeng Y.; Cheng, Kao.L.; Zavadil, Kevin R.; Guo, Jinghua G.; Prendergast, David P.
Given their natural abundance and thermodynamic stability, fluoride salts may appear as evolving components of electrochemical interfaces in Li-ion batteries and emergent multivalent ion cells. This is due to the practice of employing electrolytes with fluorine-containing species (salt, solvent, or additives) that electrochemically decompose and deposit on the electrodes. Operando X-ray absorption spectroscopy (XAS) can probe the electrode–electrolyte interface with a single-digit nanometer depth resolution and offers a wealth of insights into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure about the excited atoms. Here we explore fluorine K-edge XAS of mono- (Li, Na, and K) and di-valent (Mg, Ca, and Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the ΔSCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from the mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.
Dirac semimetals have attracted a great deal of current interests due to their potential applications in topological quantum computing, low-energy electronic devices, and single photon detection in the microwave frequency range. Herein are results from analyzing the low magnetic (B) field weak-antilocalization behaviors in a Dirac semimetal Cd3As2 thin flake device. At high temperatures, the phase coherence length lφ first increases with decreasing temperature (T) and follows a power law dependence of lφ ∝ T-0.4. Below ~3 K, lφ tends to saturate to a value of~180 nm. Another fitting parameter α, which is associated with independent transport channels, displays a logarithmic temperature dependence for T>3 K, but also tends to saturate below~3 K. The saturation value,~1.45, is very close to 1.5, indicating three independent electron transport channels, which we interpret as due to decoupling of both the top and bottom surfaces as well as the bulk. This result, to our knowledge, provides first evidence that the surfaces and bulk states can become decoupled in electronic transport in Dirac semimetal Cd3As2.
The International Atomic Energy Agency (IAEA) applies safeguards to nuclear facilities that are not operating, including those undergoing decommissioning, and the IAEA’s effort in this area is both considerable and increasing. Specifically, the IAEA Department of Safeguards’ Division of Concepts and Planning (SGCP-003: Safeguards Approaches) identified an R&D need to “Develop safeguards implementation guidelines for facilities under decommissioning and safeguards concepts for post-accident facilities under decommissioning”. Nuclear facilities undergoing decommissioning are not exempt from safeguards agreements between the IAEA and Host State, and, accordingly, the requirement for verification of no diversion of nuclear material and detection of undeclared activities at decommissioned facilities remain even after facility shutdown. However, the effort required to meet safeguards objectives diminishes as nuclear material and essential equipment are removed during the decommissioning process which shifts the emphasis from verification of ever-diminishing fissile or source material inventories to verification of changes in facility design and equipment operability.
The New Mexico Small Business Assistance Program (NMSBA) has once again paired with Optical Radio Communications Technology (ORC Tech). A New Mexico startup Limited Liability Company (LLC), with Sandia National Laboratories (SNL) Engineers at the Sensors and Textiles Innovatively Tailored for Complex, High-Efficiency Detection (STITCHED) laboratory, to aid in the development of an ultra-passive, portable, deployable wireless signal booster technology.
Artificial solid electrolyte interphases have provided a path to improved cycle life for high energy density, next-generation anodes like lithium metal. Although long cycle life is necessary for widespread implementation, understanding and mitigating the effects of aging and self-discharge are also required. In this report we investigate several coating materials and their role in calendar life aging of lithium. We find that the oxide coatings are electronically passivating whereas the LiF coating slows charge transfer kinetics. Furthermore, the Coulombic loss during self-discharge measurements improves with the oxide layers and worsens with the LiF layer. It is found that none of the coatings create a continuous conformal, electronically passivating layer on top of the deposited lithium nor are they likely to distribute evenly through a porous deposit, suggesting that none of the materials are acting as an artificial solid electrolyte interphase. Instead, they likely alter performance through modulating lithium nucleation and growth.
This project matured a new understanding (a “modern synthesis”) of the structure and evolution of science and technology. It created an understanding and framework for how Sandia National Labs, the Department of Energy, and the nation, might improve their research productivity, with significant ramifications on national security and economic competitiveness.
A zero-dimensional magnetic implosion model with a coupled equivalent circuit for the description of an imploding nested wire array or gas puff is presented. Circuit model results have been compared with data from imploding stainless steel wire arrays, and good agreement has been found. The total energy coupled to the load, , has been applied to a simple semi-analytic K-shell yield model, and excellent agreement with previously reported K-shell yields across all wire array and gas puff platforms is seen. Trade space studies in implosion radius and mass have found that most platforms operate near the predicted maximum yield. In some cases, the K-shell yield may be increased by increasing the mass or radius of the imploding array or gas puff.
The solution processability of ionogel solid electrolytes has recently garnered attention in the Li-ion battery community as a means to address the interface and fabrication issues commonly associated with most solid electrolytes. However, the trapped ionic liquid (ILE) component has hindered the electrochemical performance. In this report we present a process to tune the properties by replacing the ILE in a silica-based ionogel after fabrication with a liquid component befitting the desired application. Electrochemical cycling under various conditions showcases gels containing different liquid components incorporated into LiFePO4 (LFP)/gel/Li cells: high power (455 W kg–1 at a 1 C discharge) systems using carbonates, low temperatures (-40 °C) using ethers, or high temperatures (100 °C) using ionic liquids. Fabrication of additive-manufactured cells utilizing the exchanged carbonate-based system is demonstrated in a planar LFP/Li4Ti5O12 (LTO) system, where a marked improvement over an ionogel is found in terms of rate capability, capacity, and cycle stability (118 vs 41 mA h g–1 at C/4). This process represents a promising route to create a separator-less cell, potentially in complex architectures, where the electrolyte properties can be facilely tuned to meet the required conditions for a wide range of battery chemistries while maintaining a uniform electrolyte access throughout cast electrodes.
For 2D-temperature monitoring applications, a variant of EIT (Electrical Impedance Tomography) is evaluated computationally in this work. Literature examples of poor sensor performance in the center of the 2D domains away from the side electrodes motivated this study which seeks to overcome some of the previously noted shortcomings. In particular, the use of ‘sensing skins’ with novel tailored baseline conductivities were examined using the EIDORS package for EIT. It was found that the best approach for detecting a hot spot depends on several factors such as the current injection (stimulation) patterns, the measurement patterns, and the reconstruction algorithms. For a well-performing combination of these factors, tailored baseline conductivities were assessed and compared to the baseline uniform conductivity. It was discovered that for some EIT applications, a tailored distribution needs to be smooth and that sudden changes in the conductivity gradients should be avoided. Still, the benefits in terms of improved EIT performance were small for conditions for which the EIT measurements had been ‘optimized’ for the uniform baseline case. Within the limited scope of this study, only two specific cases showed benefits from tailored distributions. For one case, a smooth tailored distribution with increased baseline conductivity in the center provided a better separation of two centrally located hot spots. For another case, a smooth tailored distribution with reduced conductivity in the center provided better estimates of the magnitudes of two hot spots near the center of the sensing skin.
A zero-dimensional magnetic implosion model with a coupled equivalent circuit for the description of an imploding nested wire array or gas puff is presented. Circuit model results have been compared with data from imploding stainless steel wire arrays, and good agreement has been found. The total energy coupled to the load, E j × B, has been applied to a simple semi-analytic K-shell yield model, and excellent agreement with previously reported K-shell yields across all wire array and gas puff platforms is seen. Trade space studies in implosion radius and mass have found that most platforms operate near the predicted maximum yield. In some cases, the K-shell yield may be increased by increasing the mass or radius of the imploding array or gas puff.
This research effort examined the application of Nafion polymers in alcohol solvents as an anti-ice surface coating, as a mixture with hydrophilic polymers and freezing point depressant salt systems. Co-soluble systems of Nafion, polymer and salt were applied using dip coating methods to create smooth films for frost observation over a Peltier plate thermal system in ambient laboratory conditions. Cryo-DSC was applied to examine freezing events of the Nafion-surfactant mixtures, but the sensitivity of the measurement was insufficient to determine frost behavior. Collaborations with the Fog Chamber at Sandia-Albuquerque, and in environmental SAXS measurements with CINT-LANL were requested but were not able to be performed under the research duration. Since experimental characterization of these factors is difficult to achieve directly, computational modeling was used to guide the scientific basis for property improvement. Computational modeling was performed to improve understanding of the dynamic association between ionomer side groups and added molecules and deicing salts. The polyacrylic acid in water system was identified at the start of the project as a relevant system for exploring the effect of varying counterions on the properties of fully deprotonated polyacrylic acid (PAA) in the presence of water. Simulations were modeled with four different counterions, two monovalent counterions (K+ and Na+) and two divalent counterions (Ca2+ and Mg2+). The wt% of PAA in these systems was varied from ~10 to 80 wt% PAA for temperatures from 250K to 400K. In the second set of simulations, the interpenetration of water into a dry PAA film was studied for Na+ or Ca2+ counterions for temperatures between 300K and 400K. The result of this project is a sprayable Nafion film composite which resists ice nucleation at -20 °C for periods of greater than three hours. It is composed of Nafion polymer, hydrophilic polyethylene oxide polymer and CaCl2 anti-ice crosslinker. Durability and field performance properties remain to be determined.
This manual describes the use of the Xyce™ Parallel Electronic Simulator. Xyce™ has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. (2) A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. (3) Device models that are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only). (4) Object-oriented code design and implementation using modern coding practices. Xyce™ is a parallel code in the most general sense of the phrase—a message passing parallel implementation—which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel eficiency is achieved as the number of processors grows.