Publications

A stable solid electrolyte interphase (SEI) layer is key to high performing lithium ion and lithium metal batteries for metrics such as calendar and cycle life. The SEI must be mechanically robust to withstand large volumetric changes in anode materials such as lithium and silicon, so understanding the mechanical properties and behavior of the SEI is essential for the rational design of artificial SEI and anode form factors. The mechanical properties and mechanical failure of the SEI are challenging to study, because the SEI is thin at only ∼10-200 nm thick and is air sensitive. Furthermore, the SEI changes as a function of electrode material, electrolyte and additives, temperature, potential, and formation protocols. A variety of in situ and ex situ techniques have been used to study the mechanics of the SEI on a variety of lithium ion battery anode candidates; however, there has not been a succinct review of the findings thus far. Because of the difficulty of isolating the true SEI and its mechanical properties, there have been a limited number of studies that can fully de-convolute the SEI from the anode it forms on. A review of past research will be helpful for culminating current knowledge and helping to inspire new innovations to better quantify and understand the mechanical behavior of the SEI. This review will summarize the different experimental and theoretical techniques used to study the mechanics of SEI on common lithium battery anodes and their strengths and weaknesses.

Targeting host factors for anti-viral development offers several potential advantages over traditional countermeasures that include broad-spectrum activity and prevention of resistance. Characterization of host factors in animal models provides strong evidence of their involvement in disease pathogenesis, but the feasibility of performing high-throughput in vivo analyses on lists of genes is problematic. To begin addressing the challenges of screening candidate host factors in vivo, we combined advances in CRISPR-Cas9 genome editing with an immunocompromised mouse model used to study highly pathogenic viruses. Transgenic mice harboring a constitutively expressed Cas9 allele (Cas9tg/tg) with or without knockout of type I interferon receptors served to optimize in vivo delivery of CRISPR single-guide RNA (sgRNA) using Invivofectamine 3.0, a simple and easy-to-use lipid nanoparticle reagent. Invivofectamine 3.0-mediated liver-specific editing to remove activity of the critical Ebola virus host factor Niemann-Pick disease type C1 in an average of 74% of liver cells protected immunocompromised Cas9tg/tg mice from lethal surrogate Ebola virus infection. We envision that immunocompromised Cas9tg/tg mice combined with straightforward sgRNA in vivo delivery will enable efficient host factor loss-of-function screening in the liver and other organs to rapidly study their effects on viral pathogenesis and help initiate development of broad-spectrum, host-directed therapies against emerging pathogens.

This report documents the feasibility phase of the Critical Experiment Design (CED) conducted as part of integral experiment request (IER) 523. The purpose of IER-523 is to explore the effects of using 35 weight percent enriched uranium dioxide-beryllium oxide (UO₂-BeO) material on critical configurations using the Seven Percent Critical Experiment (7uPCX) at Sandia National Laboratories (Sandia). Preliminary experiment design concepts, neutronic analysis results, and proposed paths for continuing the CED process are presented.

Understanding how atoms interact with hot dense matter is essential for astrophysical and laboratory plasmas. Interactions in high-density plasmas broaden spectral lines, providing a rare window into interactions that govern, for example, radiation transport in stars. However, up to now, spectral line-shape theories employed at least one of three common approximations: second-order Taylor treatment of broadening operator, dipole-only interactions between atom and plasma, and classical treatment of perturbing electrons. In this Letter, we remove all three approximations simultaneously for the first time and test the importance for two applications: neutral hydrogen and highly ionized magnesium and oxygen. We found 15%-50% change in the spectral line widths, which are sufficient to impact applications including white-dwarf mass determination, stellar-opacity research, and laboratory plasma diagnostics.

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

Abstract not provided.

Soft-elasticity in monodomain liquid crystal elastomers (LCEs) is promising for impact-absorbing applications where strain energy is ideally absorbed at constant stress. Conventionally, compressive and impact studies on LCEs have not been performed given the notorious difficulty synthesizing sufficiently large monodomain devices. Here, we use direct-ink writing 3D printing to fabricate bulk (>cm3) monodomain LCE devices and study their compressive soft-elasticity over 8 decades of strain rate. At quasi-static rates, the monodomain soft-elastic LCE dissipated 45% of strain energy while comparator materials dissipated less than 20%. At strain rates up to 3000 s−1, our soft-elastic monodomain LCE consistently performed closest to an ideal-impact absorber. Drop testing reveals soft-elasticity as a likely mechanism for effectively reducing the severity of impacts – with soft elastic LCEs offering a Gadd Severity Index 40% lower than a comparable isotropic elastomer. Lastly, we demonstrate tailoring deformation and buckling behavior in monodomain LCEs via the printed director orientation.

Fiber optic distributed acoustic sensors (DAS) are becoming a widely used tool for seismic sensing. Here we examine recordings of two subsurface chemical explosions, DAG-1 and DAG-3, each of which was about one metric ton (TNT equivalent), that were recorded from a helical fiber installed in two boreholes 80 m away from the source location. Several clear phases including the initial P wave, a weak S wave, and a surface reflected P wave are observed on the helical DAS data. We estimate a velocity model using arrival times measured from the fiber. The DAS waveform data were compared with colocated accelerometers at specific depths in both frequency and time domains. The spectra of the DAS data matched spectra estimated from the accelerometer records. Comparisons of observed waveform shape between the accelerometer records and the fiber measurements (strain-rate) show reasonable agreement except for the data near the event depth. The DAS data and the accelerometer agreed in relative amplitudes but we had difficulties in matching absolute amplitudes, possibly due to errors in metadata. Synthetic strain-rate waveforms were calculated using a 2D wavenumber algorithm and matched the waveform shape and relative amplitudes. In general, DAS is effective at recording strong ground motions at high spatial density. Comparison of the synthetic seismograms with observed data indicate that the waveforms are not consistent with a pure isotropic explosion source and that the observed S waves originate from very near the source region.

The Kronecker product is an important matrix operation with a wide range of applications in signal processing, graph theory, quantum computing and deep learning. In this work, we introduce a generalization of the fast Johnson–Lindenstrauss projection for embedding vectors with Kronecker product structure, the Kronecker fast Johnson–Lindenstrauss transform (KFJLT). The KFJLT reduces the embedding cost by an exponential factor of the standard fast Johnson–Lindenstrauss transform’s cost when applied to vectors with Kronecker structure, by avoiding explicitly forming the full Kronecker products. We prove that this computational gain comes with only a small price in embedding power: consider a finite set of p points in a tensor product of d constituent Euclidean spaces ⊗1k=d Rnk, and let N = Πdk=1 nk. With high probability, a random KFJLT matrix of dimension m × N embeds the set of points up to multiplicative distortion (1 ± ε) provided m ≿ ε−2 log2d−1 (p) log N. We conclude by describing a direct application of the KFJLT to the efficient solution of large-scale Kronecker-structured least squares problems for fitting the CP tensor decomposition.

Combining plasmonic and magnetic properties, namely magneto-plasmonic coupling, inspires great research interest and the search for magneto-plasmonic nanostructure becomes considerably critical. Here we designed a nanopillar-in-matrix structure with core–shell alloyed nanopillars for both BaTiO3 (BTO)-Au0.5Co0.5 (AuCo) and BTO-Au0.25Cu0.25Co0.25Ni0.25 (AuCuCoNi) hybrid systems, i.e., ferromagnetic alloy cores (e.g., Co or CoNi) with plasmonic shells (e.g., Au or Au/Cu). These core–shell alloy nanopillars are uniformly embedded into a dielectric BTO matrix to form a vertically aligned nanocomposite (VAN) structure. Both hybrid systems present excellent epitaxial quality and interesting multi-functionality, e.g., high magnetic anisotropy, magneto-optical coupling response, tailorable plasmonic resonance wavelength, tunable hyperbolic properties and strong optical anisotropy. These alloyed nanopillars-in-matrix designs provide enormous potential for complex hybrid material designs with multi-functionality and demonstrate strong interface enabled magneto-plasmonic coupling along with plasmonic and magnetic performance.

Abstract not provided.

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.

Drinking water systems commonly use manual or grab sampling to monitor water quality, identify or confirm issues, and verify that corrective or emergency response actions have been effective. In this paper, the effectiveness of regulatory sampling locations for emergency response is explored. An optimization formulation based on the literature was used to identify manual sampling locations to maximize overall nodal coverage of the system. Results showed that sampling locations could be effective in confirming incidents for which they were not designed. When evaluating sampling locations optimized for emergency response against regulatory scenarios, the average performance was reduced by 3%-4%, while using optimized regulatory sampling locations for emergency response reduced performance by 7%-10%. Secondary constraints were also included in the formulation to ensure geographical and water age diversity with minimal impact on the performance. This work highlighted that regulatory sampling locations provide value in responding to an emergency for these networks.

Abstract not provided.

Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.

Particle-in-cell, direct simulation Monte Carlo simulations reveal that ion-acoustic instabilities excited in presheaths can cause significant ion heating. Ion-acoustic instabilities are excited by the ion flow toward a sheath when the neutral gas pressure is small enough and the electron temperature is large enough. A series of 1D simulations were conducted in which neutral plasma (electrons and ions) was uniformly sourced with an ion temperature of 0.026 eV and different electron temperatures (0.1 eV-50 eV). Ion heating was observed when the electron-to-ion temperature ratio exceeded the minimum value predicted by linear response theory to excite ion-acoustic instabilities at the sheath edge (T e / T i ≈ 28). When this threshold was exceeded, the temperature equilibration rate between ions and electrons rapidly increased near the sheath so that the local temperature ratio did not significantly exceed the threshold for instability. This resulted in significant ion heating near the sheath edge, which also extended back into the bulk plasma; presumably due to wave reflection from the sheath. This ion-acoustic wave heating mechanism was found to decrease for higher neutral pressures, where ion-neutral collisions damp the ion-acoustic waves and ion heating is instead dominated by inelastic collisions in the presheath.

The atomic cluster expansion is a general polynomial expansion of the atomic energy in multi-atom basis functions. Here we implement the atomic cluster expansion in the performant C++ code PACE that is suitable for use in large-scale atomistic simulations. We briefly review the atomic cluster expansion and give detailed expressions for energies and forces as well as efficient algorithms for their evaluation. We demonstrate that the atomic cluster expansion as implemented in PACE shifts a previously established Pareto front for machine learning interatomic potentials toward faster and more accurate calculations. Moreover, general purpose parameterizations are presented for copper and silicon and evaluated in detail. We show that the Cu and Si potentials significantly improve on the best available potentials for highly accurate large-scale atomistic simulations.

Lignocellulosic biomass is composed of three major biopolymers: cellulose, hemicellulose and lignin. Analytical tools capable of quickly detecting both glycan and lignin deconstruction are needed to support the development and characterization of efficient enzymes/enzyme cocktails. Previously we have described nanostructure-initiator mass spectrometry-based assays for the analysis of glycosyl hydrolase and most recently an assay for lignin modifying enzymes. Here we integrate these two assays into a single multiplexed assay against both classes of enzymes and use it to characterize crude commercial enzyme mixtures. Application of our multiplexed platform based on nanostructure-initiator mass spectrometry enabled us to characterize crude mixtures of laccase enzymes from fungi Agaricus bisporus (Ab) and Myceliopthora thermophila (Mt) revealing activity on both carbohydrate and aromatic substrates. Using time-series analysis we determined that crude laccase from Ab has the higher GH activity and that laccase from Mt has the higher activity against our lignin model compound. Inhibitor studies showed a significant reduction in Mt GH activity under low oxygen conditions and increased activities in the presence of vanillin (common GH inhibitor). Ultimately, this assay can help to discover mixtures of enzymes that could be incorporated into biomass pretreatments to deconstruct diverse components of lignocellulosic biomass.

Abstract not provided.

The friction behavior of metals is directly linked to the mechanisms that accommodate deformation. We examine the links between mechanisms of strengthening, deformation, and the wide range of friction behaviors that are exhibited by shearing metal interfaces. Specifically, the focus is on understanding the shear strength of nanocrystalline and nanostructured metals, and conditions that lead to low friction coefficients. Grain boundary sliding and the breakdown of Hall–Petch strengthening at the shearing interface are found to generally and predictably explain the low friction of these materials. While the following is meant to serve as a general discussion of the strength of metals in the context of tribological applications, one important conclusion is that tribological research methods also provide opportunities for probing the fundamental properties and deformation mechanisms of metals.

Abstract not provided.

We report on the characterization of heating rates and photoinduced electric charging on a microfabricated surface ion trap with integrated waveguides. Microfabricated surface ion traps have received considerable attention as a quantum information platform due to their scalability and manufacturability. Here, we characterize the delivery of 435-nm light through waveguides and diffractive couplers to a single ytterbium ion in a compact trap. We measure an axial heating rate at room temperature of 0.78±0.05 q/ms and see no increase due to the presence of the waveguide. Furthermore, the electric field due to charging of the exposed dielectric outcoupler settles under normal operation after an initial shift. The frequency instability after settling is measured to be 0.9 kHz.

Abstract not provided.

The study of microstructure evolution in additive manufacturing of metals would be aided by knowing the thermal history. Since temperature measurements beneath the surface are difficult, estimates are obtained from computational thermo-mechanical models calibrated against traces left in the sample revealed after etching, such as the trace of the melt pool boundary. Here we examine the question of how reliable thermal histories computed from a model that reproduces the melt pool trace are. To this end, we perform experiments in which one of two different laser beams moves with constant velocity and power over a substrate of 17-4PH SS or Ti-6Al-4V, with low enough power to avoid generating a keyhole. We find that thermal histories appear to be reliably computed provided that (a) the power density distribution of the laser beam over the substrate is well characterized, and (b) convective heat transport effects are accounted for. Poor control of the laser beam leads to potentially multiple three-dimensional melt pool shapes compatible with the melt pool trace, and therefore to multiple potential thermal histories. Ignoring convective effects leads to results that are inconsistent with experiments, even for the mild melt pools here.

Computational fluid dynamics (CFD) models were developed to simulate the impact of different ventilation scenarios on airborne exposure risks in a 72-passenger school bus. Scenarios and factors that were investigated included a moving vs. stationary bus, impacts of a heating unit within the bus, and impacts of alternative ventilation scenarios with different combinations of openings (e.g., windows, door, emergency hatch). Results of the simulations showed that when the bus was stationary, use of the heater increased receptor concentrations unless there was another opening. When the bus was moving, simulations with at least two sets of openings separated from each other in the forward and aft directions produced a through-flow condition that reduced concentrations via dilution from outside air by a factor of ten or more. A single opening in a moving bus generally increased concentrations throughout the cabin due to increased mixing with minimal ventilation. The cumulative exposure risk (time-averaged concentrations) was found to be inversely correlated to the air exchange rate. Stationary and moving-bus scenarios that yielded above ~20 air changes per hour resulted in the lowest cumulative exposures. Recommendations from this study were implemented in new safety and operating procedures by the Albuquerque Public Schools Transportation Center.

The Message Passing Interface (MPI) remains the dominant programming model for scientific applications running on today's high-performance computing (HPC) systems. This dominance stems from MPI's powerful semantics for inter-process communication that has enabled scientists to write applications for simulating important physical phenomena. MPI does not, however, specify how messages and synchronization should be carried out. Those details are typically dependent on low-level architecture details and the message characteristics of the application. Therefore, analyzing an application's MPI resource usage is critical to tuning MPI's performance on a particular platform. The result of this analysis is typically a discussion of the mean message sizes, queue search lengths and message arrival times for a workload or set of workloads. While a discussion of the arithmetic mean in MPI resource usage might be the most intuitive summary statistic, it is not always the most accurate in terms of representing the underlying data. In this paper, we analyze MPI resource usage for a number of key MPI workloads using an existing MPI trace collector and discrete-event simulator. Our analysis demonstrates that the average, while easy and efficient to calculate, is a useful metric for characterizing latency and bandwidth measurements, but may not be a good representation of application message sizes, match list search depths, or MPI inter-operation times. Additionally, we show that the median and mode are superior choices in many cases. We also observe that the arithmetic mean is not the best representation of central tendency for data that are drawn from distributions that are multi-modal or have heavy tails. The results and analysis of our work provide valuable guidance on how we, as a community, should discuss and analyze MPI resource usage data for scientific applications.

We compare the performance of 10 and 5 μm thick metal-metal waveguide terahertz quantum-cascade laser ridges operating around 2.7 THz and based on a 4-well phonon depopulation active region design. Thanks to reduced heat dissipation and lower thermal resistance, the 5 μm thick material shows an 18 K increase in continuous wave operating temperature compared to the 10 μm material, despite a lower maximum pulsed-mode operating temperature and a larger input power density. A maximum continuous wave operating temperature of 129 K is achieved using the 5 μm thick material and a 15 μm wide ridge waveguide, which lased up to 155 K in the pulsed mode. The use of thin active regions is likely to become increasingly important to address the increasing input power density of emerging 2- and 3-well active region designs that show the highest pulsed operating temperatures.

High-quality factor resonant cavities are challenging structures to model in electromagnetics owing to their large sensitivity to minute parameter changes. Therefore, uncertainty quantification (UQ) strategies are pivotal to understanding key parameters affecting the cavity response. We discuss here some of these strategies focusing on shielding effectiveness (SE) properties of a canonical slotted cylindrical cavity that will be used to develop credibility evidence in support of predictions made using computational simulations for this application.

The addition of active, nonlinear, and nonreciprocal functionalities to passive piezoelectric acoustic wave technologies could enable all-acoustic and therefore ultra-compact radiofrequency signal processors. Toward this goal, we present a heterogeneously integrated acoustoelectric material platform consisting of a 50 nm indium gallium arsenide epitaxial semiconductor film in direct contact with a 41° YX lithium niobate piezoelectric substrate. We then demonstrate three of the main components of an all-acoustic radiofrequency signal processor: passive delay line filters, amplifiers, and circulators. Heterogeneous integration allows for simultaneous, independent optimization of the piezoelectric-acoustic and electronic properties, leading to the highest performing surface acoustic wave amplifiers ever developed in terms of gain per unit length and DC power dissipation, as well as the first-ever demonstrated acoustoelectric circulator with an isolation of 46 dB with a pulsed DC bias. Finally, we describe how the remaining components of an all-acoustic radiofrequency signal processor are an extension of this work.

The isotopic composition of water in Earth’s oceans is challenging to recreate using a plausible mixture of known extraterrestrial sources such as asteroids—an additional isotopically light reservoir is required. The Sun’s solar wind could provide an answer to balance Earth’s water budget. We used atom probe tomography to directly observe an average ~1 mol% enrichment in water and hydroxyls in the solar-wind-irradiated rim of an olivine grain from the S-type asteroid Itokawa. We also experimentally confirm that H+ irradiation of silicate mineral surfaces produces water molecules. These results suggest that the Itokawa regolith could contain ~20 l m−3 of solar-wind-derived water and that such water reservoirs are probably ubiquitous on airless worlds throughout our Galaxy. The production of this isotopically light water reservoir by solar wind implantation into fine-grained silicates may have been a particularly important process in the early Solar System, potentially providing a means to recreate Earth’s current water isotope ratios.

Abstract not provided.

A resonant cavity undergoes three distinct behaviors with increasing frequency: 1) fundamental modes, localized in frequency with well defined modal distribution; 2) undermoded region, where modes are still separated, but are sufficiently perturbed by small imperfections that their spectral positions (and distributions) are statistical in nature; and 3) overmoded region, where modes overlap, field distributions follow stochastic distributions, and the slot acts as if in free space. Understanding the penetration through slots in the overmoded region is of great interest, and is the focus of this article. Since full-wave solvers may not be able to provide a timely answer for very high frequencies due to a lack of memory and/or computation resources, we develop bounding methods to estimate worst-case average and maximum fields within the cavity. After discussing the bounding formulation, we compare its results to full-wave simulations at the first, second, and third resonance supported by the slot in the case of a cylindrical cavity. Note that the bounding formulation indicates that results are nearly independent of cavity shape: only the cavity volume, frequency, and cavity quality factor affect the overmoded region, making this formulation a powerful tool to assess electromagnetic interference and electromagnetic compatibility effects within cavities.

Diffraction techniques can powerfully and nondestructively probe materials while maintaining high resolution in both space and time. Unfortunately, these characterizations have been limited and sometimes even erroneous due to the difficulty of decoding the desired material information from features of the diffractograms. Currently, these features are identified non-comprehensively via human intuition, so the resulting models can only predict a subset of the available structural information. In the present work we show (i) how to compute machine-identified features that fully summarize a diffractogram and (ii) how to employ machine learning to reliably connect these features to an expanded set of structural statistics. To exemplify this framework, we assessed virtual electron diffractograms generated from atomistic simulations of irradiated copper. When based on machine-identified features rather than human-identified features, our machine-learning model not only predicted one-point statistics (i.e. density) but also a two-point statistic (i.e. spatial distribution) of the defect population. Hence, this work demonstrates that machine-learning models that input machine-identified features significantly advance the state of the art for accurately and robustly decoding diffractograms.

The essential data for interior and thermal evolution models of the Earth and super-Earths are the density and melting of mantle silicate under extreme conditions. Here, we report an unprecedently high melting temperature of MgSiO3 at 500 GPa by direct shockwave loading of pre-synthesized dense MgSiO3 (bridgmanite) using the Z Pulsed Power Facility. We also present the first high-precision density data of crystalline MgSiO3 to 422 GPa and 7200 K and of silicate melt to 1254 GPa. The experimental density measurements support our density functional theory based molecular dynamics calculations, providing benchmarks for theoretical calculations under extreme conditions. The excellent agreement between experiment and theory provides a reliable reference density profile for super-Earth mantles. Furthermore, the observed upper bound of melting temperature, 9430 K at 500 GPa, provides a critical constraint on the accretion energy required to melt the mantle and the prospect of driving a dynamo in massive rocky planets.

In the past decade, symmetry-protected bound states in the continuum (BICs) have proven to be an important design principle for creating and enhancing devices reliant upon states with high-quality (Q) factors, such as sensors, lasers, and those for harmonic generation. However, as we show, current implementations of symmetry-protected BICs in photonic crystal slabs can only be found at the center of the Brillouin zone and below the Bragg diffraction limit, which fundamentally restricts their use to single-frequency applications. By microprinting a three-dimensional (3D) photonic crystal structure using two-photon polymerization, we demonstrate that this limitation can be overcome by altering the radiative environment surrounding the slab to be a 3D photonic crystal. This allows for the protection of a line of BICs by embedding it in a symmetry bandgap of the crystal. This concept substantially expands the design freedom available for developing next-generation devices with high-Q states.

In-silico screening of novel biofuel molecules based on chemical and fuel properties is a critical first step in the biofuel evaluation process due to the significant volumes of samples required for experimental testing, the destructive nature of engine tests, and the costs associated with bench-scale synthesis of novel fuels. Predictive models are limited by training sets of few existing measurements, often containing similar classes of molecules that represent just a subset of the potential molecular fuel space. Software tools can be used to generate every possible molecular descriptor for use as input features, but most of these features are largely irrelevant and training models on datasets with higher dimensionality than size tends to yield poor predictive performance. Feature selection has been shown to improve machine learning models, but correlation-based feature selection fails to provide scientific insight into the underlying mechanisms that determine structure–property relationships. The implementation of causal discovery in feature selection could potentially inform the biofuel design process while also improving model prediction accuracy and robustness to new data. In this study, we investigate the benefits causal-based feature selection might have on both model performance and identification of key molecular substructures. We found that causal-based feature selection performed on par with alternative filtration methods, and that a structural causal model provides valuable scientific insights into the relationships between molecular substructures and fuel properties.

The solid-state joining of oxide-dispersion-strengthened (ODS) austenitic steels was achieved using a pulsed electric current joining (PECJ) process. Microstructures of the austenitic grain structures and oxide dispersions in the joint areas were characterized using electron microscopy. Negligible grain growth was observed in austenitic grain structures, while slight coarsening of oxide dispersions occurred at a short holding time. The mechanisms of the PECJ process may involve three steps that occur simultaneously, including the sintering of mechanical alloying powders in the bonding layer, formation of oxide dispersions, and bonding of the mechanical alloying powders with the base alloy. The high hardness and irradiation resistance of ODS alloys were retained in the joint areas. This research revealed the fundamental mechanisms during the PECJ process, which is beneficial for its potential applications during the advanced manufacturing of ODS alloys.

Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.

Conventional engineering methods oftentimes have challenges in the quantification of crack nucleation processes from manufacturing defects that are relevant for engineering component lifing. We present peridynamic simulation framework for the description of the crack nucleation process from cluster of non-metallic inclusions in aluminum alloys. Our non-local simulation framework characterizes crack nucleation process as multiple micro-crack nucleation events from individual inclusions, and eventually one micro-crack dominates. We define individual stages of the crack nucleation process, i.e., nucleation, micro-crack, technical, and crack initiation, that allow a quantification and meta model development of both the individual stages and the entire crack nucleation process.

The earth texture with complex morphological geometry and compositions such as shale and carbonate rocks, is typically characterized with sparse field samples because of an expensive and time-consuming characterization process. Accordingly, generating arbitrary large size of the geological texture with similar topological structures at a low computation cost has become one of the key tasks for realistic geomaterial reconstruction and subsequent hydro-mechanical evaluation for science and engineering applications. Recently, generative adversarial neural networks (GANs) have demonstrated a potential of synthesizing input textural images and creating equiprobable geomaterial images for stochastic analysis of hydrogeological properties, for example, the feasibility of CO2 storage sites and exploration of unconventional resources. However, the texture synthesis with the GANs framework is often limited by the computational cost and scalability of the output texture size. In this study, we proposed a spatially assembled GANs (SAGANs) that can generate output images of an arbitrary large size regardless of the size of training images with computational efficiency. The performance of the SAGANs was evaluated with two and three dimensional (2D and 3D) rock image samples widely used in geostatistical reconstruction of the earth texture and Lattice-Boltzmann (LB) simulations were performed to compare pore-scale flow patterns and upscaled permeabilities of training and generated geomaterial images. We demonstrate SAGANs can generate the arbitrary large size of statistical realizations with connectivity and structural properties and flow characteristics similar to training images, and also can generate a variety of realizations even on a single training image. In addition, the computational time was significantly improved compared to standard GANs frameworks.

Layered boron compounds have attracted significant interest in applications from energy storage to electronic materials to device applications, owing in part to a diversity of surface properties tied to specific arrangements of boron atoms. Here we report the energy landscape for surface atomic configurations of MgB2 by combining first-principles calculations, global optimization, material synthesis and characterization. We demonstrate that contrary to previous assumptions, multiple disordered reconstructions are thermodynamically preferred and kinetically accessible within exposed B surfaces in MgB2 and other layered metal diborides at low boron chemical potentials. Such a dynamic environment and intrinsic disordering of the B surface atoms present new opportunities to realize a diverse set of 2D boron structures. We validated the predicted surface disorder by characterizing exfoliated boron-terminated MgB2 nanosheets. We further discuss application-relevant implications, with a particular view towards understanding the impact of boron surface heterogeneity on hydrogen storage performance.

We report carbon impurity ion incident angles and deposition rates, along with silicon erosion rates, from measurements of micro-engineered trenches on a silicon surface exposed to L-mode deuterium plasmas at the DIII-D divertor. Post exposure ex-situ analysis determined elemental maps and concentrations, carbon deposition thicknesses, and erosion of silicon surfaces. Carbon deposition profiles on the trench floor showed carbon ion shadowing that was consistent with ERO calculations of average carbon ion angle distributions (IADs) for both polar and azimuthal angles. Measured silicon net erosion rates negatively correlated with the deposited carbon concentration at different locations. Differential erosion of surfaces on two different ion-downstream trench slope structures suggested that carbon deposition rate is affected by the carbon ion incident angle and significantly suppressed the surface erosion. The results suggest the C impurity ion incident angles, determined by the IADs and surface morphology, strongly affect erosion rates as well as the main ion (D, T, He) incident angles.

Silicon-based layered nanocomposites, comprised of covalent-metal interfaces, have demonstrated elevated resistance to radiation. The amorphization of the crystalline silicon sublayer during irradiation and/or heating can provide an additional mechanism for accommodating irradiation-induced defects. In this study, we investigated the mechanical strength of irradiated Si-based nanocomposites using atomistic modeling. We first examined dose effects on the defect evolution mechanisms near silicon-gold crystalline and amorphous interfaces. Our simulations reveal the growth of an emergent amorphous interfacial layer with increasing dose, a dominant factor mitigating radiation damage. We then examined the effect of radiation on the mechanical strength of silicon-gold multilayers by constructing yield surfaces. These results demonstrate a rapid onset strength loss with dose. Nearly identical behavior is observed in bulk gold, a phenomenon that can be rooted to the formation of radiation-induced stacking fault tetrahedra which dominate the dislocation emission mechanism during mechanical loading. Taken together, these results advance our understanding of the interaction between radiation-induced point defects and metal-covalent interfaces.

Abstract not provided.

Predictions of the spatial representation of instantaneous wall-bounded flows, via coherent structure-based models, are highly sensitive to the geometry of the representative structures employed by them. In this study, we propose a methodology to extract the three-dimensional (3-D) geometry of the statistically significant eddies from multi-point wall-turbulence datasets, for direct implementation into these models to improve their predictions. The methodology is employed here for reconstructing a 3-D statistical picture of the inertial wall coherent turbulence for all canonical wall-bounded flows, across a decade of friction Reynolds number (Reτ). These structures are responsible for the Reτ-dependence of the skin-friction drag and also facilitate the inner-outer interactions, making them key targets of structure-based models. The empirical analysis brings out the geometric self-similarity of the large-scale wall-coherent motions and also suggests the hairpin packet as the representative flow structure for all wall-bounded flows, thereby aligning with the framework on which the attached eddy model (AEM) is based. The same framework is extended here to also model the very-large-scaled motions, with a consideration of their differences in internal versus external flows. Implementation of the empirically-obtained geometric scalings for these large structures into the AEM is shown to enhance the instantaneous flow predictions for all three velocity components. Finally, an active flow control system driven by the same geometric scalings is conceptualized, towards favourably altering the influence of the wall coherent motions on the skin-friction drag.

Point defects in SiC are an attractive platform for quantum information and sensing applications because they provide relatively long spin coherence times, optical spin initialization, and spin-dependent fluorescence readout in a fabrication-friendly semiconductor. The ability to precisely place these defects at the optimal location in a host material with nano-scale accuracy is desirable for integration of these quantum systems with traditional electronic and photonic structures. Here, we demonstrate the precise spatial patterning of arrays of silicon vacancy (VSi) emitters in an epitaxial 4H-SiC (0001) layer through mask-less focused ion beam implantation of Li+. We characterize these arrays with high-resolution scanning confocal fluorescence microscopy on the Si-face, observing sharp emission lines primarily coming from the V1 ′ zero-phonon line (ZPL). The implantation dose is varied over 3 orders of magnitude, leading to VSi densities from a few per implantation spot to thousands per spot, with a linear dependence between ZPL emission and implantation dose. Optically-detected magnetic resonance (ODMR) is also performed, confirming the presence of V2 VSi. Our investigation reveals scalable and reproducible defect generation.

Defects and microstructural features spanning the atomic level to the microscale play deterministic roles in the expressed properties of materials. Yet studies of material evolution in response to environmental stimuli most often correlate resulting performance with one dominant microstructural feature only. Here, the dynamic evolution of swelling in a series of Ni-based concentrated solid solution alloys under high-temperature irradiation exposure is observed using continuous, in situ measurements of thermoelastic properties in bulk specimens. Unlike traditional evaluation techniques which account only for volumetric porosity identified using electron microscopy, direct property evaluation provides an integrated response across all defect length scales. In particular, the evolution in elastic properties during swelling is found to depend significantly on the entire size spectrum of defects, from the nano- to meso-scales, some of which are not resolvable in imaging. Observed changes in thermal transport properties depend sensitively on the partitioning of electronic and lattice thermal conductivity. This emerging class of in situ experiments, which directly measure integrated performance in relevant conditions, provides unique insight into material dynamics otherwise unavailable using traditional methods.

Big Data in the area of Remote Sensing has been growing rapidly. Remote sensors are used in surveillance, security, traffic, environmental monitoring, and autonomous sensing. Real-time detection of small moving targets using a remote sensor is an ongoing, challenging problem. Since the object is located far away from the sensor, the object often appears too small. The object’s signal-to-noise-ratio (SNR) is often very low. Occurrences such as camera motion, moving backgrounds (e.g., rustling leaves), low contrast and resolution of foreground objects makes it difficult to segment out the targeted moving objects of interest. Due to the limited appearance of the target, it is tough to obtain the target’s characteristics such as its shape and texture. Without these characteristics, filtering out false detections can be a difficult task. Detecting these targets, would often require the detector to operate under a low detection threshold. However, lowering the detection threshold could lead to an increase of false alarms. In this paper, the author will introduce a new method that improves the probability to detect low SNR objects, while decreasing the number of false alarms as compared to using the traditional baseline detection technique.

Background: Mitigation of climate change requires that new routes for the production of fuels and chemicals be as oil-independent as possible. The microbial conversion of lignocellulosic feedstocks into terpene-based biofuels and bioproducts represents one such route. This work builds upon previous demonstrations that the single-celled carotenogenic basidiomycete, Rhodosporidium toruloides, is a promising host for the production of terpenes from lignocellulosic hydrolysates. Results: This study focuses on the optimization of production of the monoterpene 1,8-cineole and the sesquiterpene α-bisabolene in R. toruloides. The α-bisabolene titer attained in R. toruloides was found to be proportional to the copy number of the bisabolene synthase (BIS) expression cassette, which in turn influenced the expression level of several native mevalonate pathway genes. The addition of more copies of BIS under a stronger promoter resulted in production of α-bisabolene at 2.2 g/L from lignocellulosic hydrolysate in a 2-L fermenter. Production of 1,8-cineole was found to be limited by availability of the precursor geranylgeranyl pyrophosphate (GPP) and expression of an appropriate GPP synthase increased the monoterpene titer fourfold to 143 mg/L at bench scale. Targeted mevalonate pathway metabolite analysis suggested that 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGR), mevalonate kinase (MK) and phosphomevalonate kinase (PMK) may be pathway bottlenecks are were therefore selected as targets for overexpression. Expression of HMGR, MK, and PMK orthologs and growth in an optimized lignocellulosic hydrolysate medium increased the 1,8-cineole titer an additional tenfold to 1.4 g/L. Expression of the same mevalonate pathway genes did not have as large an impact on α-bisabolene production, although the final titer was higher at 2.6 g/L. Furthermore, mevalonate pathway intermediates accumulated in the mevalonate-engineered strains, suggesting room for further improvement. Conclusions: This work brings R. toruloides closer to being able to make industrially relevant quantities of terpene from lignocellulosic biomass.

The unprecedented quality and sampling rate of seismometer and pressure sensors of the InSight Mars mission allow us to investigate infrasound through its pressure and ground deformation signals. This study focuses on compliance effects induced by acoustic waves propagating almost horizontally close to the surface. The compliance of acoustic waves is first estimated using the compliance estimates from pressure perturbations moving at wind speed. Then, a marker of compliance events is used to detect events of ground deformation induced by pressure variations, in three frequency bands from 0.4 to 3.2 Hz, from InSight sol 180 to 690. Additional selection criteria are imposed on the detected events to focus on acoustic waves and to remove various noise sources (e.g., wind effects or seismometer artifacts). After an automated selection, the visual inspection of the records allows us to validate two infrasound candidates that cannot be related to pressure perturbations moving at wind speed nor to known noise sources. For our highest quality infrasound candidate, the relation between this event and a convective vortex occuring 10 s later is tested. The azimuth of the vortex position at the time of infrasound detection is not consistent with the arrival azimuth of the suspected infrasound inferred from the polarization of seismometer records, thus the link between these two phenomena cannot be demonstrated. Further investigations would require a better understanding of wind-related noise impacting InSight sensors and of the effects of lateral variations of subsurface mechanical properties on the ground deformations induced by atmospheric pressure variations.