Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

We report low-temperature magneto-transport measurements of an undoped Si/SiGe asymmetric double quantum well heterostructure. The density in both layers is tuned independently utilizing top and bottom gates, allowing the investigation of quantum wells at both imbalanced and matched densities. Integer quantum Hall states at total filling factor ν T = 1 and ν T = 2 are observed in both density regimes, and the evolution of their excitation gaps is reported as a function of the density. The ν T = 1 gap evolution departs from the behavior generally observed for valley splitting in the single layer regime. Furthermore, by comparing the ν T = 2 gap to the single particle tunneling energy, Δ SAS, obtained from Schrödinger-Poisson (SP) simulations, evidence for the onset of spontaneous interlayer coherence is observed for a relative filling fraction imbalance smaller than ∼ 50 %.

This report represents completion of milestone deliverable M2SF-22SN010309082 Annual Status Update for OWL, which is due on November 30, 2021 as part of the fiscal year 2022 (FY2022) work package SF-22SN01030908. This report provides an annual update on status of FY2021 activities for the work package “OWL - Inventory – SNL”. The Online Waste Library (OWL) has been designed to contain information regarding United States (U.S.) Department of Energy (DOE)-managed (as) high-level waste (DHLW), DOE-managed spent nuclear fuel (DSNF), and other wastes that are likely candidates for deep geologic disposal. Links to the current supporting documents for the data are provided when possible; however, no classified or official-use-only (OUO) data are planned to be included in OWL. There may be up to several hundred different DOE-managed wastes that are likely to require deep geologic disposal. This report contains new information on sodium-bonded spent fuel waste types and wastes forms, which are included in the next release of OWL, Version 3.0, on the Sandia National Laboratories (SNL) External Collaboration Network (ECN). The report also provides an update on the effort to include information regarding the types of vessels capable of disposing of DOE-managed waste.

In “Interfacial structures and acidity constants (pKa) of goethite from first principles molecular dynamics simulations,” authors Y. Zhang, X. Lui, J. Cheng, and X. Lu apply First Principles molecular dynamics (FPMD, also called Density Functional Theory MD, DFT/MD, or ab initio MD, AIMD), to evaluable the complete set of acidity constants (pK_a) of the hydroxyl groups on the most prominent goethite crystal facets. The pK_a of these OH and OH$^+_2$ groups are compared with available data from the multisite complexation (MUSIC) model traditionally used to estimate pK_a on mineral surfaces. The authors have presented eloquent rationale for the importance and implications of understanding goethite acidity constants in room temperature geochemistry settings. In this paper, I focus on the computational aspects, the strengths of FPMD, and its possibilities.

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid-wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R2) of ∼0.99 and a mean squared error of ∼2.9 × 10-5. Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

There is an intensive effort to control the nature of attractive interactions between ultrathin semiconductors and metals and to understand its impact on the electronic properties at the junction. Here, we present a photoelectron spectroscopy study on the interface between WS2 films and gold, with a focus on the occupied electronic states near the Brillouin zone center (i.e., the point). To delineate the spectra of WS2 supported on crystalline Au from the suspended WS2, we employ a microscopy approach and a tailored sample structure, in which the WS2/Au junction forms a semi-epitaxial relationship and is adjacent to suspended WS2 regions. The photoelectron spectra, as a function of WS2 thickness, display the expected splitting of the highest occupied states at the point. In multilayer WS2, we discovered variations in the electronic states that spatially align with the crystalline grains of underlying Au. Corroborated by density functional theory calculations, we attribute the electronic structure variations to stacking variations within the WS2 films. We propose that strong interactions exerted by Au grains cause slippage of the interfacing WS2 layer with respect to the rest of the WS2 film. Our findings illustrate that the electronic properties of transition metal dichalcogenides, and more generally 2D layered materials, are physically altered by the interactions with the interfacing materials, in addition to the electron screening and defects that have been widely considered.

Metal hydrides can store hydrogen isotopes with high volumetric density. In metal tritides, tritium beta decay can result in accumulation of helium within the solid, in some cases exceeding 10 at.% helium after only 4 years of aging. Helium is insoluble in most materials, but often does not readily escape, and instead coalesces to form nanoscale bubbles when helium concentrations are near 1 at.%. Blistering or spallation often occurs at higher concentrations. Radioactive particles shed during this process present a potential safety hazard. This study investigates the effects of high helium concentrations on erbium deuteride (ErD₂), a non-radioactive surrogate material for erbium tritide (ErT₂). To simulate tritium decay in the surrogate, high doses of 120 keV helium ions were implanted into ErD₂ films at room temperature. Scanning and transmission electron microscopy indicated spherical helium bubble formation at a critical concentration of 1.5 at.% and bubble linkage leading to nanoscale crack formation at a concentration of 7.5 at.%. Additionally, crack propagation occurred through the nanocrack region, resulting in spallation extending from the implantation peak to the surface. Electron energy loss spectroscopy was utilized to confirm the presence of high-pressure helium in the nanocracks, suggesting that helium gas plays a predominant role in deformation. This work improves the overall understanding of helium behavior in ErD₂ by using modern characterization techniques to determine: the critical helium concentration required for bubble formation, the material failure mechanism at high concentration, and the nanoscale mechanisms responsible for material failure in helium implanted ErD₂.

This report analyzes experimental data from Test Plans TP 08-02, TP 12-02, and TP 20-01 to add new log K values and Pitzer interaction parameters for Fe, Pb, Mg, Nd and B reactions to the WIPP geochemical thermodynamic database, data0.fm 1.

The objective of this report is to accept or reject the hypothesis that the experiments conducted under TP 08-02 Revision 0 (Ismail et al., 2008) were affected by CO₂(g) intrusion and sample contamination. The test of the hypothesis is accomplished by comparing the experimental data collected under the protocols of TP 08-02 Revision O and TP 20-01 Revision O (Kirkes and Zhang, 2020). The protocols of TP 20-01 Revision 0 minimize the possibilities of CO₂(g) intrusion and sample contamination. The experimental data sets obtained under both TPs will be assessed statistically to see if they are identical or not.

The Dynamic Networks Experiment 2018 (DNE18) was a collaborative effort between Los Alamos National Laboratory (LANL), Sandia National Laboratories (SNL), Lawrence Livermore National Laboratory (LLNL) and Pacific Northwest National Laboratory (PNNL) designed to evaluate methodologies for multi-modal data ingestion and processing. One component of this virtual experiment was a quantitative assessment of current capabilities for infrasound data processing, beginning with the establishment of a baseline for infrasound signal detection. To produce such baselines, SNL and LANL exploited a common dataset of infrasound data recorded across a regional network in Utah from December 2010 through February 2011. We utilize two automated signal detectors, the Adaptive F-Detector (AFD) and the Multivariate Adaptive Learning Detector (MALD) to produce automated signal detection catalogs and an analyst-produced catalog. Comparisons indicate that automatic detectors may be able to identify small amplitude, low SNR events that cannot be identified by analyst review. We document detector performance in terms of precision and recall, demonstrating that the AFD is more precise, but the MALD has higher recall. We use a synthetic dataset of signals embedded in pink noise in order to highlight shortcomings in assessing detection algorithms for low signal to noise ratio signals which are commonly of interest to the nuclear monitoring community. For comparisons utilizing the synthetic dataset, the AFD has higher recall while precision is equal for both detectors. These results indicate that both detectors perform well across a variety of background noise environments; however, both detectors fail to identify repetitive, short duration signals arriving from similar backazimuths. These failures represent specific scenarios that could be targeted for further detector development.

Compared to the classical Lanczos algorithm, the s-step Lanczos variant has the potential to improve performance by asymptotically decreasing the synchronization cost per iteration. However, this comes at a price; despite being mathematically equivalent, the s-step variant may behave quite differently in finite precision, potentially exhibiting greater loss of accuracy and slower convergence relative to the classical algorithm. It has previously been shown that the errors in the s-step version follow the same structure as the errors in the classical algorithm, but are amplified by a factor depending on the square of the condition number of the O(s)-dimensional Krylov bases computed in each outer loop. As the condition number of these s-step bases grows (in some cases very quickly) with s, this limits the s values that can be chosen and thus can limit the attainable performance. In this work, we show that if a select few computations in s-step Lanczos are performed in double the working precision, the error terms then depend only linearly on the conditioning of the s-step bases. This has the potential for drastically improving the numerical behavior of the algorithm with little impact on per-iteration performance. Our numerical experiments demonstrate the improved numerical behavior possible with the mixed precision approach, and also show that this improved behavior extends to mixed precision s-step CG. Here, we present preliminary performance results on NVIDIA V100 GPUs that show that the overhead of extra precision is minimal if one uses precisions implemented in hardware.

Characteristic modes on infinite periodic structures are studied using spectral dyadic Green’s functions. This formulation demonstrates that, in contrast to the modal analysis of finite structures, the number of radiating characteristic modes is limited by unit cell size and incident wave vector (i.e., scan angle or phase shift per unit cell). Here, the reflection tensor is decomposed into modal contributions from radiating modes, indicating that characteristic modes are a predictably sparse basis in which to study reflection phenomena.

The interplay of stress, disorder, and Coulomb screening dictating the mobility of doped cadmium oxide (CdO) is examined using Raman spectroscopy to identify the mechanisms driving dopant incorporation and scattering within this emerging infrared optical material. Specifically, multi-wavelength Raman and UV-vis spectroscopies are combined with electrical Hall measurements on a series of yttrium (X = Y) and indium (X = In) doped X:CdO thin-films. Hall measurements confirm n-type doping and establish carrier concentrations and mobilities. Spectral fitting along the low-frequency Raman combination bands, especially the TA+TO(X) mode, reveals that the evolution of strain and disorder within the lattice as a function of dopant concentration is strongly correlated with mobility. Coupling between the electronic and lattice environments was examined through analysis of first- and second-order longitudinal-optical phonon-plasmon coupled modes that monotonically decrease in energy and asymmetrically broaden with increasing dopant concentration. By fitting these trends to an impurity-induced Fröhlich model for the Raman scattering intensity, exciton-phonon and exciton-impurity coupling factors are quantified. These coupling factors indicate a continual decrease in the amount of ionized impurity scattering with increasing dopant concentration and are not as well correlated with mobility. This shows that lattice strain and disorder are the primary determining factors for mobility in donor-doped CdO. In aggregate, the study confirms previously postulated defect equilibrium arguments for dopant incorporation in CdO while at the same time identifying paths for its further refinement.

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

Borehole cement is used across the range of energy technologies to stabilize casing, to serve as a barrier to behind-casing fluid movement. Cement debonding and other flaws, both at cement interfaces and within the cement itself, can create leakage pathways that can threaten safety to personnel, and wellbore performance, with economic and regulatory consequences. A primary method to assess cement health and wellbore integrity is via acoustic methods. This project was designed with three aims: demonstrate a significant improvement in the interpretation of cement acoustic behavior, both during curing, and in interpreting effects of flaws and evolving interfaces; develop sensor technologies to improve signal-noise ratios and cement acoustic responses; and lastly, provide a borehole demonstration of at least one of these technologies. We have accomplished the first two objectives, and the third, delayed by pandemic health concerns, is proceeding as of this writing via a technology partner with the University of Texas Advanced Energy Consortium.

Elevated temperature and pressure in the earth's subsurface alters the permeability of salt formations, due to changing properties of the salt-brine interface. Molecular dynamics (MD) simulations are used to investigate the mechanisms of temperature and pressure dependence of liquid-solid interfacial tensions of NaCl, KCl, and NaCl-KCl brines in contact with (100) salt surfaces. Salt-brine dihedral angles vary between 55 and 76° across the temperature (300-450 K) and pressure range (0-150 MPa) evaluated. Temperature-dependent brine composition results in elevated dihedral angles of 65-80°, which falls above the reported salt percolation threshold of 60°. Mixed NaCl-KCl brine compositions increased this effect. Elevated temperatures excluded dissolved Na+ ions from the interface, causing the strong temperature dependence of the liquid-solid interfacial tension and the resulting dihedral angle. Therefore, at higher temperature, pressure, and brine concentrations Na-Cl systems may underpredict the dihedral angle. Higher dihedral angles in more realistic mixed brine systems maintain low permeability of salt formations due to changes in the structure and energetics of the salt-brine interface.

This project is intended to support the development of new traction drive systems that meet the targets of 100 kW/L for power electronics and 50 kW/L for electric machines with reliable operation to 300,000 miles. To meet these goals, new designs must be identified that make use of state-of-the-art and next-generation electronic materials and design methods. Designs must exploit synergies between components, for example converters designed for high-frequency switching using wide band gap devices and ceramic capacitors. This project includes: (1) a survey of available technologies; (2) the development of design tools that consider the converter volume and performance; (3) exercising the design software to evaluate performance gaps and predict the impact of certain technologies and design approaches, i.e. GaN semiconductors, ceramic capacitors, and select topologies; and (4) building and testing hardware prototypes to validate models and concepts. Early instantiations of the design tools enable co-optimization of the power module and passive elements and provide some design guidance; later instantiations will enable the co-optimization of inverter and machine. Prototype testing begins with evaluation of simpler conversion topologies (i.e. the half-bridge boost converter) and progresses with fabrication of prototype inverter drives.

Novel materials based on the aluminum oxyhydroxide boehmite phase were prepared using a glycothermal reaction in 1,4-butanediol. Under the synthesis conditions, the atomic structure of the boehmite phase is altered by the glycol solvent in place of the interlayer hydroxyl groups, creating glycoboehmite. The structure of glycoboehmite was examined in detail to determine that glycol molecules are intercalated in a bilayer structure, which would suggest that there is twice the expansion identified previously in the literature. This precursor phase enables synthesis of two new phases that incorporate either polyvinylpyrrolidone or hydroxylpropyl cellulose nonionic polymers. These new materials exhibit changes in morphology, thermal properties, and surface chemistry. All the intercalated phases were investigated using PXRD, HRSTEM, SEM, FT-IR, TGA/DSC, zeta potential titrations, and specific surface area measurement. These intercalation polymers are non-ionic and interact through wetting interactions and hydrogen bonding, rather than by chemisorption or chelation with the aluminum ions in the structure.

The global energy system is undergoing significant changes, including a shift in energy generating technologies to more renewable energy sources. However, the dependence of renewable energy sources on local environmental conditions could also increase disruptions in service through exposures to compound, extreme weather events. By fusing three diverse datasets (operations and maintenance tickets, weather data, and production data), this analysis presents a novel methodology to identify and evaluate performance impacts arising from extreme weather events across diverse geographical regions. Text analysis of maintenance tickets identified snow, hurricanes, and storms as the leading extreme weather events affecting photovoltaic plants in the United States. Statistical techniques and machine learning were then implemented to identify the magnitude and variability of these extreme weather impacts on site performance. Impacts varied between event and non-event days, with snow events causing the greatest reductions in performance (54.5%), followed by hurricanes (12.6%) and storms (1.1%). Machine learning analysis identified key features in determining if a day is categorized as low performing, such as low irradiance, geographic location, weather features, and site size. This analysis improves our understanding of compound, extreme weather event impacts on photovoltaic systems. These insights can inform planning activities, especially as renewable energy continues to expand into new geographic and climatic regions around the world.

This work studies the different types of behavior and inaccuracies that can occur when contact is not adequately accounted for in a dynamical system with freeplay, as the strength of the contact stiffness increases. The MATLAB® ode45 time integration solver, with the built-in Event Location capability, is first validated using past experimental data from a forced Duffing oscillator with freeplay. Next, numerical results utilizing event location are compared to results neglecting event location in order to highlight possible numerical errors and effects on multistable dynamical responses. Inaccuracies tend to occur in two different ways. First, neglecting event location can affect the boundaries between basins of attraction. Second, neglecting event location has little effect on the behaviors of the attractor solutions themselves besides merely resembling poorly converged solutions. Errors are less pronounced at the limits of soft or hard contact stiffness. This study shows the importance of accurately solving piecewise-smooth systems and the existing correlation between the strength of the contact force and possible numerical inaccuracies.

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices the result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

This work focuses on the space-time reduced-order modeling (ROM) method for solving large-scale uncertainty quantification (UQ) problems with multiple random coefficients. In contrast with the traditional space ROM approach, which performs dimension reduction in the spatial dimension, the space-time ROM approach performs dimension reduction on both the spatial and temporal domains, and thus enables accurate approximate solutions at a low cost. We incorporate the space-time ROM strategy with various classical stochastic UQ propagation methods such as stochastic Galerkin and Monte Carlo. Numerical results demonstrate that our methodology has significant computational advantages compared to state-of-the-art ROM approaches. By testing the approximation errors, we show that there is no obvious loss of simulation accuracy for space-time ROM given its high computational efficiency.

Heβ spectral line shapes are important for diagnosing temperature and density in many dense plasmas. This work presents Heβ line shapes measured with high spectral resolution from solid-density plasmas with minimized gradients. The line shapes show hallmark features of Stark broadening, including quantifiable redshifts and double-peaked structure with a significant dip between the peaks; these features are compared to models through a Markov chain Monte Carlo framework. Line shape theory using the dipole approximation can fit the width and peak separation of measured line shapes, but it cannot resolve an ambiguity between electron density ne and ion temperature Ti, since both parameters influence the strength of quasistatic ion microfields. Here a line shape model employing a full Coulomb interaction for the electron broadening computes self-consistent line widths and redshifts through the monopole term; redshifts have different dependence on plasma parameters and thus resolve the ne-Ti ambiguity. The measured line shapes indicate densities that are 80-100% of solid, identifying a regime of highly ionized but well-tamped plasma. This analysis also provides the first strong evidence that dense ions and electrons are not in thermal equilibrium, despite equilibration times much shorter than the duration of x-ray emission; cooler ions may arise from nonclassical thermalization rates or anomalous energy transport. The experimental platform and diagnostic technique constitute a promising new approach for studying ion-electron equilibration in dense plasmas.