SEARCH · Search NASA

Results for “mass accuracy”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Accuracy of environmental tracers and consequences for determining the Type Ia supernova magnitude step

Type Ia Supernovae (SNe Ia) are standardizable candles that allow us to measure the recent expansion rate of the Universe. Due to uncertainties in progenitor physics, potential astrophysical dependencies may bias cosmological measurements if not properly accounted for. The dependency of the intrinsic luminosity of SNe Ia with their host-galaxy environment is often used to standardize SNe Ia luminosity and is commonly parameterized as a step function. This functional form implicitly assumes two-populations of SNe Ia. In the literature, multiple environmental indicators have been considered, finding different, sometimes incompatible, step function amplitudes. We compare these indicators in the context of a two-populations model, based on their ability to distinguish the two populations. We show that local Hα-based specific star formation rate (lsSFR) and global stellar mass are better tracers than, for instance, host galaxy morphology. We show that tracer accuracy can explain the discrepancy between the observed SNe Ia step amplitudes found in the literature. Using lsSFR or global mass to identify the two populations can explain all other observations, though lsSFR is favoured. As lsSFR is strongly connected to age, our results favour a prompt and delayed population model. In any case, there exists two populations that differ in standardized magnitude by at least 0.121 ± 0.010 mag.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: III. From atomic nuclei to neutron stars

Here, we present BSkG3, the latest entry in the Brussels-Skyrme-on-a-grid series of large-scale models of nuclear structure based on an energy density functional. Compared to its predecessors, the new model offers a more realistic description of nucleonic matter at the extreme densities relevant to neutron stars. This achievement is made possible by incorporating a constraint on the infinite nuclear matter properties at high densities in the parameter adjustment, ensuring in this way that the predictions of BSkG3 for the nuclear Equation of State are compatible with the observational evidence for heavy pulsars with M > 2M ⊙ . Instead of the usual phenomenological pairing terms, we also employ a more microscopically founded treatment of nucleon pairing, resulting in extrapolations to high densities that are in line with the predictions of advanced many-body methods and are hence more suited to the study of superfluidity in neutron stars. By adopting an extended form of the Skyrme functional, we are able to reconcile the description of matter at high densities and at saturation density: the new model further refines the description of atomic nuclei offered by its predecessors. A qualitative improvement is our inclusion of ground state reflection asymmetry, in addition to the spontaneous breaking of rotational, axial, and time-reversal symmetry. Quantitatively, the model offers lowered root-mean-square deviations on 2457 masses (0.631 MeV), 810 charge radii (0.0237 fm) and an unmatched accuracy with respect to 45 primary fission barriers of actinide nuclei (0.33 MeV). Reconciling the complexity of neutron stars with those of atomic nuclei establishes BSkG3 as a tool of choice for applications to nuclear structure, the nuclear equation of state and nuclear astrophysics in general.

Physics↗

DOE OSTI2022

Measuring the electron Yukawa coupling via resonant s-channel Higgs production at FCC-ee

The Future Circular Collider (FCC-ee) offers the unique opportunity of studying the Higgs Yukawa coupling to the electron, $y_\mathrm {e}$ , via resonant s -channel production, $\mathrm {e^+e^-}\rightarrow \mathrm {H}$ , in a dedicated run at $\sqrt{s} = m_\mathrm {H}$ . The signature for direct Higgs production is a small rise in the cross sections for particular final states, consistent with Higgs decays, over the expectations for their occurrence due to Standard Model (SM) background processes involving $\mathrm {Z}^*$ , $\gamma ^*$ , or t -channel exchanges alone. Performing such a measurement is remarkably challenging for four main reasons. First, the low value of the e $^\pm $ mass leads to a tiny $y_\mathrm {e}$ coupling and correspondingly small cross section: $\sigma _\mathrm {ee\rightarrow H} \propto m_\mathrm {e}^2 = 0.57$ fb accounting for initial-state $\gamma $ radiation. Second, the $\mathrm {e^+e^-}$ beams must be monochromatized such that the spread of their centre-of-mass (c.m.) energy is commensurate with the narrow width of the SM Higgs boson, $\varGamma _\mathrm {H} = 4.1$ MeV, while keeping large beam luminosities. Third, the Higgs mass must also be known beforehand with a few-MeV accuracy in order to operate the collider at the resonance peak, $\sqrt{s} = m_\mathrm {H}$ . Last but not least, the cross sections of the background processes are many orders-of-magnitude larger than those of the Higgs decay signals. A preliminary generator-level study of 11 Higgs decay channels using a multivariate analysis, which exploits boosted decision trees to discriminate signal and background events, identifies two final states as the most promising ones in terms of statistical significance: $\mathrm {H}\rightarrow gg$ and $\mathrm {H}\rightarrow \mathrm {W}\mathrm {W}^*\!\rightarrow \ell \nu $ + 2 jets. For a benchmark monochromatization with 4.1-MeV c.m. energy spread (leading to $\sigma _\mathrm {ee\rightarrow H} = 0.28$ fb) and 10 ab $^{-1}$ of integrated luminosity, a $1.3\sigma $ signal significance can be reached, corresponding to an upper limit on the e $^\pm $ Yukawa coupling at 1.6 times the SM value: |y e | < 1.6|y e SM | at 95% confidence level, per FCC-ee interaction point per year. Directions for future improvements of the study are outlined.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Skyrme–Hartree–Fock–Bogoliubov mass models on a 3D mesh: IIb. Fission properties of BSkG2

Large-scale models of nuclear structure are currently the only way to provide consistent datasets for the many properties of thousands of exotic nuclei that are required by nucleosynthesis simulations. In [W. Ryssens et al., Eur. Phys. J. A 58, 246 (2022)], we recently presented the new BSkG2 model based on an energy density functional of the Skyrme type. Relying on a flexible three-dimensional coordinate representation of the nucleus, the model takes into account both triaxial deformation and time-reversal symmetry breaking. BSkG2 achieves a state-of-the-art global description of nuclear ground state (g.s.) properties and reproduces in particular the known masses with a root-mean-square (rms) deviation of 678 keV. Moving beyond g.s. properties, the model also reproduces all empirical values for the primary and secondary barriers as well as isomer excitation energies of actinide nuclei with rms deviations below 500 keV, i.e. with unprecedented accuracy. Here we discuss in detail the extension of our framework to the calculation of the fission barriers of 45 actinide nuclei, including odd-mass and odd-odd systems. We focus in particular on the impact of symmetry breaking which is key to the accuracy of the model: we allow systematically for axial, reflection and time-reversal symmetry breaking. The effect of the latter on the fission properties of odd-mass and odd-odd nuclei is small, but we find that allowing for shapes with triaxial or octupole deformation, as well as shapes with both, is crucial to achieving this accuracy. Here, the numerical accuracy of our coordinate space approach, the variety of nuclear configurations explored and the simultaneous successful description of fission properties and known masses makes BSkG2 the tool of choice for the large-scale study of nuclear structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

The ellipticity parameterization for an NFW profile: An overlooked angular structure in strong lens modeling

Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is a Navarro Frenk White (NFW) profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because an unmodeled azimuthal structure has been shown to be able to bias lens model results, we explored the degree to which this azimuthal structure that was introduced can affect the model accuracy. We constructed input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculated the value of the Hubble parameter H 0 one would determine from the lensing fit. We found that the fits to the NFWp input return H 0 values that are systematically biased by about 3% lower than the NFWm counterparts. We explored whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, ξ 2 . We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias was found beyond the MST due to the azimuthal structure exterior to the Einstein radius. We recommend modelers use an NFWm prescription in the future, such that the azimuthal structure can be introduced explicitly rather than implicitly.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Development and demonstration of a rolling ball viscometer for molten salts with near-minimum liquidus NaCl–KCl

The thermophysical characterization of molten chloride and fluoride salts is a key interest for the development of advanced molten salt nuclear reactors. Viscosity is one thermophysical property of interest; a high-accuracy understanding of this property is necessary for modeling the thermal energy transfer and mass flow rate for a given molten salt system. Only a few methods have consistently been used to measure the viscosity of molten chloride and fluoride salt systems, and large discrepancies have been observed in the literature between independent studies of nearly identical salt systems. As such, exploration of novel measurement methods or implementations may allow for further validation, improved measurement accuracy, and overall improved understanding about optimal approaches for measuring molten salt viscosity. Here, in this study, a novel implementation of the rolling ball viscometry method was used to measure the viscosity of a 44/56 mol % mixture of NaCl–KCl within a temperature range of 714–838 °C, with the intent to demonstrate the feasibility of this technique for molten salt viscosity measurement cross-validation. The purity of this salt was confirmed through laser-induced breakdown spectroscopy. Measured viscosities were approximately 1.3–1.0 cP within this temperature range, and the overall uncertainty was 10.8%, which was determined through Gaussian error propagation. The results agree reasonably well with measured viscosities of similar NaCl–KCl mixtures in the literature (within 2%–20% depending on the temperature and study), although the experimental error in some of these comparative studies may be significant.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Improving Dark Energy Constraints Using Low-Redshift Large-Scale Structures

The primary goal of this project was to improve constraints on dark energy measurements by improving our ability to extract cosmological information from low redshift large-scale structures. PI Clowe's project's primary aim was to reduce the bias in measurements of the masses of clusters of galaxies to a level where the evolution of the cluster mass function can be used in the Vera Rubin Observatory's Legacy Survey of Space and Time Dark Energy Science Collaboration survey to improve the accuracy of the measurement of dark energy and other cosmological parameters. Co-PI Seo's project studied observational systematics affecting large-scale clustering of galaxies, which will be used to improve dark energy constraints from the Dark Energy Spectroscopic Instrument (DESI). The cluster lensing project employed a series of simulations and observations of clusters of galaxies to test numerous potential systematic errors in cluster mass measurements using weak gravitational lensing as the accuracy of current weak lensing measurements are more than order of magnitude worse than what is required to use clusters of galaxies for accurate determination of dark energy parameters. PI Clowe and group developed and analyzed simulations to test for, and correct biases introduced in, the weak lensing measurement process. Finally, PI Clowe and group developed a method of detecting clusters using galaxy overdensities and applied the method to the BLISS and DES surveys. The success of spectroscopic dark energy mission such as the extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the Dark Energy Spectroscopic Instrument (DESI) will depend on a thorough understanding of various observational systematics in the target density fluctuations that would give rise to spurious, non-cosmological signals. PI Seo and group developed a deep learning, artificial neural network (ANN) technique that modeled and mitigated such effects, aimed at deriving more robust galaxy clustering signals not only for the baryon acoustic oscillation feature and redshift-space distortions but also for primordial non-Gaussianity constraint.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Efficiently evaluating loop integrals in the EFTofLSS using QFT integrals with massive propagators

We develop a new way to analytically calculate loop integrals in the Effective Field Theory of Large Scale-Structure. Previous available methods show severe limitations beyond the one-loop power spectrum due to analytical challenges and computational and memory costs. Our new method is based on fitting the linear power spectrum with cosmology-independent functions that resemble integer powers of quantum field theory massive propagators with complex masses. A remarkably small number of them is sufficient to reach enough accuracy. Similarly to former approaches, the cosmology dependence is encoded in the coordinate vector of the expansion of the linear power spectrum in our basis. We first produce cosmology-independent tensors where each entry is the loop integral evaluated on a given combination of basis vectors. For each cosmology, the evaluation of a loop integral amounts to contracting this tensor with the coordinate vector of the linear power spectrum. The 3-dimensional loop integrals for our basis functions can be evaluated using techniques familiar to particle physics, such as recursion relations and Feynman parametrization. We apply our formalism to evaluate the one-loop bispectrum of galaxies in redshift space. The final analytical expressions are quite simple and can be evaluated with little computational and memory cost. We show that the same expressions resolve the integration of all one-loop N-point function in the EFTofLSS. This method, which is originally presented here, has already been applied in the first one-loop bispectrum analysis of the BOSS data to constraint ΛCDM parameters and primordial non-Gaussianities [1, 2].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Method To Determine MgO and MgOHCl in Chloride Molten Salts

A commercial blend of mainly carnallite (KCl∙MgCl 2 ∙6H 2 O) is considered as a next-generation heat transfer fluid in solar thermal plants. Corrosive properties of MgCl2 hydrates must be addressed at the operating temperatures of 500-720 degrees C. For successful chemical monitoring of the carnallite heat transfer fluid, an experimental method was developed to separate and titrate for MgO and MgOHCl from solid carnallite. This new method was assessed for error and accuracy. The method's relative error for MgOHCl was -7.0% for a mass fraction of 9.0 wt % MgOHCl in the carnallite salt. The method's relative error for MgO was less than +1.0% for a mass fraction of 12.0 wt % MgO in the carnallite salt. Titration results were used to track changes in the MgOHCl concentration in carnallite salt through the carnallite's dehydration and purification.

14 SOLAR ENERGY↗

DOE OSTI2023

IsoMatchMS : Open-Source Software for Automated Annotation and Visualization of High Resolution MALDI-MS Spectra

Due to its speed, accuracy, and adaptability to various sample types, matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) has become a popular method to identify molecular isotope profiles from biological samples. Often MALDI-MS data do not include tandem MS fragmentation data, and thus the identification of compounds in samples requires external databases so that the accurate mass of detected signals can be matched to known molecular compounds. Most relevant MALDI-MS software tools developed to confirm compound identifications are focused on small molecules (e.g., metabolites, lipids) and cannot be easily adapted to protein data due to their more complex isotopic distributions. Here, we present an R package called IsoMatchMS for the automated annotation of MALDI-MS data for multiple datatypes: intact proteins, peptides, and glycans. This tool accepts already derived molecular formulas or, for proteomics applications, can derive molecular formulas from a list of input peptides or proteins including proteins with post-translational modifications. In conclusion, visualization of all matched isotopic profiles is provided in a highly accessible HTML format called a trelliscope display, which allows users to filter and sort by several parameters such as match scores and the number of peaks matched. IsoMatchMS simplifies the annotation and visualization of MALDI-MS data for downstream analyses.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Global optimization of harmonic oscillator basis in covariant density functional theory

The present investigation focuses on the improvement of the accuracy of the description of binding energies within moderately sized fermionic basis. Using the solutions corresponding to infinite fermionic basis it was shown that in the case of meson exchange (ME) covariant energy density functionals (CEDFs) the global accuracy of the description of binding energies in the finite $N_F$ = 16 - 20 bases can be drastically (by a factor ranging from ~3 up to ~9 dependent on the functional and $N_F$) improved by a global optimization of oscillator frequency of the basis. This is a consequence of the unique feature of the ME functionals in which with increasing fermionic basis size fermionic and mesonic energies approach the exact (infinite basis) solution from above and below, respectively. As a consequence, an optimal oscillator frequency $\hbar\omega_0$ of the basis can be defined which provides an accurate reproduction of exact total binding energies by the ones calculated in truncated basis. This leads to a very high accuracy of the calculations in moderately sized $N_F=20$ basis when mass dependent oscillator frequency is used: global rms differences $\delta B_{rms}$ between the binding energies calculated in infinite and truncated bases are only 0.025 MeV and 0.031 MeV for the NL5(Z) and DD-MEZ functionals, respectively. Optimized values of the oscillator frequency $\hbar\omega_0$ are provided for three major classes of CEDFs, i.e. for density dependent meson exchange functionals, nonlinear meson exchange ones and point coupling functionals.

Binding energy & masses↗

DOE OSTI2020

Detecting axionlike dark matter with linearly polarized pulsar light

QCD axions or axionlike particles are among the most popular candidates for cold dark matter (DM) in the Universe. We proposed to detect axionlike DM, using linearly polarized pulsar light as a probe. Because of birefringence effect potentially caused by an oscillating galactic axion DM background, when pulsar light travels across the Galaxy, its linear polarization angle may vary with time. With a soliton + NFW galactic DM halo profile, we show that this strategy can potentially probe an axion-photon coupling as small as ~10 –13 GeV –1 for axion mass $m_a$ ~ 10 –22 –10 –20 eV , given the current measurement accuracy. An exclusion limit stronger than CAST (~ 10 –10 GeV –1 ) and SN1987A (~ 10 –11 GeV –1 ) could be extended up to $m_a$ ~ 10 –18 eV and ~10 –19 eV, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Robust and High Throughput Tools for Characterization of Wet Cellulosic Biogas Feedstocks and Co-Products: Cooperative Research and Development Final Report, CRADA Number CRD-21-21472

This project will develop standardized analytical procedures to support understanding of the chemical composition of wet manure feedstocks, and the digestate formed following anaerobic digestion of organic waste to produce renewable natural gas (RNG). The developed methods will be validated to improve component mass closure, eliminate ambiguities around key constituents, and document the precision and accuracy of the methods. Following validation of the analytical laboratory methods for characterization of wet organic waste feedstocks and digestate, NLR will use the same samples to build and validate near-infrared (NIR) models for high throughput analysis, potentially deployable at commercial facilities, enabling increased production of low carbon intensity RNG from wet organic wastes.

09 BIOMASS FUELS↗

DOE OSTI2022

Mechanical loading prediction through accelerometry data during walking and running

ABSTRACT Currently, there is no way to assess mechanical loading variables such as peak ground reaction forces (pGRF) and peak loading rate (pLR) in clinical settings. The purpose of this study was to develop accelerometry‐based equations to predict both pGRF and pLR during walking and running. One hundred and thirty one subjects (79 females; 76.9 ± 19.6 kg) walked and ran at different speeds (2–14 km·h −1 ) on a force plate–instrumented treadmill while wearing accelerometers at their ankle, lower back and hip. Regression equations were developed to predict pGRF and pLR from accelerometry data. Leave‐one‐out cross‐validation was used to calculate prediction accuracy and Bland–Altman plots. Our pGRF prediction equation was compared with a reference equation previously published. Body mass and peak acceleration were included for pGRF prediction and body mass and peak acceleration rate for pLR prediction. All pGRF equation coefficients of determination were above 0.96, and a good agreement between actual and predicted pGRF was observed, with a mean absolute percent error (MAPE) below 7.3%. Accuracy indices from our equations were better than previously developed equations. All pLR prediction equations presented a lower accuracy compared to those developed to predict pGRF. Walking and running pGRF can be predicted with high accuracy by accelerometry‐based equations, representing an easy way to determine mechanical loading in free‐living conditions. The pLR prediction equations yielded a somewhat lower prediction accuracy compared with the pGRF equations.

Veras, Lucas↗

DOE OSTI2023

Constraining properties of asymmetric dark matter candidates from gravitational-wave observations

The accumulation of certain types of dark matter particles in neutron star cores due to accretion over long timescales can lead to the formation of a mini black hole. In this scenario, the neutron star is destabilized and implodes to form a black hole without significantly increasing its mass. When this process occurs in neutron stars in coalescing binaries, one or both stars might be converted to a black hole before they merge. Thus, in the mass range of ~1–2M ⊙ , the Universe might contain three distinct populations of compact binaries: one containing only neutron stars, the second population of only black holes, and a third, mixed population consisting of a neutron star and a black hole. However, it is unlikely to have a mixed population as the various timescales allow for both neutron stars to remain or collapse within a short timescale. In this paper, we explore the capability of future gravitational-wave detector networks, including upgrades of Advanced LIGO and Virgo, and new facilities such as the Cosmic Explorer and Einstein Telescope (XG network), to discriminate between different populations by measuring the effective tidal deformability of the binary, which is zero for binary black holes but nonzero for binary neutron stars. Furthermore, we show that observing the relative abundances of the different populations can be used to infer the timescale for neutron stars to implode into black holes, and in turn, provide constraints on the particle nature of dark matter. Furthermore, the XG network will infer the implosion timescale to within an accuracy of 0.01 Gyr at 90% credible interval and determine the dark matter mass and interaction cross section to within a factor of 2 GeV and 10 cm –2 , respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Implementing Ordinary Differential Equation Solvers in Rust Programming Language for Modeling Vehicle Powertrain Systems: Preprint

Efficient and accurate ordinary differential equation (ODE) solvers are necessary for powertrain and vehicle dynamics modeling. However, current commercial ODE solvers can be financially prohibitive, leading to a need for accessible, effective, open-source ODE solvers designed for powertrain modeling. Rust is a compiled programming language that has the potential to be used for fast and easy-to-use powertrain models, given its exceptional computational performance, robust package ecosystem, and short time required for modelers to become proficient. However, of the three commonly used (>3,000 downloads) packages in Rust with ODE solver capabilities, only one has more than four numerical methods implemented, and none are designed specifically for modeling physical systems. Therefore, the goal of the Differential Equation System Solver (DESS) was to implement accurate ODE solvers in Rust designed for the component-based problems often seen in powertrain modeling. DESS is a text-based software package that provides a flexible framework for building and solving systems of ODEs. This allows DESS to be included as a dependency for automotive powertrain models that require a variety of solvers and solver configurations. Seven explicit ODE solver methods have been implemented in DESS: Euler’s, Heun’s, midpoint, Ralston’s, classic Runge-Kutta, Bogacki-Shampine, and Cash-Karp. These represent five fixed-step methods and two adaptive-step methods. This paper shows that the solver implementations increase accuracy and computational efficiency compared to Euler's method when modeling a system of three thermal masses in Rust. DESS also includes features designed for modeling component-based physical systems. Users can define relationships between nodes in their system, which the package then translates into a system of equations, leading to simpler and more intuitive code. In the case of a three-thermal-mass system, the user can specify node thermal properties (e.g., thermal capacitance), how nodes are interconnected, and thermal conductance between nodes rather than providing a system of equations. The core contribution from this work is an open-source, text-based Rust package with ODE solvers for automotive powertrain modeling to support cost-free, fast, and accurate simulation.

ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Copacabana: a probabilistic membership assignment method for galaxy clusters

Cosmological analyses using galaxy clusters in optical/near-infrared photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for all galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89 per cent membership accuracy with a mild dependence on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass–cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities (⁠μ * ⁠) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation at low redshift for a Legacy Survey of Space and Time type survey. These results indicate the potential of Copacabana and μ * to be used in cosmological analyses of optically selected clusters in the future.

79 ASTRONOMY AND ASTROPHYSICS↗