Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “2D class averages”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

A self-supervised workflow for particle picking in cryo-EM

High-resolution single-particle cryo-EM data analysis relies on accurate particle picking. To facilitate the particle picking process, a self-supervised workflow has been developed. This includes an iterative strategy, which uses a 2D class average to improve training particles, and a progressively improved convolutional neural network for particle picking. To automate the selection of particles, a threshold is defined (%/Res) using the ratio of percentage class distribution and resolution as a cutoff. This workflow has been tested using six publicly available data sets with different particle sizes and shapes, and can automatically pick particles with minimal user input. The picked particles support high-resolution reconstructions at 3.0 Å or better. This workflow is a step towards automated single-particle cryo-EM data analysis at the stage of particle picking. It may be used in conjunction with commonly used single-particle analysis packages such as Relion , cryoSPARC , cisTEM , SPHIRE and EMAN2 .

2D class averages↗

Design of Diverse Asymmetric Pockets in De Novo Homo-oligomeric Proteins

A challenge for design of protein–small-molecule recognition is that incorporation of cavities with size, shape, and composition suitable for specific recognition can considerably destabilize protein monomers. This challenge can be overcome through binding pockets formed at homo-oligomeric interfaces between folded monomers. Interfaces surrounding the central homo-oligomer symmetry axes necessarily have the same symmetry and so may not be well suited to binding asymmetric molecules. To enable general recognition of arbitrary asymmetric substrates and small molecules, we developed an approach to designing asymmetric interfaces at off-axis sites on homo-oligomers, analogous to those found in native homo-oligomeric proteins such as glutamine synthetase. We symmetrically dock curved helical repeat proteins such that they form pockets at the asymmetric interface of the oligomer with sizes ranging from several angstroms, appropriate for binding a single ion, to up to more than 20 Å across. Of the 133 proteins tested, 84 had soluble expression in E. coli, 47 had correct oligomeric states in solution, 35 had small-angle X-ray scattering (SAXS) data largely consistent with design models, and 8 had negative-stain electron microscopy (nsEM) 2D class averages showing the structures coming together as designed. Both an X-ray crystal structure and a cryogenic electron microscopy (cryoEM) structure are close to the computational design models. The nature of these proteins as homo-oligomers allows them to be readily built into higher-order structures such as nanocages, and the asymmetric pockets of these structures open rich possibilities for small-molecule binder design free from the constraints associated with monomer destabilization.

59 BASIC BIOLOGICAL SCIENCES↗

Massive all-atom analysis of 2D materials with quantum properties (Final report)

Improvements in microscopy have enabled the acquisition of data at a scale that is difficult to process manually, making automated machine learning approaches to analyzing experimental images essential. In this project, we developed and applied machine learning (ML) workflows for atomic resolution scanning transmission electron microscopy (STEM) images. This development included improving both methodology as well as generating user-friendly codes. We developed machine learning architectures which, after training, automatically identify the location and types of defects throughout a material. We used these data to produce class-averaged images of 2D atomic coordinates with up to 0.3 pm precision, uncovering the structure and oscillations of long-range strain fields around point defects in WSe 2-2x Te 2x . We also resolved a long-standing problem in this field in the training of ML models, a lack of labeled experimental data, by developing a cycle-GAN that transformed simulated-generated labeled data into labeled data indistinguishable from experiment and therefore suitable for training. This removed the remaining parts of the ML data processing workflow where human intervention was still critical and therefore a bottleneck to working at scale. Codes have been developed and released for this full machine learning workflow. ML approaches to partially automate STEM acquisition were also developed. Finally we applied ML and other advanced data processing methods to several materials science problems in two-dimensional materials, including studying the evolution of hyperuniformity with defect concentration in WSe2, understanding phase transformations in transition metal dichalcogenides during in-situ heating in the STEM, and exploring how 2D interfaces transform from twisted into aligned structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Deep Learning Enabled Strain Mapping of Single-Atom Defects in Two-Dimensional Transition Metal Dichalcogenides with Sub-Picometer Precision

Two-dimensional (2D) materials offer an ideal platform to study the strain fields induced by individual atomic defects, yet challenges associated with radiation damage have so far limited electron microscopy methods to probe these atomic-scale strain fields. In this work, we demonstrate an approach to probe single-atom defects with sub-picometer precision in a monolayer 2D transition metal dichalcogenide, WSe 2–2x Te 2x . We utilize deep learning to mine large data sets of aberration-corrected scanning transmission electron microscopy images to locate and classify point defects. By combining hundreds of images of nominally identical defects, we generate high signal-to-noise class averages which allow us to measure 2D atomic spacings with up to 0.2 pm precision. Our methods reveal that Se vacancies introduce complex, oscillating strain fields in the WSe 2–2x Te 2x lattice that correspond to alternating rings of lattice expansion and contraction. These results indicate the potential impact of computer vision for the development of high-precision electron microscopy methods for beam-sensitive materials.

2D materials↗

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗

2D dilaton-gravity, deformations of the minimal string, and matrix models

A large class of two-dimensional dilaton-gravity theories in asymptotically AdS2 spacetimes are holographically dual to a matrix integral, interpreted as an ensemble average over Hamiltonians. Viewing these theories as Jackiw–Teitelboim gravity with a gas of defects, we extend this duality to a broader class of dilaton potentials compared to previous work by including conical defects with small deficit angles. In order to do this we show that these theories are equal to the large p limit of a natural deformation of the (2, p) minimal string theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A highly efficient atomically thin curved PdIr bimetallene electrocatalyst

Abstract The multi-metallene with an ultrahigh surface area has great potential in precise tuning of surface heterogeneous d-electronic correlation by surface strain effect for the distinctive surface electronic structure, which is a brand new class of promising 2D electrocatalyst for sustainable energy device application. However, achieving such an atomically thin multi-metallene still presents a great challenge. Herein, we present a new synthetic method for an atomic-level palladium-iridium (PdIr) bimetallene with an average thickness of only ∼1.0 nm for achieving superior catalysis in the hydrogen evolution reaction (HER) and the formic acid oxidation reaction (FAOR). The curved PdIr bimetallene presents a top-ranked high electrochemical active area of 127.5 ± 10.8 m2 gPd+Ir−1 in the reported noble alloy materials, and exhibits a very low overpotential, ultrahigh activity and improved stability for HER and FAOR. DFT calculation reveals that the PdIr bimetallene herein has a unique lattice tangential strain, which can induce surface distortion while concurrently creating a variety of concave-convex featured micro-active regions formed by variously coordinated Pd sites agglomeration. Such a strong strain effect correlates the abnormal on-site active 4d10-t2g-orbital Coulomb correlation potential and directly elevates orbital-electronegativity exposure within these active regions, resulting in a preeminent barrier-free energetic path for significant enhancement of FAOR and HER catalytic performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self‐Lubricating Tribo‐Catalytic Activity of 2D High Entropy Alloy Nanoflakes

High Entropy Alloys (HEAs) have garnered attention due to their remarkable tribological attributes. Predominantly, failure mechanisms in HEAs emanate from stress‐induced dislocations, culminating in crack propagation and film delamination. In this study, we report on the synthesis of 2D HEA of (MoWNbTaV) 0.2 S 2 which facilitates shear‐induced energy dissipation at sliding interfaces. The ball‐on‐disk tribological investigations demonstrate unprecedentedly low average coefficients of friction (0.076) and wear rates (10 −9 mm 3 (N∙m) −1 ) under high contact pressures (0.936 GPa) within ambient conditions. Employing multi‐scale characterizations alongside molecular dynamic simulations, we elucidate that the presence of the HEA triggers tribocatalytic activity under high contact pressures emerging as a pivotal factor in extending lubricant lifespan during tribological tests. The resilient lubriciousness coupled with the facile spray coating methodology of (MoWNbTaV) 0.2 S 2 in ambient environments paves the way for the development of a new class of solid lubricants based on 2D HEA.

2D solid lubricant↗

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗

LAVA Simulations for the 3rd AIAA CFD High Lift Prediction Workshop with Body Fitted Grids

In response to the 3rd AIAA CFD High Lift Prediction Workshop, the workshop cases were analyzed using Reynolds-averaged Navier-Stokes flow solvers within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework. For the workshop cases the advantages and limitations of both overset-structured an unstructured polyhedral meshes were assessed. The workshop included 3 cases: a 2D airfoil validation case, a mesh convergence study using the High Lift Common Research Model, and a nacelle/pylon integration study using the JAXA (Japan Aerospace Exploration Agency) Standard Model. The 2D airfoil case from the workshop is used to verify the implementation of the Spalart-Allmaras turbulence model along with some of its variants within the solver. The High Lift Common Research Model case is used to assess solver performance and accuracy at varying mesh resolutions, as well as identify the minimum mesh fidelity required for LAVA on this class of problem. The JAXA Standard Model case is used to assess the solver's sensitivity to the turbulence model and to compare the structured and unstructured mesh paradigms. These workshop cases have helped establish best practices for high lift flow configurations for the LAVA solver.

LAVA↗

Classical combinatorial optimization scaling for random Ising models on 2D heavy-hex graphs

Motivated by near term quantum computing hardware limitations, combinatorial optimization problems that can be addressed by current quantum algorithms and noisy hardware with little or no overhead are used to probe capabilities of quantum algorithms such as the quantum approximate optimization algorithm. In this study, a specific class of near term quantum computing hardware defined combinatorial optimization problems, Ising models on heavy-hex graphs both with and without geometrically local cubic terms, are examined for their classical computational hardness via empirical computation time scaling quantification. Specifically the time-to-solution (TTS) metric using the classical heuristic simulated annealing is measured for finding optimal variable assignments (ground states), as well as the time required for the optimization software Gurobi to find an optimal variable assignment. Because of the sparsity of these Ising models, the classical algorithms are able to find optimal solutions efficiently even for large instances (i.e. 100 000 spin variables). The Ising models both with and without geometrically local cubic terms exhibit average-case linear-time or weakly quadratic scaling when solved exactly using Gurobi, and the Ising models with no cubic terms show evidence of exponential-time TTS scaling when sampled using simulated annealing. These findings point to the necessity of developing and testing more complex, namely more densely connected, optimization problems in order for quantum computing to ever have a practical advantage over classical computing. Our results are another illustration that different classical algorithms can indeed have exponentially different running times, thus making the identification of the best practical classical technique important in any quantum computing vs. classical computing comparison.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗

AEPF: Attention-Enabled Point Fusion for 3D Object Detection

Current state-of-the-art (SOTA) LiDAR-only detectors perform well for 3D object detection tasks, but point cloud data are typically sparse and lacks semantic information. Detailed semantic information obtained from camera images can be added with existing LiDAR-based detectors to create a robust 3D detection pipeline. With two different data types, a major challenge in developing multi-modal sensor fusion networks is to achieve effective data fusion while managing computational resources. With separate 2D and 3D feature extraction backbones, feature fusion can become more challenging as these modes generate different gradients, leading to gradient conflicts and suboptimal convergence during network optimization. To this end, we propose a 3D object detection method, Attention-Enabled Point Fusion (AEPF). AEPF uses images and voxelized point cloud data as inputs and estimates the 3D bounding boxes of object locations as outputs. An attention mechanism is introduced to an existing feature fusion strategy to improve 3D detection accuracy and two variants are proposed. These two variants, AEPF-Small and AEPF-Large, address different needs. AEPF-Small, with a lightweight attention module and fewer parameters, offers fast inference. AEPF-Large, with a more complex attention module and increased parameters, provides higher accuracy than baseline models. Experimental results on the KITTI validation set show that AEPF-Small maintains SOTA 3D detection accuracy while inferencing at higher speeds. AEPF-Large achieves mean average precision scores of 91.13, 79.06, and 76.15 for the car class’s easy, medium, and hard targets, respectively, in the KITTI validation set. Results from ablation experiments are also presented to support the choice of model architecture.

Chemistry↗

Characterizing Detailed Grain Shape and Size Distribution Properties of Lunar Regolith

Introduction: As the nation prepares to return to the Moon, there is an increasing need for testing tools, instruments, and equipment in simulated environments on Earth to ensure successful operations during lunar missions. Regolith will affect all aspects of future lunar missions, from plume interactions during landing to space suit and tool design [1]. Because of this, it is important to understand the grain shape and size properties of lunar regolith and how those influence regolith behavior in order to prepare for these missions. This knowledge is also vital to create more accurate lunar regolith simulants for testing equipment in a lunar environment. While particle size analyses have been performed on most Apollo soils using simple sieving, shape has only been crudely addressed [2]. This work analyzes 4 lunar regolith samples to provide a better understanding of these size and shape parameters and will provide a more accurate baseline of data to create high-fidelity lunar regolith simulants. Methods: New technologies exist today that are capable of measuring size and shape simultaneously for hundreds of thousands of particles in a single measurement. We conducted a rigorous analysis of the particle size distribution (PSD) as well as the size-dependent 2D and 3D shape parameters of lunar regolith samples of different compositions and maturity levels. This analysis was done using a Microtrac SYNC which provides a unique combination of tri-laser diffraction and Dynamic Image Analysis (DIA). Sample Selection. Four regolith samples were selected for analysis based on maturity level and lunar terrain type: • 10084 – Mature high-Ti mare regolith • 15601 – Immature low-Ti mare regolith • 64501 – Mature highland regolith • 67461 – Immature highland regolith Sample Analysis. After receiving the samples, each sample was imaged with an optical microscope (Figure 1). To obtain 2D and 3D particle size and shape, 0.1 g of each sample was then added to the SYNC for analysis, which outputs more than 30 size and shape parameters for each individual grain as well as the complete size distribution from 0.01-2000 μm by blending laser diffraction and DIA together. For the 0.1 g sample masses requested, ~100,000 grains per sample were captured by DIA. Results: From the PSD analysis, it can be seen that the average particle size of 10084 is ~24.5 µm which is much smaller than the average particle size of the Apollo sample collection (~72 µm) [3], while the average particle size of samples 15601, 64501, and 67461 are larger than the Apollo sample average at ~106 µm, ~103 µm, and ~118.5 µm respectively. The size and shape measurements for the samples output ~30 parameters, four of which were focused on for this study: sphericity, aspect ratio, roundness, and concavity (Figure 3). However, after analyzing the plots it was found that only sphericity and aspect ratio showed differences between the samples. These two parameters are measured on a scale from 0 to 1, with 1 being a perfect sphere with equal dimensions. The results show that the sphericity values are slightly high-er for the mature samples 10084 and 64501 (~0.96) than they are for samples 15601 and 67461 (~0.95) (Figure 4). The aspect ratio values are slightly lower for samples 64501 (~0.7) and 67461 (~0.74) than they are for samples 10084 (~0.76) and 15601 (~0.8) (Figure 5). Discussion: Mature regoliths are those that have been exposed to micrometeorite and solar wind bombardment for long periods of time, breaking up particles and causing them to become more rounded [4]. Therefore, the smaller particle sizes and higher sphericity values for samples 10084 and 64501 are expected due to their higher maturity level compared to 15601 and 67461. However, the aspect ratio values are not dependent on maturity and are instead dependent on terrain type. The lower aspect ratio values for the high-land samples is potentially due to the higher plagioclase content, which occurs in elongated particles and does not break down as easily as the pyroxenes or olivines that are present in the mare. Conclusion and Future Work: The results of this work provide a baseline of high-quality data that will contribute to the creation of high-fidelity lunar simulants and will greatly benefit NASA’s efforts of establishing a human presence on the Moon. Future work includes performing unsupervised image classification on the ~105 particle images per sample in order to identify different classes of grains. These grain classes can then be linked to detailed shape properties, and the relative abundance of each class in the samples can be compared. Acknowledgments: We would like to thank the Extraterrestrial Materials Analysis Group (ExMAG) and the Astromaterials Allocation Review Board (AARB) for allocating lunar regolith samples 10084,9010, 15601,365, 64501,249, and 67461,171 for this work. References: [1] Taylor, L. A. et al. (2005) AIAA #2510. [2] Katagiri et al. (2015) ASCE. [3] Carrier III, W. D. (2005) Lunar Geotech Institute Tech Report. [4] McKay, D. S. et al. (1991) The Lunar Sourcebook, Chapter 7.

S R Deitrick↗