Publications

Cyber-physical systems have behaviour that crosses domain boundaries during events such as planned operational changes and malicious disturbances. Traditionally, the cyber and physical systems are monitored separately and use very different toolsets and analysis paradigms. The security and privacy of these cyber-physical systems requires improved understanding of the combined cyber-physical system behaviour and methods for holistic analysis. Therefore, the authors propose leveraging clustering techniques on cyber-physical data from smart grid systems to analyse differences and similarities in behaviour during cyber-, physical-, and cyber-physical disturbances. Since clustering methods are commonly used in data science to examine statistical similarities in order to sort large datasets, these algorithms can assist in identifying useful relationships in cyber-physical systems. Through this analysis, deeper insights can be shared with decision-makers on what cyber and physical components are strongly or weakly linked, what cyber-physical pathways are most traversed, and the criticality of certain cyber-physical nodes or edges. This paper presents several types of clustering methods for cyber-physical graphs of smart grid systems and their application in assessing different types of disturbances for informing cyber-physical situational awareness. The collection of these clustering techniques provide a foundational basis for cyber-physical graph interdependency analysis.

Solid–water interfaces are crucial for clean water, conventional and renewable energy, and effective nuclear waste management. However, reflecting the complexity of reactive interfaces in continuum-scale models is a challenge, leading to oversimplified representations that often fail to predict real-world behavior. This is because these models use fixed parameters derived by averaging across a wide physicochemical range observed at the molecular scale. Recent studies have revealed the stochastic nature of molecular-level surface sites that define a variety of reaction mechanisms, rates, and products even across a single surface. To bridge the molecular knowledge and predictive continuum-scale models, we propose to represent surface properties with probability distributions rather than with discrete constant values derived by averaging across a heterogeneous surface. This conceptual shift in continuum-scale modeling requires exponentially rising computational power. By incorporating our molecular-scale understanding of solid–water interfaces into continuum-scale models we can pave the way for next generation critical technologies and novel environmental solutions.

Additive manufacturing (AM) technology, specifically 3D printing, holds great promise for in-orbit manufacturing. In-space printing can significantly reduce the mass, cost, and risk of long-term space exploration by enabling replacement parts to be made as needed and reducing dependence on Earth. However, printing in a zero-gravity environment poses challenges due to the absence of a rigid ground for the print platform, which can result in vibrational and rotational forces that may impact printing integrity. To address this issue, this paper proposes a novel linear magnetic position tracking algorithm, named Navigation Integrating Magnets By Linear Estimation (NIMBLE), for dynamic vibration compensation during 3D printing of truss structures in space. Compared to the most commonly used nonlinear optimization method, the NIMBLE algorithm is more than two orders of magnitude faster. With only a single 3-axis magnet sensor and a small NdFeB magnet, the NIMBLE algorithm provides a simple and easily implemented tracking solution for in-orbit 3D printing.

Information security and computing, two critical technological challenges for post-digital computation, pose opposing requirements – security (encryption) requires a source of unpredictability, while computing generally requires predictability. Each of these contrasting requirements presently necessitates distinct conventional Si-based hardware units with power-hungry overheads. This work demonstrates Cu0.3Te0.7/HfO2 (‘CuTeHO’) ion-migration-driven memristors that satisfy the contrasting requirements. Under specific operating biases, CuTeHO memristors generate truly random and physically unclonable functions, while under other biases, they perform universal Boolean logic. Using these computing primitives, this work experimentally demonstrates a single system that performs cryptographic key generation, universal Boolean logic operations, and encryption/decryption. Circuit-based calculations reveal the energy and latency advantages of the CuTeHO memristors in these operations. This work illustrates the functional flexibility of memristors in implementing operations with varying component-level requirements.

Material Testing 2.0 (MT2.0) is a paradigm that advocates for the use of rich, full-field data, such as from digital image correlation and infrared thermography, for material identification. By employing heterogeneous, multi-axial data in conjunction with sophisticated inverse calibration techniques such as finite element model updating and the virtual fields method, MT2.0 aims to reduce the number of specimens needed for material identification and to increase confidence in the calibration results. To support continued development, improvement, and validation of such inverse methods—specifically for rate-dependent, temperature-dependent, and anisotropic metal plasticity models—we provide here a thorough experimental data set for 304L stainless steel sheet metal. The data set includes full-field displacement, strain, and temperature data for seven unique specimen geometries tested at different strain rates and in different material orientations. Commensurate extensometer strain data from tensile dog bones is provided as well for comparison. We believe this complete data set will be a valuable contribution to the experimental and computational mechanics communities, supporting continued advances in material identification methods.

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. The CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

The use of photon active interrogation to detect special nuclear material has held significant theoretical promise, as the interrogating source particles, photons, are fundamentally different from one of the main signatures of special nuclear material: neutrons produced in nuclear fission. However, neutrons produced by photonuclear reactions in the accelerator target, collimator, and environment can obscure the fission neutron signal. These (γ,n) neutrons could be discriminated from fission neutrons by their energy spectrum, but common detectors sensitive to the neutron spectrum, like organic scintillators, are typically hampered by the intense photon background characteristic of photon-based active interrogation. In contrast, high-pressure 4He-based scintillation detectors are well -suited to photon active interrogation, as they are similarly sensitive to fast neutrons and can measure their spectrum, but show little response to gamma rays. In this work, a photon active interrogation system utilizing a 4He scintillation detector and a 9 MeV linac-bremsstrahlung x-ray source was experimentally evaluated. The detector was shown to be capable of operating in intense gamma-ray environments and detecting photofission neutrons from 238U when interrogated by this x-ray source. The photofission neutrons show clear spectral separation from (γ,n) neutrons produced in lead, a common shielding material.

The bulk-boundary correspondence in topological crystalline insulators (TCIs) links the topological properties of the bulk to robust observables on the edges, e.g., the existence of robust edge modes or fractional charge. In one dimension, TCIs protected by reflection symmetry have been realized in a variety of systems in which each unit cell has spatially distributed degrees of freedom (SDOF). However, these realizations exhibit sensitivity of the resulting edge modes to variations in edge termination and to the local breaking of the protective spatial symmetries by inhomogeneity. Here we demonstrate topologically protected edge states in a monoatomic, orbital-based TCI that mitigates both of these issues. By collapsing all SDOF within the unit cell to a singular point in space, we eliminate the ambiguity in unit-cell definition and hence remove a prominent source of boundary termination variability. The topological observables are also more tolerant to disorder in the orbital energies. To validate this concept, we experimentally realize a lattice of mechanical resonators where each resonator acts as an "atom"that harbors two key orbital degrees of freedom having opposite reflection parity. Our measurements of this system provide direct visualization of the sp-hybridization between orbital modes that leads to a nontrivial band inversion in the bulk.

In this study we present a replication method to determine surface roughness and to identify surface features when a sample cannot be directly analyzed by conventional techniques. As a demonstration, this method was applied to an unused spent nuclear fuel dry storage canister to determine variation across different surface features. In this study, an initial material down-selection was performed to determine the best molding agent and determined that non-modified Polytek PlatSil23-75 provided the most accurate representation of the surface while providing good usability. Other materials that were considered include Polygel Brush-On 35 polyurethane rubber (with and without Pol-ease 2300 release agent), Polytek PlatSil73-25 silicone rubber (with and without PlatThix thickening agent and Pol-ease 2300 release agent), and Express STD vinylpolysiloxane impression putty. The ability of PlatSil73-25 to create an accurate surface replica was evaluated by creating surface molds of several locations on surface roughness standards representing ISO grade surfaces N3, N5, N7, and N8. Overall, the molds were able to accurately reproduce the expected roughness average (Ra) values, but systematically over-estimated the peak-valley maximum roughness (Rz) values. Using a 3D printed sample cell, several locations across the stainless steel spent nuclear fuel canister were sampled to determine the surface roughness. These measurements provided information regarding variability in normal surface roughness across the canister as well as a detailed evaluation on specific surface features (e.g., welds, grind marks, etc.). The results of these measurements can support development of dry storage canister ageing management programs, as surface roughness is an important factor for surface dust deposition and accumulation. This method can be applied more broadly to different surfaces beyond stainless steel to provide rapid, accurate surface replications for analytical evaluation by profilometry.

Scalable quantum processors require high-fidelity universal quantum logic operations in a manufacturable physical platform. Donors in silicon provide atomic size, excellent quantum coherence and compatibility with standard semiconductor processing, but no entanglement between donor-bound electron spins has been demonstrated to date. Here we present the experimental demonstration and tomography of universal one- and two-qubit gates in a system of two weakly exchange-coupled electrons, bound to single phosphorus donors introduced in silicon by ion implantation. We observe that the exchange interaction has no effect on the qubit coherence. We quantify the fidelity of the quantum operations using gate set tomography (GST), and we use the universal gate set to create entangled Bell states of the electrons spins, with fidelity 91.3 ± 3.0%, and concurrence 0.87 ± 0.05. These results form the necessary basis for scaling up donor-based quantum computers.

In the machine learning problem of multilabel classification, the objective is to determine for each test instance which classes the instance belongs to. In this work, we consider an extension of multilabel classification, called multilabel proportion prediction, in the context of radioisotope identification (RIID) using gamma spectra data. We aim to not only predict radioisotope proportions, but also identify out-of-distribution (OOD) spectra. We achieve this goal by viewing gamma spectra as discrete probability distributions, and based on this perspective, we develop a custom semi-supervised loss function that combines a traditional supervised loss with an unsupervised reconstruction error function. Our approach was motivated by its application to the analysis of short-lived fission products from spent nuclear fuel. In particular, we demonstrate that a neural network model trained with our loss function can successfully predict the relative proportions of 37 radioisotopes simultaneously. The model trained with synthetic data was then applied to measurements taken by Pacific Northwest National Laboratory (PNNL) to conduct analysis typically done by subject-matter experts. We also extend our approach to successfully identify when measurements are OOD, and thus should not be trusted, whether due to the presence of a novel source or novel proportions.

Modern lens designs are capable of resolving greater than 10 gigapixels, while advances in camera frame-rate and hyperspectral imaging have made data acquisition rates of Terapixel/second a real possibility. The main bottlenecks preventing such high data-rate systems are power consumption and data storage. In this work, we show that analog photonic encoders could address this challenge, enabling high-speed image compression using orders-of-magnitude lower power than digital electronics. Our approach relies on a silicon-photonics front-end to compress raw image data, foregoing energy-intensive image conditioning and reducing data storage requirements. The compression scheme uses a passive disordered photonic structure to perform kernel-type random projections of the raw image data with minimal power consumption and low latency. A back-end neural network can then reconstruct the original images with structural similarity exceeding 90%. This scheme has the potential to process data streams exceeding Terapixel/second using less than 100 fJ/pixel, providing a path to ultra-high-resolution data and image acquisition systems.

Composite materials with different microstructural material symmetries are common in engineering applications where grain structure, alloying and particle/fiber packing are optimized via controlled manufacturing. In fact these microstructural tunings can be done throughout a part to achieve functional gradation and optimization at a structural level. To predict the performance of particular microstructural configuration and thereby overall performance, constitutive models of materials with microstructure are needed. In this work we provide neural network architectures that provide effective homogenization models of materials with anisotropic components. These models satisfy equivariance and material symmetry principles inherently through a combination of equivariant and tensor basis operations. We demonstrate them on datasets of stochastic volume elements with different textures and phases where the material undergoes elastic and plastic deformation, and show that the these network architectures provide significant performance improvements.

The current present in a galvanic couple can define its resistance or susceptibility to corrosion. However, as the current is dependent upon environmental, material, and geometrical parameters it is experimentally costly to measure. To reduce these costs, Finite Element (FE) simulations can be used to assess the cathodic current but also require experimental inputs to define boundary conditions. Due to these challenges, it is crucial to accelerate predictions and accurately predict the current output for different environments and geometries representative of in-service conditions. Machine learned surrogate models provides a means to accelerate corrosion predictions. However, a one-time cost is incurred in procuring the simulation and experimental dataset necessary to calibrate the surrogate model. Therefore, an active learning protocol is developed through calibration of a low-cost surrogate model for the cathodic current of an exemplar galvanic couple (AA7075-SS304) as a function of environmental and geometric parameters. The surrogate model is calibrated on a dataset of FE simulations, and calculates an acquisition function that identifies specific additional inputs with the maximum potential to improve the current predictions. This is accomplished through a staggered workflow that not only improves and refines prediction, but identifies the points at which the most information is gained, thus enabling expansion to a larger parameter space. The protocols developed and demonstrated in this work provide a powerful tool for screening various forms of corrosion under in-service conditions.

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. A bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

Downhole logging tools are commonly used to characterize multi-thousand-foot geothermal wells. The elevated temperatures, pressures, and harsh chemical environments present significant challenges for the long-term operation of these tools, especially when real-time data transmission to the surface is required via data cable lines. Teflon-based single or multi-conductor cables with grease-filled cable heads are typically used for downhole tools. However, over extended periods of operation, the grease used to seal the conductors can slowly dissolve into the well fluid, creating electrical shorts and disabling data transmission. Additionally, when temperatures exceed 260 °C, Teflon can soften, potentially allowing parallel conductors to make contact and cause shorts. Between 2009 and 2015, Draka Cableteq USA, now part of the Prysmian Group, developed a multi-conductor/fiber cable and a four-conductor cable capable of operating above 300 °C. While a full study was conducted on the conductor/fiber cable, the evaluation of the four-conductor cable remained incomplete. With the increasing need for long-term high-temperature (HT) operation of logging tools, Sandia National Laboratories is now completing the evaluation of the four-conductor cable. The four-conductor cable has two major novel aspects. Firstly, its glass braid insulation can operate above 300 °C, eliminating the potential for shorts. Secondly, the insulated conductors are encased in metal tubing along the full length of the cable, creating a high-pressure seal between the cable and the tool. This metal tubing eliminates the need for a grease seal, a major limiting factor in the operation time of common cable lines. Sandia National Laboratories will conduct multiple tests to characterize the cable at temperatures above 300 °C and pressures up to 5,000 psi. This cable would enable tools to operate continuously at elevated temperatures, pressures, and in harsh fluids for extended periods, potentially lasting months.

Abstract not provided.

Herein, we report on the ultrafast photodissociation of nickel tetracarbonyl─a prototypical metal-ligand model system─at 197 nm. Using mid-infrared transient absorption spectroscopy to probe the bound C≡O stretching modes, we find evidence for the picosecond time scale production of highly vibronically excited nickel dicarbonyl and nickel monocarbonyl, in marked contrast with a prior investigation at 193 nm. Further spectral evolution with a 50 ps time constant suggests an additional dissociation step; the absence of any corresponding growth in signal strongly indicates the production of bare Ni, a heretofore unreported product from single-photon excitation of nickel tetracarbonyl. Thus, by probing the deep UV-induced photodynamics of a prototypical metal carbonyl, this Letter adds time-resolved spectroscopic signatures of these dynamics to the sparse literature at high excitation energies.

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. We then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

Carbon dots have attracted widespread interest for sensing applications based on their low cost, ease of synthesis, and robust optical properties. We investigate structure-function evolution on multiemitter fluorescence patterns for model carbon-nitride dots (CNDs) and their implications on trace-level sensing. Hydrothermally synthesized CNDs with different reaction times were used to determine how specific functionalities and their corresponding fluorescence signatures respond upon the addition of trace-level analytes. Archetype explosives molecules were chosen as a testbed due to similarities in substituent groups or inductive properties (i.e., electron withdrawing), and solution-based assays were performed using ratiometric fluorescence excitation-emission mapping (EEM). Analyte-specific quenching and enhancement responses were observed in EEM landscapes that varied with the CND reaction time. We then used self-organizing map models to examine EEM feature clustering with specific analytes. The results reveal that interactions between carbon-nitride frameworks and molecular-like species dictate response characteristics that may be harnessed to tailor sensor development for specific applications.

There is growing interest in material candidates with properties that can be engineered beyond traditional design limits. Compositionally complex oxides (CCO), often called high entropy oxides, are excellent candidates, wherein a lattice site shares more than four cations, forming single-phase solid solutions with unique properties. However, the nature of compositional complexity in dictating properties remains unclear, with characteristics that are difficult to calculate from first principles. Here, compositional complexity is demonstrated as a tunable parameter in a spin-transition oxide semiconductor La1− x(Nd, Sm, Gd, Y)x/4CoO3, by varying the population x of rare earth cations over 0.00≤ x≤ 0.80. Across the series, increasing complexity is revealed to systematically improve crystallinity, increase the amount of electron versus hole carriers, and tune the spin transition temperature and on-off ratio. At high a population (x = 0.8), Seebeck measurements indicate a crossover from hole-majority to electron-majority conduction without the introduction of conventional electron donors, and tunable complexity is proposed as new method to dope semiconductors. First principles calculations combined with angle resolved photoemission reveal an unconventional doping mechanism of lattice distortions leading to asymmetric hole localization over electrons. Thus, tunable complexity is demonstrated as a facile knob to improve crystallinity, tune electronic transitions, and to dope semiconductors beyond traditional means.

Harmonic and subharmonic RF injection locking is demonstrated in a terahertz (THz) quantum-cascade vertical-external-cavity surface-emitting laser (QC-VECSEL). By tuning the RF injection frequency around integer multiples and submultiples of the cavity round-trip frequency, different harmonic and subharmonic orders can be excited in the same device. Modulation-dependent behavior of the device has been studied with recorded lasing spectral broadening and locking bandwidths in each case. In particular, harmonic injection locking results in the observation of harmonic spectra with bandwidths over 200 GHz. A semiclassical Maxwell-density matrix formalism has been applied to interpret QC-VECSEL dynamics, which aligns well with experimental observations.

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm² system that included dehydrated and hydrated regions in the same interlayer.

The stochastic weighted particle method (SWPM) is a generalization of the Direct Simulation Monte Carlo (DSMC) method where particle weights are variable and dynamic. SWPM is backed by a strong theoretical foundation but has not been critically evaluated for problems of practical interest. A thorough assessment of SWPM for boundary-driven flows reveals significant numerical artifacts near the boundary, notably a diverging heat flux. To correct the boundary heat flux, two modifications to SWPM are proposed: separated grouping and a spatially-dependent weight transfer function. To gauge the relative efficiency of SWPM in comparison to DSMC, a high-Mach-number wheel flow which forms a strong density gradient is also simulated.

Researchers are exploring adding wave energy converters to existing oceanographic buoys to provide a predictable source of renewable power. A ”pitch resonator” power take-off system has been developed that generates power using a geared flywheel system designed to match resonance with the pitching motion of the buoy. However, the novelty of the concept leaves researchers uncertain about various design aspects of the system. This work presents a novel design study of a pitch resonator to inform design decisions for an upcoming deployment of the system. The assessment uses control co-design via WecOptTool to optimize control trajectories for maximal electrical power production while varying five design parameters of the pitch resonator. Given the large search space of the problem, the control trajectories are optimized within a Monte Carlo analysis to identify optimal designs, followed by parameter sweeps around the optimum to identify trends between the design parameters. The gear ratio between the pitch resonator spring and flywheel are found to be the most sensitive design variables to power performance. The assessment also finds similar power generation for various sizes of resonator components, suggesting that correctly designing for optimal control trajectories at resonance is more critical to the design than component sizing.