Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cardinal”

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 73 records · Page 4

High-dimensional isotomics, part 2: Observations of over 100 constraints on methionine's isotome

The abundances of different isotopic forms of a compound, or isotopologues, will vary based on its physical and chemical history. The number of isotopologues increases combinatorically with the size of a molecule, and even small molecules such as amino acids have thousands of potentially observable isotopic variants. However, due to the analytical challenges of separating and observing isotopologues, only a few dimensions of isotopic diversity are routinely measured. Overcoming these challenges requires both an experimental method to observe many isotopic properties and a theoretical framework for interpreting these experiments. In Part 1, we presented such a theoretical framework; here, we demonstrate an experimental method, which we apply to methionine. Our approach uses a Q Exactive HF Orbitrap to perform several “M + N experiments”, where a sample is ionized, a subset of its isotopologues with cardinal mass N daltons greater than the unsubstituted isotopologue is selected and fragmented, and the proportions of all detectable isotopic forms of those fragment ions are quantified. We perform M + 1, M + 2, M + 3, and M + 4 experiments of a methionine sample and standard where the sample has a 100 ‰ enrichment of 13C at the methyl carbon relative to the natural 13C abundance at that position in the standard, and is otherwise identical to the standard. We observe isotopic forms of 8 fragment ion species for each version of the M + N experiment. With the assistance of a forward model of expected mass spectra, we identify isotopic peaks for each fragment ion based on observed mass and abundance, screen these for data quality, and quantify abundances for 146 unique isotopic peaks at precisions of ≈ 0.3–3 ‰. We present our direct observations and use them to reconstruct the concentrations of 19 individual singly, doubly, and triply-substituted isotopologues; doing so gives fewer constraints and broader error bars than working with the direct observations, but may be more interpretable for some applications. We also examine possibilities for measuring additional peaks, which are primarily limited by the detection limit of the Orbitrap-IRMS method. We then suggest some possible uses of our direct measurements for chemical forensics and hypothesis testing. Furthermore, our results demonstrate the diversity of isotopic constraints currently observable and interpretable for organic molecules.

58 GEOSCIENCES↗

Stochastic economic dispatch of wind power under uncertainty using clustering-based extreme scenarios

Operation of power systems with high penetrations of renewable energy sources requires tools for robust decision making under uncertainty. Stochastic economic dispatch and stochastic unit commitment are effective techniques for planning and operation under uncertainty, whose effectiveness depends on the cardinality and quality of the scenario set. Here, this article proposes a machine learning method using -means clustering for capturing relevant physical information from a large population of analog scenarios. Extreme scenario samples drawn from the clusters are used in a two-stage stochastic economic dispatch computation. The effectiveness of the proposed approach is assessed on a synthetic 200-bus system with a geographic footprint over Illinois, USA for four months from each season of WIND Toolkit data. The combination of -means clustering with importance sampling is shown to reduce the total operational cost by over 43% compared to sampling from populations based on heuristic clustering-based methods. Additionally, the variability in the mean cost is about 56% lower than the variability using Monte Carlo sampling. Moreover, the operational cost with the presented approach is shown to be close to the cost calculated based on a hindsight exact wind profile, signifying a highly accurate quantification of wind uncertainty by the presented -means clustering based sampling method.

17 WIND ENERGY↗

LuGo: An enhanced quantum phase estimation implementation

Quantum Phase Estimation (QPE) is a cardinal algorithm in quantum computing that plays a crucial role in various applications, including cryptography, molecular simulation, and solving systems of linear equations. However, the standard implementation of QPE faces challenges related to time complexity and circuit depth, which limit its practicality for large-scale computations. We introduce LuGo, a novel framework designed to enhance the performance of QPE by reducing circuit duplication, as well as using parallelization techniques to achieve faster generation of the QPE circuit and gate reduction. We validate the effectiveness of our framework by generating quantum linear solver circuits, which require both QPE and inverse QPE, to solve linear systems of equations. LuGo achieves significant improvements in both computational efficiency and hardware requirements without compromising on accuracy. Compared to a standard QPE implementation, LuGo reduces time consumption to generate a circuit that solves a 2 6 × 2 6 system matrix by a factor of 50.68 and over 31× reduction of quantum gates and circuit depth, with no fidelity loss on an ideal quantum simulator. Furthermore, we demonstrated the versatility and scalability of LuGo enabled HHL algorithm by simulating a canonical Hele-Shaw fluid problem using a quantum simulator. With these advantages, LuGo paves the way for more efficient implementations of QPE, enabling broader applications across several quantum computing domains.

Quantum algorithm↗

14C and U-series disequilibria age constraints from recent eruptions at Sete Cidades volcano, Azores

Accurate dating of young (<10 ka) volcanic rocks poses numerous challenges, yet determining eruption ages is critical to understanding magmatic evolution and hazards associated with active volcanic centers. 14C dating is the preferred technique for dating young geologic materials, but a few studies have suggested that volcanic degassing can result in erroneously old 14C dates (e.g. Pasquier-Cardin et al., 1999). Finding an isotopic system with which to accurately constrain eruptive recurrence intervals at recently active volcanic centers is critical to assessing future volcanic hazards. In this study, we have obtained 14C data from modern terrestrial gastropods on Sa~o Miguel island, Azores to assess the current impact of volcanic degassing on 14C ages near active and extinct volcanic centers, and for one active volcanic center, we have compared 14C of paleosols with 226Ra-230Th maximum eruptive ages for pumices immediately overlying the respective paleosols. 14C ages of the modern terrestrial gastropod Cornu aspersum demonstrate the effects of dilution from a range of current outgassing levels at three active trachytic stratovolcanoes on Sa~o Miguel island: Sete Cidades, Fogo, Furnas, as well as the extinct Povoac¸a~o volcano. A modern gastropod sampled at Furnas, where there is evidence of extensive present-day degassing, has a 14C age of 394 ±97 y cal BP. At Fogo volcano, also a location with active degassing, the 14C age of a modern gastropod is 1324 ±60 y cal BP. Both of these ages clearly show effects of dilution of atmospheric 14C due to volcanic degassing. In contrast, there is currently no apparent outgassing occurring at Sete Cidades or Povoac¸a~o, and modern gastropod ages from these volcanoes are post-bomb. We also obtained 14C data from paleosols and 226Ra-230Th ages associated pumices from the youngest eruptive sequence of Sete Cidades volcano, the Pepom P1-P17 deposits, the ages of which are poorly constrained but thought to be <5,000 y (Queiroz, 1997). 226Ra-230Th disequilibrium in a P1 glass separate indicates a maximum eruptive age of 4440 ±195 y. 14C ages of paleosols below P1 and P3 are older than those derived from 226Ra-230Th dating, at 7293 ±141 y cal BP and 7698 ± 172 y cal BP, respectively. We interpret the 14C ages as being erroneously old due to abundant volcanic degassing at the time of eruption. Similarly, but to a lesser extent, the 226Ra-230Th maximum eruptive age of glass separated from the P8 eruption is 2659 ±55 y, whereas the age of the paleosol sampled directly beneathP8 is older (3054 ±194 y BP). In contrast, paleosols sampled at multiple field localities from below the most recent eruptive deposit (P17) are within error of each other (476 ±56 cal y BP and 585 ±68 cal y BP), while 226Ra-230Th maximum eruptive ages in glass separates from two different P17 sampling sites are 762 ±14 y and 730 y ±15 y. This is consistent with the observation that Sete Cidades volcano is experiencing limited hydrothermal activity at the present time and likely had similar low levels of degassing just prior to the last eruption. The average eruptive recurrence interval for the P1-P17 deposits, based on 226Ra-230Th maximum eruptive ages from P1 and P8 and 14C ages from P17, is ~220 years. This recurrence interval is similar to the ~200 y eruptive recurrence interval at Furnas volcano and much shorter than the ~1200 y eruptive interval at Fogo volcano. This has important implications with respect to volcanic hazards and should be considered in any future volcanic hazard planning.

U-series disequilibria, Carbon-14, Cornu aspersum,↗

Variation in microbial community structure correlates with heavy-metal contamination in soils decades after mining ceased

Microorganisms play vital roles in Earth's biogeochemical cycles. Identifying disturbances in microbial communities due to anthropogenic contamination can provide insights into the health of ecosystems. Picher, Oklahoma, was the site of large-scale mining operations for Pb, Zn, and other heavy metals until the mid-1950s, operating within the Tri-State Mining District (TSMD) of Missouri, Kansas and Oklahoma. Although mining ceased decades ago, high concentrations of heavy metals (>1000 ppm) remain in area soil and water systems. Previously, we mapped metal concentrations on samples collected from mine tailings in Picher and along cardinal-direction transects within an 8.05-km radius of the town. To elucidate changes in microbial community structure due to regional metal contamination, 16S rRNA gene sequences and qPCR calculations of total Bacteria and Archaea were analyzed against these metal concentrations. Bacteria were negatively and significantly correlated with Pb, Cd, Zn, and Mg; however, Archaea was only significantly and positively correlated with pH. Illumina sequencing of 16S rRNA genes showed significant differences in microbial communities of chat and west transect samples. Comparison of soil chemistry with community structure indicated that Al, Pb, Cd, and Zn significantly impacted community composition and distribution of individual OTUs. Mapping the distribution of heavy-metal contamination and microbial communities in these soils represents the first step in understanding effects of long-term, heavy-metal contamination at a basic trophic level.

59 BASIC BIOLOGICAL SCIENCES↗

ATP/P2X7-NLRP3 axis of dendritic cells participates in the regulation of airway inflammation and hyper-responsiveness in asthma by mediating HMGB1 expression and secretion

Highlights: • ATP/P2X7 axis induces Th2, 17 cell inflammatory response and AHR in allergic asthma. • Inhibiting NLRP3 blocks Th2, 17 cell inflammatory response and AHR in allergic asthma. • Activating ATP/P2X7 axis increases NLRP3 inflammasome complex expression. • Activating ATP/P2X7 axis induces HMGB1 expresssion and release from DCs. • Inhibiting NLRP3 reduces HMGB1 expresssion and release from DCs. The ATP/P2X7 axis of dendritic cells (DCs) mediates the activation of NLRP3 inflammasome and promotes secretion of interleukin (IL)−1β and IL-18 to induce T helper (Th) 2, Th17 differentiation in the pathogenesis of asthma. NLRP3 inflammasome also regulates high mobility protein 1 (HMGB1) release in DCs. Recent studies demonstrated the correlation between HMGB1 expression and airway inflammation and hyper-responsiveness (AHR) in asthma. However, the relationship between the ATP/P2X7-NLRP3 axis and HMGB1 in DCs in asthma is still unclear. ATP, apyrase, Brilliant Blue G, BzATP, glibenclamide, and Z-YVAD-FMK were administered to ovalbumin (OVA)-induced murine asthmatic model. For in vitro studies, bone marrow-derived mononuclear cells (BMDCs) were primed with LPS and stimulated with the same reagents. Activation of the ATP/P2X7 axis aggravated airway inflammation and AHR in the lung and induced Th2, Th17 polarization in asthmatic mice. Inhibition of NLRP3 inflammasome weakened cardinal features of asthma and blocked Th2, Th17 polarization. In vitro and vivo, ATP/P2X7 axis activated NLRP3 inflammasome and induced HMGB1 expression and release from DCs. Inhibition of NLRP3 inflammasome reduced HMGB1 expression and release. The ATP/P2X7-NLRP3 axis of DCs participates in mediating airway inflammation, AHR, and promoting Th2, Th17 inflammatory responses in asthmatic mice by inducing HMGB1 expression and secretion.

60 APPLIED LIFE SCIENCES↗

Approximation algorithms in combinatorial scientific computing

We survey recent work on approximation algorithms for computing degree-constrained subgraphs in graphs and their applications in combinatorial scientific computing. The problems we consider include maximization versions of cardinality matching, edge-weighted matching, vertex-weighted matching and edge-weighted $b$ -matching, and minimization versions of weighted edge cover and $b$ -edge cover. Exact algorithms for these problems are impractical for massive graphs with several millions of edges. For each problem we discuss theoretical foundations, the design of several linear or near-linear time approximation algorithms, their implementations on serial and parallel computers, and applications. Our focus is on practical algorithms that yield good performance on modern computer architectures with multiple threads and interconnected processors. We also include information about the software available for these problems.

Mathematics↗

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗

Toward a Balanced Description of Ground and Excited States with Transcorrelated F12 Methods

By correlating only the 1-particle states occupied in the reference determinant, the conventional design for the single-reference R12/F12 explicitly correlated methods biases them toward the ground-state description, thereby making the treatment of response properties of the ground state, and energies and other properties of excited states less robust. While the use of multireference methods and/or extensions of the standard SP-projected geminals can achieve a more balanced description of ground and excited states, here we show that the same goals can be achieved by extending the action of F12 correlators to the occupied and valence unoccupied 1-particle states only. This design choice reflects the strong dependence of the optimal correlation length scale of the F12 ansatz on the orbital energies/structure, and helps to avoid the unphysical raising of the ground-state energy if the F12 geminals are used to correlate pairs of all 1-particle states. The improved F12 geminal design is incorporated into the unitary transcorrelation framework to produce a unitary 2-body Hamiltonian that incorporates the short-range dynamical correlation physics for ground and low-energy excited states in a balanced manner. This explicitly correlated effective Hamiltonian reduces the basis set requirement on the correlation-consistent basis cardinal number by 1 or more over the uncorrelated counterpart for the description of ground-state coupled-cluster singles and doubles (CCSD) energies, vertical excitation energies, and harmonic vibrational frequencies of equation-of-motion CCSD low-energy excited states.

Hamiltonians↗

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗

First Analysis of Earth's Stratospheric γ ‐Ray Emission in Geographical Coordinates With Fermi LAT

Abstract From low‐Earth orbit, the Earth is an extremely bright γ ‐ray source compared to distant astrophysical objects. The high‐energy γ ‐ray emission, from the Earth's upper atmosphere, is produced by the interactions between cosmic rays (CRs) and air molecules. As CRs are predominantly charged particles, the geomagnetic field acts as a spectrometer, forbidding CRs below a certain energy to arrive and interact with the atmosphere at a given geomagnetic location from a given direction. We use photon data between 0.2 and 20 GeV as observed by the Fermi Large Area Telescope (LAT) to create, for the first time, stratospheric γ ‐ray emission maps of the Earth in the geographical coordinate system, for different viewing directions from LAT altitude. These maps confirm a strong association between measured γ ‐ray fluxes and geomagnetic cutoffs for CRs. We also report the spectrum of the Earth's stratospheric γ ‐ray emission in cardinal directions and use the geomagnetic spectrometer effect to infer the atmosphere's yield function of γ rays per incoming CR proton as a function of CR rigidity at the altitude of the Fermi spacecraft. This work is relevant to the study of near‐Earth radiation exposure for astronauts and spacecraft, modeling the geomagnetic field, and CR‐air interactions.

Madlee, S.↗

Precision measurement of electron-electron scattering in GaAs/AlGaAs using transverse magnetic focusing

Electron-electron (e-e) interactions assume a cardinal role in solid-state physics. Quantifying the e-e scattering length is hence critical. In this paper we show that the mesoscopic phenomenon of transverse magnetic focusing (TMF) in two-dimensional electron systems forms a precise and sensitive technique to measure this length scale. Conversely we quantitatively demonstrate that e-e scattering is the predominant effect limiting TMF amplitudes in high-mobility materials. Using high-resolution kinetic simulations, we show that the TMF amplitude at a maximum decays exponentially as a function of the e-e scattering length, which leads to a ready approach to extract this length from the measured TMF amplitudes. The approach is applied to measure the temperature-dependent e-e scattering length in high-mobility GaAs/AlGaAs heterostructures. The simulations further reveal current vortices that accompany the cyclotron orbits - a collective phenomenon counterintuitive to the ballistic transport underlying a TMF setting.

36 MATERIALS SCIENCE↗

Dark Energy Survey Year 6 results: Clustering redshifts and importance sampling of self-organized-maps 𝑛⁡(𝑧) realizations for 3 × 2 ⁢pt samples

This work is part of a series establishing the redshift framework for the 3 × 2 ⁢pt analysis of the Dark Energy Survey Year 6 (DES Y6). For DES Y6, photometric redshift distributions are estimated using self-organizing maps (SOMs), calibrated with spectroscopic and many-band photometric data. To overcome limitations from color-redshift degeneracies and incomplete spectroscopic coverage, we enhance this approach by incorporating clustering-based redshift constraints (clustering-z, or WZ) from angular cross-correlations with BOSS and eBOSS galaxies and eBOSS quasar samples. We define a WZ likelihood and apply importance sampling to a large ensemble of SOM-derived 𝑛⁡(𝑧) realizations, selecting those consistent with the clustering measurements to produce a posterior sample for each lens and source bin. The analysis uses angular scales corresponding to 1.5–5 Mpc to optimize signal-to-noise ratio while mitigating modeling uncertainties and marginalizes over redshift-dependent galaxy bias and other systematics informed by the N-body simulation CARDINAL . While a sparser spectroscopic reference sample limits WZ constraining power at 𝑧 >1.1, particularly for source bins, we demonstrate that combining SOM with WZ improves redshift accuracy and enhances the overall cosmological constraining power of DES Y6. As a result, we estimate an improvement in 𝑆 8 of approximately 10% for cosmic shear and 3 ×2⁢pt analysis, primarily due to the WZ calibration of the source samples.

Cosmological parameters↗

Association between optically identified galaxy clusters and the underlying dark matter halos

Clusters of galaxies trace massive dark matter halos in the Universe, but they can include multiple halos projected along lines of sight. Here, we study the halos contributing to clusters using the Cardinal simulation, which mimics the Dark Energy Survey data. We use the red-sequence-based cluster finding algorithm redMaPPer as a case study. For each cluster, we identify the halos hosting its member galaxies, and we define the main halo as the one contributing the most to the cluster's richness ($λ$, the estimated number of member galaxies). At $z=0.3$, for clusters with $λ> 60$, the main halo typically contributes to $92\%$ of the richness, and this fraction drops to $67\%$ for $λ\approx 20$. Defining "clean" clusters as those with $\geq50\%$ of the richness contributed by the main halo, we find that $100\%$ of the $λ> 60$ clusters are clean, while $73\%$ of the $λ\approx 20$ clusters are clean. Three halos can usually account for more than $80\%$ of the richness of a cluster. The main halos associated with redMaPPer clusters have a completeness ranging from $98\%$ at virial mass $10^{14.6}~h^{-1}M_{\odot}$ to $64\%$ at $10^{14}~h^{-1}M_{\odot}$. In addition, we compare the inferred cluster centers with true halo centers, finding that $30\%$ of the clusters are miscentered with a mean offset $40\%$ of the cluster radii, in agreement with recent X-ray studies. These systematics worsen as redshift increases, but we expect that upcoming surveys extending to longer wavelengths will improve the cluster finding at high redshifts. Our results affirm the robustness of the redMaPPer algorithm and provide a framework for benchmarking other cluster-finding strategies.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark energy survey: Modeling strategy for multiprobe cluster cosmology and validation for the full six-year dataset

Here, we introduce an updated To&Krause2021 model for joint analyses of cluster abundances and large-scale two-point correlations of weak lensing and galaxy and cluster clustering (termed CL+3×2 pt analysis) and validate that this model meets the systematic accuracy requirements of analyses with the statistical precision of the final Dark Energy Survey (DES) Year 6 (Y6) dataset. The validation program consists of two distinct approaches, (i) identification of modeling and parametrization choices and impact studies using simulated analyses with each possible model misspecification and (ii) end-to-end validation using mock catalogs from customized Cardinal simulations that incorporate realistic galaxy populations and DES-Y6-specific galaxy and cluster selection and photometric redshift modeling, which are the key observational systematics. In combination, these validation tests indicate that the model presented here meets the accuracy requirements of DES-Y6 for CL+3×2 pt based on a large list of tests for known systematics. In addition, we also validate that the model is sufficient for several other data combinations: the CL+GC subset of this data vector (excluding galaxy–galaxy lensing and cosmic shear two-point statistics) and the CL+3×2 pt+BAO+SN (combination of CL+3×2 pt with the previously published Y6 DES baryonic acoustic oscillation and Y5 supernovae data).

79 ASTRONOMY AND ASTROPHYSICS↗

AMG Preconditioners based on parallel hybrid coarsening and multi-objective graph matching

We describe preliminary results from a multi-objective graph matching algorithm, in the coarsening step of an aggregation-based Algebraic MultiGrid (AMG) preconditioner, for solving large and sparse linear systems of equations on high-end parallel computers. We have two objectives. First, we wish to improve the convergence behavior of the AMG method when applied to highly anisotropic problems. Second, we wish to extend the parallel package \texttt{PSCToolkit} to exploit multi-threaded parallelism at the node level on multi-core processors. Our matching proposal balances the need to simultaneously compute high weights and large cardinalities by a new formulation of the weighted matching problem combining both these objectives using a parameter $\lambda$. We compute the matching by a parallel $2/3-\varepsilon$-approximation algorithm for maximum weight matchings. Results with the new matching algorithm show that for a suitable choice of the parameter $\lambda$ we compute effective preconditioners in the presence of anisotropy, i.e., smaller solve times, setup times, iterations counts, and operator complexity.

D'Ambra, Pasqua↗