Publications

Abstract not provided.

We consider estimating the discretization error in a nonlinear functional J(u) in the setting of an abstract variational problem: find u∈V such that B(u,φ)=L(φ)∀φ∈V, as approximated by a Galerkin finite element method. Here, V is a Hilbert space, B(⋅,⋅) is a bilinear form, and L(⋅) is a linear functional. We consider well-known error estimates η of the form J(u)−J(uh)≈η=L(z)−B(uh,z), where uh denotes a finite element approximation to u, and z denotes the solution to an auxiliary adjoint variational problem. We show that there exist nonlinear functionals for which error estimates of this form are not reliable, even in the presence of an exact adjoint solution z. An estimate η is said to be reliable if there exists a constant C∈R>0 independent of uh such that |J(u)−J(uh)|≤C|η|. We present several example pairs of bilinear forms and nonlinear functionals where reliability of η is not achieved.

Binary Code Similarity Analysis (BCSA) has a wide spectrum of applications, including plagiarism detection, vulnerability discovery, and malware analysis, thus drawing significant attention from the security community. However, conventional techniques often face challenges in balancing both accuracy and scalability simultaneously. To overcome these existing problems, a surge of deep learning-based work has been recently proposed. Unfortunately, many researchers still find it extremely difficult to conduct relevant studies or extend existing approaches. First, prior work typically relies on proprietary benchmark without making the entire dataset publicly accessible. Consequently, a large-scale, well-labeled dataset for binary code similarity analysis remains precious and scarce. Moreover, previous work has primarily focused on comparing at the function level, rather than exploring other finer granularities. Therefore, we argue that the lack of a fine-grained dataset for BCSA leaves a critical gap in current research. To address these challenges, we construct a benchmark dataset for fine-grained binary code similarity analysis called BinSimDB, which contains equivalent pairs of smaller binary code snippets, such as basic blocks. Specifically, we propose BMerge and BPair algorithms to bridge the discrepancies between two binary code snippets caused by different optimization levels or platforms. Furthermore, we empirically study the properties of our dataset and evaluate its effectiveness for the BCSA research. The experimental results demonstrate that BinSimDB significantly improves the performance of binary code similarity comparison.

Scripta Materialia

Journal of Computational and Applied Mathematics

Although synthetic biology can produce valuable chemicals in a renewable manner, its progress is still hindered by a lack of predictive capabilities. Media optimization is a critical, and often overlooked, process which is essential to obtain the titers, rates and yields needed for commercial viability. Here, we present a molecule- and host-agnostic active learning process for media optimization that is enabled by a fast and highly repeatable semi-automated pipeline. Its application yielded 60% and 70% increases in titer, and 350% increase in process yield in three different campaigns for flaviolin production in Pseudomonas putida KT2440. Explainable Artificial Intelligence techniques pinpointed that, surprisingly, common salt (NaCl) is the most important component influencing production. The optimal salt concentration is very high, comparable to seawater and close to the limits that P. putida can tolerate. The availability of fast Design-Build-Test-Learn (DBTL) cycles allowed us to show that performance improvements for active learning are rarely monotonous. This work illustrates how machine learning and automation can change the paradigm of current synthetic biology research to make it more effective and informative, and suggests a cost-effective and underexploited strategy to facilitate the high titers, rates and yields essential for commercial viability.

Material design and accessible manufacturing are often at odds with each other, calling for creative solutions to adapt high-performance materials to available processes. This challenge is represented well by in-mold electronics, an innovative approach to the manufacture of 3D circuitry and electronic components that offers game-changing advantages. In-mold electronics relies on vacuum forming processes, which are historically limited to thermoplastics. Extending these methods to include thermosets would enable manufacturing of robust components with desirable properties. Here, we provide a solution to make thermoset materials amenable to vacuum forming. Specifically, an ambient polymerization is used to transition a liquid monomeric solution to an elastomeric gel. These free-standing gels can then be vacuum formed, and the reaction can be completed via frontal polymerization. Thermoset materials produced with this method have properties that provide benefits over traditionally employed thermoplastic substrates and enable 3D device integration into environmentally demanding architectural, automotive, and extraterrestrial structures.

Magnetic anisotropy is essential for many applications of ferromagnetic/ferrimagnetic materials, including permanent magnets and magnetic recording media. Attempts have been made recently to build up 3-D nanoparticle and quantum dot assemblies, however, it is not understood yet if a nanoparticle assembly can possess high magnetic anisotropy with low anisotropic materials. In this article, we report our discovery of high magnetic anisotropy resulted from Fe3O4 nanoparticle chains. We started with closely-packed nanoparticle assemblies of spherical Fe3O4 nanoparticles that exhibit low magnetocrystalline anisotropy and shape anisotropy, and corresponding negligible coercivity. When the nanoparticle assemblies are compressed under pressure, they form bundles or arrays that consist of Fe3O4 chains with a length scale of several hundred nanometers. Magnetic measurements show that these Fe3O4 chain arrays possess a high uniaxial magnetic anisotropy (Keff ~ 2.9×105J/m³) and significant magnetic coercivity. Our simulations reveal that interparticle magnetic dipolar interactions contribute to this type of superstructure magnetic anisotropy. This study demonstrates the feasibility and approaches to create “patterned” high magnetic anisotropy in nanoparticle superstructures/assemblies.

A simple co-design example in a reduced parameter space is presented for an oscillating flap device. Initially, the WEC geometry and mass properties are considered along with drivetrain gear ratio, inertia, motor constant and stiffness under both PI and optimal control. This parameter space is reduced to those to which performance is most sensitive for a fixed geometry. The gear ratio, drivetrain stiffness, and flap mass are found to be the most impactful design criteria as they can create orders of magnitude variations in power performance. The performance of the optimized system is compared with several sub-optimal variants in terms of electrical and mechanical power capture, transmission coefficients, and transducer power gain. Notably, though substantial power capture improvements are demonstrated when an optimal controller is employed, this power capture remains sensitive to appropriate selections of drivetrain and flap design parameters, implying that control co-design procedures remain necessary for high-performing WECs. A number of practical caveats and extensions to the presented co-design methodology are suggested, including the characterization of system static friction, especially in the presence of high gear ratios.

Viral hemorrhagic Lassa fever (LF), caused by Lassa virus (LASV), is a significant public health concern endemic in West Africa with high morbidity and mortality rates, limited treatment options, and potential for international spread. Despite advances in interrogating its epidemiology and clinical manifestations, the molecular mechanisms driving pathogenesis of LASV and other arenaviruses remain incompletely understood. This review synthesizes current knowledge regarding the role of LASV host-virus interactions in mediating the pathogenesis of LF, with emphasis on interactions between viral and host proteins. Through investigation of these critical protein–protein interactions, we identify potential therapeutic targets and discuss their implications for development of medical countermeasures including antiviral drugs. This review provides an update in recent literature of significant LASV host-virus interactions important in informing the development of targeted therapies and improving clinical outcomes for LF patients. Knowledge gaps are highlighted as opportunities for future research efforts that would advance the field of LASV and arenavirus pathogenesis.

Recent trends within computational and data sciences show an increasing recognition and adoption of computational workflows as tools for productivity and reproducibility that also democratize access to platforms and processing know-how. As digital objects to be shared, discovered, and reused, computational workflows benefit from the FAIR principles, which stand for Findable, Accessible, Interoperable, and Reusable. The Workflows Community Initiative’s FAIR Workflows Working Group (WCI-FW), a global and open community of researchers and developers working with computational workflows across disciplines and domains, has systematically addressed the application of both FAIR data and software principles to computational workflows. We present recommendations with commentary that reflects our discussions and justifies our choices and adaptations. These are offered to workflow users and authors, workflow management system developers, and providers of workflow services as guidelines for adoption and fodder for discussion. The FAIR recommendations for workflows that we propose in this paper will maximize their value as research assets and facilitate their adoption by the wider community.

Ship tracks, long thin artificial cloud features formed from the pollutants in ship exhaust, are satellite-observable examples of aerosol-cloud interactions (ACI) that can lead to increased cloud albedo and thus increased solar reflectivity, phenomena of interest in solar radiation management. In addition to ship tracks being of interest to meteorologists and policy makers, their observed cloud perturbations provide benchmark evidence of ACI that remain poorly captured by climate models. To broadly analyze the effects of ship tracks, high-resolution satellite imagery data highlighting their presence are required. To support this, we provide a hand labelled dataset to serve as a benchmark for a variety of subsequent analyses. Established from a previous dataset that identified ship track presence using NASA’s MODIS Aqua satellite imager, our first-of-its-kind dataset is comprised of image masks: capturing full ship track regions, including their contours, emission points and dispersive patterns. In total, 300 images, or around 2,500 masked ship tracks, observed under varying conditions are provided, and may facilitate training of machine learning algorithms to automate extraction.

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width deq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: deq decreases for θ>30°, deq remains unchanged for θ=30°, and deq increases for θ<30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. These findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Nuclear Materials and Energy

Abstract not provided.

International Journal of Engineering Sciences

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of −11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of −35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. New experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

A predictive analytical model is presented for the electrical conductivity of multi-principal-element alloys (MPEAs), including those containing aluminum, transition metals, and refractory metals. Given that the lattice parameter of the Wigner-Seitz cell of an MPEA is similarly variable to a bulk metallic glass, it is postulated that electron scattering can be approximated by a series of two-level systems. The resulting reduced-order model enabled an accurate determination of electrical resistivity and electron thermal conductivity based on the scattering of electrons in a two-level system across a Bloch-potential-based virtual crystal approximation. Model results are compared to experimental four-point probe electrical resistivity measurements between 300 K and 700 K for Al0.3CoCrCuFeNi, CoCrFeMnNi, (CoCrFeMnNi)0.98W0.02, (CoCrFeMnNi)0.95W0.05, and Nb4Ta4V3Ti, for model validation.

Background: Non-invasive simulations of coronary hemodynamics have improved clinical risk stratification and treatment outcomes for coronary artery disease, compared to relying on anatomical imaging alone. However, simulations typically use empirical approaches to distribute total coronary flow amongst the arteries in the coronary tree, which ignores patient variability, the presence of disease, and other clinical factors. Further, uncertainty in the clinical data often remains unaccounted for in the modeling pipeline. Objective: We present an end-to-end uncertainty-aware pipeline to (1) personalize coronary flow simulations by incorporating vessel-specific coronary flows as well as cardiac function; and (2) predict clinical and biomechanical quantities of interest with improved precision, while accounting for uncertainty in the clinical data. Methods: We assimilate patient-specific measurements of myocardial blood flow from clinical CT myocardial perfusion imaging to estimate branch-specific coronary artery flows. Simulated noise in the clinical data is used to estimate the joint posterior distributions of the model parameters using adaptive Markov Chain Monte Carlo sampling. Additionally, the posterior predictive distribution for the relevant quantities of interest is determined using a new approach combining multi-fidelity Monte Carlo estimation with non-linear, data-driven dimensionality reduction. This leads to improved correlations between high- and low-fidelity model outputs. Results: Our framework accurately recapitulates clinically measured cardiac function as well as branch-specific coronary flows under measurement noise uncertainty. We observe substantial reductions in confidence intervals for estimated quantities of interest compared to single-fidelity Monte Carlo estimation and state-of-the-art multi-fidelity Monte Carlo methods. This holds especially true for quantities of interest that showed limited correlation between the low- and high-fidelity model predictions. In addition, the proposed multi-fidelity Monte Carlo estimators are significantly cheaper to compute than traditional estimators, under a specified confidence level or variance. Conclusions: The proposed pipeline for personalized and uncertainty-aware predictions of coronary hemodynamics is based on routine clinical measurements and recently developed techniques for CT myocardial perfusion imaging. The proposed pipeline offers significant improvements in precision and reduction in computational cost.

This paper introduces the topics that need consideration to allow an effective feasibility analysis of electron beams and/or X-rays technologies prior to potential investment and implementation of these technologies for industrial applications. It also highlights considerations for planning, developing and construction of an electron beam or X-ray facility The focus of this paper is on business topics and it provides details of where additional technical information can be found. The paper also describes current applications for electron beam and X-ray technologies to demonstrate their potential.

The demand for low-cost, low-energy, and highly selective gas capture and separations is an ongoing driver of porous material development. Porous liquids have been identified as a promising gas separation material by creating permanent porosity in inorganic solvents through inclusion of nanoporous materials that sterically exclude solvent from their internal porosity. Among the nanoporous materials that can be used to form porous liquids, porous-organic cages (POCs) have been one of the most popular due to the inherent tunability of POCs. “Scrambled” POCs with varying functionalities on the POC vertices have been developed and incorporated into porous liquid compositions, increasing their gas adsorption capacity. An unexplored avenue to tailor the properties of porous liquids is through scrambling the functionality of the core of the POC. Therefore, we have synthesized a new POC, a CC3-OH derivative with scrambled hydroxides on the core and evaluated the impact on the CO2 uptake capacity in silicon oil-based porous liquids. Core scrambling of the POC resulted in a twofold increase CO2 adsorption capacity in the porous liquid, an emergent property that is a dramatic increase beyond a linear combination of the gas adsorption capacity of the neat solvent and the POC. Density functional theory modeling of the CC3 POC and its hydroxide-based derivatives identified that free rotation of the linker hydroxide allowed for forced interaction between the CO2 molecule and the hydroxide in the pore window. Solvation of the POC may release scrambled core hydroxides from intramolecular bonding with a neighboring imine, allowing for increased gas uptake in the porous liquid over the neat POC. These results identify a key structural relationship of POCs that enables emergent properties in porous liquids and can guide future development of liquid phase gas capture and separation materials for environmental and industrial applications.

Distributed acoustic sensing (DAS) systems, which use existing telecommunication fibers, offer high-resolution capabilities ideal for recording anthropogenic sources. However, the complexity of urban environments and the large amount of data recorded by DAS require automated methods to efficiently detect and categorize anthropogenic sources. We evaluate how well three machine learning models (k-nearest neighbor [k-NN], convolutional neural networks, and recurrent-convolutional neural networks) can identify various anthropogenic sources recorded by DAS. Our findings reveal that both k-NN and neural network methods perform well in high signal-to-noise ratio (SNR) settings. However, their accuracy decreases at SNRs < 4. We also use Kalman filtering, a form of Bayesian inference, on backprojected locations of these sources to recover locations that generally fall within standard smartphone Global Positioning System errors. By combining machine learning and Kalman filter results, we calculate a multidimensional model of moving anthropogenic sources. These results demonstrate the potential of DAS data in urban seismology for accurately identifying and locating such sources. Depending on the research objectives, these sources can be further studied or filtered out to improve the quality of seismic data for earthquake studies. Such methods provide a valuable tool for urban seismology and seismic hazard analysis.

Remote coastal communities, which could be early adopters of wave energy projects, have concerns over costs, conflicts, and potential risks of development. Designers and developers are challenged to address these community concerns as they continue to develop wave energy technologies. One potential means of reducing costs, conflicts, and risks, especially for demonstration and pilot-scale projects, could be planning a deployment that operates for only a portion of the year—a seasonal deployment. In this paper we examine the impacts of a seasonal deployment in terms of cost, electricity production, operations and maintenance, environmental impacts, and community benefits. We take a holistic, comparative approach to feasibility that can be replicated for other comparative studies. We estimate electricity production using a point absorber WEC modeled near Sitka, AK, USA and optimized for the given sea conditions. We determine that, for remote community sized projects, seasonal deployments could result in small cost savings (less than 10 %), but larger decreases in annual energy production (around 30 % for our case study area). Seasonal deployments could be preferable in places with seasonal energy needs, if failures and device access become a major hindrance to wave energy technology development, or as a cautionary approach to introducing new technology to the oceans. We also determine that a highly seasonal wave resource is not necessarily a requirement for seasonal deployments to be considered. Seasonal deployments are an alternative to year-round deployments that can be considered in places where marine spatial conflict is a seasonal concern.

Current methods for modeling non-adiabatic molecular dynamics face fundamental limitations when treating geometric phase effects, quantum mechanical phenomena where nuclear wavepackets acquire phase shifts when encircling conical intersections. Existing approaches either neglect these effects entirely or rely on potential energy surfaces arising from the Born-Oppenheimer approximation, which introduce artificial singularities and can overestimate geometric phase contributions. This project deve

Abstract not provided.