Publications

Systems engineering today faces a wide array of challenges, ranging from new operational environments to disruptive technological — necessitating approaches to improve research and development (R&D) efforts. Yet, emphasizing the Aristotelian argument that the “whole is greater than the sum of its parts” seems to offer a conceptual foundation creating new R&D solutions. Invoking systems theoretic concepts of emergence and hierarchy and analytic characteristics of traceability, rigor, and comprehensiveness is potentially beneficial for guiding R&D strategy and development to bridge the gap between theoretical problem spaces and engineering-based solutions. In response, this article describes systems–theoretic process analysis (STPA) as an example of one such approach to aid in early-systems R&D discussions. STPA—a ‘top-down’ process that abstracts real complex system operations into hierarchical control structures, functional control loops, and control actions—uses control loop logic to analyze how control actions (designed for desired system behaviors) may become violated and drive the complex system toward states of higher risk. By analyzing how needed controls are not provided (or out of sequence or stopped too soon) and unneeded controls are provided (or engaged too long), STPA can help early-system R&D discussions by exploring how requirements and desired actions interact to either mitigate or potentially increase states of risk that can lead to unacceptable losses. This article will demonstrate STPA's benefit for early-system R&D strategy and development discussion by describing such diverse use cases as cyber security, nuclear fuel transportation, and US electric grid performance. Together, the traceability, rigor, and comprehensiveness of STPA serve as useful tools for improving R&D strategy and development discussions. In conclusion, leveraging STPA as well as related systems engineering techniques can be helpful in early R&D planning and strategy development to better triangulate deeper theoretical meaning or evaluate empirical results to better inform systems engineering solutions.

Magnetic properties of more than 20 Cantor alloy samples of varying composition were investigated over a temperature range of 5 K to 300 K and in fields of up to 70 kOe using magnetometry and muon spin relaxation. Two transitions are identified: a spin-glass-like transition that appears between 55K and 190K, depending on composition, and a ferrimagnetic transition that occurs at approximately 43K in multiple samples with widely varying compositions. The magnetic signatures at 43K are remarkably insensitive to chemical composition. A modified Curie-Weiss model was used to fit the susceptibility data and to extract the net effective magnetic moment for each sample. The resulting values for the net effective moment were either diminished with increasing Cr or Mn concentrations or enhanced with decreasing Fe, Co, or Ni concentrations. Beyond a sufficiently large effective moment, the magnetic ground state transitions from ferrimagnetism to ferromagnetism. The effective magnetic moments, together with the corresponding compositions, are used in a global linear regression analysis to extract element-specific effective magnetic moments, which are compared to the values obtained by ab initio based density functional theory calculations. Finally, these moments provide the information necessary to controllably tune the magnetic properties of Cantor alloy variants.

A new capability for modeling graded density reactive flow materials in the shock physics hydrocode, CTH, is demonstrated here. Previously, materials could be inserted in CTH with graded material properties, but the sensitivity of the material was not adjusted based on these properties. Of particular interest are materials that are graded in density, sometimes due to pressing or other assembly operations. The sensitivity of explosives to both density and temperature has been well demonstrated in the literature, but to-date the material parameters for use in a simulation were fit to a single condition and applied to the entire material, or the material had to be inserted in sections and each section assigned a condition. The reactive flow model xHVRB has been extended to shift explosive sensitivity with initial density, so that sensitivity is also graded in the material. This capability is demonstrated for use in three examples. The first models detonation transfer in a graded density pellet of HNS, the second is a shaped charge with density gradients in the explosive, and the third is an explosively formed projectile.

Maximizing the production of heterologous biomolecules is a complex problem that can be addressed with a systems-level understanding of cellular metabolism and regulation. Specifically, growth-coupling approaches can increase product titers and yields and also enhance production rates. However, implementing these methods for non-canonical carbon streams is challenging due to gaps in metabolic models. Over four design-build-test-learn cycles, we rewire Pseudomonas putida KT2440 for growth-coupled production of indigoidine from para-coumarate. We explore 4,114 potential growth-coupling solutions and refine one design through laboratory evolution and ensemble data-driven methods. The final growth-coupled strain produces 7.3 g/L indigoidine at 77% maximum theoretical yield in para-coumarate minimal medium. The iterative use of growth-coupling designs and functional genomics with experimental validation was highly effective and agnostic to specific hosts, carbon streams, and final products and thus generalizable across many systems.

For reactive burn models in hydrocodes, an equilibrium closure assumption is typically made between the unreacted and product equations of state. In the CTH [1] (not an acronym) hydrocode the assumption of density and temperature equilibrium is made by default, while other codes make a pressure and temperature equilibrium assumption. The main reason for this difference is the computational efficiency in making the density and temperature assumption over the pressure and temperature one. With fitting to data, both assumptions can accurately predict reactive flow response using the various models, but the model parameters from one code cannot necessarily be used directly in a different code with a different closure assumption. A new framework is intro-duced in CTH to allow this assumption to be changed independently for each reactive material. Comparisons of the response and computational cost of the History Variable Reactive Burn (HVRB) reactive flow model with the different equilibrium assumptions are presented.

Process variations within Field Programmable Gate Arrays (FPGAs) provide a rich source of entropy and are therefore well-suited for the implementation of Physical Unclonable Functions (PUFs). However, careful considerations must be given to the design of the PUF architecture as a means of avoiding undesirable localized bias effects that adversely impact randomness, an important statistical quality characteristic of a PUF. Here in this paper, we investigate a ring-oscillator (RO) PUF that leverages localized entropy from individual look-up table (LUT) primitives. A novel RO construction is presented that enables the individual paths through the LUT primitive to be measured and isolated at high precision, and an analysis is presented that demonstrates significant levels of localized design bias. The analysis demonstrates that delay-based PUFs that utilize LUTs as a source of entropy should avoid using FPGA primitives that are localized to specific regions of the FPGA, and instead, a more robust PUF architecture can be constructed by distributing path delay components over a wider region of the FPGA fabric. Compact RO PUF architectures that utilize multiple configurations within a small group of LUTs are particularly susceptible to these types of design-level bias effects. The analysis is carried out on data collected from a set of identically designed, hard macro instantiations of the RO implemented on 30 copies of a Zynq 7010 SoC.

Organic co-crystals have emerged as a promising class of semiconductors for next-generation optoelectronic devices due to their unique photophysical properties. This paper presents a joint experimental-theoretical study comparing the crystal structure, spectroscopy, and electronic structure of two charge transfer co-crystals. Reported herein is a novel co-crystal Npe:TCNQ, formed from 4-(1-naphthylvinyl)pyridine (Npe) and 7,7,8,8-tetracyanoquinodimethane (TCNQ) via molecular self-assembly. This work also presents a revised study of the co-crystal composed of Npe and 1,2,4,5-tetracyanobenzene (TCNB) molecules, Npe:TCNB, herein reported with a higher-symmetry (monoclinic) crystal structure than previously published. Npe:TCNB and Npe:TCNQ dimer clusters are used as theoretical model systems for the co-crystals; the geometries of the dimers are compared to geometries of the extended solids, which are computed with periodic boundary conditions density functional theory. UV-Vis absorption spectra of the dimers are computed with time-dependent density functional theory and compared to experimental UV-Vis diffuse reflectance spectra. Both Npe:TCNB and Npe:TCNQ are found to exhibit neutral character in the S0 state and ionic character in the S1 state. The high degree of charge transfer in the S1 state of both Npe:TCNB and Npe:TCNQ is rationalized by analyzing the changes in orbital localization associated with the S1 transitions.

Abstract not provided.

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10⁶ times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents fiscal year 2023 (FY23) updates on the Advanced Drying Cycle Simulator (ADCS). This apparatus was built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. These waterproof heaters are the next generation design to heater rods developed and tested at Sandia National Laboratories in FY20. In FY23, a series of four simulated commercial drying tests was completed. This report presents the temperature and pressure histories of the drying tests as well as axial temperature profiles that can be compared to data from the Electric Power Research Institute (EPRI) High Burnup Demonstration TN-32B cask. Water content measurements and dew point calculations from a Hiden Analytical HPR-30 mass spectrometer are also presented in this report. Due to familiarization with this first-of-a-kind system, refinements to equipment calibration and test procedures have been identified to better match commercial drying cycles for future simulated tests. However, the presented data demonstrate the successful construction and operation of a viable research platform for quantifying residual water content closely approaching that expected in dry storage canisters during commercial drying procedures.

Work evaluating spent nuclear fuel (SNF) dry storage canister surface environments and canister corrosion progressed significantly in FY23, with the goal of developing a scientific understanding of the processes controlling initiation and growth of stress corrosion cracking (SCC) cracks in stainless steel canisters in relevant storage environments. The results of the work performed at Sandia National Laboratories (SNL) will guide future work and will contribute to the development of better tools for predicting potential canister penetration by SCC.

In this study, we demonstrate the ability of polarity inversion of sputtered aluminum scandium nitride thin films through post-fabrication processes with domain widths as small as 220 nm at a periodicity of 440 nm. An approach using photo- and electron-beam lithography to generate sub-quarter micrometer feature size with adjustable duty cycle through a lift-off process is presented. The film with a coercive field Ec+ of 5.35 MV/cm was exercised first with a 1 kHz triangular double bipolar wave and ultimately poled with a 0.5 kHz double monopolar wave using a Radiant Precision Premier II tester. The metal polar (M-polar) and nitrogen polar (N-polar) domains were identified and characterized through potassium hydroxide wet etching as well as piezoresponse force microscopy (PFM). Well-distinguished boundaries between the oppositely polarized domain regions were confirmed through the phase diagram of the PFM results. The relationship between the electrode width, poling voltage, and domain growth was experimentally studied and statistically analyzed, where 7.96 nm/V domain width broadening vs escalating poling voltage was observed. This method produces extremely high domain spatial resolution in III-nitride materials via poling and is transferable to a CMOS-compatible photolithography process. The spatial resolution of the periodically poled Al0.68Sc0.32N is suitable for second-harmonic generation of deep ultraviolet through quasi-phase-matching and RF MEMS operating in the X-Band spectrum.

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. Comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

Tritium permeability in zirconium-based tritium getter critically impacts tritium storage and environmental safety during operation of tritium-producing burnable absorber rods (TPBARs). Previous experiments indicated that during irradiation operation, the hydrogen equilibrium pressured is increased. Further experimental and modeling studies suggested that the enhanced tritium release observed for reactor scale assemblies might be related to a thermal diffusion known as the Soret effect. A direct measurement of the Soret factor, however, has not been performed. To improve TPBAR and other nuclear applications, here we have applied two non-equilibrium molecular dynamics methods to study thermal diffusion of hydrogen isotopes in low-concentration zirconium hydrides. One of the methods produces sufficiently converged results to distinguish crystal orientation, isotope type, and concentration effects. In conclusion, with this method, crystal orientation, isotope type, and concentration effects are discussed.

Local crystallographic features negatively affect quantum spin defects by changing the local electrostatic environment, often resulting in degraded or varied qubit optical and coherence properties. Few tools exist that enable the deterministic synthesis and study of such intricate systems on the nano-scale, making defect-to-defect strain environment quantification difficult. In this paper, we highlight state-of-the-art capabilities from the U.S. Department of Energy’s Nanoscale Science Research Centers that directly address these shortcomings. Specifically, we demonstrate how complementary capabilities of nano-implantation and nano-diffraction can be used to demonstrate the quantum relevant, spatially deterministic creation of neutral divacancy centers in 4H silicon carbide, while investigating and characterizing these systems on the ≤ 25 nm scale with strain sensitivities on the order of 1 × 10 − 6 , relevant to defect formation dynamics. This work lays the foundation for ongoing studies into the dynamics and deterministic formation of low strain homogeneous quantum relevant spin defects in the solid state.

Microneedle sensors could enable minimally-invasive, continuous molecular monitoring – informing on disease status and treatment in real-time. Wearable sensors for pharmaceuticals, for example, would create opportunities for treatments personalized to individual pharmacokinetics. Here, we demonstrate a commercial-off-the-shelf (COTS) approach for microneedle sensing using an electrochemical aptamer-based sensor that detects the high-toxicity antibiotic, vancomycin. Wearable monitoring of vancomycin could improve patient care by allowing targeted drug dosing within its narrow clinical window of safety and efficacy. To produce sensors, we miniaturize the electrochemical aptamer-based sensors to a microelectrode format, and embed them within stainless steel microneedles (sourced from commercial insulin pen needles). The microneedle sensors achieve quantitative measurements in body-temperature undiluted blood. Further, the sensors effectively maintain electrochemical signal within porcine skin. This COTS approach requires no cleanroom fabrication or specialized equipment, and produces individually-addressable, sterilizable microneedle sensors capable of easily penetrating the skin. In the future, this approach could be adapted for multiplexed detection, enabling real-time monitoring of a range of biomarkers.

The Sibling Pin test campaign is a Department of Energy (DOE) research activity within the Spent Fuel and Waste Science and Technology (SFWST) program that is tasked with characterization of high burnup (HBU) fuel in support of the High Burnup Spent Fuel Data Project. Of the 25 fuel rods in the Sibling Pin inventory, approximately 9 rod lengths have been consumed during the first phase (Phase I) of the test campaign leaving approximately 16 rod lengths for the second phase (Phase II) of testing. This plan outlines the Phase II testing and the motivations for performing these tests. Priorities for Phase II testing are based on previously identified knowledge gaps, lessons-learned from Phase I work, the original objectives of the High Burnup Spent Fuel Data Project and the Sibling Pin test campaign, and input from external stakeholders. The priorities for Phase II testing are to obtain data to characterize the effects of annealing on cladding mechanical properties and fuel rod performance, to quantify the creep behavior of cladding materials and fuel rods and the effects of creep deformations on the performance of cladding and fuel rods, and to gather data to support the final closure of the hydride reorientation and radial hydride induced embrittlement gap for HBU fuel rods.

Translating the surging interest in neuromorphic electronic components, such as those based on nonlinearities near Mott transitions, into large-scale commercial deployment faces steep challenges in the current lack of means to identify and design key material parameters. These issues are exemplified by the difficulties in connecting measurable material properties to device behavior via circuit element models. Here, the principle of local activity is used to build a model of VO2/SiN Mott threshold switches by sequentially accounting for constraints from a minimal set of quasistatic and dynamic electrical and high-spatial-resolution thermal data obtained via in situ thermoreflectance mapping. By combining independent data sets for devices with varying dimensions, the model is distilled to measurable material properties, and device scaling laws are established. The model can accurately predict electrical and thermal conductivities and capacitances and locally active dynamics (especially persistent spiking self-oscillations). The systematic procedure by which this model is developed has been a missing link in predictively connecting neuromorphic device behavior with their underlying material properties, and should enable rapid screening of material candidates before employing expensive manufacturing processes and testing procedures.

Motivated by recent experimental results we calculate from first-principles the lifetime of low-energy quasiparticles in bilayer graphene (BLG). Here, we take into account the scattering rate arising from electron-electron interactions within the GW approximation for the electron self-energy and consider several p-type doping levels ranging from 0 to ρ ≈ 2.4 × 10¹² holes/cm². In the undoped case we find that the average inverse lifetime scales linearly with energy away from the charge neutrality point, with values in good agreement with experiments. The decay rate is approximately three times larger than in monolayer graphene, a consequence of the enhanced screening in BLG. In the doped case, the dependence of the inverse lifetime on quasiparticle energy acquires a non-linear component due to the opening of an additional decay channel mediated by acoustic plasmons.

Here, we study the impact of deep-level defects on trap-assisted Auger–Meitner recombination in c-plane InGaN/GaN LEDs using a small-signal electroluminescence (SSEL) method and deep-level optical spectroscopy (DLOS). Carrier dynamics information, including carrier lifetime, recombination rate, and carrier density, is obtained from SSEL, while DLOS is used to obtain the deep-level defect density. Through fitting the nonradiative recombination rates of wafers with different deep-level defect densities, we obtain the Shockley–Read–Hall (SRH) and trap-assisted Auger–Meitner recombination (TAAR) coefficients. We show that defect-related nonradiative recombination, including both SRH and TAAR, accounts for a relatively small fraction of the total nonradiative recombination, which is dominated by intrinsic Auger–Meitner recombination. The interplay between carrier localization and Coulomb enhancement has a different impact on radiative and intrinsic Auger–Meitner recombination. Evidence is presented that the imbalance between the change of radiative and intrinsic Auger–Meitner recombination is the primary cause of the efficiency droop at high carrier densities in the samples studied.

We study the impact of deep-level defects on trap-assisted Auger-Meitner recombination in c-plane InGaN/GaN LEDs using a small-signal electroluminescence (SSEL) method and deep-level optical spectroscopy (DLOS). Carrier dynamics information, including carrier lifetime, recombination rate, and carrier density, is obtained from SSEL, while DLOS is used to obtain the deep-level defect density. Through fitting the nonradiative recombination rates of wafers with different deep-level defect densities, we obtain the Shockley-Read-Hall (SRH) and trap-assisted Auger-Meitner recombination (TAAR) coefficients. We show that defect-related nonradiative recombination, including both SRH and TAAR, accounts for a relatively small fraction of the total nonradiative recombination, which is dominated by intrinsic Auger-Meitner recombination. The interplay between carrier localization and Coulomb enhancement has a different impact on radiative and intrinsic Auger-Meitner recombination. Evidence is presented that the imbalance between the change of radiative and intrinsic Auger-Meitner recombination is the primary cause of the efficiency droop at high carrier densities in the samples studied.

Cast Monel alloys are used in applications requiring a combination of good mechanical properties and excellent resistance to corrosion. Despite prevalent industrial use, relatively few studies have been conducted to investigate the relationships between composition, solidification behavior, and microstructure. Given that these alloys are used in the cast and welded conditions, these factors have a significant influence over the material properties. Here, in this work, microstructural characterization, electron probe microanalysis, X-ray diffraction, and differential scanning calorimetry were used to study how changes in Si and Nb concentrations affected the solidification path and microstructure of Monel alloys. It was found that increasing Nb concentration stabilized higher amounts of MC carbides and suppressed graphite formation during solidification. It was also found that the high nominal concentration and segregation of Si to the liquid led to the formation of Ni₃₁Si₁₂ and other silicides via terminal eutectic reactions at the end of solidification. A pseudo-binary solidification diagram was constructed using experimental data and was applied to predict the mass fraction of solidified eutectic as a function of composition. The modeled microstructures were found to be in good agreement with experimentally measured phase fractions.

Here, we perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

Alkene ozonolysis generates short-lived Criegee intermediates that are a significant source of hydroxyl (OH) radicals. This study demonstrates that roaming of the separating OH radicals can yield alternate hydroxycarbonyl products, thereby reducing the OH yield. Specifically, hydroxybutanone has been detected as a stable product arising from roaming in the unimolecular decay of the methyl-ethyl-substituted Criegee intermediate (MECI) under thermal flow cell conditions. The dynamical features of this novel multistage dissociation plus a roaming unimolecular decay process have also been examined with ab initio kinetics calculations. Experimentally, hydroxybutanone isomers are distinguished from the isomeric MECI by their higher ionization threshold and distinctive photoionization spectra. Moreover, the exponential rise of the hydroxybutanone kinetic time profile matches that for the unimolecular decay of MECI. A weaker methyl vinyl ketone (MVK) photoionization signal is also attributed to OH roaming. Complementary multireference electronic structure calculations have been utilized to map the unimolecular decay pathways for MECI, starting with 1,4 H atom transfer from a methyl or methylene group to the terminal oxygen, followed by roaming of the separating OH and butanonyl radicals in the long-range region of the potential. Roaming via reorientation and the addition of OH to the vinyl group of butanonyl is shown to yield hydroxybutanone, and subsequent C-O elongation and H-transfer can lead to MVK. A comprehensive theoretical kinetic analysis has been conducted to evaluate rate constants and branching yields (ca. 10-11%) for thermal unimolecular decay of MECI to conventional and roaming products under laboratory and atmospheric conditions, consistent with the estimated experimental yield (ca. 7%).

As a part of NASA's efforts in space, options are being examined for an Artemis moon base project to be deployed. This project requires a system of interconnected, but separate, DC microgrids for habitation, mining, and fuel processing. This in-place use of power resources is called in-situ resource utilization (ISRU). These microgrids are to be separated by 9-12 km and each contains a photovoltaic (PV) source, energy storage systems (ESS), and a variety of loads, separated by level of criticality in operation. The separate microgrids need to be able to transfer power between themselves in cases where there are generation shortfall, faults, or other failures in order to keep more critical loads running and ensure safety of personnel and the success of mission goals. In this work, a 2 grid microgrid system is analyzed involving a habitation unit and a mining unit separated by a tie line. A set of optimal controls that has been developed, including power flow controls on the tie line, dispatch of PV generation, and dispatch of non-critical loads, is analyzed, and validated in hardware on the Secure Scalable Microgrid Testbed (SSMTB). This testbed includes hardware emulators for a variety of energy sources, energy storage devices, pulsed loads, and other loads.

This report describes specific activities in the Fiscal Year (FY) 2023 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package funded by the Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), Office of Spent Fuel and Waste Disposition (SFWD).

Experimental studies and ab initio quantum chemistry calculations were combined to investigate the process by which a Fenton reaction breaks down polystyrene sulfonate. The experimental results show that both molecular weight reduction and loss of aromaticity occur nearly simultaneously, a finding that is supported by the calculations. The results show that more than half of the material is broken down to low molecular weight compounds (< 500 g/mol) with two molar equivalents of H2O2 per styrene monomer. The calculations provide insights into the reaction pathways and indicate that at least two hydroxyl radicals are required to cleave backbone C–C bonds or to eliminate aromaticity. The calculations also show that, of the aromatic carbons, hydroxyl radical is most likely to add to the carbon bonded to sulfur. This finding explains the loss of hydrogen sulfite anion early in the process and also the efficient reduction of Fe(III) to Fe(II) through semiquinone formation. Taken together the experimental and computational results indicate that the reaction is very efficient and that very little H2O2 is lost to unproductive reactions. This high efficiency is attributed to the close association of Fe atoms with the sulfonate group such that hydroxyl radicals are generated near the polymer chains.

Abstract not provided.

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Reτ=505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

Abstract not provided.

Update to prior 5.14 user manual. I think updates are minor and mostly in the Johnson-cook section. In there those updates are more writing and less on technical changes.

Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. In this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.

Calibrated measurements of lightning optical emissions are critical for both quantifying the impacts of lightning in our atmosphere and devising detection instruments with sufficient dynamic range capable of yielding close to 100% detection efficiency. However, to date, there is only a limited number of investigations that have attempted to take such calibrated measurements. In this work, we report the power radiated by lightning in both visible and infrared bands, assuming isotropic emission, and accounting for atmospheric absorption. More precisely, we report peak radiated power and total radiated energy in the combined visible plus near-infrared range (VNIR, 0.34–1.1 μm), around the Hα line (652–667 nm), and for the 2–2.5 μm infrared band. The estimated peak power and total energy radiated by negative cloud-to-ground return strokes in the VNIR range is 130 MW and 20 kJ, respectively. Additionally, we detected peak radiated powers of 12 and 0.19 MW in the Hα and infrared bands, respectively. We cross-reference the optical data set with peak current reported by a lightning detection network. The resulting trend is that optical power emitted around the Hα line scales with peak return stroke current according to a power law with exponent equal to 1.25. This trend, which should be approximately true across the entire visible spectrum, can be attributed to the plasma negative differential resistance of the lightning return stroke channel. We conclude by discussing the challenges in performing calibrated measurements of lightning optical power in different bands and comparing the results with previously-collected data with different experimental setups, observation conditions, and calibration methods.

We propose a novel statistical inference methodology for multiway count data that is corrupted by false zeros that are indistinguishable from true zero counts. Our approach consists of zero-truncating the Poisson distribution to neglect all zero values. This simple truncated approach dispenses with the need to distinguish between true and false zero counts and reduces the amount of data to be processed. Inference is accomplished via tensor completion that imposes low-rank tensor structure on the Poisson parameter space. Our main result shows that an N-way rank-R parametric tensor M ∈ (0, ∞)I×.....×I generating Poisson observations can be accurately estimated by zero-truncated Poisson regression from approximately IR2 log22(I) non-zero counts under the nonnegative canonical polyadic decomposition. Our result also quantifies the error made by zero-truncating the Poisson distribution when the parameter is uniformly bounded from below. Therefore, under a low-rank multiparameter model, we propose an implementable approach guaranteed to achieve accurate regression in under-determined scenarios with substantial corruption by false zeros. Several numerical experiments are presented to explore the theoretical results.

Previous studies have shown that selective laser flash sintering (SLFS) can be initiated in dielectrics that exhibit ionic or electronic conduction at high temperature. These materials required high laser powers to reach the temperatures where electrical conduction is sufficient to initiate SLFS. In this study, SLFS in lanthanum chromite (LC), an intrinsic electronic conductor with high conductivity, and lanthanum strontium chromite (LSC), which is doped to further increase electronic conductivity, were investigated with a focus on understanding the initiation mechanisms. Results show that the initiation of SLFS in LC and LSC occurs when electronic charge carriers are activated and flow to the electrode where the current is measured. A combination of carriers produced at the electrode, temperature-activated intrinsic charge carriers, and extrinsic charge carriers present in LSC due to doping are responsible for the facile initiation of SLFS.

Abstract not provided.

Abstract not provided.

Achieving brain-like efficiency in computing requires a co-design between the development of neural algorithms, brain-inspired circuit design, and careful consideration of how to use emerging devices. The recognition that leveraging device-level noise as a source of controlled stochasticity represents an exciting prospect of achieving brain-like capabilities in probabilistic neural algorithms, but the reality of integrating stochastic devices with deterministic devices in an already-challenging neuromorphic circuit design process is formidable. Here, we explore how the brain combines different signaling modalities into its neural circuits as well as consider the implications of more tightly integrated stochastic, analog, and digital circuits. Further, by acknowledging that a fully CMOS implementation is the appropriate baseline, we conclude that if mixing modalities is going to be successful for neuromorphic computing, it will be critical that device choices consider strengths and limitations at the overall circuit level.

Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the combined-field integral equation. We demonstrate the effectiveness of these approaches for cases with and without coding errors.

Global climate change has prompted many national plans for rapid emissions reductions. For example, the United States recently committed to transitioning to 100% carbon-free electricity by 2035 and net-zero emissions economy-wide by 2050. Parallel to conversations surrounding emissions reductions is the call for energy justice, or the demand for more equitable distribution of energy-related burdens and benefits among communities. To date, energy justice has evolved as a mostly academic conversation, which may limit its utility to praxis. In response, we offer an interdisciplinary framework that aims to organize existing knowledge and lessons learned from energy development. Specifically, we developed the Meaningful Marine Renewable Energy (MRE) Development Framework and conducted a literature review using MRE as a case study. MRE was chosen because it is a nascent renewable energy technology in the US with projects mostly in demonstration stages and no commercial deployment, making it a useful case study to apply lessons learned from other energy sectors and other countries. Discussion of current resources being developed among the MRE community and their implications for furthering energy justice priorities are also explored. We conclude the review with a compiled list of questions meant to support stakeholders in translating theoretical concepts of Meaningful MRE Development to practice. Although the Meaningful MRE framework was developed using MRE as a use case, our interdisciplinary theoretical framework can be applied beyond MRE to other sustainable and renewable energy projects.

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

Abstract not provided.

This study presents a deep learning based methodology for both remote sensing and design of acoustic scatterers. The ability to determine the shape of a scatterer, either in the context of material design or sensing, plays a critical role in many practical engineering problems. This class of inverse problems is extremely challenging due to their high-dimensional, nonlinear, and ill-posed nature. To overcome these technical hurdles, we introduce a geometric regularization approach for deep neural networks (DNN) based on non-uniform rational B-splines (NURBS) and capable of predicting complex 2D scatterer geometries in a parsimonious dimensional representation. Then, this geometric regularization is combined with physics-embedded learning and integrated within a robust convolutional autoencoder (CAE) architecture to accurately predict the shape of 2D scatterers in the context of identification and inverse design problems. An extensive numerical study is presented in order to showcase the remarkable ability of this approach to handle complex scatterer geometries while generating physically-consistent acoustic fields. The study also assesses and contrasts the role played by the (weakly) embedded physics in the convergence of the DNN predictions to a physically consistent inverse design.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material's mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. To construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer's glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

Salt offers an optimal medium for the permanent isolation of heat-producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal-hydrological-mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high-resolution fiber-optic-based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.

Abstract not provided.

Deflection and stress calculated from an experimentally validated, high-fidelity finite element model (FEM) of a photovoltaic module experiencing mechanical load was compared to results from a simplified FEM treating the module laminate as a homogenized composite using a rule of mixtures approach, and further compared to analytical calculations treating the module as a Kirchoff-Love flat plate. The goal of this study was to determine the error incurred by analyzing module mechanics with varying levels of simplification, since resolving the aspect ratios of a module is computationally expensive. Homogenized FEMs were found to underpredict peak deflection under a 1.0 kPa load by between 13 and 19% for lower and upper bound application of the rule of mixtures. However, module shape was captured, implying that a useful replication of a resolved model could be achieved with a reduced, calibrated material stiffness. Homogenized stress results captured glass layer tensile stress components to within 46 to 52% at a sample location of interest, though agreement was poor through the remainder of the laminate due to the lack of material resolution. For plate theory, deflection was overpredicted by 45 to 67% for upper and lower bound homogenizations, and frame-adjacent module shapes were not adequately replicated. Stress results mirrored FEM trends but magnitudes were not well correlated to resolved model values. These results support the use of homogenized laminate models for module shape derivation, though resolved models remain necessary for stress analyses. The accuracy of plate theory was found to be inadequate for most applications.

Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas — subsurface transport, turbulence, and anomalous materials — in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.

Abstract not provided.

Abstract not provided.

Efficient operation of battery energy storage systems requires that battery temperature remains within a specific range. Current techno-economic models neglect the parasitic loads heating and cooling operations have on these devices, assuming they operate at constant temperature. In this work, these effects are investigated considering the optimal sizing of battery energy storage systems when deployed in cold environments. A peak shaving application is presented as a linear programming problem which is then formulated in the PYOMO optimization programming language. The building energy simulation software EnergyPlus is used to model the heating, ventilation, and air conditioning load of the battery energy storage system enclosure. Case studies are conducted for eight locations in the United States considering a nickel manganese cobalt oxide lithium ion battery type and whether the power conversion system is inside or outside the enclosure. The results show an increase of 42% to 300% in energy capacity size, 43% to 217% in power rating, and 43% to 296% increase in capital cost dependent on location. This analysis shows that the heating, ventilation, and air conditioning load can have a large impact on the optimal sizes and cost of a battery energy storage system and merit consideration in techno-economic studies.

The carbon phase diagram is rich with polymorphs which possess very different physical and optical properties ideal for different scientific and engineering applications. An understanding of the dynamically driven phase transitions in carbon is particularly important for applications in inertial confinement fusion, as well as planetary and meteorite impact histories. Experiments on the Z Pulsed Power Facility at Sandia National Laboratories generate dynamically compressed high-pressure states of matter with exceptional uniformity, duration, and size that are ideal for investigations of fundamental material properties. X-ray diffraction (XRD) is an important material physics measurement because it enables direct observation of the strain and compression of the crystal lattice, and it enables the detection and identification of phase transitions. Several unique challenges of dynamic compression experiments on Z prevent using XRD systems typically utilized at other dynamic compression facilities, so novel XRD diagnostics have been designed and implemented. We performed experiments on Z to shock compress carbon (pyrolytic graphite) samples to pressures of 150–320 GPa. The Z-Beamlet Laser generated Mn-Heα (6.2 keV) X-rays to probe the shock-compressed carbon sample, and the new XRD diagnostics measured changes in the diffraction pattern as the carbon transformed into its high-pressure phases. Quantitative analysis of the dynamic XRD patterns in combination with continuum velocimetry information constrained the stability fields and melting of high-pressure carbon polymorphs.

The nitrogen-vacancy (NV) color center in diamond has demonstrated great promise in a wide range of quantum sensing. Recently, there have been a series of proposals and experiments using NV centers to detect spin noise of quantum materials near the diamond surface. This is a rich complex area of study with novel nano-magnetism and electronic behavior, that the NV center would be ideal for sensing. However, due to the electronic properties of the NV itself and its host material, getting high quality NV centers within nanometers of such systems is challenging. Band bending caused by space charges formed at the metal-semiconductor interface force the NV center into its insensitive charge states. Here, we investigate optimizing this interface by depositing thin metal films and thin insulating layers on a series of NV ensembles at different depths to characterize the impact of metal films on different ensemble depths. We find an improvement of coherence and dephasing times we attribute to ionization of other paramagnetic defects. The insulating layer of alumina between the metal and diamond provide improved photoluminescence and higher sensitivity in all modes of sensing as compared to direct contact with the metal, providing as much as a factor of 2 increase in sensitivity, decrease of integration time by a factor of 4, for NV T 1 relaxometry measurements.

Abstract not provided.

Abstract not provided.

Quantum point contacts (QPC) are the building blocks of quantum dot qubits and semiconducting quantum electrical metrology circuits. QPCs also make highly sensitive electrical amplifiers with the potential to operate in the quantum-limited regime. Though the inherent operational bandwidth of QPCs can eclipse the THz regime, the impedance mismatch with the external circuitry limits the operational frequency to a few kHz. Lumped-element impedance-matching circuits are successful only up to a few hundreds of MHz in frequency. QPCs are characterised by a complex impedance consisting of quantized resistance, capacitance, and inductance elements. Characterising the complex admittance at higher frequencies and understanding the coupling of QPC to other circuit elements and electromagnetic environments will provide valuable insight into its sensing and backaction properties. In this work, we couple a QPC galvanically to a superconducting stub tuner impedance matching circuit realised in a coplanar waveguide architecture to enhance the operation frequency into the GHz regime and investigate the electrical amplification and complex admittance characteristics. The device, operating at ~1.96 $GHz$ exhibits a conductance sensitivity of 2.92 X 10^-5(e²/h)/$\sqrt{Hz}$ with a bandwidth of 13 $MHz$. Besides, the RF reflected power unambiguously reveals the complex admittance characteristics of the QPC, shining more light on the behaviour of quantum tunnel junctions at higher operational frequencies.

The ion velocity distribution functions of thermonuclear plasmas generated by spherical laser direct drive implosions are studied using deuterium-tritium (DT) and deuterium-deuterium (DD) fusion neutron energy spectrum measurements. A hydrodynamic Maxwellian plasma model accurately describes measurements made from lower temperature (<10 keV), hydrodynamiclike plasmas, but is insufficient to describe measurements made from higher temperature more kineticlike plasmas. The high temperature measurements are more consistent with Vlasov-Fokker-Planck (VFP) simulation results which predict the presence of a bimodal plasma ion velocity distribution near peak neutron production. These measurements provide direct experimental evidence of non-Maxwellian ion velocity distributions in spherical shock driven implosions and provide useful data for benchmarking kinetic VFP simulations.

The goal of quantum tracking control is to identify shaped fields to steer observable expectation values along designated time-dependent tracks. The fields are determined via an iteration-free procedure, which is based on inverting the underlying dynamical equations governing the controlled observables. In this paper, we generalize the ideas in [Phys. Rev. A 98, 043429 (2018)2469-992610.1103/PhysRevA.98.043429] to the task of orienting symmetric top molecules in three dimensions. To this end, we derive equations for the control fields capable of directly tracking the expected value of the three-dimensional dipole orientation vector along a desired path in time. We show this framework can be utilized for tracking the orientation of linear molecules as well, and present numerical illustrations of these principles for symmetric-top tracking control problems.

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. This evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

Abstract not provided.

Abstract not provided.

Although visualizations are a useful tool for helping people to understand information, they can also have unintended effects on human cognition. This is especially true for uncertain information, which is difficult for people to understand. Prior work has found that different methods of visualizing uncertain information can produce different patterns of decision making from users. However, uncertainty can also be represented via text or numerical information, and few studies have systematically compared these types of representations to visualizations of uncertainty. We present two experiments that compared visual representations of risk (icon arrays) to numerical representations (natural frequencies) in a wildfire evacuation task. Like prior studies, we found that different types of visual cues led to different patterns of decision making. In addition, our comparison of visual and numerical representations of risk found that people were more likely to evacuate when they saw visualizations than when they saw numerical representations. These experiments reinforce the idea that design choices are not neutral: seemingly minor differences in how information is represented can have important impacts on human risk perception and decision making.

Abstract not provided.

Abstract not provided.

The injection and storage of anthropogenic CO2 in the subsurface is being deployed as a climate change mitigation tool; however, diagenetic-paragenetic heterogeneity in sandstone reservoirs often contributes to interval specific chemomechanical changes that affect injection and can increase leakage risk. Here, we address reservoir heterogeneities’ impact on chemomechanical changes in a macroporous-dominated lithofacies of Morrow B sandstone, a formation containing several diagenetically-distinct hydraulic facies while undergoing enhanced oil recovery (EOR) and carbon dioxide (CO2) sequestration. We performed three flow-through experiments using a CO2-charged or uncharged formation water combined with four indirect tensile strength tests per post-test sample. We then used the microstructure and paragenetic sequence to understand chemomechanical weakening with key observations as follows: dissolution of carbonates and feldspars changed porosity; increased permeability led to reclassifying each sample in a different hydraulic flow unit; decreased ultrasonic velocity; and did not lead to a loss of tensile strength. Tensile strength maintenance occurred due to the low abundance and minor dissolution of siderite, the stability of quartz, and the relative position of diagenetic ankerite within feldspar. This macroporous-dominated lithofacies is the primary reservoir for the Morrow B Sandstone, and is analogous to other porous sandstone reservoirs. It represents an end-member of a chemomechanically low-risk siliceous CO2 sequestration and CO2-EOR reservoir.

Living document updating the 5.14 version, SAND2023-03603 of May 2023.

Frontal polymerization involves the propagation of a thermally driven polymerization wave through a monomer solution to rapidly generate high-performance polymeric materials with little energy input. The balance between latent catalyst activation and sufficient reactivity to sustain a front can be difficult to achieve and often results in systems with poor storage lives. This is of particular concern for frontal ring-opening metathesis polymerization (FROMP) where gelation occurs within a single day of resin preparation due to the highly reactive nature of Grubbs-type catalysts. In this report we demonstrate the use of encapsulated catalysts to provide remarkable latency to frontal polymerization systems, specifically using the highly active dicyclopentadiene monomer system. Negligible differences were observed in the frontal velocities or thermomechanical properties of the resulting polymeric materials. FROMP systems with encapsulated catalyst particles are shown with storage lives exceeding 12 months and front rates that increase over a well-characterized 2 month period. Moreover, the modularity of this encapsulation method is demonstrated by encapsulating a platinum catalyst for the frontal polymerization of silicones by using hydrosilylation chemistry.

We explore the use of reduced physics models for efficient kinetic particle simulations of space charge limited (SCL) emission in inner magnetically insulated transmission lines (inner MITLs), with application to Sandia National Laboratories' Z machine. We propose a drift kinetic (guiding center) model of electron motion in place of a fully kinetic model and electrostatic-magnetostatic fields in place of electromagnetic fields. The validity of these approximations is suggested by the operational parameters of the Z machine, namely, current pulse lengths of order 100 ns compared with Larmor periods typically smaller than 10-11 s, typical Larmor radii of a few (tens) of microns (magnetic fields of tens to hundreds of Tesla) compared with MITL dimensions of a few centimeters, and transient time of light waves along the inner MITL of order a fraction of a nanosecond. Guiding center orbits eliminate the fast electron gyromotion, which enables the use of tens to hundreds of times larger time steps in the numerical particle advance. Electrostatic-magnetostatic fields eliminate the Courant-Friedrichs-Lewy (CFL) numerical stability limit on the time step and allow the use of higher grid resolutions or, alternatively, larger time steps in the fields advance. Overall, potential computational cost savings of tens to hundreds of times exists. The applicability of the reduced physics models is examined on two problems. First, in the simulation of space charge limited emission of electrons from the cathode surface due to high electric fields in a radial inner MITL geometry with a short load. In particular, it is shown that a drift kinetic-based particle-in-cell (PIC) model with electrostatic-magnetostatic fields is able to accurately reproduce well-known physics of electron vortex formation, spatially and temporally. Second, deeper understanding is gained of the mechanism behind vortex formation in this MITL geometry by considering an exemplar problem of an electron block of charge. This simpler setup reveals that the main mechanism of vortex formation can be attributed to pure drift motion of the electrons, that is, the (fully kinetic) gyromotion of the electrons is inessential to the process. This exemplar problem also suggests a correlation of the spatial dimensions of vortices to the thickness of the electron layer, as observed in SCL simulations. It also confirms that the electromagnetic nature of the fields does not play an essential role. Finally, an improved hybrid fully kinetic and drift kinetic model for electron motion is proposed, as means of capturing finite Larmor radius (FLR) effects; the particular FLR physics that is missed by the drift kinetic model is the particle-wall interaction. By initializing SCL emitted electrons as fully kinetic and later transitioning them to drift kinetic, according to simple criteria, the accuracy of SCL simulations can be improved, while preserving the potential for computational efficiency.

Abstract not provided.

The wave energy control competition established a benchmark problem which was offered as an open challenge to the wave energy system control community. The competition had two stages: In the first stage, competitors used a standard wave energy simulation platform (WEC-Sim) to evaluate their controllers while, in the second stage, competitors were invited to test their controllers in a real-time implementation on a prototype system in a wave tank. The performance function used was based on converted energy across a range of standard sea states, but also included aspects related to economic performance, such as peak/average power, peak force, etc. This paper compares simulated and experimental results and, in particular, examines if the results obtained in a linear system simulation are borne out in reality. Overall, within the scope of the device tested, the range of sea states employed, and the performance metric used, the conclusion is that high-performance WEC controllers work well in practice, with good carry-over from simulation to experimentation. However, the availability of a good WEC mathematical model is deemed to be crucial.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We demonstrate the discrimination of ground-state hyperfine manifolds of a cesium atom in an optical tweezer using a simple probe beam with Formula Presented% detection fidelity and 0.9(2)% detection-driven loss of bright-state atoms. Our detection infidelity of Formula Presented% is an order of magnitude better than previously published low-loss readout results for alkali-metal atoms in optical tweezers. We achieve these results by identifying and mitigating an extra depumping mechanism due to stimulated Raman transitions induced by trap light in the presence of probe light. In this work, complex optical systems and stringent vacuum pressures are not required, enabling straightforward adoption of our techniques on contemporary experiments.

Abstract not provided.

This project applies methods in Bayesian inference and modern statistical methods to quantify the value of new experimental data, in the form of new or modified diagnostic configurations and/or experiment designs. We demonstrate experiment design methods that can be used to identify the highest priority diagnostic improvements or experimental data to obtain in order to reduce uncertainties on critical inferred experimental quantities and select the best course of action to distinguish between competing physical models. Bayesian statistics and information theory provide the foundation for developing the necessary metrics, using two high impact experimental platforms on Z as exemplars to develop and illustrate the technique. We emphasize that the general methodology is extensible to new diagnostics (provided synthetic models are available), as well as additional platforms. We also discuss initial scoping of additional applications that began development in the last year of this LDRD.

This paper describes a process for forming a buried field shield in GaN by an etch-and-regrowth process, which is intended to protect the gate dielectric from high fields in the blocking state. GaN trench MOSFETs made at Sandia serve as the baseline to show the limitations in making a trench gated device without a method to protect the gate dielectric. Device data coupled with simulations show device failure at 30% of theoretical breakdown for devices made without a field shield. Implementation of a field shield reduces the simulated electric field in the dielectric to below 4 MV/cm at breakdown, which eliminates the requirement to derate the device in order to protect the dielectric. For realistic lithography tolerances, however, a shield-to-channel distance of 0.4 μm limits the field in the gate dielectric to 5 MV/cm and requires a small margin of device derating to safeguard a long-term reliability and lifetime of the dielectric.

Abstract not provided.

We combine crossed-beam velocity map imaging with high-level ab initio/transition state theory modeling of the reaction of S(3P) with 1,3-butadiene and isoprene under single collision conditions. For the butadiene reaction, we detect both H and H2 loss from the initial adduct, and from reaction with isoprene, we see both H loss and methyl loss. Theoretical calculations confirm these arise following intersystem crossing to the singlet surface forming long-lived intermediates. For the butadiene reaction, these lose H2 to form thiophene as the dominant channel, H to form the detected 2H-thiophenyl radical, or ethene, giving thioketene. For isoprene, additional reaction products are suggested by theory, including the observed H and methyl loss radicals, but also methyl thiophene, thioformaldehyde, and thioketene. The results for S(3P) + 1,3-butadiene, showing direct cyclization to the aromatic product and yielding few bimolecular product channels, are in striking contrast to those for the analogous O(3P) reaction.

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have 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 to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

This paper explores the utility of organizational system modeling frameworks to provide valuable insight into information flows within organizations and subsequently the opportunities for increasing resilience against disinformation campaigns targeting the system's ability to utilize information within its decision making. Disinformation is a growing challenge for many organizations and in recent years has created delay in decision making. Here the paper has utilized the viable systems model (VSM) to characterize organizational systems and used this approach to outline potential subsystem requirements to promote resilience of the system. The results of this paper can support the development of simulations and models considering the human elements within the system as well as support the development of quantitative measures of resilience.

Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains (~1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

This paper is concerned with goal-oriented a posteriori error estimation for nonlinear functionals in the context of nonlinear variational problems solved with continuous Galerkin finite element discretizations. A two-level, or discrete, adjoint-based approach for error estimation is considered. The traditional method to derive an error estimate in this context requires linearizing both the nonlinear variational form and the nonlinear functional of interest which introduces linearization errors into the error estimate. In this paper, we investigate these linearization errors. In particular, we develop a novel discrete goal-oriented error estimate that accounts for traditionally neglected nonlinear terms at the expense of greater computational cost. We demonstrate how this error estimate can be used to drive mesh adaptivity. Here, we show that accounting for linearization errors in the error estimate can improve its effectivity for several nonlinear model problems and quantities of interest. We also demonstrate that an adaptive strategy based on the newly proposed estimate can lead to more accurate approximations of the nonlinear functional with fewer degrees of freedom when compared to uniform refinement and traditional adjoint-based approaches.

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using higher backfill pressures in the canister, up to approximately 800 kPa, compared to their horizontal counterparts. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO₂, were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to add to previous testing that characterized SCCs under well-controlled boundary conditions through the inclusion of testing improvements that establish initial conditions in a more consistent way. While the engineered microchannel has dimensions similar to actual SCCs, it does not reproduce the tortuous path the aerosol laden flow would have to traverse for eventual transmission. SCCs can be rapidly grown in a laboratory setting given the right conditions, and initial characterization and clean-flow testing has begun on lab grown crack samples provided to Sandia National Laboratories (SNL). Many such samples are required to produce statistically relevant transmission results, and SNL is developing a procedure to produce samples in welded steel plates. These ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

Grid-scale battery energy storage systems (BESSs) are vulnerable to false data injection attacks (FDIAs), which could be used to disrupt state of charge (SoC) estimation. Inaccurate SoC estimation has negative impacts on system availability, reliability, safety, and the cost of operation. In this article a combination of a Cumulative Sum (CUSUM) algorithm and an improved input noise-aware extended Kalman filter (INAEKF) is proposed for the detection and identification of FDIAs in the voltage and current sensors of a battery stack. The series-connected stack is represented by equivalent circuit models, the SoC is modeled with a charge reservoir model and the states are estimated using the INAEKF. Further, the root mean squared error of the states’ estimation by the modified INAEKF was found to be superior to the traditional EKF. By employing the INAEKF, this article addresses the research gap that many state estimators make asymmetrical assumptions about the noise corrupting the system. Additionally, the INAEKF estimates the input allowing for the identification of FDIA, which many alternative methods are unable to achieve. The proposed algorithm was able to detect attacks in the voltage and current sensors in 99.16% of test cases, with no false positives. Utilizing the INAEKF compared to the standard EKF allowed for the identification of FDIA in the input of the system in 98.43% of test cases.

Intimately intertwined atomic and electronic structures of point defects govern diffusion-limited corrosion and underpin the operation of optoelectronic devices. For some materials, complex energy landscapes containing metastable defect configurations challenge first-principles modeling efforts. Here, we thoroughly reevaluate native point defect geometries for the illustrative case of α-Al2O3 by comparing three methods for sampling candidate geometries in density functional theory calculations: displacing atoms near a naively placed defect, initializing interstitials at high-symmetry points of a Voronoi decomposition, and Bayesian optimization. We find symmetry-breaking distortions for oxygen vacancies in some charge states, and we identify several distinct oxygen split-interstitial geometries that help explain literature discrepancies involving this defect. We also report a surprising and, to our knowledge, previously unknown trigonal geometry favored by aluminum interstitials in some charge states. These new configurations may have transformative impacts on our understanding of defect migration pathways in aluminum-oxide scales protecting metal alloys from corrosion. Overall, the Voronoi scheme appears most effective for sampling candidate interstitial sites because it always succeeded in finding the lowest-energy geometry identified in this study, although no approach found every metastable configuration. Finally, we show that the position of defect levels within the band gap can depend strongly on the defect geometry, underscoring the need to conduct careful searches for ground-state geometries in defect calculations.

Using compressive mechanical forces, such as pressure, to induce crystallographic phase transitions and mesostructural changes while modulating material properties in nanoparticles (NPs) is a unique way to discover new phase behaviors, create novel nanostructures, and study emerging properties that are difficult to achieve under conventional conditions. In recent decades, NPs of a plethora of chemical compositions, sizes, shapes, surface ligands, and self-assembled mesostructures have been studied under pressure by in-situ scattering and/or spectroscopy techniques. As a result, the fundamental knowledge of pressure-structure-property relationships has been significantly improved, leading to a better understanding of the design guidelines for nanomaterial synthesis. In the present review, we discuss experimental progress in NP high-pressure research conducted primarily over roughly the past four years on semiconductor NPs, metal and metal oxide NPs, and perovskite NPs. We focus on the pressure-induced behaviors of NPs at both the atomic- and mesoscales, inorganic NP property changes upon compression, and the structural and property transitions of perovskite NPs under pressure. We further discuss in depth progress on molecular modeling, including simulations of ligand behavior, phase-change chalcogenides, layered transition metal dichalcogenides, boron nitride, and inorganic and hybrid organic-inorganic perovskites NPs. These models now provide both mechanistic explanations of experimental observations and predictive guidelines for future experimental design. We conclude with a summary and our insights on future directions for exploration of nanomaterial phase transition, coupling, growth, and nanoelectronic and photonic properties.

Pultrusion manufacturing of fiber reinforced polymers has been shown to yield some of the highest mechanical properties for unidirectional composites, having a high degree of fiber alignment with consistent performance. Pultrusions offer a low-cost manufacturing approach for producing unidirectional composites with a constant cross-section and are used in many applications, including spar caps of wind turbine blades. However, as an intermediate processing step for wind blades, the additional cost of manufacturing pultrusions must be accompanied by sufficient increases in mechanical performance and system benefits. Wind turbine blades are manufactured using vacuum-assisted resin transfer molding with infused unidirectional fiberglass or carbon pultrusions for the spar cap. Infused fiberglass composites are among the most cost-effective structural materials available and replacing this material in the cost-driven wind industry has proven challenging, where infused fiberglass spar caps are still the predominant material system in use. To evaluate alternative material systems in a pultruded composite form, it is necessary to understand the costs for this additional manufacturing step which are shown to add 33%–55% on top of the material costs. A pultrusion cost model has been developed and used to quantify cost sensitivities to various processing parameters. The mechanical performance for pultruded composites is improved versus resin-infusion manufacturing with a 17% increase in design strength at a constant fiber volume fraction, but also enables higher achievable fiber volume fractions. The cost-specific mechanical performance is compared as a function of processing parameters for pultruded composites to identify the opportunities for alternative material and manufacturing approaches for wind turbine spar caps. Finally, four materials are compared in a representative wind turbine blade model to assess the performance of pultruded carbon fiber systems and pultruded fiberglass relative to infused fiberglass, where the pultruded systems produce lower weight blades with various cost distinctions.

Interlocking metasurfaces (ILMs) are a new class of mechanical metasurfaces built from architected arrays of interlocking features that can serve as a nonpermanent, robust joining technology. An ILM's strength is governed by the structural material, orientation, and topology of its latching unit cells. The presented work optimized the topologies of ILM unit cells to maximize strength in tensile and shear loading using gradient-based parametric optimization and genetic algorithms. Experimental validation confirmed that the optimized designs achieved considerable strength increases compared to a human intuitive design. In several cases, the optimized designs were approximately double the effective interfacial strength compared to that achieved via expert intuition alone. The strength improvement was seen for isolated unit cells and arrays of interacting unit cells (metasurfaces). An analysis of the topologies of the optimized designs showed that tall dendritic geometric features with large contact surfaces yield robust solutions in tension, while short and broad geometric features with large contact surfaces yield better results in shear loading. This study revealed the importance of shape optimization to maximize ILM effectiveness under single- and multi-objective scenarios.

The objective of the crystalline disposal work packages is to advance our understanding of long-term disposal of used fuel in crystalline rocks and to develop necessary experimental and computational capabilities to evaluate various disposal concepts in such media.

Barium titanate (BTO) is a widely researched ferroelectric useful for energy storage. While BTO’s surface chemistry is commonly studied using density functional theory, little has been published on the TiO₂ surface. Here, we determined that BTO’s surface response can be decoupled from the ferroelectric response by using a pre-optimized ferroelectric slab and allowing only the top three atomic z-layers to respond to ligand binding. Multiple favorable binding modes were identified for hydrogen, hydroxyl, water, and tert-butyl phosphonic acid on BTO’s TiO₂ surface. Of these ligands, tBuPA dominates surface binding with binding energies as low as -2.61 eV for its nine configurations.

The hydrogen sorption properties of single-phase bcc (TiVNb)100-xCrx alloys (x = 0-35) are reported. All alloys absorb hydrogen quickly at 25 °C, forming fcc hydrides with storage capacity depending on the Cr content. A thermodynamic destabilization of the fcc hydride is observed with increasing Cr concentration, which agrees well with previous compositional machine learning models for metal hydride thermodynamics. The steric effect or repulsive interactions between Cr-H might be responsible for this behavior. The cycling performances of the TiVNbCr alloy show an initial decrease in capacity, which cannot be explained by a structural change. Pair distribution function analysis of the total X-ray scattering on the first and last cycled hydrides demonstrated an average random fcc structure without lattice distortion at short-range order. If the as-cast alloy contains a very low density of defects, the first hydrogen absorption introduces dislocations and vacancies that cumulate into small vacancy clusters, as revealed by positron annihilation spectroscopy. Finally, the main reason for the capacity drop seems to be due to dislocations formed during cycling, while the presence of vacancy clusters might be related to the lattice relaxation. Having identified the major contribution to the capacity loss, compositional modifications to the TiVNbCr system can now be explored that minimize defect formation and maximize material cycling performance.

Immune checkpoint immunotherapy (ICI) can re-activate immune reactions against neoantigens, leading to remarkable remission in cancer patients. Nevertheless, only a minority of patients are responsive to ICI, and approaches for prediction of responsiveness are needed to improve the success of cancer treatments. While the tumor mutational burden (TMB) correlates positively with responsiveness and survival of patients undergoing ICI, the influence of the subcellular localizations of the neoantigens remains unclear. Here, we demonstrate in both a mouse melanoma model and human clinical datasets of 1,722 ICI-treated patients that a high proportion of membrane-localized neoantigens, particularly at the plasma membrane, correlate with responsiveness to ICI therapy and improved overall survival across multiple cancer types. We further show that combining membrane localization and TMB analyses can enhance the predictability of cancer patient response to ICI. Our results may have important implications for establishing future clinical guidelines to direct the choice of treatment toward ICI.

This FY2023 report is the second update to the Disposal Research (DR) Research and Development (R&D) 5-year plan for the Spent Fuel and Waste Science and Technology (SFWST) Campaign DR R&D activities. In the planning for FY2020 in the U.S. Department of Energy (DOE) NE-81 SFWST Campaign, the DOE requested development of a high-level summary plan for activities in the DR R&D program for the next five (5)-year period, with periodic updates to this summary plan. The DR R&D 5-year plan was provided to the DOE based initially on the FY2020 priorities and program structure (initial 2020 version of this 5-year plan) and provides a strategic summary guide to the work within the DR R&D technical areas (Control Accounts, CA), focusing on the highest priority technical thrusts. This 5-year plan is a living document (planned to be updated periodically) that provides review of SFWST R&D accomplishments (as seen on the 2021 revision of this 5-year plan), describes changes to technical R&D prioritization based on (a) progress in each technical area (including external technical understanding) with specific accomplishments and (b) any changes in SFWST Campaign objectives and/or funding levels (i.e., Program Direction). Updates to this 5-year plan include the DR R&D adjustments to high-priority knowledge gaps to be investigated in the near-term, as well as the updated longer-term DR R&D directions for the program activities. This plan fulfills the Milestone M2SF23SN010304083 in DR Work Package (WP) SF-23SN01030408 (GDSA - Framework Development – SNL).

Recent advances in the growth of aluminum scandium nitride films on silicon suggest that this material platform could be applied for quantum electromechanical applications. Here, we model, fabricate, and characterize microwave frequency silicon phononic delay lines with transducers formed in an adjacent aluminum scandium nitride layer to evaluate aluminum scandium nitride films, at 32% scandium, on silicon interdigital transducers for piezoelectric transduction into suspended silicon membranes. We achieve an electromechanical coupling coefficient of 2.7% for the extensional symmetric-like Lamb mode supported in the suspended material stack and show how this coupling coefficient could be increased to at least 8.5%, which would further boost transduction efficiency and reduce the device footprint. The one-sided transduction efficiency, which quantifies the efficiency at which the source of microwave photons is converted to microwave phonons in the silicon membrane, is 10% at 5 GHz at room temperature and, as we discuss, there is a path to increase this toward near-unity efficiency based on a combination of modified device design and operation at cryogenic temperatures.