Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase problem”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

A Comprehensive Approach towards Multi-Objective EMI Filter Design Optimization in High-Frequency SiC-Based Motor Drives

With rapid penetration of high-frequency, highdensity power electronics into the industry, controlling the system generated EMI becomes one of the major design challenges. Designing practical optimized EMI filters require simultaneous consideration of multiple aspects - including electromagnetic couplings, magnetic materials, practical component parasitics and their impacts at high-frequencies, and also the overall manufacturing costs. The multi-dimensional nature of the problem often results in overcompensated designs that hurt the system efficiency and power density. This paper focuses on a comprehensive design platform (developed in MATLAB) towards designing high-density, optimized, highly efficient EMI filters by analyzing the generated noise spectrum and evaluating a vast array of potential solutions before outputting an optimized filter solution. The proposed tool includes database-based component selection, optimal topology selection, multistage filter design, overall optimization for volume, mass, cost and total loss. The proposed filter has been thoroughly evaluated in a PSIM simulation environment emulating the hardware in CRD300DA12EXM3, a 300kW three-phase inverter from Wolfspeed with an RL load representing a motor drive application. Furthermore, a 20kHz switching frequency is chosen and impacts of the motor high-frequency characteristics on the generated EMI noise spectrum of the overall system have been analyzed.

42 ENGINEERING↗
DOE OSTI2024

Dartmouth Theory Group: The Physics of the Universe

1. What were the major goals of the project? Under this project, the senior investigator Robert Caldwell (PI) conducted research that addresses outstanding problems in cosmology. This project built on the PI’s ongoing research program in dark energy and cosmic acceleration, cosmology as a probe of fundamental physics, and new gravitational phenomena. Other faculty conducting research under this project were: Stephon Alexander, Marcelo Gleiser, and Devin Walker. Specific results have included predictions of new physics that informs and helps guide Cosmic Frontier observational probes of dark energy, dark matter, and CMB probes pursuing B-mode polarization. 2. What was accomplished under these goals? The project accomplished the primary goals: to extend the state of knowledge of dark energy, inflation, and the phases of the early Universe; to identify new observations and measurements that may help identify the underlying physics of the cosmos; to develop new analysis techniques that enhance the value of exist

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal

Highlights: • We look at the role played by generic scale invariance (GSI) in the phase transition between a U(1) spin liquid and a Fermi liquid. • We identify dangerously irrelevant operators in the renormalization group transformation, involved with the GSI. • We highlight novel mean-field dynamics of this Mott quantum criticality, resulting from the GSI. We investigate the role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal, where there exist massless degrees of freedom in addition to quantum critical fluctuations. Here, the Mott quantum criticality is described by critical charge fluctuations, and additional gapless excitations are U(1) gauge-field fluctuations coupled to a spinon Fermi surface in the spin-liquid state, which turn out to play a central role in the Mott transition. An interesting feature of this problem is that the scaling dimension of effective leading local interactions between critical charge fluctuations differs from that of the coupling constant between U(1) gauge fields and matter-field fluctuations in the presence of a Fermi surface. As a result, there appear dangerously irrelevant operators, which can cause conceptual difficulty in the implementation of renormalization group (RG) transformations. Indeed, we find that the curvature term along the angular direction of the spinon Fermi surface is dangerously irrelevant at this spin-liquid Mott quantum criticality, responsible for divergence of the self-energy correction term in U(1) gauge-field fluctuations. Performing the RG analysis in the one-loop level based on the dimensional regularization method, we reveal that such extremely overdamped dynamics of U(1) gauge-field fluctuations, which originates from the emergent one-dimensional dynamics of spinons, does not cause any renormalization effects to the effective dynamics of both critical charge fluctuations and spinon excitations. However, it turns out that the coupling between U(1) gauge-field fluctuations and both matter-field excitations still persists at this Mott transition, which results in novel mean-field dynamics to explain the nature of the spin-liquid Mott quantum criticality. We discuss physical implications of effective one-dimensional spin dynamics and extremely overdamped gauge dynamics at the Mott quantum criticality.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electronic spectroscopy and ionization potentials for YbOH and YbOCH 3

The polyatomic molecules YbOH and YbOCH 3 have been recognized as being of potential value for spectroscopic experiments that explore charge-parity and time reversal symmetry violation effects. These measurements require very high precision, which, in turn, will necessitate that the molecules be manipulated at ultra-cold temperatures. Both YbOH and YbOCH 3 have electronic transitions that appear suitable for laser cooling ($\tilde{A}$ 2 Π 1/2 -$\tilde{X}$ 2 Σ + and $\tilde{A}$2E 1/2 -$\tilde{X}$ 2 A 1 , respectively) but the currently available spectroscopic data is not sufficient to determine the extent to which population leaks may compromise the optical cooling processes. A further complication is that the quantum states of interest for these measurements will need to be selectively populated. The $\tilde{A}$-$\tilde{X}$ band systems of both YbOH and YbOCH 3 show evidence of vibronic perturbations, such that there are unassigned vibronic features at energies that are just above the origin bands. Here, in the present study we have recorded spectra for the $\tilde{A}$ 2 Π 1/2 -$\tilde{X}$ 2 Σ + transition of jet-cooled YbOD to facilitate the vibronic assignments. In addition, spectra for the $\tilde{B}$ 2 Σ + - $\tilde{X}$ 2 Σ + transition of YbOH were recorded, establishing the origin band at 20473.8 cm -1 . Previously, the reaction of Yb with CH 3 OH has been used to generate gas-phase YbOCH 3 . As this reaction also yields YbOH, there have been complications in spectroscopic studies of YbOCH 3 due to overlap of the $\tilde{A} -\tilde{X}$ band systems. To identify specific regions of overlap, resonantly enhanced two-photon ionization spectra were recorded using mass-resolved detection of the YbOH + and YbOCH 3 + ions. These data confirmed the overlap of vibronic bands near 17640 and 17680 cm -1 . Two-photon ionization spectroscopy also provided accurate ionization energies of IE(YbOH)=45788(10) and IE(YbOCH 3 )=45283(10) cm -1 . The IE for YbOH is relevant to problems encountered in previous attempts to determine the bond dissociation energy of YbOH + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Self-reported health impacts of do-it-yourself air cleaner use in a smoke-impacted community

Smoke exposure from wildfires or residential wood burning for heat is a public health problem for many communities. Do-It-Yourself (DIY) portable air cleaners (PACs) are promoted as affordable alternatives to commercial PACs, but evidence of their effect on health outcomes is limited. Pilot test an evaluation of the effect of DIY PAC usage on self-reported symptoms, and investigate barriers and facilitators of PAC use, among members of a tribal community that routinely experiences elevated concentrations of fine particulate matter (PM 2.5 ) from smoke. We conducted studies in Fall 2021 (“wildfire study”; N = 10) and Winter 2022 (“wood stove study”; N = 17). Each study included four sequential one-to-two-week phases: 1) initial, 2) DIY PAC usage ≥8 h/day, 3) commercial PAC usage ≥8 h/day, and 4) air sensor with visual display and optional PAC use. We continuously monitored PAC usage and indoor/outdoor PM 2.5 concentrations in homes. Concluding each phase, we conducted phone surveys about participants’ symptoms, perceptions, and behaviors. We analyzed symptoms associated with PAC usage and conducted an analysis of indoor PM 2.5 concentrations as a mediating pathway using mixed effects multivariate linear regression. We categorized perceptions related to PACs into barriers and facilitators of use. No association was observed between PAC usage and symptoms, and the mediation analysis did not indicate that small observed trends were attributable to changes in indoor PM 2.5 concentrations. Small sample sizes hindered the ability to draw conclusions regarding the presence or absence of causal associations. DIY PAC usage was low; loud operating noise was a barrier to use. This research is novel in studying health effects of DIY PACs during wildfire and wood smoke exposures. Such research is needed to inform public health guidance. Recommendations for future studies on PAC use during smoke exposure include building flexibility of intervention timing into the study design.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Analysis and Optimization of Seismo-Acoustic Monitoring Networks with Bayesian Optimal Experimental Design

The Bayesian optimal experimental design (OED) problem seeks to identify data, sensor configurations, or experiments which can optimally reduce uncertainty. The goal of OED is to find an experiment that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge about expected data. Therefore, within the context of seismic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types, and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize the ability to locate seismic events from arrival time data of detected seismic phases. In order to do utilize Bayesian OED we must develop four elements:1. A likelihood function that describes the uncertainty of detection and travel times; 2. A Bayesian solver that takes a prior and likelihood to identify the posterior; 3. An algorithm to compute EIG; and, 4. An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we can explore many relevant questions to monitoring such as: how and what multiphenomenology data can be used to optimally reduce uncertainty, how to trade off sensor fidelity and earth model uncertainty, and how sensor types, number, and locations influence uncertainty

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

$t_1$-Noise eliminated dipolar heteronuclear multiple-quantum coherence solid-state NMR spectroscopy

Heteronuclear correlation (HETCOR) spectroscopy is one of the key tools in the arsenal of the solid-state NMR spectroscopist to probe chemical and spatial proximities between two different nuclei and enhance spectral resolution. Dipolar heteronuclear multiple-quantum coherence (D-HMQC) is a powerful technique that can be potentially utilized to obtain 1 H detected 2D HETCOR solid-state NMR spectra of any NMR active nucleus. A long-standing problem in 1 H detected D-HMQC solid-state NMR experiments is the presence of t 1 -noise which reduces sensitivity and impedes spectral interpretation. In this contribution, we describe novel pulse sequences, termed t 1 -noise eliminated (TONE) D-HMQC, that minimize t 1 -noise and can provide higher sensitivity and resolution than conventional D-HMQC. Monte-Carlo and numerical simulations confirm that t 1 -noise in conventional D-HMQC primarily occurs because random MAS frequency fluctuations cause variations in the NMR signal amplitude from scan to scan, leading to imperfect cancellation of uncorrelated signals by phase cycling. The TONE D-HMQC sequence uses 1 H π-pulses to refocus the evolution of 1 H CSA across each SR4 2 1 recoupling block, improving the stability of the pulse sequence to random MAS frequency fluctuations. The 1 H refocusing pulses also restore the orthogonality of in-phase and anti-phase magnetization for all crystallite orientations at the end of each recoupling block, enabling the use of 90° flip-back or LG spin-lock trim pulses to reduce the intensity of uncorrelated signals. We demonstrate the application of these methods to acquire 1 H detected 2D 1 H{ 35 Cl} and 1 H{ 13 C} HETCOR spectra of histidine·HCl·H 2 O with reduced t 1 -noise. To show generality, we also apply these methods to obtain 2D 1 H{ 17 O} spectra of 20%- 17 O fmoc-alanine and for the first time at natural abundance, 2D 1 H{ 25 Mg} HETCOR spectra of magnesium hydroxide. The TONE D-HMQC sequences are also used to probe 1 H– 25 Mg and 1 H– 27 Al proximities in Mg–Al layered double hydroxides and confirm the even mixing of Mg and Al in these materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Reaching the equilibrium state of the frustrated triangular Ising magnet Ca 3 Co 2 O 6

Ca 3 Co 2 O 6 is a frustrated magnet consisting of a triangular arrangement of chains of Ising spins. It shows regular magnetization steps versus magnetic field every 1.2 T that are metastable with very slow dynamics. This has puzzled the community for many years and given rise to numerous potential theories. Here we approach the problem by seeking the elusive magnetic equilibrium state at T=2 K. To this end, we explore two approaches: (1) bypassing the slow dynamics produced by changing fields by instead field-cooling directly to the target temperature and (2) quantum annealing in transverse magnetic fields. While we observe no measurable effect of the quantum annealing in fields up to 7 T, which is likely due to the large Ising anisotropy of Co spins in this material, we find that for the field cooling in longitudinal fields we achieve the predicted equilibrium 1/3 magnetization. We perform Monte Carlo simulations of the ground state phase diagram and we also simulate the quantum annealing process and find good agreement between experiment and theory. Thus we present an investigation of the elusive ground state properties of the canonical frustrated triangular system Ca 3 Co 2 O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2024

Chapter 7: Learning Stable Local Volt/Var Controllers in Distribution Grids

This chapter describes a framework to synthesize provably stable local Volt/Var controllers for distributed energy resources (DERs) in power distribution grids (DGs). The goal is to control the reactive power injections of DERs to improve the system performance as quantified by a generic optimal reactive power flow (ORPF) problem. To achieve this, we jointly design for each DER the control function, which prescribes the reactive power update rule, and the equilibrium function, which approximates the ORPF solutions from local measurements of voltages and powers. We provide conditions on the equilibrium functions and the control parameters ensuring the stability of the closed-loop system. In particular, we discuss the trade-offs between each set of conditions accounting for practical considerations, like fully exploiting the DERs' generation capabilities and reducing the optimality gap. These conditions are then translated into learning constraints on the neural networks' parameters that are enforced in the training phase. We validate our framework with numerical simulations on the IEEE 37-bus network and through a comparison with an optimized version of standard piece wise linear control rules.

closed-loop asymptotic stability↗
DOE OSTI2020

Discrete ordinates analysis of the forced-flight variance reduction technique in Monte Carlo neutral particle transport simulations

This paper presents mathematical formulations and methods to predict the effect of forced-flight variance reduction on Monte Carlo tally variance and calculation time. This includes deducing biasing operators that are then used to construct a history-score probability density function (HSPDF), which represents all possible Monte Carlo random walks and gives the probability of a Monte Carlo history scoring in a tally from a particular phase-space position. The history-score moment equations (HSMEs), the statistical moments of the HSPDF, are then derived to calculate the statistical behavior of the Monte Carlo tally when forced-flight variance reduction is applied. In addition, the future-time equation (FTE) is derived to predict the Monte Carlo computational time as a result of applying forced-flight variance reduction. The solutions of the HSMEs and FTE can be used to predict Monte Carlo computational cost. Furthermore, this work also describes a discrete ordinates method to solve the forced-flight HSMEs and FTE. Several 1-D and 2-D test problems verify that the derivations are performed and implemented correctly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Newton versus the machine: solving the chaotic three-body problem using deep neural networks

ABSTRACT Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system’s chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.

Breen, Philip G.↗
DOE OSTI2021

Iterated projected position algorithm for constructing exponentially localized generalized Wannier functions for periodic and nonperiodic insulators in two dimensions and higher

Localized bases play an important role in understanding electronic structure. In periodic insulators, a natural choice of localized basis is given by the Wannier functions which depend on a choice of unitary transform known as a gauge transformation. Over the past few decades, there have been many works that have focused on optimizing the choice of the gauge so that the corresponding Wannier functions are maximally localized or reflect some symmetry of the underlying system. In this work, we consider fully nonperiodic materials where the usual Wannier functions are not well defined and gauge optimization is impractical. To tackle the problem of calculating exponentially localized generalized Wannier functions in both periodic and nonperiodic systems, we discuss the ‘iterated projected position (IPP)” algorithm. Here, the IPP algorithm is based on matrix diagonalization and therefore unlike optimization-based approaches, it does not require initialization and cannot get stuck at a local minimum. Furthermore, the IPP algorithm is guaranteed by a rigorous analysis to produce exponentially localized functions under certain mild assumptions. We numerically demonstrate that the IPP algorithm can be used to calculate exponentially localized bases for the Haldane model, the Kane-Mele model (in both Z 2 invariant even and Z 2 invariant odd phases), and the p x + ip y model on a quasicrystal lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Cryogenic electron tomography reveals novel structures in the apical complex of Plasmodium falciparum

Intracellular infectious agents, like the malaria parasite, Plasmodium falciparum, face the daunting challenge of how to invade a host cell. This problem may be even harder when the host cell in question is the enucleated red blood cell, which lacks the host machinery co-opted by many pathogens for internalization. Evolution has provided P. falciparum and related single-celled parasites within the phylum Apicomplexa with a collection of organelles at their apical end that mediate invasion. This apical complex includes at least two sets of secretory organelles, micronemes and rhoptries, and several structural features like apical rings and a putative pore through which proteins may be introduced into the host cell during invasion. We perform cryogenic electron tomography (cryo-ET) equipped with Volta Phase Plate on isolated and vitrified merozoites to visualize the apical machinery. Through tomographic reconstruction of cellular compartments, we see new details of known structures like the rhoptry tip interacting directly with a rosette resembling the recently described rhoptry secretory apparatus (RSA), or with an apical vesicle docked beneath the RSA. Subtomogram averaging reveals that the apical rings have a fixed number of repeating units, each of which is similar in overall size and shape to the units in the apical rings of tachyzoites of Toxoplasma gondii. Comparison of these polar rings in Plasmodium and Toxoplasma parasites also reveals them to have a structurally conserved assembly pattern. These results provide new insight into the essential and structurally conserved features of this remarkable machinery used by apicomplexan parasites to invade their respective host cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Fast Butterfly-Compressed Hadamard–Babich Integrator for High-Frequency Helmholtz Equations in Inhomogeneous Media with Arbitrary Sources

Here we present a butterfly-compressed representation of the Hadamard-Babich (HB) ansatz for the Green's function of the high-frequency Helmholtz equation in smooth inhomogeneous media. For a computational domain discretized with Nv discretization cells, the proposed algorithm first solves and tabulates the phase and HB coefficients via eikonal and transport equations with observation points and point sources located at the Chebyshev nodes using a set of much coarser computation grids, and then butterfly compresses the resulting HB interactions from all Nv cell centers to each other. The overall CPU time and memory requirement scale as O(Nv log2 Nv) for any bounded two-dimensional (2D) domains with arbitrary excitation sources. A direct extension of this scheme to bounded 3D domains yields an O(Nv4/3) CPU complexity, which can be further reduced to quasi-linear complexities with proposed remedies. The scheme can also efficiently handle scattering problems involving inclusions in inhomogeneous media. Although the current construction of our HB integrator does not accommodate caustics, the resulting HB integrator itself can be applied to certain sources, such as concave-shaped sources, to produce caustic effects. Compared to finite-difference frequency domain methods, the proposed HB integrator is free of numerical dispersion and requires fewer discretization points per wavelength. As a result, it can solve wave propagation problems well beyond the capability of existing solvers. Remarkably, the proposed scheme can accurately model wave propagation in 2D domains with 640 wavelengths per direction and in 3D domains with 54 wavelengths per direction on a state-of-the-art supercomputer at Lawrence Berkeley National Laboratory.

Hadamard--Babich ansatz↗
DOE OSTI2020

Crystal Structure Prediction of Binary Alloys via Deep Potential

Predicting crystal structure has been a challenging problem in physics and materials science for a long time. A reliable energy calculation engine combined with an efficient global search algorithm, such as particle swarm optimization algorithm or genetic algorithm, is needed to conduct crystal structure prediction. In recent years, machine learning-based interatomic potential energy surface models have been proposed, potentially allowing us to perform crystal structure prediction for systems with the accuracy of density functional theory (DFT) and the speed of empirical force fields. In this paper, we employ a previously developed Deep Potential model to predict the intermetallic compound of the aluminum–magnesium system, and find six meta-stable phases with negative or nearly zero formation energy. In particular, Mg 12 Al 8 shows excellent ductility and Mg 5 Al 27 has a high Young's modulus. Based on our benchmark results, we propose a relatively robust structure screening criterion that selects potentially stable structures from the Deep Potential-based convex hull and performs DFT refinement. By using this criterion, the computational cost needed to construct the convex hull with ab initio accuracy can be dramatically reduced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Measuring Thread Timing to Assess the Feasibility of Early-Bird Message Delivery Across Systems and Scales

Early-bird communication is a communication/computation overlap technique that leverages fine-grained communication to improve application run-time. Communication is divided such that each individual thread can initiate transmission of its portion of the data upon completion rather than waiting for a dedicated communication phase. The benefit of early-bird communication depends on the completion timing of the individual threads: On the one hand, if all threads are complete at nearly the same time, the overheads of sending multiple messages will accumulate, leading to performance that is worse than if a single message had been sent. On the other hand, if thread completions are spread out in time, those that complete earlier can send data while others continue working, leading to performance that is better than if a single message had been sent. The challenge is that the completion times are currently unknown and can vary based on application, problem size, system software, and underlying hardware. In this paper, we address this lacuna by measuring and evaluating the potential overlap afforded by early-bird communication for a selection of proxy applications. These measurements help us understand whether a given application could benefit from early-bird communication. Here, we present our technique for gathering this data and evaluate data collected from three proxy applications: MiniFE, MiniMD, and MiniQMC. Each application is run on three systems with distinct CPU architectures and strong scales across three run sizes. To characterize the behavior of these workloads, we study the trends of thread timings at both a macro level, across all threads across all runs of an application, and a micro level, that is, within a single process of a single run. We observe that our tested applications exhibit significantly different thread arrival distributions. The machine used had a significant impact, with the window of potential overlap varying by as much as an order of magnitude.

97 MATHEMATICS AND COMPUTING↗