Publications

Abstract not provided.

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O-H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

CFD (Computational Fluid Dynamic) simulation of aerosol-laden natural convective flow and particle deposition in a spent fuel storage canister with 37 assemblies is currently computationally prohibitive. PWR (Pressurized Water Reactor) assemblies have up to 289 pins or tubes with several spacer grids to align the pins. Spacer grids with mixing vanes induce swirling during operation to increase heat transfer. Each spacer grid contains hundreds of small structures such as retaining clips, channel walls, and openings. The largest canisters store 37 PWR assemblies thus, there are numerous pins, tubes, and spacer grids for which the flow region between and around these structures needs to be determined along with the movement and deposition of aerosol particles. Because of the complicated geometry, modeling the intricate flow even for just one assembly is currently impractical. Nonetheless, we are developing techniques for a practical model to assess the natural aerosol particle deposition process in a canister in the event that a release occurs from one or more fuel pins. In the previous work it was demonstrated that CFD can model the flow through a PWR spacer grid with mixing vanes, including particle deposition, in a reasonable amount of time on a personal computer. In this work, the analysis is extended to include the bypass region between an assembly and the canister basket walls. It is shown that the flow velocity in the bypass region is about three times that of the interstitial region between the pins. The lengths before and after the spacer grid are also extended to determine when the flow becomes fully developed. In addition, the approach of computationally “stitching together” segments of an assembly is demonstrated with the plan to ultimately model a full assembly. The fraction of particles that are deposited in a segment with a spacer grid is determined as a function of particle size and flow velocity.

The paper deals with a new effective numerical technique on unfitted Cartesian meshes for simulations of heterogeneous elastic materials. Here, we develop the optimal local truncation error method (OLTEM) with 27- point stencils (similar to those for linear finite elements) for the 3-D time-independent elasticity equations with irregular interfaces. Only displacement unknowns at each internal Cartesian grid point are used. The interface conditions are added to the expression for the local truncation error and do not change the width of the stencils. The unknown stencil coefficients are calculated by the minimization of the local truncation error of the stencil equations and yield the optimal second order of accuracy for OLTEM with the 27-point stencils on unfitted Cartesian meshes. A new post-processing procedure for accurate stress calculations has been developed. Similar to basic computations it uses OLTEM with the 27-point stencils and the elasticity equations. The post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing numerical techniques (e.g., with finite elements). Numerical experiments show that at an accuracy of 0.1% for stresses, OLTEM with the new post-processing procedure significantly (by 10⁵-10⁹ times) reduces the number of degrees of freedom compared to linear finite elements. OLTEM with the 27-point stencils yields even more accurate results than high-order finite elements with wider stencils.

Junctions are discontinuities in flat grain boundaries that arise in all polycrystalline materials and are thought to play important roles in the response of a grain boundary network to thermal and mechanical loads. A key open question concerns the mechanisms by which solute segregation to junctions impacts properties of the grain boundary. Here, in this work, we investigate the influence of grain boundary facet junctions on solute embrittlement, and we present an analytical model that uses the hydrostatic stress field contributed by dislocations at multiple junctions to describe these effects. Specifically, we study junctions between {112} facets of various lengths in Au $\langle111\rangle$ Σ3 tilt grain boundaries. Copper and silver solutes are employed to determine if the effect of junctions on solute segregation and embrittlement is dependent on size relative to the host. Combined, atomistic simulation data and the analytical model show that Cu and Ag have opposite segregation responses to junctions due to the sign of the hydrostatic stress field induced by junctions. However, a positive shift in the embrittling potency is computed near junctions regardless of solute type or the stress state of the segregation site. Hence, for the conditions studied, junctions consistently shift the energetic landscape towards embrittlement.

Bonding diamond to the back side of gallium nitride (GaN) electronics has been shown to improve thermal management in lateral devices; however, engineering challenges remain with the bonding process and characterizing the bond quality for vertical device architectures. Here, integration of these two materials is achieved by room-temperature compression bonding centimeter-scale GaN and a diamond die via an intermetallic bonding layer of Ti/Au. Recent attempts at GaN/diamond bonding have utilized a modified surface activation bonding (SAB) method, which requires Ar fast atom bombardment immediately followed by bonding within the same tool under ultrahigh vacuum (UHV) conditions. The method presented here does not require a dedicated SAB tool yet still achieves bonding via a room-temperature metal-metal compression process. Imaging of the buried interface and the total bonding area is achieved via transmission electron microscopy (TEM) and confocal acoustic scanning microscopy (C-SAM), respectively. The thermal transport quality of the bond is extracted from spatially resolved frequency-domain thermoreflectance (FDTR) with the bonded areas boasting a thermal boundary conductance of >100 MW/m2·K. Additionally, Raman maps of GaN near the GaN-diamond interface reveal a low level of compressive stress, <80 MPa, in well-bonded regions. FDTR and Raman were coutilized to map these buried interfaces and revealed some poor thermally bonded areas bordered by high-stress regions, highlighting the importance of spatial sampling for a complete picture of bond quality. Overall, this work demonstrates a novel method for thermal management in vertical GaN devices that maintains low intrinsic stresses while boasting high thermal boundary conductances.

Testing of a compact Bremsstrahlung diode at the High Energy Radiation Megavolt Electron Source III (HERMES-III) was performed at Sandia National Laboratories in November, 2023. The compact diode described here is the first prototype diode in a campaign to optimize a Bremsstrahlung diode in terms of size and dose production. The goal was to test the diode at 13MV, and the experiment realized between 10-12MV at the diode. Modeling and simulation of this geometry was performed after the test, shedding insight into several phenomena seen by experimental diagnostics. Modifications to the diode and experiment are proposed for future experiments to help explain the phenomena and approach a better final design.

Ionic assemblies, or clusters, determine the structure and dynamics of ionizable polymers and enable their many applications. Fundamental to attaining well-defined materials is controlling the balance between the van der Waals interactions that govern the backbone behavior and the forces that drive the formation of ionic clusters. Here, using small-angle neutron scattering and fully atomistic molecular dynamics simulations, the structure of a model ionomer, sulfonated polystyrene in toluene solutions, was investigated as the cluster cohesion was tweaked by the addition of ethanol. The static structure factor was measured by both techniques and correlated with the size of the ionic clusters as the polymer concentration was varied. The conjunction of SANS results and molecular insight from MD simulations enabled the determination of the structure of these inhomogeneous networks on multiple length scales. We find that across the entire concentration range studied, a network driven by the formation of ionic clusters was formed, where the size of the clusters drives the inhomogeneity of these systems. Tweaking the ionic clusters through the addition of ethanol impacts the packing of the sulfonated groups, their shape, and their size distribution, which, in turn, affects the structure of these networks.

Coarse-grained molecular dynamics simulations are used to study the diffusion of thin nanorods in entangled polymer melts for varying nanorod length and roughness. While prior studies observed a nanorod parallel diffusion constant scaling inversely with rod length D_∥ ~ l^–1, here, we show that this scaling is not universal and depends sensitively on the nanorod surface roughness. We observe D_∥ ~ l^–k, where k < 1 and decreases with decreasing surface roughness. The weaker scaling is driven by the non-Gaussian diffusion of nanorods due to the emergence of an intermittent hopping process that becomes more pronounced with decreasing roughness at the monomer scale. Analysis shows that the mean hop size grows for smoother rods but shows little to no variation with rod length. The mean hopping frequency shows no dependence on either rod length or roughness, suggesting it originates from the polymer melt environment. Further, our results show that the small-scale features of the nanorod surface strongly influence the large-scale and long-time transport of nanorods in polymer matrices, creating new material design opportunities for precisely engineered nanocomposites.

Characterizing and identifying cells in multicellular in vitro models remain a substantial challenge. Here, we utilize hyperspectral confocal Raman microscopy and principal component analysis coupled with linear discriminant analysis to form a label-free, noninvasive approach for classifying bone cells and osteosarcoma cells. Through the development of a library of hyperspectral Raman images of the K7M2-wt osteosarcoma cell lines, 7F2 osteoblast cell lines, RAW 264.7 macrophage cell line, and osteoclasts induced from RAW 264.7 macrophages, we built a linear discriminant model capable of correctly identifying each of these cell types. The model was cross-validated using a k-fold cross validation scheme. The results show a minimum of 72% accuracy in predicting cell type. We also utilize the model to reconstruct the spectra of K7M2 and 7F2 to determine whether osteosarcoma cancer cells and normal osteoblasts have any prominent differences that can be captured by Raman. We find that the main differences between these two cell types are the prominence of the β-sheet protein secondary structure in K7M2 versus the α-helix protein secondary structure in 7F2. Additionally, differences in the CH2 deformation Raman feature highlight that the membrane lipid structure is different between these cells, which may affect the overall signaling and functional contrasts. Overall, we show that hyperspectral confocal Raman microscopy can serve as an effective tool for label-free, nondestructive cellular classification and that the spectral reconstructions can be used to gain deeper insight into the differences that drive different functional outcomes of different cells.

The First Tri-Lab Workshop on Formal Verification was held in Santa Fe, New Mexico, on December 5th, 2023. This workshop gathered staff from Sandia, Los Alamos, and Lawrence Livermore National Laboratories and NASA’s Jet Propulsion Laboratory. This report summarizes and expands on the presentations given and discussion had at this workshop. In this report, we describe the current capabilities and research needs related to formal methods at the NNSA labs. In particular, we identify medium-term and long-term research gaps in programming languages, formalization efforts of complex systems, embedded systems verification, hardware verification, cybersecurity, formal methods usability, workflows, numerical methods, the use of formal methods for artificial intelligence (and its converse, artificial intelligence for formal methods), and collaboration opportunities and considerations on these topics. We conclude with a small number of exemplar research problems related to these topics.

Superconducting topological systems formed by a strong 3D topological insulator (TI) in proximity to a conventional s-wave superconductor (SC) have been intensely studied, as they may host Majorana zero modes. However, there are limited experimental realizations of TI-SC systems in which robust superconducting pairing is induced on the surface states of the TI and a topological superconducting state is established. Here, we fabricate a TI-SC system by depositing, via a focused ion beam, tungsten (W) nanoscale clusters on the surface of TI Bi0.91Sb0.09. We find that the resulting heterostructure supports phase-slip lines (PSLs) that act as effective Josephson junctions (JJs). We probe the response of the system to microwave radiation. We find that for some ac frequencies, and powers, the resulting Shapiro steps’ structure of the voltage-current characteristic exhibits a missing first step and an unexpectedly wide second Shapiro step. The theoretical analysis of the measurements shows that the unusual Shapiro response arises from the interplay between a static JJ and a dynamic one and allows us to identify the conditions under which the missing first step can be attributed to the topological nature of the JJs formed by the PSLs. Our results suggest an approach to induce superconductivity in a TI, a route to realizing highly transparent topological JJs, and show how the response of superconducting systems to microwave radiation can be used to infer the dynamics of PSLs. Highly transparent topological junctions are promising candidates to realize vector field sensors with very high sensitivity. In addition, due to the nontrivial Berry phase of the TI’s surface states such junctions can be in a topological state which is ideal to create topologically protected qubits.

Testing is pivotal for early identification of disease and subsequent infection control. Pathogens’ nucleic acid sequence can change due to naturally-occurring genetic drift or intentional modification. Because of the reliance on molecular assays for human, animal, and plant disease diagnosis, we must understand how nucleotide mutations affect test accuracy. Primers designed against original lineages of a pathogen may be less efficient at detecting variants with genetic changes in priming regions. Here, we made single- and multi-point mutations in priming regions of a model SARS-CoV-2 template that was used as input for a loop-mediated isothermal amplification (LAMP) assay. We found that many of the modifications impacted assay sensitivity, amplification speed, or both. Further research exploring mutations at every position in each of the eight priming regions should be conducted to evaluate trends and determine generalizability.