Publications

Nonlocal operators of fractional type are a popular modeling choice for applications that do not adhere to classical diffusive behavior; however, one major challenge in nonlocal simulations is the selection of model parameters. In this study we propose an optimization-based approach to parameter identification for fractional models with an optional truncation radius. We formulate the inference problem as an optimal control problem where the objective is to minimize the discrepancy between observed data and an approximate solution of the model, and the control variables are the fractional order and the truncation length. For the numerical solution of the minimization problem we propose a gradient-based approach, where we enhance the numerical performance by an approximation of the bilinear form of the state equation and its derivative with respect to the fractional order. Several numerical tests in one and two dimensions illustrate the theoretical results and show the robustness and applicability of our method.

On May 26, 2021, Sandia National Laboratories (SNL) convened a diverse group of experts spanning private industry, academia, the United States military and federal government, and the national laboratories, and hosted a series of panels to gain their insight on critical emergent research and capability development needs to support national cyber strategy objectives. Two panelists of experts presented their prepared remarks, followed by open discussion from over 250 audience members.

A demonstration of 2D hole gases in GeSn/Ge heterostructures with a mobility as high as 20 000 cm² V^–1 s^–1 is given. Both the Shubnikov–de Haas oscillations and integer quantum Hall effect are observed, indicating high sample quality. The Rashba spin-orbit coupling (SOC) is investigated via magneto-transport. Further, a transition from weak localization to weak anti-localization is observed, which shows the tunability of the SOC strength by gating. The magneto-transport data are fitted to the Hikami–Larkin–Nagaoka formula. The phase-coherence and spin-relaxation times, as well as spin-splitting energy and Rashba coefficient of the k-cubic term, are extracted. Furthermore, the analysis reveals that the effects of strain and confinement potential at a high fraction of Sn suppress the Rashba SOC caused by the GeSn/Ge heterostructures.

We report a detailed study of the tunnel barriers within a single-hole GaAs/AlGaAs double quantum dot device (DQD). For quantum information applications as well as fundamental studies, careful tuning and reliable measurements of the barriers are important requirements. In order to tune a DQD device adequately into the single-hole electric dipole spin resonance regime, one has to employ a variety of techniques to cover the extended range of tunnel couplings. In this work, we demonstrate four separate techniques, based upon charge sensing, quantum transport, time-resolved pulsing, and electron dipole spin resonance spectroscopy to determine the couplings as a function of relevant gate voltages and magnetic field. Measurements were performed under conditions of both symmetric and asymmetric tunnel couplings to the leads. Good agreement was observed between different techniques when measured under the same conditions. The results indicate that even in this relatively simple circuit, the requirement to tune multiple gates and the consequences of real potential profiles result in non-intuitive dependencies of the couplings as a function of the plunger gate voltage and the magnetic field.

With evolving landscape of DC power transmission and distribution, a reliable and fast protection against faults is critical, especially for medium- and high-voltage applications. Thus, solid-state circuit breakers (SSCB), consisting of cascaded silicon carbide (SiC) junction field-effect transistors (JFET), utilize the intrinsic normally-ON characteristic along with their low ON-resistance. This approach provides an efficient and robust protection solution from detrimental short-circuit events. However, for applications that require high-voltage blocking capability, a proper number of JFETs need be connected in series to achieve the desired blocking voltage rating. Ensuring equal voltage balancing across the JFETs during the switching transitions as well as the blocking stage is critical and hence, this paper presents a novel passive balancing network for series connected JFETs for DC SSCB applications. The dynamic voltage balancing network to synchronize both the turn ON and OFF intervals is described analytically. Moreover, the static voltage balancing network is implemented to establish equal sharing of the total blocking voltage across the series connection of JFETs. The proposed dynamic and steady-state balancing networks are validated by SPICE simulation and experimental results.

Ultrafast all-optical switching using Mie resonant metasurfaces requires both on-demand tunability of the wavefront of the light and ultrafast time response. However, devising a switching mechanism that has a high contrast between its "on"and "off"states without compromising speed is challenging. Here, we report the design of a tunable Mie resonant metasurface that achieves this behavior. Our approach utilizes a diffractive array of semiconductor resonators that support both dipolar and quadrupolar Mie resonances. By balancing the strengths of the dipole and quadrupole resonances, we can suppress radiation into the first diffraction order, thus creating a clearly delineated "off"-state at the operating wavelength. Then, we use optical injection of free- carriers to spectrally shift the multipoles and rebalance the multipole strengths, thereby enabling radiation into the diffraction order - all on an ultrafast timescale. We demonstrate ultrafast off-to-on switching with Ion/Ioff ≈ 5 modulation of the diffracted intensity and ultrafast on-to-off switching with Ion/Ioff ≈ 9 modulation. Both switches exhibit a fast τtr ≈ 2.7 ps relaxation time at 215 μJ cm-2 pump fluence. Further, we show that for higher fluences, the temporal response of the metasurface is governed by thermo-optic effects. This combination of multipole engineering with lattice diffraction opens design pathways for tunable metasurface-based integrated devices.

We adapt the robust phase estimation algorithm to the evaluation of energy differences between two eigenstates using a quantum computer. This approach does not require controlled unitaries between auxiliary and system registers or even a single auxiliary qubit. As a proof of concept, we calculate the energies of the ground state and low-lying electronic excitations of a hydrogen molecule in a minimal basis on a cloud quantum computer. The denominative robustness of our approach is then quantified in terms of a high tolerance to coherent errors in the state preparation and measurement. Conceptually, we note that all quantum phase estimation algorithms ultimately evaluate eigenvalue differences.

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both systematic and statistical errors below ≈0.05 eV. The ground state exhibits a fixed-node error of only 1.3(2)% of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. We have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments.

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

The first-principles computation of the surfaces of metals is typically accomplished through slab calculations of finite thickness. The extraction of a convergent surface formation energy from slab calculations is dependent upon defining an appropriate bulk reference energy. I describe a method for an independently computed, slab-consistent bulk reference that leads to convergent surface formation energies from slab calculations that also provides realistic uncertainties for the magnitude of unavoidable nonlinear divergence in the surface formation energy with slab thickness. The accuracy is demonstrated on relaxed, unreconstructed low-index aluminum surfaces with slabs with up to 35 layers.

The stability of low-index platinum surfaces and their electronic properties is investigated with density functional theory, toward the goal of understanding the surface structure and electron emission, and identifying precursors to electrical breakdown, on nonideal platinum surfaces. Propensity for electron emission can be related to a local work function, which, in turn, is intimately dependent on the local surface structure. The (1×N) missing row reconstruction of the Pt(110) surface is systematically examined. The (1×3) missing row reconstruction is found to be the lowest in energy, with the (1×2) and (1×4) slightly less stable. In the limit of large (1×N) with wider (111) nanoterraces, the energy accurately approaches the asymptotic limit of the infinite Pt(111) surface. This suggests a local energetic stability of narrow (111) nanoterraces on free Pt surfaces that could be a common structural feature in the complex surface morphologies, leading to work functions consistent with those on thermally grown Pt substrates.

Lawrence Pratt’s career following completion of his Ph.D. at the University of Illinois Urbana Champaign has taken him from Harvard University, to the University of California, Berkeley, and Los Alamos National Laboratory. Most recently, he joined the faculty of the Department of Chemical and Biomolecular Engineering at Tulane University in 2008. Over his career Lawrence has been a leader in theoretical physical chemistry, making influential contributions to a number of areas including the theory of the hydrophobic effect, the development of transition path sampling, contributions to orbital free density functional theory, and the theory of liquids and solutions.

High-energy-density physics is the field of physics concerned with studying matter at extremely high temperatures and densities. Such conditions produce highly nonlinear plasmas, in which several phenomena that can normally be treated independently of one another become strongly coupled. The study of these plasmas is important for our understanding of astrophysics, nuclear fusion and fundamental physics—however, the nonlinearities and strong couplings present in these extreme physical systems makes them very difficult to understand theoretically or to optimize experimentally. Here we argue that machine learning models and data-driven methods are in the process of reshaping our exploration of these extreme systems that have hitherto proved far too nonlinear for human researchers. From a fundamental perspective, our understanding can be improved by the way in which machine learning models can rapidly discover complex interactions in large datasets. From a practical point of view, the newest generation of extreme physics facilities can perform experiments multiple times a second (as opposed to approximately daily), thus moving away from human-based control towards automatic control based on real-time interpretation of diagnostic data and updates of the physics model. To make the most of these emerging opportunities, we suggest proposals for the community in terms of research design, training, best practice and support for synthetic diagnostics and data analysis.

Multivariate time series are used in many science and engineering domains, including health-care, astronomy, and high-performance computing. A recent trend is to use machine learning (ML) to process this complex data and these ML-based frameworks are starting to play a critical role for a variety of applications. However, barriers such as user distrust or difficulty of debugging need to be overcome to enable widespread adoption of such frameworks in production systems. To address this challenge, we propose a novel explainability technique, CoMTE, that provides counterfactual explanations for supervised machine learning frameworks on multivariate time series data. Using various machine learning frameworks and data sets, we compare CoMTE with several state-of-the-art explainability methods and show that we outperform existing methods in comprehensibility and robustness. We also show how CoMTE can be used to debug machine learning frameworks and gain a better understanding of the underlying multivariate time series data.

High-performance computing (HPC) researchers have long envisioned scenarios where application workflows could be improved through the use of programmable processing elements embedded in the network fabric. Recently, vendors have introduced programmable Smart Network Interface Cards (SmartNICs) that enable computations to be offloaded to the edge of the network. There is great interest in both the HPC and high-performance data analytics (HPDA) communities in understanding the roles these devices may play in the data paths of upcoming systems. This paper focuses on characterizing both the networking and computing aspects of NVIDIA’s new BlueField-2 SmartNIC when used in a 100Gb/s Ethernet environment. For the networking evaluation we conducted multiple transfer experiments between processors located at the host, the SmartNIC, and a remote host. These tests illuminate how much effort is required to saturate the network and help estimate the processing headroom available on the SmartNIC during transfers. For the computing evaluation we used the stress-ng benchmark to compare the BlueField-2 to other servers and place realistic bounds on the types of offload operations that are appropriate for the hardware. Our findings from this work indicate that while the BlueField-2 provides a flexible means of processing data at the network’s edge, great care must be taken to not overwhelm the hardware. While the host can easily saturate the network link, the SmartNIC’s embedded processors may not have enough computing resources to sustain more than half the expected bandwidth when using kernel-space packet processing. From a computational perspective, encryption operations, memory operations under contention, and on-card IPC operations on the SmartNIC perform significantly better than the general-purpose servers used for comparisons in our experiments. Therefore, applications that mainly focus on these operations may be good candidates for offloading to the SmartNIC.

We analyze the dislocation content of grain boundary (GB) phase junctions, i.e., line defects separating two different GB phases coexisting on the same GB plane. While regular GB disconnections have been characterized for a variety of interfaces, GB phase junctions formed by GBs with different structures and different numbers of excess atoms have not been previously studied. We apply a general Burgers circuit analysis to calculate the Burgers vectors b of junctions in two ς5 Cu boundaries previously simulated with molecular dynamics. The Burgers vectors of these junctions cannot be described by the displacement shift complete (DSC) lattice alone. We show that, in general, the normal component of b is not equal to the difference in the GB excess volumes but contains another contribution from the numbers of GB atoms per unit area ΔN∗ required to transform one GB phase into another. In the boundaries studied, the latter component dominates and even changes the sign of b. We derive expressions for the normal and tangential components of b in terms of the DSC lattice vectors and the non-DSC part due to ΔN∗ and additional GB excess properties, including excess volume and shears. These expressions provide a connection between GB phase transformations driven by the GB free energy difference and the motion of GB junctions under applied normal and shear stresses. The proposed analysis quantifies b and therefore makes it possible to calculate the elastic part of the energy of these defects, evaluate their contribution to the nucleation barrier during GB phase transformations, and treat elastic interactions with other defects.

Plasmas formed in microscale gaps at DC and plasmas formed at radiofrequency (RF) both deviate in behavior compared to the classical Paschen curve, requiring lower voltage to achieve breakdown due to unique processes and dynamics, such as field emission and controlled rates of electron/ion interactions. Both regimes have been investigated independently, using high precision electrode positioning systems for microscale gaps or large, bulky emitters for RF. However, no comprehensive study of the synergistic phenomenon between the two exists. The behavior in such a combined system has the potential to reach sub-10 V breakdown, which combined with the unique electrical properties of microscale plasmas could enable a new class of RF switches, limiters and tuners.This work describes the design and fabrication of novel on-wafer microplasma devices with gaps as small as 100 nm to be operated at GHz frequencies. We used a dual-sacrificial layer process to create devices with microplasma gaps integrated into RF compatible 50 Ω coplanar waveguide transmission lines, which will allow this coupled behaviour to be studied for the first time. These devices are modelled using conventional RF simulations as well as the Sandia code, EMPIRE, which is capable of modelling the breakdown and formation of plasma in microscale gaps driven by high frequencies. Synchronous evaluation of the modelled electrical and breakdown behaviour is used to define device structures, predict behaviour and corroborate results. We further report preliminary independent testing of the microscale gap and RF behaviour. DC testing shows modified-Paschen curve behaviour for plasma gaps at and below four microns, demonstrating decreased breakdown voltage with reduced gap size. Additionally, preliminary S-parameter measurements of as-prepared and connectorized devices have elucidated RF device behaviour. Together, these results provide baseline data that enables future experiments as well as discussion of projected performance and applications for these unique devices.

Studies of size effects on thermal conductivity typically necessitate the fabrication of a comprehensive film thickness series. In this Letter, we demonstrate how material fabricated in a wedged geometry can enable similar, yet higher-throughput measurements to accelerate experimental analysis. Frequency domain thermoreflectance (FDTR) is used to simultaneously determine the thermal conductivity and thickness of a wedged silicon film for thicknesses between 100 nm and 17 μm by considering these features as fitting parameters in a thermal model. FDTR-deduced thicknesses are compared to values obtained from cross-sectional scanning electron microscopy, and corresponding thermal conductivity measurements are compared against several thickness-dependent analytical models based upon solutions to the Boltzmann transport equation. Our results demonstrate how the insight gained from a series of thin films can be obtained via fabrication of a single sample.

Using the local moment counter charge (LMCC) method to accurately represent the asymptotic electrostatic boundary conditions within density functional theory supercell calculations, we present a comprehensive analysis of the atomic structure and energy levels of point defects in cubic silicon carbide (3C-SiC). Finding that the classical long-range dielectric screening outside the supercell induced by a charged defect is a significant contributor to the total energy. we describe and validate a modified Jost screening model to evaluate this polarization energy. This leads to bulk-converged defect levels in finite size supercells. With the LMCC boundary conditions and a standard Perdew-Burke-Ernzerhof (PBE) exchange correlation functional, the computed defect level spectrum exhibits no band gap problem: the range of defect levels spans ∼2.4eV, an effective defect band gap that agrees with the experimental band gap. Comparing with previous literature, our LMCC-PBE defect results are in consistent agreement with the hybrid-exchange functional results of Oda et al. [J. Chem. Phys. 139, 124707 (2013)JCPSA60021-960610.1063/1.4821937] rather than their PBE results. The difference with their PBE results is attributed to their use of a conventional jellium approximation rather than the more rigorous LMCC approach for handling charged supercell boundary conditions. The difference between standard dft and hybrid functional results for defect levels lies not in a band gap problem but rather in solving a boundary condition problem. The LMCC-PBE entirely mitigates the effect of the band gap problem on defect levels. The more computationally economical PBE enables a systematic exploration of 3C-SiC defects, where, most notably, we find that the silicon vacancy undergoes Jahn-Teller-induced distortions from the previously assumed Td symmetry, and that the divacancy, like the silicon vacancy, exhibits a site-shift bistability in p-type conditions.

Inspired by the parallelism and efficiency of the brain, several candidates for artificial synapse devices have been developed for neuromorphic computing, yet a nonlinear and asymmetric synaptic response curve precludes their use for backpropagation, the foundation of modern supervised learning. Spintronic devices - which benefit from high endurance, low power consumption, low latency, and CMOS compatibility - are a promising technology for memory, and domain-wall magnetic tunnel junction (DW-MTJ) devices have been shown to implement synaptic functions such as long-term potentiation and spike-timing dependent plasticity. In this work, we propose a notched DW-MTJ synapse as a candidate for supervised learning. Using micromagnetic simulations at room temperature, we show that notched synapses ensure the non-volatility of the synaptic weight and allow for highly linear, symmetric, and reproducible weight updates using either spin transfer torque (STT) or spin-orbit torque (SOT) mechanisms of DW propagation. We use lookup tables constructed from micromagnetics simulations to model the training of neural networks built with DW-MTJ synapses on both the MNIST and Fashion-MNIST image classification tasks. Accounting for thermal noise and realistic process variations, the DW-MTJ devices achieve classification accuracy close to ideal floating-point updates using both STT and SOT devices at room temperature and at 400 K. Our work establishes the basis for a magnetic artificial synapse that can eventually lead to hardware neural networks with fully spintronic matrix operations implementing machine learning.

The rapidly increasing use of electronics in high-radiation environments and the continued evolution in transistor architectures and materials demand improved methods to characterize the potential damaging effects of radiation on device performance. Here, electron-beam-induced current is used to map hot-carrier transport in model metal-oxide semiconductor field-effect transistors irradiated with a 300 KeV focused He+ beam as a localized line spanning across the gate and bulk Si. By correlating the damage to the electronic properties and combining these results with simulations, the contribution of spatially localized radiation damage on the device characteristics is obtained. This identified damage, caused by the He+ beam, is attributed to localized interfacial Pb centers and delocalized positive fixed-charges, as surmised from simulations. Comprehension of the long-term interaction and mobility of radiation-induced damage are key for future design of rad-hard devices.

This project is focused on developing advanced combustion strategies for mixing-controlled compression ignition (MCCI, i.e., diesel-cycle) engines that are synergistic with renewable and/or unconventional fuels in a manner that enhances domestic energy security, economic competitiveness, and environmental quality. During this reporting period, the two focus areas were ducted fuel injection (DFI) and surrogate diesel fuels.

Gallium is incorporated into the strain-balanced In(Ga)As/InAsSb superlattice system to achieve the same mid-wave infrared cutoff tunability as conventional Ga-free InAs/InAsSb type-II superlattices, but with an additional degree of design freedom to enable optimization of absorption and transport properties. Time-resolved photoluminescence measurements of InGaAs/InAsSb superlattice characterization- and doped device structures are reported from 77 to 300 K and compared to InAs/InAsSb. The low-injection photoluminescence decay yields the minority carrier lifetime, which is analyzed with a recombination rate model, enabling the determination of the temperature-dependent Shockley-Read-Hall, radiative, and Auger recombination lifetimes and extraction of defect energy levels and capture cross section defect concentration products. The Shockley-Read-Hall-limited lifetime of undoped InGaAs/InAsSb is marginally reduced from 2.3 to 1.4 μs due to the inclusion of Ga; however, given that Ga improves the vertical hole mobility by a factor of >10×, a diffusion-limited InGaAs/InAsSb superlattice nBn could expect a lower bound of 2.5× improvement in diffusion length with significant impact on photodetector quantum efficiency and radiation hardness. At temperatures below 120 K, the doped device structures are Shockley-Read-Hall limited at 0.5 μs, which shows promise for detector applications.

Pulsed laser irradiation is used to investigate the local initiation of rapid, self-propagating formation reactions in Al/Pt multilayers. The single pulse direct laser ignition of these 1.6 μm thick freestanding foils was characterized over 10 decades of pulse duration (10 ms to 150 fs). Finite element, reactive heat transport modeling of the near-threshold conditions has identified three distinct ignition pathways. For milli- to microsecond pulses, ignition occurs following sufficient absorption of laser energy to enable diffusion of Al and Pt between layers such that the heat released from the corresponding exothermic reaction overcomes conductive losses outside the laser-irradiated zone. When pulse duration is decreased into the nanosecond regime, heat is concentrated near the surface such that the Al locally melts, and a portion of the top-most bilayers react initially. The favorable kinetics and additional heat enable ignition. Further reducing pulse duration to hundreds of femtoseconds leads to a third ignition pathway. While much of the energy from these pulses is lost to ablation, the remaining heat beneath the crater can be sufficiently concentrated to drive a transverse self-propagating reaction, wherein the heat released from mixing at each interface occurs under kinetic conditions capable of igniting the subsequent layer.

Multisensor networks deployed at nuclear facilities can be leveraged to collect data used as inputs to machine learning models predicting nuclear safeguard relevant information. This work demonstrates an application of this idea by regressing nuclear reactor power levels, a key indicator for nuclear safeguard verification, at the McClellan Nuclear Research Center using data collected by five Merlyn multisensor platforms with LASSO and LSTM models. This work also demonstrates the use of Leave One Node Out to measure the importance of each multisensor for this regression problem providing insight into model explainability and allowing inferential hypotheses about the nuclear facility to be made. This work can be used as a starting point for future development of methods for regression on reactor power levels at nuclear facilities using multisensor network data.

Following the 1986 disastrous accident at Chernobyl, gaining independence following the fall of the former Soviet Union, two revolutions, and on-going Russian intervention, Ukrainea is seeking to expand its use of commercial nuclear energy and to further reduce its dependence on Russia for energy. Nuclear energy in Ukraine has made a significant contribution toward achieving sustainable development and social goals. Ukraine is very vulnerable to the effects of climate change and nuclear is a cornerstone of the country’s effort to combat the effects. With broad support from the U.S. government, Ukraine is cooperating with the U.S. on fresh nuclear fuel supply, civilian nuclear security, outage optimization, spent fuel storage, and future nuclear technology. And is looking to tap U.S. expertise on plant management, and electric grid operations and expansion. NEIc hosted this virtual event to provide an update on nuclear energy projects and plans in Ukraine. Featured speakers included Yaroslav Demchenkov, Deputy Minister of Energy, Taras Kachka, Deputy Minister – Trade Representative, Petro Kotin, CEO, Energoatom,d and Ann K. Ganzer, Senior Bureau Official, Bureau of International Security and Nonproliferation, U.S. Dept. of State. Maria Korsnick, NEI CEO, moderated the discussion. About forty individuals attended representing the Ukraine and U.S. governments, and U.S. suppliers.

This paper describes the verification and validation (V&V) framework developed for the stochastic Particle-in-Cell, Direct Simulation Monte Carlo code Aleph. An ideal framework for V&V from the viewpoint of the authors is described where a physics problem is defined, and relevant physics models and parameters to the defined problem are assessed and captured in a Phenomena Identification and Ranking Table (PIRT). Furthermore, numerous V&V examples guided by the PIRT for a simple gas discharge are shown to demonstrate the V&V process applied to a real-world simulation tool with the overall goal to demonstrably increase the confidence in the results for the simulation tool and its predictive capability. Although many examples are provided here to demonstrate elements of the framework, the primary goal of this work is to introduce this framework and not to provide a fully complete implementation, which would be a much longer document. Comparisons and contrasts are made to more usual approaches to V&V, and techniques new to the low-temperature plasma community are introduced. Specific challenges relating to the sufficiency of available data (e.g., cross sections), the limits of ad hoc validation approaches, the additional difficulty of utilizing a stochastic simulation tool, and the extreme cost of formal validation are discussed.

Reactive gas formation in pores of metal–organic frameworks (MOFs) is a known mechanism of framework destruction; understanding those mechanisms for future durability design is key to next generation adsorbents. Herein, an extensive set of ab initio molecular dynamics (AIMD) simulations are used for the first time to predict competitive adsorption of mixed acid gases (NO2 and H2O) and the in-pore reaction mechanisms for a series of rare earth (RE)-DOBDC MOFs. Spontaneous formation of nitrous acid (HONO) is identified as a result of deprotonation of the MOF organic linker, DOBDC. The unique DOBDC coordination to the metal clusters allows for proton transfer from the linker to the NO2 without the presence of H2O and may be a factor in DOBDC MOF durability. This is a previously unreported mechanisms of HONO formation in MOFs. With the presented methodology, prediction of future gas interactions in new nanoporous materials can be achieved.

Plastics, with their ubiquitous presence in our daily lives and environment, pose an uncomfortable conundrum. Producers and consumers are aware of the value of these organic ingredients in material flow, yet their persistence and disruption to the ecological milieu desperately stipulate a shift in the status quo. Biodegradable plastics-as the name suggests-has its appeal in ensuring the safe return of carbon to ecosystems by complete assimilation of the degraded product as a food source for soil or aquatic microorganisms. However, despite more than a decade of commercial presence, these plastics are still far from replacing the demand for fossil-fuel-based commodity plastics. We discuss this apparent disconnect herein through a material value chain perspective. We review the current state of commercial biodegradable plastics and contrast it against the desired state of the zero-waste-focused circular economy. To close the gap, we suggest critical research needs concerning the structure and properties of biodegradable plastics, testing standards, application development, and waste management. The ultimate success in displacing conventional plastics with biodegradable alternatives will be predicated on collaboration between all stakeholders across the product value chain.

This paper serves as the Interface Control Document (ICD) for the Seascape automated test harness developed at Sandia National Laboratories. The primary purposes of the Seascape system are: (1) provide a place for accruing large, curated, labeled data sets useful for developing and evaluating detection and classification algorithms (including, but not limited to, supervised machine learning applications) (2) provide an automated structure for specifying, running and generating reports on algorithm performance. Seascape uses GitLab, Nexus, Solr, and Banana, open source codes, together with code written in the Python language, to automatically provision and configure computational nodes, queue up jobs to accomplish algorithms test runs against the stored data sets, gather the results and generate reports which are then stored in the Nexus artifact server.

The use of evidence theory and associated cumulative plausibility functions (CPFs), cumulative belief functions (CBFs), cumulative distribution functions (CDFs), complementary cumulative plausibility functions (CCPFs), complementary cumulative belief functions (CCBFs), and complementary cumulative distribution functions (CCDFs) in the analysis of loss of assured safety (LOAS) for weak link (WL)/strong link (SL) systems is introduced and illustrated. Article content includes cumulative and complementary cumulative belief, plausibility, and probability for (i) time at which LOAS occurs for a one WL/two SL system, (ii) time at which a two-link system fails, (iii) temperature at which a two-link system fails, and (iv) temperature at which LOAS occurs for a one WL/two SL system. The presented results can be generalized to systems with more than one WL and two SLs.

The use of evidence theory and associated cumulative plausibility functions (CPFs), cumulative belief functions (CBFs), cumulative distribution functions (CDFs), complementary cumulative plausibility functions (CCPFs), complementary cumulative belief functions (CCBFs), and complementary cumulative distribution functions (CCDFs) in the analysis of time and temperature margins associated with loss of assured safety (LOAS) for one weak link (WL)/two strong link (SL) systems is illustrated. Article content includes cumulative and complementary cumulative belief, plausibility, and probability for (i) SL/WL failure time margins defined by (time at which SL failure potentially causes LOAS) – (time at which WL failure potentially prevents LOAS), (ii) SL/WL failure temperature margins defined by (the temperature at which SL failure potentially causes LOAS) – (the temperature at which WL failure potentially prevents LOAS), and (iii) SL/SL failure temperature margins defined by (the temperature at which SL failure potentially causes LOAS) – (the temperature of SL whose failure potentially causes LOAS at the time at which WL failure potentially prevents LOAS).

Liquid organic hydrogen carriers such as alcohols and polyols are a high-capacity means of transporting and reversibly storing hydrogen that demands effective catalysts to drive the (de)hydrogenation reactions under mild conditions. We employed a combined theory/experiment approach to develop MOF-74 catalysts for alcohol dehydrogenation and examine the performance of the open metal sites (OMS), which have properties analogous to the active sites in high-performance single-site catalysts and homogeneous catalysts. Methanol dehydrogenation was used as a model reaction system for assessing the performance of five monometallic M-MOF-74 variants (M = Co, Cu, Mg, Mn, Ni). Co-MOF-74 and Ni-MOF-74 give the highest H2 productivity. However, Ni-MOF-74 is unstable under reaction conditions and forms metallic nickel particles. To improve catalyst activity and stability, bimetallic (NixMg1-x)-MOF-74 catalysts were developed that stabilize the Ni OMS and promote the dehydrogenation reaction. An optimal composition exists at (Ni0.32Mg0.68)-MOF-74 that gives the greatest H2 productivity, up to 203 mL gcat-1 min-1 at 300 °C, and maintains 100% selectivity to CO and H2 between 225-275 °C. The optimized catalyst is also active for the dehydrogenation of other alcohols. DFT calculations reveal that synergistic interactions between the open metal site and the organic linker lead to lower reaction barriers in the MOF catalysts compared to the open metal site alone. This work expands the suite of hydrogen-related reactions catalyzed by MOF-74 which includes recent work on hydroformulation and our earlier reports of aryl-ether hydrogenolysis. Moreover, it highlights the use of bimetallic frameworks as an effective strategy for stabilizing a high density of catalytically active open metal sites. This journal is

This interim report is an update of ongoing experimental and modeling work on bentonite material described in Jové Colón et al. (2019, 2020) from past international collaboration activities. As noted in Jové Colón et al. (2020), work on international repository science activities such as FEBEX-DP and DECOVALEX19 is either no longer continuing by the international partners. Nevertheless, research activities on the collected sample materials and field data are still ongoing. Descriptions of these underground research laboratory (URL) R&D activities are described elsewhere (Birkholzer et al. 2019; Jové Colón et al. 2020) but will be explained here when needed. The current reports recent reactive-transport modeling on the leaching of sedimentary rock.

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫ C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. Here, we find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

Channelrhodopsins (ChR) are light-sensitive cation channels used in optogenetics, a technique that applies light to control cells (e.g., neurons) that have been modified genetically to express those channels. Although mutations are known to affect pore kinetics, little is known about how mutations induce changes at the molecular scale. To address this issue, we first measured channel opening and closing rates of a ChR chimera (C1C2) and selected variants (N297D, N297V, and V125L). Then, we used atomistic simulations to correlate those rates with changes in pore structure, hydration, and chemical interactions among key gating residues of C1C2 in both closed and open states. Overall, the experimental results show that C1C2 and its mutants do not behave like ChR2 or its analogous variants, except V125L, making C1C2 a unique channel. Our atomistic simulations confirmed that opening of the channel and initial hydration of the gating regions between helices I, II, III, and VII of the channel occurs with 1) the presence of 13-cis retinal; 2) deprotonation of a glutamic acid gating residue, E129; and 3) subsequent weakening of the central gate hydrogen bond between the same glutamic acid E129 and asparagine N297 in the central region of the pore. Also, an aspartate (D292) is the unambiguous primary proton acceptor for the retinal Schiff base in the hydrated channel.

Alternative candidate precursors to [Hf(BH4)4] for low-temperature chemical vapor deposition of hafnium diboride (HfB2) films were identified using density functional theory simulations of molecules with the composition [Hf(BH4)2L2], where L = -OH, -OMe, -O-t-Bu, -NH2, -N═C═O, -N(Me)2, and -N(CH2)5NH2 (1-piperidin-2-amine referred to as Pip2A). Disassociation energies (ED), potential energy surface (PES) scans, ionization potentials, and electron affinities were all calculated to identify the strength of the Hf-L bond and the potential reactivity of the candidate precursor. Ultimately, the low ED (2.07 eV) of the BH4 ligand removal from the Hf atom in [Hf(BH4)4] was partially attributed to an intermediate state where [Hf(BH4)3(H)] and BH3 is formed. Of the candidate precursors investigated, three exhibited a similar mechanism, but only -Pip2A had a PES scan that indicated binding competitive with [Hf(BH4)4], making it a viable candidate for further study.

Two-dimensional (2D) metal-boride-derived nanostructures have been a focus of intense research for the past decade, with an emphasis on new synthetic approaches, as well as on the exploration of possible applications in next-generation advanced materials and devices. Their unusual mechanical, electronic, optical, and chemical properties, arising from low dimensionality, present a new paradigm to the science of metal borides that has traditionally focused on their bulk properties. This Perspective discusses the current state of research on metal-boride-derived 2D nanostructures, highlights challenges that must be overcome, and identifies future opportunities to fully utilize their potential.

In this work, experimental results from a study on the evolution of gas jets ejected through the orifices of a pre-chamber in a heavy-duty optical engine are presented. The work examines conditions without fuel inside the main-chamber, which helps to describe the dynamics of the ejected gas jets without the interference of subsequent combustion in the main-chamber. Experimental diagnostics consist of high-speed visible intensified imaging and low-speed infrared imaging. Additionally a one-dimensional gas jet model is used to characterize the spatial distribution of the ejected flow, including parameters such as tip penetration, which are then validated based on experimental results. Different stages in the ejection of pre-chamber jets are identified, with chemical activity restricted to a maximum distance of 5 to 10 orifice diameters downstream of the orifice as indicated by the recorded visible radiation. Sensitivity of cycle-to-cycle variations in pre-chamber jet development to the air-to-fuel ratio in the pre-chamber observed in the experiments is in most part attributed to the variations in the timing of combustion initiation in the pre-chamber. The influence of the ejection flow on the penetration of the gas jet on a cycle-to-cycle basis is presented using the one-dimensional model. The one-dimensional model also indicates that the local flow exhibits highest sensitivity to operating conditions during the start of ejection until the timing when maximum flow is attained. Differences that exist during the decreasing mass-flow ejection time-period tend to smear out in part due to the transient slowdown of the ejection process.

A high-speed temperature diagnostic based on spontaneous Raman scattering (SRS) was demonstrated using a pulse-burst laser. The technique was first benchmarked in near-adiabatic H2-air flames at a data-acquisition rate of 5 kHz using an integrated pulse energy of 1.0 J per realization. Both the measurement precision and accuracy in the flame were within 3% of adiabatic predictions. This technique was then evaluated in a challenging free-piston shock tube environment operated at a shock Mach number of 3.5. SRS thermometry resolved the temperature in post-incident and post-reflected shock flows at a repetition rate of 3 kHz and clearly showed cooling associated with driver expansion waves. Collectively, this Letter represents a major advancement for SRS in impulsive facilities, which had previously been limited to steady state regions or single-shot acquisition.

While using a prosthesis, transtibial amputees can experience pain and discomfort brought on by large pressure gradients at the interface between the residual limb and the prosthetic socket. Current prosthetic interface solutions attempt to alleviate these pressure gradients using soft homogenous liners to reduce and distribute pressures. This research investigates an additively manufactured metamaterial inlay with a tailored mechanical response to reduce peak pressure gradients around the limb. The inlay uses a hyperelastic behaving metamaterial (US10244818) comprised of triangular pattern unit cells, 3D printed with walls of various thicknesses controlled by draft angles. The hyperelastic material properties are modeled using a Yeoh third-order model. The third-order coefficients can be adjusted and optimized, which corresponds to a change in the unit cell wall thickness to create an inlay that can meet the unique offloading needs of an amputee. Finite element analysis simulations evaluated the pressure gradient reduction from (1) a standard homogenous silicone liner, (2) a prosthetist's inlay prescription that utilizes three variations of the metamaterial, and (3) a metamaterial solution with optimized Yeoh third-order coefficients. Compared to a traditional homogenous silicone liner for two unique limb loading scenarios, the prosthetist prescribed inlay and the optimized material inlay can achieve equal or greater pressure gradient reduction capabilities. These preliminary results show the potential feasibility of implementing this metamaterial as a method of personalized medicine for transtibial amputees by creating a customizable interface solution to meet the unique performance needs of an individual patient.

The potential for control design to dramatically improve the economic viability of wave energy has generated a great deal of interest and excitement. However, for a number of reasons, the promised benefits from better control designs have yet to be widely realized by wave energy devices and wave energy remains a relatively nascent technology. This brief paper summarizes a simple, yet powerful approach to wave energy dynamics modeling, and subsequent control design based on impedance matching. Our approach leverages the same concepts that are exploited by a simple FM radio to achieve a feedback controller for wave energy devices that approaches optimal power absorption. If fully utilized, this approach can deliver immediate and consequential reductions to the cost of wave energy. Additionally, this approach provides the necessary framework for control co-design of a wave energy converter, in which an understanding of the control logic allows for synchronous design of the device control system and hardware.

Ferritic/martensitic steels, such as HT-9, are known for their complex microstructural features and mechanical properties. In this paper, in-situ micro-tensile tests and traditional fractography methods were utilized to study the fracture behavior of proton-irradiated HT-9 steels. First, to evaluate the viability of micro-tensile tests for nuclear material qualification process, meso‑tensile tests on as-received HT-9 steels were performed. Fracture mechanisms of unirradiated HT-9 steels at both length scales were compared and underlying mechanisms discussed. The direct comparison of micro- and meso‑scale data shows a distinctive size effect demonstrated by the increase in yield stress (YS). Upon completion of initial assessment, specimens were irradiated with 4 MeV+ protons to three fluences, all of which were lower than 0.01 displacements per atom (dpa). As expected, the YS increases with irradiation. However, at 7 × 10−3 dpa, the reversal of the trend was observed, and the YS exhibited sharp decline. We demonstrate that at lower length scales, grain structure has a more profound impact on the mechanical properties of irradiated materials, which provides information needed to fill in the gap in current understanding of the HT-9 fracture at different length scales.

This work seeks to advance the state of the art in HPC I/O performance analysis and interpretation. In particular, we demonstrate effective techniques to: (1) model output performance in the presence of I/O interference from production loads; (2) build features from write patterns and key parameters of the system architecture and configurations; (3) employ suitable machine learning algorithms to improve model accuracy. We train models with five popular regression algorithms and conduct experiments on two distinct production HPC platforms. We find that the lasso and random forest models predict output performance with high accuracy on both of the target systems. We also explore use of the models to guide adaptation in I/O middleware systems, and show potential for improvements of at least 15% from model-guided adaptation on 70% of samples, and improvements up to 10 × on some samples for both of the target systems.

We evaluate the sensitivity of neuromorphic inference accelerators based on silicon-oxide-nitride-oxide-silicon (SONOS) charge trap memory arrays to total ionizing dose (TID) effects. Data retention statistics were collected for 16 Mbit of 40-nm SONOS digital memory exposed to ionizing radiation from a Co-60 source, showing good retention of the bits up to the maximum dose of 500 krad(Si). Using this data, we formulate a rate-equation-based model for the TID response of trapped charge carriers in the ONO stack and predict the effect of TID on intermediate device states between 'program' and 'erase.' This model is then used to simulate arrays of low-power, analog SONOS devices that store 8-bit neural network weights and support in situ matrix-vector multiplication. We evaluate the accuracy of the irradiated SONOS-based inference accelerator on two image recognition tasks - CIFAR-10 and the challenging ImageNet data set - using state-of-the-art convolutional neural networks, such as ResNet-50. We find that across the data sets and neural networks evaluated, the accelerator tolerates a maximum TID between 10 and 100 krad(Si), with deeper networks being more susceptible to accuracy losses due to TID.

We evaluate the resilience of CoFeB/MgO/CoFeB magnetic tunnel junctions (MTJs) with perpendicular magnetic anisotropy (PMA) to displacement damage induced by heavy-ion irradiation. MTJs were exposed to 3-MeV Ta2+ ions at different levels of ion beam fluence spanning five orders of magnitude. The devices remained insensitive to beam fluences up to $10^{11}$ ions/cm2, beyond which a gradual degradation in the device magnetoresistance, coercive magnetic field, and spin-transfer-torque (STT) switching voltage were observed, ending with a complete loss of magnetoresistance at very high levels of displacement damage (>0.035 displacements per atom). The loss of magnetoresistance is attributed to structural damage at the MgO interfaces, which allows electrons to scatter among the propagating modes within the tunnel barrier and reduces the net spin polarization. Ion-induced damage to the interface also reduces the PMA. This study clarifies the displacement damage thresholds that lead to significant irreversible changes in the characteristics of STT magnetic random access memory (STT-MRAM) and elucidates the physical mechanisms underlying the deterioration in device properties.

The Zenifim Formation is being considered as a potential disposal formation for a deep borehole nuclear repository concept in Israel. Site selection and repository construction are intended to ensure that waste is separated from circulating groundwater, but long-term deformation of the wellbore could potentially create fluid flow pathways. To understand how time-dependent rock strength could affect wellbore stability, we conducted creep tests under low to moderate confining pressures on retrieved core from the Zenifim formation. During creep, samples strain slowly as gradual damage accumulation progressively weakens the samples. Failure eventually occurred through the near-instantaneous formation of a shear fracture. Experimental results were used to calibrate a continuum damage poro-elastic model for sandstones. The calibrated damage-poro-elastic model successfully simulates different types of loading experiments including quasi-static and creep. The state of strain in experiments is close to yield during loading as the yield cap continuously evolves with damage accumulation. For creep tests, most damage occurs during triaxial loading. Minor damage accumulation occurs under constant load until the final stage of creep, where damage accelerates and promotes unstable fracturing.

Additive Manufacturing (AM) of porous polymeric materials, such as foams, recently became a topic of intensive research due their unique combination of low density, impressive mechanical properties, and stress dissipation capabilities. Conventional methods for fabricating foams rely on complex and stochastic processes, making it challenging to achieve precise architectural control of structured porosity. In contrast, AM provides access to a wide range of printable materials, where precise spatial control over structured porosity can be modulated during the fabrication process enabling the production of foam replacement structures (FRS). Current approaches for designing FRS are based on intuitive understanding of their properties or an extensive number of finite element method (FEM) simulations. These approaches, however, are computationally expensive and time consuming. Therefore, in this work, we present a novel methodology for determining the mechanical compression response of direct ink write (DIW) 3D printed FRS using a simple cross-sectional image. By obtaining measurement data for a relatively small number of samples, an artificial neural network (ANN) was trained, and a computer vision algorithm was used to make inferences about foam compression characteristics from a single cross-sectional image. Finally, a genetic algorithm (GA) was used to solve the inverse design problem, generating the AM printing parameters that an engineer should use to achieve a desired compression response from a DIW printed FRS. The methods developed herein present an avenue for entirely autonomous design and analysis of additively manufactured structures using artificial intelligence.

Grain-scale microstructure evolution during additive manufacturing is a complex physical process. As with traditional solidification methods of material processing (e.g. casting and welding), microstructural properties are highly dependent on the solidification conditions involved. Additive manufacturing processes however, incorporate additional complexity such as remelting, and solid-state evolution caused by subsequent heat source passes and by holding the entire build at moderately high temperatures during a build. We present a three-dimensional model that simulates both solidification and solid-state evolution phenomena using stochastic Monte Carlo and Potts Monte Carlo methods. The model also incorporates a finite-difference based thermal conduction solver to create a fully integrated microstructural prediction tool. The three modeling methods and their coupling are described and demonstrated for a model study of laser powder-bed fusion of 300-series stainless steel. The investigation demonstrates a novel correlation between the mean number of remelting cycles experienced during a build, and the resulting columnar grain sizes.

Traditional probabilistic methods for the simulation of advection-diffusion equations (ADEs) often overlook the entropic contribution of the discretization, e.g., the number of particles, within associated numerical methods. Many times, the gain in accuracy of a highly discretized numerical model is outweighed by its associated computational costs or the noise within the data. We address the question of how many particles are needed in a simulation to best approximate and estimate parameters in one-dimensional advective-diffusive transport. To do so, we use the well-known Akaike Information Criterion (AIC) and a recently-developed correction called the Computational Information Criterion (COMIC) to guide the model selection process. Random-walk and mass-transfer particle tracking methods are employed to solve the model equations at various levels of discretization. Numerical results demonstrate that the COMIC provides an optimal number of particles that can describe a more efficient model in terms of parameter estimation and model prediction compared to the model selected by the AIC even when the data is sparse or noisy, the sampling volume is not uniform throughout the physical domain, or the error distribution of the data is non-IID Gaussian.

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

There is rapidly expanding interest in exploiting the spin of valence-band holes rather than conduction-band electrons for spin qubit semiconductor circuits composed of coupled quantum dots. The hole platform offers stronger spin-orbit interaction (SOI), large difference between in-dot-plane and out-of-dot-plane g-factors, i.e. g-factor anisotropy, and a significantly reduced hyperfine coupling to nuclei in the host material. These attributes collectively can deliver fast all-electric coherent spin manipulation, efficient spin-flip inter-dot tunneling channels, a voltage tunable effective g-factor, a g-factor adjustable to nearly zero in an appropriately oriented external magnetic field, and long spin relaxation and coherence times. Here, we review our recent work on the physics of heavy holes confined in a planar GaAs/AlGaAs double quantum dot system with strong SOI. For a single-hole, we have performed resonant tunneling magneto-spectroscopy to extract spin-flip and spin-conserving tunneling strengths, implemented spin-flip Landau-Zener-Stückelberg-Majorana (LZSM) interferometry, determined the spin relaxation time T 1 as a function of magnetic field using a fast single-shot latched charge technique, electrically tuned the effective g-factor revealed by electric dipole spin resonance, and found signatures of the hyperfine interaction and dynamic nuclear polarization with holes. For two-holes, we have measured the energy spectrum in the presence of strong SOI (and so not limited by Pauli spin blockade), quantified the heavy-hole (HH) g-factor anisotropy on tilting the magnetic field, described a scheme to employ HHs whose g-factor is tunable to nearly zero for an in-plane magnetic field for a coherent photon-to-spin interface, and observed a well-defined LZSM interference pattern at small magnetic fields on pulsing through the singlet-triplet anti-crossing.

Turbulent signals are intermittent with large instantaneous fluctuations. Such large fluctuations lead to small Kolmogorov scales that are hard to resolve in numerical simulations [P. K. Yeung, K. R. Sreenivasan, and S. B. Pope, Effects of finite spatial and temporal resolution in direct numerical simulations of incompressible isotropic turbulence, Phys. Rev. Fluids 3, 064603 (2018)2469-990X10.1103/PhysRevFluids.3.064603]. The present paper follows the above basic logic, but instead of dissipation events in isotropic turbulence, we study wall-shear stress events in plane channel flow. Wall-shear stress fluctuations are increasingly more intermittent as the Reynolds number increases. Hence, one has to employ higher grid resolutions as the Reynolds number increases in order to resolve a given percentage of wall-shear stress events. The objective of this paper is to quantify effects of the grid resolutions on the rare and high intensity wall-shear stress events. We find that the standard grid resolution resolves about 99% of the wall-shear stress events at Reτ=180. A slightly higher grid resolution has to be employed in order to resolve 99% of the wall-shear stress events at higher Reynolds numbers, and if the standard grid resolution is used for, e.g., a Reτ=10000 channel flow, one resolves about 90%-95% wall-shear stress events.

Stark polarization spectroscopy is used to investigate the temporal evolution of the electric field distribution in the cathode region of a nanosecond pulsed discharge in helium at 120 Torr. The measurements are performed on the He I transition at 492.19 nm, during the early stages of the discharge formation. The experimental results are compared with the predictions of a 1D fluid model. Time-resolved ICCD images show that the discharge develops as a diffuse, cathode-directed ionization wave with a Townsend-like feature before transitioning into a glow-like structure. Near anode instabilities characterized by filament formation were observed near the high voltage electrode. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. The cathode electric field corresponding to this sheath with a thickness of 250 μm is about 40 kV cm-1. A subsequent steep increase of the discharge current leads to a further reduction of the sheath width. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field.

Integration-technology feature shrink increases computing-system susceptibility to single-event effects (SEE). While modeling SEE faults will be critical, an integrated processor's scope makes physically correct modeling computationally intractable. Without useful models, presilicon evaluation of fault-tolerance approaches becomes impossible. To incorporate accurate transistor-level effects at a system scope, we present a multiscale simulation framework. Charge collection at the 1) device level determines 2) circuit-level transient duration and state-upset likelihood. Circuit effects, in turn, impact 3) register-transfer-level architecture-state corruption visible at 4) the system level. Thus, the physically accurate effects of SEEs in large-scale systems, executed on a high-performance computing (HPC) simulator, could be used to drive cross-layer radiation hardening by design. We demonstrate the capabilities of this model with two case studies. First, we determine a D flip-flop's sensitivity at the transistor level on 14-nm FinFet technology, validating the model against published cross sections. Second, we track and estimate faults in a microprocessor without interlocked pipelined stages (MIPS) processor for Adams 90% worst case environment in an isotropic space environment.

Energy storage using lithium-ion cells dominates consumer electronics and is rapidly becoming predominant in electric vehicles and grid-scale energy storage, but the high energy densities attained lead to the potential for release of this stored chemical energy. This article introduces some of the paths by which this energy might be unintentionally released, relating cell material properties to the physical processes associated with this potential release. The selected paths focus on the anode–electrolyte and cathode–electrolyte interactions that are of typical concern for current and near-future systems. Relevant material processes include bulk phase transformations, bulk diffusion, surface reactions, transport limitations across insulating passivation layers, and the potential for more complex material structures to enhance safety. We also discuss the development, parameterization, and application of predictive models for this energy release and give examples of the application of these models to gain further insight into the development of safer energy storage systems.

Nanocrystalline Al thin films have been strained in situ in a transmission electron microscope using two separate nanomechanical techniques involving a push-to-pull device and a microelectromechanical system (MEMS) device. Deformation-induced grain growth was observed to occur via stress-assisted grain boundary migration with extensive grain growth occurring in the necked region, indicating that the increase in local stress drives the boundary migration. Under applied tensile stresses close to the ultimate tensile strength of 450 MPa for a nanocrystalline Al specimen, measured boundary migration speeds are 0.2 – 0.7 nm s−1 for grains outside necked region and increases to 2.5 nm s−1 for grains within the necked region where the local estimated tensile stresses are elevated to around 630 MPa. By tracking grain boundary motion over time, molecular dynamics simulations showed qualitative agreement in terms of pronounced grain boundary migration with the experimental observations. The combined in situ observation and molecular dynamics simulation results underscore the important role of stress-driven grain growth in plastically deforming nanocrystalline metals, leading to intergranular fracture through predominant grain boundary sliding in regions with large localized deformation.

Bayesian inference is used to calibrate a bottom-up home PLC network model with unknown loads and wires at frequencies up to 30 MHz. A network topology with over 50 parameters is calibrated using global sensitivity analysis and transitional Markov Chain Monte Carlo (TMCMC). The sensitivity-informed Bayesian inference computes Sobol indices for each network parameter and applies TMCMC to calibrate the most sensitive parameters for a given network topology. A greedy random search with TMCMC is used to refine the discrete random variables of the network. This results in a model that can accurately compute the transfer function despite noisy training data and a high dimensional parameter space. The model is able to infer some parameters of the network used to produce the training data, and accurately computes the transfer function under extrapolative scenarios.

Ferroelectric phase stability in hafnium oxide is reported to be influenced by factors that include composition, biaxial stress, crystallite size, and oxygen vacancies. In the present work, the ferroelectric performance of atomic layer deposited Hf0.5Zr0.5O2 (HZO) prepared between TaN electrodes that are processed under conditions to induce variable biaxial stresses is evaluated. The post-processing stress states of the HZO films reveal no dependence on the as-deposited stress of the adjacent TaN electrodes. All HZO films maintain tensile biaxial stress following processing, the magnitude of which is not observed to strongly influence the polarization response. Subsequent composition measurements of stress-varied TaN electrodes reveal changes in stoichiometry related to the different preparation conditions. HZO films in contact with Ta-rich TaN electrodes exhibit higher remanent polarizations and increased ferroelectric phase fractions compared to those in contact with N-rich TaN electrodes. HZO films in contact with Ta-rich TaN electrodes also have higher oxygen vacancy concentrations, indicating that a chemical interaction between the TaN and HZO layers ultimately impacts the ferroelectric orthorhombic phase stability and polarization performance. The results of this work demonstrate a necessity to carefully consider the role of electrode processing and chemistry on performance of ferroelectric hafnia films.

Abstract not provided.

Abstract not provided.

This work describes the diagnostic implementation and image processing methods to quantitatively measure diesel spray mixing injected into a high-pressure, high-temperature environment. We used a high-repetition-rate pulse-burst laser developed in-house, a high-speed CMOS camera, and optimized the optical configuration to capture Rayleigh scattering images of the vaporized fuel jets inside a constant volume chamber. The experimental installation was modified to reduce reflections and flare levels to maximize the images’ signal-to-noise ratios by anti-reflection coatings on windows and surfaces, as well as series of optical baffles. Because of the specificities of the high-speed system, several image processing techniques had to be developed and implemented to provide quantitative fuel concentration measurements. These methods involve various correction procedures such as camera linearity, laser intensity fluctuation, dynamic background flare, as well as beam-steering effects. Image inpainting was also applied to correct the Rayleigh scattering signal from large scatterers (e.g. particulates). The experiments demonstrate that applying planar laser Rayleigh scattering at high repetition rate to quantitatively resolve the mixing of fuel and ambient gases in diesel jets is challenging, but possible. The thorough analysis of the experimental uncertainty and comparisons to past data prove that such measurements can be accurate, whilst providing valuable information about the mixing processes of high-pressure diesel jets.

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.

This report will summarize the group's work to provide recommendations to secure nuclear facilities before, during and after an extreme weather event. Extreme weather events can have drastic impacts to nuclear facilities as seen by the earthquake and subsequent tsunami at the Fukushima Daiichi Nuclear Power Plant in 2011. Recent hurricanes in the United States including Hurricane Harvey demonstrate the devastating effects these storms can have on infrastructure and the surrounding communities. The group is attempting to identify the gaps that potential small modular reactor (SMR) facilities will need to address in order to provide adequate site security before, during and after extreme weather events. This effort proceeded in three parts to provide insights and recommendations to secure Small Modular Reactor facilities for extreme weather events:(1) a literature review of academic articles as well as relevant documents including the existing regulatory framework and recommendations from the IAEA, NRC, and DOE, (2) subject matter expert interviews from a wide variety of security backgrounds, and (3) modeling and simulation on a hypothetical SMR facility. Special attention was paid to the interactions between stakeholders and nuclear facility design considerations, particularly the topics of safety and security. Engineering design issues from safety and security perspectives were discussed and included in simulation. Each step informed the proceeding, with the result including full tabletop scenarios of EWE impacts to security system effectiveness on the hypothetical model. This systems-level analysis provides results to inform recommendations to secure SMR facilities.

Abstract not provided.

Abstract not provided.

In this work, we investigate the critical coupling of a single gold disk antenna with a focused beam by evaluating its absorption and scattering using spectral interferometry microcopy.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Mechanical properties on alluvium blocks and core samples were determined to support the Source Physics Experiment Dry Alluvium Geology experimental series. Because material was not available directly from the experimental location, the alluvium blocks and core samples are intended to serve as surrogate material . P - and S - wave velocity was measured on cubes cut from the alluvium blocks and core with the intention to study variation from water content and measured direction (material anisotropy). Indirect tensile tests were conducted dry and with moisture ranging from 6 to 9.1%. For the range of water content tested, increasing moisture level resulted in slower P - and S - wave velocities. P - and S - wave variability is less influenced by material heterogeneity than moisture content. P - wave velocity ranges from 629 m/s to 2599 m/s and S - wave velocity ranges from 288 m/s to 1200 m/s. Counter to the velocity measurement findings, material variability on indirect tensile strength has a greater effect than moisture content. Compared to dry strength and at moisture levels from 6 to 9% the block's tensile strength was lowered by at least a factor of 5. Indirect tensile strength for the first block averaged 0.35 MPa and 0.25 MPa for dry and 8.9% moisture respectively. For the second block indirect tensile strength averaged 0.05 MPa for both dry and 6.4% moisture.

Abstract not provided.

This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users' Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users' Guide.

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers; A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models; Device models that are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only); Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase—a message passing parallel implementation—which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

This manual describes the installation and use of the Xyce™ XDM Netlist Translator. XDM simplifies the translation of netlists generated by commercial circuit simulator tools into Xyce-compatible netlists. XDM currently supports translation from PSpice, HSPICE, and Spectre netlists into Xyce™ netlists.

Abstract not provided.

Abstract not provided.

Abstract not provided.