Search NASA⌕ Search

SEARCH · Search NASA

Results for “irregular”

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.

At least 361 records · Page 20

The parallelism motifs of genomic data analysis

Genomic datasets are growing dramatically as the cost of sequencing continues to decline and small sequencing devices become available. Enormous community databases store and share these data with the research community, but some of these genomic data analysis problems require large-scale computational platforms to meet both the memory and computational requirements. These applications differ from scientific simulations that dominate the workload on high-end parallel systems today and place different requirements on programming support, software libraries and parallel architectural design. For example, they involve irregular communication patterns such as asynchronous updates to shared data structures. We consider several problems in high-performance genomics analysis, including alignment, profiling, clustering and assembly for both single genomes and metagenomes. We identify some of the common computational patterns or ‘motifs’ that help inform parallelization strategies and compare our motifs to some of the established lists, arguing that at least two key patterns, sorting and hashing, are missing. This article is part of a discussion meeting issue ‘Numerical algorithms for high-performance computational science’.

59 BASIC BIOLOGICAL SCIENCES↗

Gaussian processes meet NeuralODEs: a Bayesian framework for learning the dynamics of partially observed systems from scarce and noisy data

We present a machine learning framework (GP-NODE) for Bayesian model discovery from partial, noisy and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling and Gaussian Process priors over the observed system states. This allows us to exploit temporal correlations in the observed data, and efficiently infer posterior distributions over plausible models with quantified uncertainty. The use of the Finnish Horseshoe as a sparsity-promoting prior for free model parameters also enables the discovery of parsimonious representations for the latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed GP-NODE method including predator–prey systems, systems biology and a 50-dimensional human motion dynamical system. This article is part of the theme issue ‘Data-driven prediction in dynamical systems’.

Science & Technology - Other Topics↗

Formation of niobium hydride precipitates in superconducting qubits

We report evidence for the formation of niobium hydride phase within niobium films on silicon substrates in superconducting qubits fabricated at Rigetti Computing. For this study, we combined complementary techniques—including room-temperature and cryogenic atomic force microscopy (AFM), synchrotron x-ray diffraction, and time-of-flight secondary ion mass spectroscopy (ToF-SIMS)—to directly reveal the existence of niobium hydride precipitates on the surface of superconducting qubits. Upon cryogenic cooling, we observed variation in the size and morphology of the hydrides, ranging from small (∼5 nm) irregular shapes to large (∼10–100 nm) domains within the Nb grains, which were fully converted to niobium hydrides. Since niobium hydrides are nonsuperconducting and can easily change in size and location upon different cooldowns to cryogenic temperature, our finding highlights a previously unknown source of decoherence in superconducting qubits. This contributes to quasiparticle losses, offering a potential explanation for changes in qubit performance upon cooldowns. Finally, by leveraging the RF performance of a 3D bulk Nb resonator, we quantify RF dissipation in a superconducting qubit caused by hydrogen concentration variation, and propose a practical engineering pathway to mitigate the formation of Nb hydrides for superconducting qubit applications.

Sung, Zuhawn [Fermi National Accelerator Laborator↗

Nonequilibrium effects and transverse spherocity in ultrarelativistic proton-nucleus collisions

We investigate the effects of nonequilibrium dynamics in small colliding systems by comparing a nonequilibrium transport approach, Parton-Hadron-String Dynamics (PHSD), with a (2+1)-dimensional viscous hydrodynamic model, VISHNew. Focusing on p + Pb collisions at energies available at the Large Hadron Collider, we extract the initial conditions for the hydrodynamic model from PHSD, in order to reduce the impact of the early out-of-equilibrium dynamics and focus on the traces of nonequilibrium in the ensuing medium evolution. We find that in the transport approach quantities like energy density and bulk viscous pressure are highly inhomogeneous on the transverse plane during the whole evolution, whereas the hydrodynamic simulations dissolve more efficiently the initial spatial irregularities, still keeping a high degree of inhomogeneity due to the smaller size and lifetime of the medium produced in p + Pb collisions with respect to heavy-ion reactions. As a first step that will help to identify the impact of these nonequilibrium effects on final observables in proton-nucleus collisions, we perform an analysis of the transverse spherocity, an event-shape observable able to distinguish between jetty and isotropic configurations of transverse momenta. We found that the spherocity distribution in PHSD is slightly shifted towards the isotropic limit with respect to the hydrodynamic result. Even though this dissimilarity is partially due to the difference in the final charged particle production, it mainly comes from the different description within the two frameworks of the medium produced in small colliding systems. This finding supports the idea that multidifferential measurements, such as those based on event categorization according to multiplicity and spherocity, are useful to study final-state observables in ultrarelativistic proton-nucleus collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for beyond-mean-field signatures in heavy-ion fusion reactions

Examination of high-resolution, experimental fusion excitation functions for 16,17,18 O + 12 C reveals a remarkable irregular behavior that is rooted in the structure of both the colliding nuclei and the quasi-molecular composite system. The impact of the l-dependent fusion barriers is assessed using a time-dependent Hartree-Fock model, viewed as a mean-field reference. Barrier penetrabilities, taken directly from a density-constrained calculation, provide a significantly improved description of the experimental data as compared to the standard Hill-Wheeler approach. Furthermore, the remaining deviations between the parameter-free theoretical mean-field predictions and experimental fusion cross sections are exposed and discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-precision mass measurement of 103 Sn restores smoothness of the mass surface

As a step towards the ultimate goal of a high-precision mass measurement of doubly magic 100 Sn, the mass of 103 Sn was measured at the Low Energy Beam and Ion Trap (LEBIT) located at the Facility for Rare Isotope Beams (FRIB). Utilizing the time-of-flight ion cyclotron resonance technique, a mass uncertainty of 3.7 keV was achieved, an improvement by more than an order of magnitude compared to a recent measurement performed in 2023 at the Cooler Storage Ring (CSRe) in Lanzhou. Although the LEBIT and CSRe mass measurements of 103 Sn are in agreement, they diverge from the experimental mass value reported in the 2016 version of the Atomic Mass Evaluation (AME2016), which was derived from the measured 𝑄 𝛽 + value and the mass of 103 In. In AME2020, this indirectly measured 103 Sn mass was classified as a “seriously irregular mass” and replaced with an extrapolated value, which aligns with the most recent measured values from CSRe and LEBIT. As such, the smoothness of the mass surface is confidently reestablished for 103 Sn. Here, LEBIT's mass measurement of 103 Sn enabled a significant reduction in the mass uncertainties of five parent isotopes which are now dominated by uncertainties in their respective 𝑄 values.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Denoising diffusion models with geometry adaptation for high fidelity calorimeter simulation

Simulation is crucial for all aspects of collider data analysis, but the available computing budget in the High Luminosity LHC era will be severely constrained. Generative machine learning models may act as surrogates to replace physics-based full simulation of particle detectors, and diffusion models have recently emerged as the state of the art for other generative tasks. Here, we introduce CaloDiffusion, a denoising diffusion model trained on the public CaloChallenge datasets to generate calorimeter showers. Our algorithm employs 3D cylindrical convolutions, which take advantage of symmetries of the underlying data representation. To handle irregular detector geometries, we augment the diffusion model with a new geometry latent mapping (GLaM) layer to learn forward and reverse transformations to a regular geometry that is suitable for cylindrical convolutions. The showers generated by our approach are nearly indistinguishable from the full simulation, as measured by several different metrics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Radio measurements of the depth of air-shower maximum at the Pierre Auger Observatory

The (AERA), part of the Pierre Auger Observatory, is currently the largest array of radio antenna stations deployed for the detection of cosmic rays, spanning an area of 17 km 2 with 153 radio stations. It detects the radio emission of extensive air showers produced by cosmic rays in the 30–80 MHz band. Here, we report the AERA measurements of the ( X max ), a probe for mass composition, at cosmic-ray energies between 10 17.5 and 10 18.8 eV , which show agreement with earlier measurements with the fluorescence technique at the Pierre Auger Observatory. We show advancements in the method for radio X max reconstruction by comparison to dedicated sets of / air-shower simulations, including steps of reconstruction-bias identification and correction, which is of particular importance for irregular or sparse radio arrays. Using the largest set of radio air-shower measurements to date, we show the radio X max resolution as a function of energy, reaching a resolution better than 15 g cm − 2 at the highest energies, demonstrating that radio X max measurements are competitive with the established high-precision fluorescence technique. In addition, we developed a procedure for performing an extensive data-driven study of systematic uncertainties, including the effects of acceptance bias, reconstruction bias, and the investigation of possible residual biases. These results have been cross-checked with air showers measured independently with both the radio and fluorescence techniques, a setup unique to the Pierre Auger Observatory. Published by the American Physical Society 2024

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Patterning of multicomponent elastic shells by gaussian curvature

Recent findings suggest that shell protein distribution and the morphology of bacterial microcompartments regulate the chemical fluxes facilitating reactions which dictate their biological function. Here, we explore how the morphology and component patterning are coupled through the competition of mean and gaussian bending energies in multicomponent elastic shells that form three-component irregular polyhedra. We observe two softer components with lower bending rigidities allocated on the edges and vertices while the harder component occupies the faces. When subjected to a nonzero interfacial line tension, the two softer components further separate and pattern into subdomains that are mediated by the gaussian curvature. We find that this degree of fractionation is maximized when there is a weaker line tension and when the ratio of bending rigidities between the two softer domains ≈2. Our results reveal a patterning mechanism in multicomponent shells that can capture the observed morphologies of bacterial microcompartments, and moreover, can be realized in synthetic vesicles.

Monte Carlo methods↗

Shearless effective barriers to chaotic transport induced by even twin islands in nontwist systems

For several decades now it has been known that systems with shearless invariant tori, nontwist Hamiltonian systems, possess barriers to chaotic transport. These barriers are resilient to breakage under perturbation and therefore regions where they occur are natural places to look for barriers to transport. Here we describe a kind of effective barrier that persists after the shearless torus is broken. Because phenomena are generic, for convenience we study the standard nontwist map (SNM), an area-preserving map that violates the twist condition locally in the phase space. The barrier occurs in nontwist systems when twin even period islands are present, which happens for a broad range of parameter values in the SNM. With a phase space composed of regular and irregular orbits, the movement of chaotic trajectories is hampered by the existence of shearless curves, total barriers, and a network of partial barriers formed by the stable and unstable manifolds of the hyperbolic points. Being a degenerate system, the SNM has twin islands and, consequently, twin hyperbolic points. We show that the structures formed by the manifolds intrinsically depend on period parity of the twin islands. For this even scenario the structure that we call a torus free barrier occurs because the manifolds of different hyperbolic points form an intricate chain atop a dipole configuration and the transport of chaotic trajectories through the chain becomes a rare event. This structure impacts the emergence of transport, the escape basin for chaotic trajectories, the transport mechanism, and the chaotic saddle. The case of odd periodic orbits is different: we find for this case the emergence of transport immediately after the breakup of the last invariant curve, and this leads to a scenario of higher transport, with intricate escape basin boundary and a chaotic saddle with nonuniformly distributed points.

classical mechanics↗

Defect Spirograph: Dynamical Behavior of Defects in Spatially Patterned Active Nematics

We report that topological defects in active liquid crystals can be confined by introducing gradients of activity. Here, we examine the dynamical behavior of two defects confined by a sharp gradient of activity that separates an active circular region and a surrounding passive nematic material. Continuum simulations are used to explain how the interplay among energy injection into the system, hydrodynamic interactions, and frictional forces governs the dynamics of topologically required self-propelling + 1/2 defects. Our findings are rationalized in terms of a phase diagram for the dynamical response of defects in terms of activity and frictional damping strength. Different regions of the underlying phase diagram correspond to distinct dynamical modes, namely immobile defects, steady rotation of defects, bouncing defects, bouncing-cruising defects, dancing defects, and multiple defects with irregular dynamics. These dynamic states raise the prospect of generating synchronized defect arrays for microfluidic applications.

36 MATERIALS SCIENCE↗

Discovery of Charge Order and Corresponding Edge State in Kagome Magnet FeGe

Kagome materials often host exotic quantum phases, including spin liquids, Chern gap, charge density wave, and superconductivity. Existing scanning microscopy studies of the kagome charge order have been limited to nonkagome surface layers. Here, we tunnel into the kagome lattice of FeGe to uncover features of the charge order. Our spectroscopic imaging identifies a 2×2 charge order in the magnetic kagome lattice, resembling that discovered in kagome superconductors. Spin mapping across steps of unit cell height demonstrates the existence of spin-polarized electrons with an antiferromagnetic stacking order. We further uncover the correlation between antiferromagnetism and charge order anisotropy, highlighting the unusual magnetic coupling of the charge order. Finally, we detect a pronounced edge state within the charge order energy gap, which is robust against the irregular shape fluctuations of the kagome lattice edges. Here, we discuss our results with the theoretically considered topological features of the kagome charge order including unconventional magnetism and bulk-boundary correspondence.

36 MATERIALS SCIENCE↗

Machine Learning to Select Experiments Driven by Fundamental Science and Applications for Targeted Nuclear Data Improvement

This work describes a blueprint for a process that accelerates progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application-driven experiments to maximally reduce pertinent data uncertainties? Answering this question entails solving a high-dimensional and complex optimization problem that is best solved with advanced statistic techniques often classified as machine learning. We apply this process within the framework of nuclear data with the aim to select an experiment combination that will reduce uncertainties in 239 Pu nuclear data for neutron energies between 1 and 600 keV. In this field, fundamental-physics driven data, called differential, look at one nuclear physics observable at a time. They are contrasted to application-driven, integral, data where one or few resulting values inform a broad set of nuclear data across several nuclides and energies. The candidates for integral experiments are criticality measurements that were refined by a genetic algorithm to be maximally sensitive to 239 Pu fission cross sections in the desired energy range. Twenty-three candidate differential experiments were investigated and span multiple nuclear physics observables (e.g., total, capture cross sections) for isotopes appearing in the integral experiments. The optimal combination among these candidate experiments was investigated via generalized least squares fitting, augmented with Gaussian processes to ameliorate statistical irregularities in data, and the D-optimality criterion. The latter evaluates for each pair of candidates the joint reduction in uncertainties of all 12200 nuclear data appearing in the integral experiments compared to the knowledge we have from 168 past experiments, theory, and nuclear data. We chose as differential measurements those that investigate 63 Cu and 239 Pu total cross sections, based on D-optimality rank and feasibility constraints. Two integral (criticality) experiments were selected: An experiment with Al 2 ⁢O 3 and graphite interleaved with Pu and a thick Cu reflector explores 1–30 keV, while we target the 30–600 keV range with an experiment that swaps boron in place of graphite with a different geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

FMX – the Frontier Microfocusing Macromolecular Crystallography Beamline at the National Synchrotron Light Source II

Two new macromolecular crystallography (MX) beamlines at the National Synchrotron Light Source II, FMX and AMX, opened for general user operation in February 2017 [Schneider et al. (2013). J. Phys. Conf. Ser. 425 , 012003; Fuchs et al. (2014). J. Phys. Conf. Ser. 493 , 012021; Fuchs et al. (2016). AIP Conf. Proc. SRI2015 , 1741 , 030006]. FMX, the micro-focusing Frontier MX beamline in sector 17-ID-2 at NSLS-II, covers a 5–30 keV photon energy range and delivers a flux of 4.0 × 10 12 photons s −1 at 1 Å into a 1 µm × 1.5 µm to 10 µm × 10 µm (V × H) variable focus, expected to reach 5 × 10 12 photons s −1 at final storage-ring current. This flux density surpasses most MX beamlines by nearly two orders of magnitude. The high brightness and microbeam capability of FMX are focused on solving difficult crystallographic challenges. The beamline's flexible design supports a wide range of structure determination methods – serial crystallography on micrometre-sized crystals, raster optimization of diffraction from inhomogeneous crystals, high-resolution data collection from large-unit-cell crystals, room-temperature data collection for crystals that are difficult to freeze and for studying conformational dynamics, and fully automated data collection for sample-screening and ligand-binding studies. FMX's high dose rate reduces data collection times for applications like serial crystallography to minutes rather than hours. With associated sample lifetimes as short as a few milliseconds, new rapid sample-delivery methods have been implemented, such as an ultra-high-speed high-precision piezo scanner goniometer [Gao et al. (2018). J. Synchrotron Rad. 25 , 1362–1370], new microcrystal-optimized micromesh well sample holders [Guo et al. (2018). IUCrJ , 5 , 238–246] and highly viscous media injectors [Weierstall et al. (2014). Nat. Commun. 5 , 3309]. The new beamline pushes the frontier of synchrotron crystallography and enables users to determine structures from difficult-to-crystallize targets like membrane proteins, using previously intractable crystals of a few micrometres in size, and to obtain quality structures from irregular larger crystals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deconvoluting thermomechanical effects in X-ray diffraction data using machine learning

X-ray diffraction is ideal for probing the sub-surface state during complex or rapid thermomechanical loading of crystalline materials. However, challenges arise as the size of diffraction volumes increases due to spatial broadening and because of the inability to deconvolute the effects of different lattice deformation mechanisms. Here, we present a novel approach that uses combinations of physics-based modeling and machine learning to deconvolve thermal and mechanical elastic strains for diffraction data analysis. The method builds on a previous effort to extract thermal strain distribution information from diffraction data. The new approach is applied to extract the evolution of the thermomechanical state during laser melting of an Inconel 625 wall specimen which produces significant residual stress upon cooling. A combination of heat transfer and fluid flow, elasto-plasticity and X-ray diffraction simulations is used to generate training data for machine-learning (Gaussian process regression, GPR) models that map diffracted intensity distributions to underlying thermomechanical strain fields. First-principles density functional theory is used to determine accurate temperature-dependent thermal expansion and elastic stiffness used for elasto-plasticity modeling. The trained GPR models are found to be capable of deconvoluting the effects of thermal and mechanical strains, in addition to providing information about underlying strain distributions, even from complex diffraction patterns with irregularly shaped peaks.

36 MATERIALS SCIENCE↗

Effective Missing Value Imputation Methods for Building Monitoring Data

To understand behaviors of natural and man-made events, such as energy consumption of buildings, which accounts for 40% of energy uses in the US, we deploy automated monitoring devices to record periodic observations. However, such experimental and observation data often contains problems and irregularities that have to be cleaned up before analyses. Due to various conditions affecting sensor operations, the communication channels, recording steps, or the recording media, the recorded data might have missing values, errors, or anomalous values. An effective way to clean up these problems is to replace these missing values, errors and anomalous values with expected values, a process generally known as imputation. In this work, we survey commonly used missing value imputation techniques and compare their performance on a set of building monitoring data. To compare the different types of sensor measurements with widely varying characteristics, we use normalized root mean squared error (NRMSE) as the key metric for the effectiveness of the imputation methods. We additionally consider periodicity and run time when considering comparing methods. Through extensive testing, we find that for small gap sizes, up to 8 consecutive missing values, linear interpolation performs the best; for larger gaps stretching up to 48 consecutive missing values, K-nearest neighbors provides the most accurate imputations; for even larger gaps, more computational intensive methods, such as matrix factorization, achieve the smallest NRMSE. Additionally, we observe that these computationally intensive algorithms not only provide accurate imputations for large gaps, but are also more robust across all types of sensors.

Cho, B↗

MassiveGNN: Efficient Training via Prefetching for Massively Connected Distributed Graphs

Graph Neural Networks (GNN) are indispensable in learning from graph-structured data, yet their rising computational costs, especially on massively connected graphs, pose significant challenges in terms of execution performance. To tackle this, distributed-memory solutions such as partitioning the graph to concurrently train multiple replicas of GNNs are in practice. However, approaches requiring a partitioned graph usually suffer from communication overhead and load imbalance, even under optimal partitioning and communication strategies due to irregularities in the neighborhood minibatch sampling. This paper proposes practical trade-offs for improving the sampling and communication overheads for representation learn- ing on distributed graphs (using popular GraphSAGE architecture) by developing a parameterized prefetch and eviction scheme on top of the state-of-the-art Amazon DistDGL distributed GNN framework, demonstrating about 15–40% improvement in end-to-end training performance on the NERSC Perlmutter supercomputer for various OGB datasets.

Machine Leanring, high performance comptuing, grap↗

Signal Processing Based Method for Real-Time Anomaly Detection in High-Performance Computing

Performance anomalies can manifest as irregular execution times or abnormal execution events for many reasons, including network congestion and resource contention. Detecting such anomalies in real-time by analyzing the details of performance traces at scale is impractical due to the sheer volume of data High-Performance Computing (HPC) applications produce. In this paper, we propose formulating HPC performance anomaly detection as a signal-processing problem where anomalies can be treated as noise. We evaluate our proposed method in comparison with two other commonly used anomaly detection techniques of varying complexity based on their detection accuracy and scalability. Since real-time in-situ anomaly detection at a large scale requires lightweight methods that can handle a large volume of streaming data, we find that our proposed method provides the best trade-off. We then implement the proposed method in Chimbuko, the first online, distributed, and scalable workflow-level performance trace analysis framework. We compare our proposed signal-based anomaly detection algorithm with two other methods using a function of their accuracy, F1 score, and detection overhead. Our experiments demonstrate that our proposed approach achieves a 99% improvement for the benchmark datasets and a 93% improvement with Chimbuko traces.

99 GENERAL AND MISCELLANEOUS↗