Publications

As the demand for higher-performance batteries has increased, so has the body of research on theoretical high-capacity anode materials. However, the research has been hindered because the high-capacity anode material properties and interactions are not well understood, largely due to the difficulty of observing cycling in situ. Using electrochemical scanning transmission electron microscopy (ec-STEM), we report the real-time observation and electrochemical analysis of pristine tin (Sn) and titanium dioxide-coated Sn (TiO₂@Sn) electrodes during lithiation/delithiation. As expected, we observed a volume expansion of the pristine Sn electrodes during lithiation, but we further observed that the expansion was followed by Sn detachment from the current collector. Remarkably, although the TiO₂@Sn electrodes also exhibited similar volume expansion during lithiation, they showed no evidence of Sn detachment. We found that the TiO₂ surface layer acted as an electrochemically activated artificial solid-electrolyte interphase that serves to conduct Li ions. As a physical coating, it mechanically prevented Sn detachment following volume changes during cycling, providing significant degradation resistance and 80% Coulombic efficiency for a complete lithiation/delithiation cycle. Interestingly, upon delithiation, TiO₂@Sn electrode displayed a self-healing mechanism of small pore formation in the Sn particle followed by agglomeration into several larger pores as delithiation continued.

Automatic mechanism generation is used to determine mechanisms for the CO₂ hydrogenation on Ni(111) in a two-stage process while considering the correlated uncertainty in DFT-based energetic parameters systematically. In a coarse stage, all the possible chemistry is explored with gas-phase products down to the ppb level, while a refined stage discovers the core methanation submechanism. Five thousand unique mechanisms were generated, which contain minor perturbations in all parameters. Global uncertainty assessment, global sensitivity analysis, and degree of rate control analysis are performed to study the effect of this parametric uncertainty on the microkinetic model predictions. Comparison of the model predictions with experimental data on a Ni/SiO₂ catalyst find a feasible set of microkinetic mechanisms within the correlated uncertainty space that are in quantitative agreement with the measured data, without relying on explicit parameter optimization. Global uncertainty and sensitivity analyses provide tools to determine the pathways and key factors that control the methanation activity within the parameter space. Together, these methods reveal that the degree of rate control approach can be misleading if parametric uncertainty is not considered. The procedure of considering uncertainties in the automated mechanism generation is not unique to CO₂ methanation and can be easily extended to other challenging heterogeneously catalyzed reactions.

Nanostructures with a high density of interfaces, such as in nanoporous materials and nanowires, resist radiation damage by promoting the annihilation and migration of defects. This study details the size effect and origins of the radiation damage mechanisms in nanowires and nanoporous structures in model face-centered (gold) and body-centered (niobium) cubic nanostructures using accelerated multi-cascade atomistic simulations and in-situ ion irradiation experiments. Our results reveal three different size-dependent mechanisms of damage accumulation in irradiated nanowires and nanoporous structures: sputtering for very small nanowires and ligaments, the formation and accumulation of point defects and dislocation loops in larger nanowires, and a face-centered-cubic to hexagonal-close-packed phase transformation for a narrow range of wire diameters in the case of gold nanowires. Smaller nanowires and ligaments have a net effect of lowering the radiation damage as compared to larger wires that can be traced back to the fact that smaller nanowires transition from a rapid accumulation of defects to a saturation and annihilation mechanism at a lower dose than larger nanowires. These irradiation damage mechanisms are accompanied with radiation-induced surface roughening resulting from defect-surface interactions. Comparisons between nanowires and nanoporous structures show that the various mechanisms seen in nanowires provide adequate bounds for the defect accumulation mechanisms in nanoporous structures with the difference attributed to the role of nodes connecting ligaments in nanoporous structures. Taken together, our results shed light on the compounded, size-dependent mechanisms leading to the radiation resistance of nanowires and nanoporous structures.

We describe a method to automatically generate an ion implantation recipe, a set of energies and fluences, to produce a desired defect density profile in a solid using the fewest required energies. We simulate defect density profiles for a range of ion energies, fit them with an appropriate function, and interpolate to yield defect density profiles at arbitrary ion energies. Given N energies, we then optimize a set of N energy-fluence pairs to match a given target defect density profile. Finally, we find the minimum N such that the error between the target defect density profile and the defect density profile generated by the N energy-fluence pairs is less than a given threshold. Inspired by quantum sensing applications with nitrogen-vacancy centers in diamond, we apply our technique to calculate optimal ion implantation recipes to create uniform-density 1 μm surface layers of 15N or vacancies (using 4He).

The trapping of paraffins is beneficial compared to selective olefin adsorption for adsorptive olefin purification from a process engineering point of view. Here we demonstrate the use of a series of Zn2(X-bdc)2(dabco) (where X-bdc2−is bdc2−= 1,4-benzenedicarboxylate with substituting groups X, DM-bdc2−= 2,5-dimethyl-1,4-benzenedicarboxylate or TM-bdc2−= 2,3,5,6-tetramethyl-1,4-benzenedicarboxylate and dabco = diazabicyclo[2.2.2.]octane) metal-organic frameworks (MOFs) for the adsorptive removal of ethane from ethylene streams. The best performing material from this series is Zn2(TM-bdc)2(dabco) (DMOF-TM), which shows a high ethane uptake of 5.31 mmol g−1at 110 kPa, with a good IAST selectivity of 1.88 towards ethane over ethylene. Through breakthrough measurements a high productivity of 13.1 L kg−1per breakthrough is revealed with good reproducibility over five consecutive cycles. Molecular simulations show that the methyl groups of DMOF-TM are forming a van der Waals trap with the methylene groups from dabco, snuggly fitting the ethane. Further, rarely used high pressure coadsorption measurements, in pressure regimes that most scientific studies on hydrocarbon separation on MOFs ignore, reveal an increase in ethane capacity and selectivity for binary mixtures with increased pressures. The coadsorption measurements reveal good selectivity of 1.96 at 1000 kPa, which is verified also through IAST calculations up to 3000 kPa. This study overall showcases the opportunities that pore engineering by alkyl group incorporation and pressure increase offer to improve hydrocarbon separation in reticular materials.

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Fuel Cycle Technology (FCT) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). Two high priorities for SFWST disposal R&D are design concept development and disposal system modeling. These priorities are directly addressed in the SFWST Geologic Disposal Safety Assessment (GDSA) control account, which is charged with developing a geologic repository system modeling and analysis capability, and the associated software, GDSA Framework, for evaluating disposal system performance for nuclear waste in geologic media. GDSA Framework is supported by SFWST Campaign and its predecessor the Used Fuel Disposition (UFD) campaign. This report fulfills the GDSA Uncertainty and Sensitivity Analysis Methods work package (SF-21SN01030404) level 3 milestone, Uncertainty and Sensitivity Analysis Methods and Applications in GDSA Framework (FY2021) (M3SF-21SN010304042). It presents high level objectives and strategy for development of uncertainty and sensitivity analysis tools, demonstrates uncertainty quantification (UQ) and sensitivity analysis (SA) tools in GDSA Framework in FY21, and describes additional UQ/SA tools whose future implementation would enhance the UQ/SA capability of GDSA Framework. This work was closely coordinated with the other Sandia National Laboratory GDSA work packages: the GDSA Framework Development work package (SF-21SN01030405), the GDSA Repository Systems Analysis work package (SF-21SN01030406), and the GDSA PFLOTRAN Development work package (SF-21SN01030407). This report builds on developments reported in previous GDSA Framework milestones, particularly M3SF 20SN010304032.

Organic materials are an attractive choice for structural components due to their light weight and versatility. However, because they decompose at low temperatures relative to tradiational materials they pose a safety risk due to fire and loss of structural integrity. To quantify this risk, analysts use chemical kinetics models to describe the material pyrolysis and oxidation using thermogravimetric analysis. This process requires the calibration of many model parameters to closely match experimental data. Previous efforts in this field have largely been limited to finding a single best-fit set of parameters even though the experimental data may be very noisy. Furthermore the chemical kinetics models are often simplified representations of the true de- composition process. The simplification induces model-form errors that the fitting process cannot capture. In this work we propose a methodology for calibrating decomposition models to thermogravimetric analysis data that accounts for uncertainty in the model-form and experimental data simultaneously. The methodology is applied to the decomposition of a carbon fiber epoxy composite with a three-stage reaction network and Arrhenius kinetics. The results show a good overlap between the model predictions and thermogravimetric analysis data. Uncertainty bounds capture devia- tions of the model from the data. The calibrated parameter distributions are also presented. In conclusion, the distributions may be used in forward propagation of uncertainty in models that leverage this material.

Designing polymers with controlled nanoscale morphologies and scalable synthesis is of great interest in the development of fluorine-free materials for proton-exchange membranes in fuel cells. This study focuses on a precision polyethylene with phenylsulfonic acid branches at every fifth carbon, p5PhSA, with a high ion-exchange capacity (4.2 mmol/g). The polymers self-assemble into hydrophilic and hydrophobic co-continuous nanoscale domains. In the hydrated state, the hydrophilic domain, composed of polar sulfonic acid moieties and water, serves as a pathway for efficient mesoscopic proton conductivity. The morphology and proton transport of p5PhSA are evaluated under hydrated conditions using in situ X-ray scattering and electrochemical impedance spectroscopy techniques. At 40 °C and 95% relative humidity, the proton conductivity of p5PhSA is 0.28 S/cm, which is four times greater than Nafion 117 under the same conditions. Atomistic molecular dynamics (MD) simulations are also used to elucidate the interplay between the structure and the water dynamics. The MD simulations show strong nanophase separation between the percolated hydrophilic and hydrophobic domains over a wide range of water contents. The percolated hydrophilic nanoscale domain facilitates the rapid proton transport in p5PhSA and demonstrates the potential of precise hydrocarbon-based polymers as processible and effective proton-exchange membranes.

Continued creation of harmful emissions such as NO_x and soot from compression-ignition engines utilizing mixing-controlled combustion systems (i.e., diesel engines) remains a problem and is the subject of on-going research. The inherently high efficiency, relatively low cost, and numerous other desirable attributes of such engines, coupled with a widely supported infrastructure, motivates their continued advancement. Recently, a scientifically distinct and mechanically simple technology called ducted fuel injection (DFI) has shown a robust ability to allow such engines to operate with simultaneously low engine-out soot and NO_x emissions when it is employed with simulated exhaust-gas recirculation. To better understand the property ranges of sustainable, oxygenated-fuel blending stocks that will most improve engine performance, two oxygenated blendstocks were separately blended with a commercial diesel base fuel and tested within a heavy-duty diesel optical engine equipped with a four-duct DFI configuration. Conventional and crank-angle-resolved optical diagnostics were used to elucidate the effects of fuel ignition quality, oxygenate molecular structure, and overall oxygen content on engine performance.

The representation of material heterogeneity (also referred to as "spatial variation") plays a key role in the material failure simulation method used in ALEGRA. ALEGRA is an arbitrary Lagrangian-Eulerian shock and multiphysics code developed at Sandia National Laboratories and contains several methods for incorporating spatial variation into simulations. A desirable property of a spatial variation method is that it should produce consistent stochastic behavior regardless of the mesh used (a property referred to as "mesh independence"). However, mesh dependence has been reported using the Weibull distribution with ALEGRA's spatial variation method. This report describes efforts towards providing additional insight into both the theory and numerical experiments investigating such mesh dependence. In particular, we have implemented a discrete minimum order statistic model with properties that are theoretically mesh independent.

Metallic enclosures are commonly used to protect electronic circuits against unwanted electromagnetic (EM) interactions. However, these enclosures may be sealed with imperfect mechanical seams or joints. These joints form narrow slots that allow external EM energy to couple into the cavity and then to the internal circuits. This coupled EM energy can severely affect circuit operations, particularly at the cavity resonance frequencies when the cavity has a high Q factor. To model these slots and the corresponding EM coupling, a thin-slot sub-cell model [1] , developed for slots in infinite ground plane and extended to numerical modeling of cavity-backed apertures, was successfully implemented in Sandia's electromagnetic code EIGER [2] and its next-generation counterpart Gemma [3]. However, this thin-slot model only considers resonances along the length of the slot. At sufficiently high frequencies, the resonances due to the slot depth must also be considered. Currently, slots must be explicitly meshed to capture these depth resonances, which can lead to low-frequency instability (due to electrically small mesh elements). Therefore, a slot sub-cell model that considers resonances in both length and depth is needed to efficiently and accurately capture the slot coupling.

This paper describe our team's experience using minimega, a network emulation system using node and network virtualization, to support evaluation of a set of networked and distributed systems for topology discovery, traffic classification and engineering in the DARPA Searchlight program [18]. We present the methodology we developed to encode network and traffic definitions into an experiment description model, and how our tools compile this model onto the underlying minimega API. We then present three cases studies which demonstrate the ability of our EDM to support experiments with diverse network topologies, diverse traffic mixes, and networks with specialized layer-2 connectivity requirements. We conclude with the overall takeaways from using minimega to support our evaluation process.

Heteroepitaxy of β-phase gallium oxide (β-Ga₂O₃) thin films on foreign substrates shows promise for the development of next-generation deep ultraviolet solar blind photodetectors and power electronic devices. In this work, the influences of the film thickness and crystallinity on the thermal conductivity of ($\bar{2}01$)-oriented β-Ga₂O₃ heteroepitaxial thin films were investigated. Unintentionally doped β-Ga₂O₃ thin films were grown on c-plane sapphire substrates with off-axis angles of 0° and 6° toward $\langle$$11\bar{2}0$$\rangle$ via metal–organic vapor phase epitaxy (MOVPE) and low-pressure chemical vapor deposition. The surface morphology and crystal quality of the β-Ga₂O₃ thin films were characterized using scanning electron microscopy, X-ray diffraction, and Raman spectroscopy. The thermal conductivities of the β-Ga₂O₃ films were measured via time-domain thermoreflectance. The interface quality was studied using scanning transmission electron microscopy. The measured thermal conductivities of the submicron-thick β-Ga₂O₃ thin films were relatively low as compared to the intrinsic bulk value. The measured thin film thermal conductivities were compared with the Debye–Callaway model incorporating phononic parameters derived from first-principles calculations. The comparison suggests that the reduction in the thin film thermal conductivity can be partially attributed to the enhanced phonon-boundary scattering when the film thickness decreases. They were found to be a strong function of not only the layer thickness but also the film quality, resulting from growth on substrates with different offcut angles. Growth of β-Ga₂O₃ films on 6° offcut sapphire substrates was found to result in higher crystallinity and thermal conductivity than films grown on on-axis c-plane sapphire. However, the β-Ga₂O₃ films grown on 6° offcut sapphire exhibit a lower thermal boundary conductance at the β-Ga₂O₃/sapphire heterointerface. In addition, the thermal conductivity of MOVPE-grown ($\bar{2}01$)-oriented β-(Al_xGa_1–x)₂O₃ thin films with Al compositions ranging from 2% to 43% was characterized. Because of phonon-alloy disorder scattering, the β-(Al_xGa_1–x)₂O₃ films exhibit lower thermal conductivities (2.8–4.7 W/m∙K) than the β-Ga₂O₃ thin films. The dominance of the alloy disorder scattering in β-(Al_xGa_1–x)₂O₃ is further evidenced by the weak temperature dependence of the thermal conductivity. This work provides fundamental insight into the physical interactions that govern phonon transport within heteroepitaxially grown β-phase Ga₂O₃ and (Al_xGa_1–x)₂O₃ thin films and lays the groundwork for the thermal modeling and design of β-Ga₂O₃ electronic and optoelectronic devices.

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (Tg) and chain length and thus have the potential for phase separation. One crosslinker has a high Tg corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low Tg. A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

We study the deformation of tantalum under extreme loading conditions. Experimental velocity data are drawn from both ramp loading experiments on Sandia's Z-machine and gas gun compression experiments. The drive conditions enable the study of materials under pressures greater than 100 GPa. We provide a detailed forward model of the experiments including a model of the magnetic drive for the Z-machine. Utilizing these experiments, we simultaneously infer several different types of physically motivated parameters describing equation of state, plasticity, and anelasticity via the computational device of Bayesian model calibration. Characteristics of the resulting calculated posterior distributions illustrate relationships among the parameters of interest via the degree of cross correlation. The calibrated velocity traces display good agreement with the experiments up to experimental uncertainty as well as improvement over previous calibrations. Examining the Z-shots and gun-shots together and separately reveals a trade-off between accuracy and transferability across different experimental conditions. Implications for model calibration, limitations from model form, and suggestions for improvements are discussed.

Electric vehicles (EVs) represent an important socio-economic development opportunity for islands and remote locations because they can lead to reduced fuel imports, electricity storage, grid services, and environmental and health benefits. This paper presents an overview of opportunities, challenges, and examples of EVs in islands and remote power systems, and is meant to provide background to researchers, utilities, energy offices, and other stakeholders interested in the impacts of electrification of transportation. The impact of uncontrolled EV charging on the electric grid operation is discussed, as well as several mitigation strategies. Of particular importance in many islands and remote systems is taking advantage of local resources by combining renewable energy and EV charging. Policy and economic issues are presented, with emphasis on the need for an overarching energy policy to guide the strategies for EVs growth. The key conclusion of this paper is that an orderly transition to EVs, one that maximizes benefits while addressing the challenges, requires careful analysis and comprehensive planning.

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

Superionic phases of bulk anhydrous salts based on large cluster-like polyhedral (carba)borate anions are generally stable only well above room temperature, rendering them unsuitable as solid-state electrolytes in energy-storage devices that typically operate at close to room temperature. To unlock their technological potential, strategies are needed to stabilize these superionic properties down to subambient temperatures. One such strategy involves altering the bulk properties by confinement within nanoporous insulators. In the current study, the unique structural and ion dynamical properties of an exemplary salt, NaCB11H12, nanodispersed within porous, high-surface-area silica via salt-solution infiltration were studied by differential scanning calorimetry, X-ray powder diffraction, neutron vibrational spectroscopy, nuclear magnetic resonance, quasielastic neutron scattering, and impedance spectroscopy. Combined results hint at the formation of a nanoconfined phase that is reminiscent of the high-temperature superionic phase of bulk NaCB11H12, with dynamically disordered CB11H12- anions exhibiting liquid-like reorientational mobilities. However, in contrast to this high-temperature bulk phase, the nanoconfined NaCB11H12 phase with rotationally fluid anions persists down to cryogenic temperatures. Moreover, the high anion mobilities promoted fast-cation diffusion, yielding Na+ superionic conductivities of ∼0.3 mS/cm at room temperature, with higher values likely attainable via future optimization. It is expected that this successful strategy for conductivity enhancement could be applied as well to other related polyhedral (carba)borate-based salts. Thus, these results present a new route to effectively utilize these types of superionic salts as solid-state electrolytes in future battery applications.

Poly(carbon monofluoride), or (CF)n, is a layered fluorinated graphite material consisting of nanosized platelets. Here, we present experimental multidimensional solid-state NMR spectra of (CF)n, supported by density functional theory (DFT) calculations of NMR parameters, which overhauls our understanding of structure and bonding in the material by elucidating many ways in which disorder manifests. We observe strong 19F NMR signals conventionally assigned to elongated or "semi-ionic"C-F bonds and find that these signals are in fact due to domains where the framework locally adopts boat-like cyclohexane conformations. We calculate that C-F bonds are weakened but are not elongated by this conformational disorder. Exchange NMR suggests that conformational disorder avoids platelet edges. We also use a new J-resolved NMR method for disordered solids, which provides molecular-level resolution of highly fluorinated edge states. The strings of consecutive difluoromethylene groups at edges are relatively mobile. Topologically distinct edge features, including zigzag edges, crenellated edges, and coves, are resolved in our samples by solid-state NMR. Disorder should be controllable in a manner dependent on synthesis, affording new opportunities for tuning the properties of graphite fluorides.

Finding dense regions of graphs is fundamental in graph mining. We focus on the computation of dense hierarchies and regions with graph nuclei - -a generalization of k-cores and trusses. Static computation of nuclei, namely through variants of 'peeling', are easy to understand and implement. However, many practically important graphs undergo continuous change. Dynamic algorithms, maintaining nucleus computations on dynamic graph streams, are nuanced and require significant effort to port between nuclei, e.g., from k-cores to trusses. We propose a unifying framework to maintain nuclei in dynamic graph streams. First, we show no dynamic algorithm can asymptotically beat re-computation, highlighting the need to experimentally understand variability. Next, we prove equivalence between k-cores on a special hypergraph and nuclei. Our algorithm splits the problem into maintaining the special hypergraph and maintaining k-cores on it. We implement our algorithm and experimentally demonstrate improvements up to 108 x over re-computation. We show algorithmic improvements on k-cores apply to trusses and outperform truss-specific implementations.

Costs to permit Marine Energy projects are poorly understood. In this paper we examine environmental compliance and permitting costs for 19 projects in the U.S., covering the last 2 decades. Guided discussions were conducted with developers over a 3-year period to obtain historical and ongoing project cost data relative to environmental studies (e.g., baseline or pre-project site characterization as well as post-installation effects monitoring), stakeholder outreach, and mitigation, as well as qualitative experience of the permitting process. Data are organized in categories of technology type, permitted capacity, pre-and post-installation, geographic location, and funding types. We also compare our findings with earlier logic models created for the Department of Energy (i.e., Reference Models). Environmental studies most commonly performed were for Fish and Fisheries, Noise, Marine Habitat/Benthic Studies and Marine Mammals. Studies for tidal projects were more expensive than those performed for wave projects and the range of reported project costs tended to be wider than ranges predicted by logic models. For eight projects reporting full project costs, from project start to FERC or USACE permit, the average amount for environmental permitting compliance was 14.6%.

Recently, lithium nitride (Li3N) has been proposed as a chemical warfare agent (CWA) neutralization reagent for its ability to produce nucleophilic ammonia molecules and hydroxide ions in aqueous solution. Quantum chemical calculations can provide insight into the Li3N neutralization process that has been studied experimentally. Here, we calculate reaction-free energies associated with the Li3N-based neutralization of the CWA VX using quantum chemical density functional theory and ab initio methods. We find that alkaline hydrolysis is more favorable to either ammonolysis or neutral hydrolysis for initial P-S and P-O bond cleavages. Reaction-free energies of subsequent reactions are calculated to determine the full reaction pathway. Notably, products predicted from favorable reactions have been identified in previous experiments.

This document is a guide to setting the system and processing configuration for the Geophysical Monitoring System (GMS) Station State-of-Health (SOH) Monitoring and Interactive Analysis (IAN) applications.

This document is a user’s manual for the Geophysical Monitoring System (GMS) Station State of Health (SOH) Monitoring User Interface.

This document is a user’s manual for the Geophysical Monitoring System (GMS) Interactive Analysis (IAN) User Interface.

High-throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. We develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community-level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.

Coupled Mode Theory (CMT) is a classic model that addresses many diverse problems in physics and electrical engineering. Although CMT is well-established in areas such as directional coupler design, its significance in antenna engineering is less well-recognized. Recently, Characteristic Mode Analysis (CMA) has shown how CMT quantitatively models a wide variety of multi-mode patch structures. In this paper, we extend these ideas to a double-tuned slot antenna utilizing some unique features of the method of moments code FEKO. The modeled slot antenna has a stable radiation pattern throughout its 50% impedance bandwidth. CMA decomposition indicates the antenna consists of two resonators: the slot and a lumped LC resonator that represents the feed probe. The two are electromagnetically coupled to produce broad impedance bandwidth. An equivalent circuit whose values are derived directly from the CMA results is given.

Degradation and partial shading impact the long-term reliability and power production of photovoltaic (PV) modules and power plants. Time-series power (Pmp) and current–voltage (I-V) curve datastreams from PV modules enable a remote diagnostic approach to quantify active degradation mechanisms and identify partial shading. We study three to nine years of these datastreams, including 3.6 million I-V curves and 36 million Pmp values, from eight PV modules, four each of double-glass and glass-backsheet module architectures, located in three distinctly different Köppen-Geiger climate zones, to determine the module's performance loss rates (PLR), identify active degradation mechanisms and power loss modes, along with partial shading by local objects. Considering both module architectures, PLR results indicate that the BSh climate zone is the most aggressive for module degradation, while the Alpine ET zone is the mildest climate. PLR of double-glass modules located in BWh and BSh climate zones are different due to the significantly greater uniform current loss (ΔPIsc) for double-glass modules in BSh, at a 5% significance level. Power loss for four out of five modules located in the BWh and BSh climates are dominated by uniform current degradation. Statistical analysis of multistep I-V curves detects partial shading experienced by three studied modules with details of the shading profile, the shading Poynting vector diagram for the obstacle's relative position, shading scenarios, and duration. This work demonstrates how remote monitoring and diagnosis of Pmp & I-V time-series of modules can provide quantitative operations and maintenance insights into system performance, degradation mechanisms, and shading.

Zirconium-based metal-organic frameworks, including UiO-66 and related frameworks, have become the focus of considerable research in the area of chemical warfare agent (CWA) decontamination. However, little work has been reported exploring these metal-organic frameworks (MOFs) for CWA sensing applications. For many sensing approaches, the growth of high-quality thin films of the active material is required, and thin film growth methods must be compatible with complex device architectures. Several approaches to synthesize thin films of UiO-66 have been described but many of these existing methods are complex or time consuming. We describe the development of a simple and rapid microwave assisted synthesis of oriented UiO-66 thin films on unmodified silicon (Si) and gold (Au) substrates. Thin films of UiO-66 and UiO-66-NH2 can be grown in as little as 2 min on gold substrates and 30 min on Si substrates. The film morphology and orientation are characterized and the effects of reaction time and temperature on thin film growth on Au are investigated. Both reaction time and temperature impact the overgrowth of protruding discrete crystallites in the thin film layer but, surprisingly, no strong correlation is observed between film thickness and reaction time or temperature. We also briefly describe the synthesis of Zr/Ce solid solution thin films of UiO-66 on Au and report the first synthesis of a solid solution thin film MOF. Finally, we demonstrate the utility of the microwave method for the facile functionalization of two sensor architectures, plasmonic nanohole arrays and microresonators, with UiO-66 thin films.

Understanding the role of physical processes contributing to breakdown is critical for many applications in which breakdown is undesirable, such as capacitors, and applications in which controlled breakdown is intended, such as plasma medicine, lightning protection, and materials processing. The electron emission from the cathode is a critical source of electrons which then undergo impact ionization to produce electrical breakdown. In this study, the role of secondary electron yields due to photons (γ ph) and ions (γ i) in direct current breakdown is investigated using a particle-in-cell direct simulation Monte Carlo model. The plasma studied is a one-dimensional discharge in 50 Torr of pure helium with a platinum cathode, gap size of 1.15 cm, and voltages of 1.2-1.8 kV. The current traces are compared with experimental measurements. Larger values of γ ph generally result in a faster breakdown, while larger values of γ i result in a larger maximum current. The 58.4 nm photons emitted from He(21P) are the primary source of electrons at the cathode before the cathode fall is developed. Of the values of γ ph and γ i investigated, those which provide the best agreement with the experimental current measurements are γ ph = 0.005 and γ i = 0.01. These values are significantly lower than those in the literature for pristine platinum or for a graphitic carbon film which we speculate may cover the platinum. This difference is in part due to the limitations of a one-dimensional model but may also indicate surface conditions and exposure to a plasma can have a significant effect on the secondary electron yields. The effects of applied voltage and the current produced by a UV diode which was used to initiate the discharge, are also discussed.

Concerns about cyber threats to space systems are increasing. Researchers are developing intrusion detection and protection systems to mitigate these threats, but sparsity of cyber threat data poses a significant challenge to these efforts. Development of credible threat data sets are needed to overcome this challenge. This paper describes the extension/development of three data generation algorithms (generative adversarial networks, variational auto-encoders, and generative algorithm for multi-variate timeseries) to generate cyber threat data for space systems. The algorithms are applied to a use case that leverages the NASA Operational Simulation for Small Satellites (NOS$^{3})$ platform. Qualitative and quantitative measures are applied to evaluate the generated data. Strengths and weaknesses of each algorithm are presented, and suggested improvements are provided. For this use case, generative algorithm for multi-variate timeseries performed best according to both qualitative and quantitative measures.

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. This introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

Concerns about cyber threats to space systems are increasing. Researchers are developing intrusion detection and protection systems to mitigate these threats, but sparsity of cyber threat data poses a significant challenge to these efforts. Development of credible threat data sets are needed to overcome this challenge. This paper describes the extension/development of three data generation algorithms (generative adversarial networks, variational auto-encoders, and generative algorithm for multi-variate timeseries) to generate cyber threat data for space systems. The algorithms are applied to a use case that leverages the NASA Operational Simulation for Small Satellites (NOS$^{3})$ platform. Qualitative and quantitative measures are applied to evaluate the generated data. Strengths and weaknesses of each algorithm are presented, and suggested improvements are provided. For this use case, generative algorithm for multi-variate timeseries performed best according to both qualitative and quantitative measures.

In each of our brains, 86 billion neurons work in parallel, processing inputs from senses and memories to produce the many feats of human cognition. The brains of other creatures are less broadly capable, but those animals often exhibit innate aptitudes for particular tasks, abilities honed by millions of years of evolution.

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. The model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

Abstract: Uptake of noble gases into heterogeneous geologic core samples was measured using a piezometric methodology. In addition to measuring accessible porosity—as with gas pycnometry—by monitoring the rate of pressure decay, this method can also be used to estimate the gas effective diffusivity in the sample. In contrast to previous applications of this method, where milligram quantities of fractured grains are characterized, here approximately kilogram core samples were left intact when tested. In doing so, a more representative sample of the heterogeneous field geology is provided. Additionally, alteration of the pore structure and connectivity during sample preparation is avoided. To scale the piezometric method from milligrams to kilograms, the system was designed to operate at medium vacuum (1 to 100 Pa) to restrict transport in pores less than approximately 60 µm to large Knudsen numbers. To test the system performance, two samples of interest were selected: a rhyolitic welded tuff from Blue Canyon Dome at the Energetic Materials Research and Testing Center and a zeolitized non-welded rhyolitic tuff from the Nevada National Security Site. Three noble gases were utilized in this test series; Argon and xenon as they are of direct interest to nuclear monitoring efforts and helium as it is a weakly adsorbing reference standard. Additionally, mercury intrusion porosimetry measurements were made on subsamples of the core to compare the observed porosity by the two methods and to discuss gas transport rates in the context of the measured pore distribution. Article Highlights: The piezometric method was extended to measure transport in intact geologic core samples between 800 and 1400 g.Transport in the pores spaces was restricted to Knudsen flow using medium vacuum, enabling a closed-form solution.Argon and xenon in a zeolitized tuff core exhibited significant adsorption and enhanced transport relative to helium.

In situ TEM/SEM microscopy to investigate the structural evolution that occurs due to various extreme environments

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. The coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier-Stokes/eddy dissipation concept-based methodology.

We discuss thinned InAsSb resonant infrared detectors that are designed to enable high quantum efficiency by using interleaved nanoantennas to read out two wavelengths from each pixel simultaneously.

In previous studies using the University of Nevada, Reno's (UNR's) high-impedance Zebra Marx generator (1.9 ω, 1.7 MA, 100 ns), Double Planar Wire Arrays (DPWAs) proved to be excellent radiators, and Double Planar Foil Liners (DPFLs) proved useful for future inertial confinement fusion applications. This article presents the results of joint UNR/UM (University of Michigan) experiments with aluminum (Al) DPWAs, Al DPFLs, and tungsten (W) DPWAs using UM's Michigan Accelerator for Inductive Z-Pinch Experiments (MAIZE) generator, a low-impedance Linear Transformer Driver (LTD) (0.1 ω, 0.5-1 MA, and 100-250 ns). The main goals of this study were twofold: the first was a pioneering effort to test whether a relatively heavy Al DPFL could successfully be imploded on a low-impedance university-scale LTD like the MAIZE generator, and, if so, to analyze the results and make comparisons to the optimized, lighter DPWA configurations that have been previously studied. The DPWAs consisted of two planes of micrometer-scale diameter Al or W wires, while the DPFLs consisted of two planes of micrometer-scale thickness Al foils. Diagnostics include filtered Si-diodes, an absolutely calibrated filtered PCD, x-ray pinhole cameras, spectrometers, and gated optical self-emission imaging. The implosion dynamics and radiative properties of Al DPWAs and DPFLs and W DPWAs on the MAIZE LTD are discussed and compared. Time-dependent load inductance calculations derived from measurements of the load current and a MAIZE circuit model provide a relative measurement of pinch strength. In experiments on MAIZE, W planar wire arrays exhibited a higher peak load inductance throughout the pinch than Al DPWAs and DPFLs, while x-ray pulses from Al DPFLs had the longest emission duration.

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.