SEARCH · Engineering Papers

Results for “galaxies: general”

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 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Star forming regions in gas-rich SO galaxies

The first results of an H alpha imaging survey of HI rich SO galaxies, which were searched for HII regions and other sources of emission, are presented. The charge coupled device H alpha interference filter images were made of 16 galaxies. Eight of these galaxies show evidence for on-going star formation, one has nuclear emission but no HII regions, and the remaining seven have no emissions detected within well defined upper limits. With the exception of one notably peculiar galaxy in which the emission from HII regions appears pervasive, the HII regions are either organized into inner-disk rings or randomly distributed throughout the disk. A few of these galaxies are found to be clearly not SO's; or peculiar objects atypical of the SO class. Using simple models star formation rates (SFRs) and gas depletion times from the observed H alpha fluxes were estimated. In general, the derived SFRs are much lower than those found in isolated field spiral galaxies and the corresponding gas depletion time scales are also longer.

Pogge, Richard W.↗

DOE OSTI2026

The DESI DR1 peculiar velocity survey: Growth rate measurements from the galaxy power spectrum

The large-scale structure of the Universe and its evolution encapsulate a wealth of cosmological information. A powerful means of unlocking this knowledge lies in measuring the auto-power spectrum and/or the cross-power spectrum of the galaxy density and momentum fields, followed by the estimation of cosmological parameters based on these spectrum measurements. In this study, we generalize the cross-power spectrum model to accommodate scenarios in which the density and momentum fields are derived from distinct galaxy surveys. The growth rate of the large-scale structures of the Universe, commonly represented as fσ 8 , was extracted by jointly fitting the monopole and quadrupole moments of the auto-density power spectrum, the monopole of the auto-momentum power spectrum, and the dipole of the cross-power spectrum. Our estimators, theoretical models, and parameter-fitting framework were tested using mocks, confirming their robustness and accuracy in retrieving the fiducial growth rate from simulation. These techniques were then applied to analyse the power spectrum of the DESI Bright Galaxy Survey and Peculiar Velocity Survey. The fit result of the growth rate is fσ8 = 0.440$^{+0.080}_{−0.096}$ at effective redshift zeff = 0.07. By synthesizing the fitting outcomes from correlation functions, maximum likelihood estimation, and the power spectrum, a consensus value is yielded of fσ 8 (z eff = 0.07) = 0.450$^{+0.055}_{−0.055}$, and correspondingly we obtain γ = 0.580$^{+0.110}_{−0.110}$, Ω m = 0.301$^{+0.011}_{−0.011}$, and σ 8 = 0.834$^{+0.032}_{−0.032}$. The measured fσ8 and γ are consistent with the prediction of the Λ cold dark matter model and general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

Imaging the host galaxies of high-redshift radio-quiet QSOs

We present new deep K-band and optical images of four radio-quiet QSOs at z approximately = 1 and six radio-quiet QSOs at z approximately = 2.5, as well as optical images only of six more at z approximately = 2.5. We have examined the images carefully for evidence of extended 'fuzz' from any putative QSO host galaxy. None of the z approximately = 2.5 QSOs shows any extended emission, and only two of the z approximately = 1 QSOs show marginal evidence for extended emission. Our 3 sigma detection limits in the K images, m(sub K) approximately = 21 for an isolated source, would correspond approximately to an unevolved L(sup star) elliptical galaxy at z = 2.5 or 2-3 mag fainter than an L(sup star) elliptical at z = 1, although our limits on host galaxy light are weaker than this due to the difficulty of separating galaxy light from QSO light. We simulate simple models of disk and elliptical host galaxies, and find that the marginal emission around the two z approximately = 1 QSOs can be explained by disks or bulges that are approximately 1-2 mag brighter than an unevolved L(sup star) galaxy in one case and approximately 1.5-2.5 mag brighter than L(sub star) in the other. For two other z approximately = 1 QSOs, we have only upper limits (L approximately = L(sup star)). The hosts of the high-redshift sample must be no brighter than about 3 mag above an unevolved L(sup star) galaxy, and are at least 1 magnitude fainter than the hosts of radio-loud QSOs at the same redshift. If the easily detected K-band light surrounding a previous sample of otherwise similar but radio-loud QSOs is starlight, then it must evolve on timescales of greater than or approximately equal to 10(exp 8) yr (e.g., Chambers & Charlot 1990); therefore our non-detection of host galaxy fuzz around radio-quiet QSOs supports the view that high-redshift radio-quiet and radio-loud QSOs inhabit different host objects, rather than being single types of objects that turn their radio emission on and off over short timescales. This is consistent with the general trend at low redshifts that radio-loud QSOs are found in giant elliptical galaxies while radio-quiet QSOs are found in less luminous disk galaxies. It also suggests that the processes responsible for the spectacular properties of radio-loud AGNs at high redshifts might not be generally relevent to the (far more numerous) radio-quiet population.

Lowenthal, James D.↗

DOE OSTI2022

Deblurring galaxy images with Tikhonov regularization on magnitude domain

Abstract We propose a regularization-based deblurring method that works efficiently for galaxy images. The spatial resolution of a ground-based telescope is generally limited by seeing conditions and is much worse than space-based telescopes. This circumstance has generated considerable research interest in the restoration of spatial resolution. Since image deblurring is a typical inverse problem and often ill-posed, solutions tend to be unstable. To obtain a stable solution, much research has adopted regularization-based methods for image deblurring, but the regularization term is not necessarily appropriate for galaxy images. Although galaxies have an exponential or Sérsic profile, the conventional regularization assumes the image profiles to behave linearly in space. The significant deviation between the assumption and real situations leads to blurring of the images and smoothing out the detailed structures. Clearly, regularization on logarithmic domain, i.e., magnitude domain, should provide a more appropriate assumption, which we explore in this study. We formulate a problem of deblurring galaxy images by an objective function with a Tikhonov regularization term on a magnitude domain. We introduce an iterative algorithm minimizing the objective function with a primal–dual splitting method. We investigate the feasibility of the proposed method using simulation and observation images. In the simulation, we blur galaxy images with a realistic point spread function and add both Gaussian and Poisson noise. For the evaluation with the observed images, we use galaxy images taken by the Subaru HSC-SSP. Both of these evaluations show that our method successfully recovers the spatial resolution of the deblurred images and significantly outperforms the conventional methods. The code is publicly available from the GitHub 〈https://github.com/kzmurata-astro/PSFdeconv_amag〉.

Murata, Kazumi↗

NASA NTRS1994

Testing the gravitational instability hypothesis?

We challenge a widely accepted assumption of observational cosmology: that successful reconstruction of observed galaxy density fields from measured galaxy velocity fields (or vice versa), using the methods of gravitational instability theory, implies that the observed large-scale structures and large-scale flows were produced by the action of gravity. This assumption is false, in that there exist nongravitational theories that pass the reconstruction tests and gravitational theories with certain forms of biased galaxy formation that fail them. Gravitational instability theory predicts specific correlations between large-scale velocity and mass density fields, but the same correlations arise in any model where (a) structures in the galaxy distribution grow from homogeneous initial conditions in a way that satisfies the continuity equation, and (b) the present-day velocity field is irrotational and proportional to the time-averaged velocity field. We demonstrate these assertions using analytical arguments and N-body simulations. If large-scale structure is formed by gravitational instability, then the ratio of the galaxy density contrast to the divergence of the velocity field yields an estimate of the density parameter Omega (or, more generally, an estimate of beta identically equal to Omega(exp 0.6)/b, where b is an assumed constant of proportionality between galaxy and mass density fluctuations. In nongravitational scenarios, the values of Omega or beta estimated in this way may fail to represent the true cosmological values. However, even if nongravitational forces initiate and shape the growth of structure, gravitationally induced accelerations can dominate the velocity field at late times, long after the action of any nongravitational impulses. The estimated beta approaches the true value in such cases, and in our numerical simulations the estimated beta values are reasonably accurate for both gravitational and nongravitational models. Reconstruction tests that show correlations between galaxy density and velocity fields can rule out some physically interesting models of large-scale structure. In particular, successful reconstructions constrain the nature of any bias between the galaxy and mass distributions, since processes that modulate the efficiency of galaxy formation on large scales in a way that violates the continuity equation also produce a mismatch between the observed galaxy density and the density inferred from the peculiar velocity field. We obtain successful reconstructions for a gravitational model with peaks biasing, but we also show examples of gravitational and nongravitational models that fail reconstruction tests because of more complicated modulations of galaxy formation.

Babul, Arif↗

NASA NTRS1990

Evidence against a simple two-component model for the far-infrared emission from galaxies

Two of the first Infrared Astronomy Satellite (IRAS) results were that galaxies have a wide range of values for the ratio of 60 micron to 100 micron flux density (0.2 less than or equal to S sub 60/S sub 100 less than or equal to 1.0) and that this ratio is correlated with L sub fir, L sub b, L sub fir being the total far-infrared luminosity and L sub b being the luminosity at visible wavelengths (de Jong et al. 1984; Soifer et al. 1984). From these results arose the following simple model for the far-infrared emission from galaxies (de Jong et al. 1984), which has remained the standard model ever since. In this model, the far-infrared emission comes from two dust components: warm dust (T approx. equals 50 K) intermingled with, and heated by, young massive OB stars in molecular clouds and HII regions, and colder dust (T approx. equals 20 K) associated with the diffuse atomic hydrogen in the interstellar medium and heated by the general interstellar radiation field. As the number of young stars in a galaxy increases, S sub 60/S sub 100 increases, because there is a greater proportion of warm dust, and so does L sub fir/L sub b, because most of the radiation from the young stars is absorbed by the dust, leading to a swifter increase in far-infrared emission than in visible light. Although this model explains the basic IRAS results, it is inelegant. It uses two free parameters to fit two data (the 60 and 100 micron flux densities)-and there are now several observations that contradict it. Despite these major problems with the two-component model, it is not clear what should be put in its place. When considering possible models for the far-infrared emission from galaxies, the observational evidence for our own galaxy must be considered. Researchers suspect that the study by Boulanger and Perault (1988) of the far-infrared properties of the local interstellar medium may be particularly relevant. They showed that molecular clouds are leaky - that most of the light from OB stars in molecular clouds does not heat the dust in the clouds, but instead leaks out. The consequence of this is that that while most of the far-infrared emission from the solar neighborhood is from dust associated with diffuse HI, this dust is mostly heated by young stars.

Eales, Stephen A.↗

NASA NTRS1977

Cosmic-ray propagation and containment

The cosmic rays are an active gaseous component of the disk of the galaxy, and their propagation and containment is a part of the general dynamics of the disk. The sources of cosmic rays are a matter of speculation. The disk is inflated by the cosmic-ray gas pressure, P, comparable to the magnetic pressure B super 2/ 8 pi, but the rate of inflation is unknown. The time spent by the individual cosmic-ray particles in the disk is inversely proportional to the cosmic-ray production rate and may be anything from 100,000 to more than 10 million years. It is evident from the decay of Be(10) that the cosmic rays circulate through a volume of space perhaps ten times the thickness of the gaseous disk, suggesting a magnetic halo extending out approximately 1 kpc from either face of the disk. The cosmic rays may be responsible for the halo by inflating the magnetic fields of the disk. Extension of the fields to 1 kpc would imply a high production rate and short life of cosmic rays in the dense gaseous disk of the galaxy.

Parker, E. N.↗

DOE OSTI2025

Bridging Scales in Black Hole Accretion and Feedback: Relativistic Jet Linking the Horizon to the Host Galaxy

Simulating black hole (BH) accretion and feedback from the BH horizon to galactic scales is extremely challenging, as it involves a vast range of scales. Recently, our multizone method has successfully achieved global dynamical steady states of hot accretion flows in 3D general relativistic magnetohydrodynamic simulations by tracking the bidirectional interaction between a nonspinning BH and its host galaxy. In this paper, we present technical improvements to the method and apply it to spin a * = 0.9 BHs, which power relativistic jets. We first test the new multizone setup with a smaller Bondi radius, R B ≈ 400 r g , where r g is the gravitational radius. The strongly magnetized accretion launches a relativistic jet with an intermediate feedback efficiency η ∼ 30%, in between that of a prograde (η ∼ 100%) and retrograde (η ∼ 10%) torus. Interestingly, both prograde and retrograde simulations also eventually converge to the same intermediate efficiency when evolved long enough, as accumulated magnetic fields remove gas rotation. We then extend strongly magnetized simulations to larger Bondi radii, R B ≈ 2 × 10 3 , 2 × 10 4 , 2 × 10 5 r g . We find that the BH accretion rate $\dot{M}$ is suppressed with respect to the Bondi rate as $\dot{M}_{\textrm{B}}$ as $\dot{M}/\dot{M}_{\textrm{B}} ∝ R_{\textrm{B}}^{-1/2}$. However, despite some variability, the time-averaged feedback efficiency remains at η ∼ 30%, independent of R B . This suggests that BH feedback efficiency in hot accretion flows is mainly governed by the BH spin (a * ) rather than by the galactic properties (R B ). From these first-principles simulations, we provide a feedback subgrid prescription for cosmological simulations: $\dot{E}_{\textrm{fb}} = 2$ x $10^{-3}[R_{\textrm{B}}/(2$ x $10^5 r_g)]^{-1/2}$ $\dot{M}_{\textrm{B}}c^2$ for BH spin a * = 0.9.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1990

NGC 3312: A victim of ram pressure sweeping

Researchers are undertaking a volume limited survey of the Hydra I cluster in neutral hydrogen using the National Radio Astronomy Observatory's Very Large Array (VLA). The main purpose is to study the effects of a dense environment on the gaseous component of the galaxies. Observational evidence has been accumulating recently that ram pressure sweeping does occur in the centers of clusters, but it is possible that tidal interactions play a role as well. Results of high resolution HI imaging of NGC 3312, the large peculiar spiral near the cluster center are presented. Hydra I (= A1060) is the nearest rich cluster beyond Virgo and, as such, presents a unique opportunity to do a complete survey of a cluster. It is similar to the Virgo cluster in many of its general physical characteristics, such as size, x ray luminosity, velocity dispersion, and galaxy content (high spiral fraction). However, Hydra I appears to be more regular and relaxed. This is evident in the x ray distribution in its central region, which is radially symmetric and centered on the dominant galaxy, NGC 3311, a cD-like elliptical. The observed x ray luminosity implies a central gas density of 4.5 x 10 to the 3rd power cm(-3). Gallagher (1978) argued from optical images of NGC 3312 that this galaxy might be an ideal candidate to directly study effects of the ram pressure process; it might currently be undergoing stripping of its interstellar medium. The researchers' data are consistent with this suggestion, but other origins of the peculiar appearance cannot yet be ruled out.

Mcmahon, P. M.↗

NASA NTRS1994

The motions of clusters and group of galaxies

The distributions of peculiar velocities of rich clusters and of groups of galaxies are investigated for different cosmological models and are compared with observations. Four cosmological models are studied: standard cold dark matter (CDM) (omega = 1); low-density CDM (omega = 0.3); hot dark matter (HDM) (omega = 1); and primeval baryonic isocurvature (PBI) (omega = 0.3). All models are normalized to the microwave background fluctuations observed by Cosmic Background Explorer (COBE). We find that rich clusters of galaxies exhibit a Maxwellian distribution of peculiar velocities in all models, as expected from a Gaussian initial density field. The clusters appear to be fundamental and efficient tracers of the large-scale velocity field. The cluster three-dimensional velocity distribution typically peaks at v approximately 600 km/s and extends to high cluster velocities of v approximately 2000 km/s. The low-density CDM model exhibits somewhat lower velocities: it peaks at approximately 400 km/s and extends to approximately 1200 km/s. Approximately 10% (approximately 1% for low-density CDM) of all model rich clusters move with high peculiar velocities of V greater than or = 10(exp 3) km/s. The highest velocity clusters frequently originate in dense superclusters. The model velocity distributions of rich clusters are compared with the model velocity distributions of small groups of galaxies, and of the total matter. The group velocity distribution is, in general, similar to the velocity distribution of the rich clusters. The matter velocity distribution is similar to that of the rich clusters for the omega = 0.3 models; these models exhibit Maxwellian velocity distributions for clusters, for groups, and for matter that are all similar to one another. The mass distribution in the omega = 1 models, however, exhibits a longer tail of high velocities than do the clusters. This high-velocity tail originates mostly from the high velocities that exist within rich clusters. The model velocity distributions of groups and clusters of galaxies are compared with observations. The data exhibit a larger high-velocity tail, to v(sub r) greater than or = 2000 km/s, than is seen in the model simulations (except for HDM). Because of the large observational uncertainties, however, the data are consistent at a approximately 1-3 sigma level with the model predictions and with a Gaussian initial density field. Accurate observations of cluster peculiar velocities, especially at the high-velocity tail, shoud provide powerful constraints on the cosmological models.

Bahcall, Neta A.↗

NASA NTRS1991

Far-infrared properties of NGC 6822

IRAS FIR maps of NGC 6822, one of the few noninteracting Magellanic irregular galaxies resolved in the FIR by IRAS, are examined and compared to optical and H I data. The diffuse FIR emission at large radii from the center of the galaxy is found to be due to Galactic cirrus. There is a general correspondence between H-alpha and FIR emission in that FIR peaks are associated with major H II complexes, but the ratio of the FIR to the H-alpha emission varies between these complexes by about a factor of five. H I mass-to-FIR luminosity ratios cover a similar range for the three major H I cloud complexes in NGC 6822. More than about half of the FIR flux is produced in regions near H II region complexes. The remaining flux comes from diffuse FIR emission that extends over the entire optically visible galaxy. A model for the galaxy is discussed in which the galaxy is fairly transparent and the FIR is strongest where there is a strong UV source and higher optical depth due to dust in gas clouds.

Gallagher, John S., III↗

NASA NTRS1994

Ultraviolet imaging of the AGN+starburst galaxy NGC 1068

Images of the Seyfert 2 galaxy NGC 1068 were obtained at two ultraviolet wavelengths by the Ultraviolet Imaging Telescope (UIT). These data represent the first detailed UV imagery of a composite (active galactic nucleus + starburst) disk galaxy. NGC 1068 cotains multiple components at UV wavelengths: the central active galactic nucleus; a population of very luminous starburst knots; a bright oval inner disk; and a fainter, more circular halo. The most luminous knot, which is located approximately 750 pc from the nucleus at PA 315 deg, is approximately 80 times the luminosity of 30 Doradus and gives NGC 1068 a 'double nucleus' appearance in the UV. Significant extended emission is observed throughout the disk, unlike other disk galaxies so far observed in the UV. The radial brightness profile in both UV bandpasses generally follows an exponential decline to approximately 5 kpc. A faint halo extending to approximately 13 kpc is likely to be a galaxian-sized reflection nebula where ambient dust scatters the intense UV continuum from the inner galaxy. UV colors show a striking asymmetric morphology, which is correlated with the observed molecular CO emission.

Neff, Susan G.↗

DOE OSTI2026

Dwarf Galaxy Number Counts within 25 Mpc: Predictions from Local Group Analogs in TNG50

The modern generation of wide-field galaxy surveys, such as LSST, Euclid, and Roman, will enable studies of dwarf galaxies (10 6 ≤ M * /M ⊙ ≤ 10 9 ) beyond the Local Group (LG) in unprecedented detail. Improved theoretical understanding of this population is necessary to guide these observations, since predictions in this regime are generally limited to specific environments like the LG. We present predictions for the population of dwarf galaxies from the TNG50 run of the IllustrisTNG suite of cosmological hydrodynamical simulations, focusing on the environments within 1 < D/Mpc < 25 of LG analogs at z = 0. In the simulated sample, there are ∼1000 and ∼12,000 dwarf galaxies within 10 and 25 Mpc, respectively. We compare our results with the 50 Mpc Galaxy Catalog and estimate that current observations are highly incomplete at low masses: for 10 6 ≤ M * /M ⊙ ≤ 10 7 (−13 ≲ M r ≲ −10), we find completeness fractions of ∼23% within 10 Mpc and ∼4% within 25 Mpc. The simulated galaxies below the completeness limits of the observations exist in a range of environments, with notable populations of field dwarfs at all distances and satellites around centrals with masses 10 8 ≲ M * /M ⊙ ≲ 10 11 within 10–25 Mpc. We find that there are ∼8 times more quiescent dwarf galaxies in the TNG50 sample than are currently cataloged. Our results suggest that upcoming observations should uncover a substantial population of dwarf galaxies, and that ≳15% of these will be red, currently quenched galaxies in the field.

Shread, Evangela E. [California Institute of Techn↗

NASA NTRS1979

A cosmological redshift-distance square law

This paper directly examines the claims of Segal (1976) that the (m,z) Hubble diagram is fitted best by a square law z = Kr-squared rather than by the traditional Hubble law z = Hr in the low-redshift range, z no more than about 0.01, corresponding to galaxies brighter than 14th mag. Segal attempts to fit a distance relation to the (m,z) scatter diagram in which each individual galaxy is plotted. The exact relation between the mean redshift for all galaxies in a small magnitude interval and the apparent magnitude is calculated. This relation is independent of luminosity function and peculiar velocity distribution about the general expansion, and is not affected by sample incompleteness as a function of apparent magnitude or the clustering of galaxies in the sample. Segal's method is affected by all of these and requires a highly sophisticated statistical analysis to deal with the non-Gaussian pointwise scatter. The present analysis favors the Hubble law and conclusively rules out the square law for the small redshift region.

Soneira, R. M.↗

NASA NTRS1978

X-ray and gamma-ray line production by nonthermal ions

X-ray production at about 6.8 keV by the 2p to 1s transition in fast hydrogen- and helium-like iron ions is calculated, following both electron capture to excited levels and collisional excitation. A refinement of the Oppenheimer-Brinkman-Kramers (1930) approximation is used to obtain an improved charge-exchange cross-section. This, and the corresponding ionization cross section, were used to determine equilibrium charge fractions for iron ions as functions of their energy. The effective X-ray line production cross section was found to be sharply peaked in energy at about 8 to 12 MeV per amu. Since fast ions of similar energies can also excite nuclear levels, the ratio of selected strong gamma-rays line emissivities to the X-ray line emissivity is calculated. These calculations are employed to set limits on the intensity of gamma-rays line emission from the galactic center and the radio galaxy Cen A, and it is found that these limits are generally lower than those reported in the literature.

Bussard, R. W.↗

NASA NTRS1995

Beryllium and Boron abundances in population II stars

The scientific focus of this program was to undertake UV spectroscopic abundance analyses of extremely metal poor stars with attention to determining abundances of light elements such as beryllium and boron. The abundances are likely to reflect primordial abundances within the early galaxy and help to constrain models for early galactic nucleosynthesis. The general metal abundances of these stars are also important for understanding stellar evolution.

Source record↗

NASA NTRS1999

Interpreting Broad Double-Peaked Emission Lines in Active Galactic Nuclei

The principal objectives of this project were to probe the inner regions of active galactic nuclei and to test general relativity in the strong-field limit. The approach takes advantage of broad atomic line emission observed from material deep in the potential well of an active galactic nucleus which contains key information as to the physics of the system. Line profiles in a wide range of wavebands from optical to X-ray have provided compelling evidence of the existence of a relativistic accretion disk around a supermassive black hole in a number of galaxies. The simplest model posits a geometrically thin disk in Keplerian orbit, with general relativistic effects in evidence. This model is the point of departure for the proposed work. We developed a high-performance numerical code to calculate photon trajectories in a Schwarzschild or Kerr metric and implemented it on parallel supercomputers. This code includes a general purpose ray tracer that calculates line profiles, light curves, and other observable quantities for a wide variety of emitter configurations. The versatility comes from the fact that the ray tracing algorithm does not depend on any symmetries regarding emitter locations. The speed comes from parallel implementation which enables us to sample hitherto unattainable volumes of disk model parameter space. During the period 1 March 1997 through 28 February 1998, two papers, supported in whole or in part by this grant, were published in refereed journals. They are reproduced in their entirety in the next two sections of this report.

Halpern, Jules↗

DOE OSTI2025

Reducing model error using optimized galaxy selection: weak lensing cluster mass estimation

Galaxy clusters are one of the most powerful probes to study extensions of General Relativity and the Standard Cosmological Model. Upcoming surveys like the Vera Rubin Observatory’s Legacy Survey of Space and Time are expected to revolutionise the field, by enabling the analysis of cluster samples of unprecedented size and quality. To reach this era of high-precision cluster cosmology, the mitigation of sources of systematic error is crucial. A particularly important challenge is bias in cluster mass measurements induced by the mismodelling of photometric redshift estimates of source galaxies. This work proposes a method to optimise the source sample selection in cluster weak lensing analyses drawn from wide-field survey lensing catalogs to reduce the bias on reconstructed cluster masses. We use a combinatorial optimisation scheme and methods from variational inference to select galaxies in latent space to produce a probabilistic galaxy source sample catalog for highly accurate cluster mass estimation. We show that our method reduces the critical surface mass density Σ crit relative modelling bias on the 60-70% level, while maintaining up to 90% of galaxies. We highlight that our methodology has applications beyond cluster mass estimation as an approach to jointly combine galaxy selection and model inference under sources of systematics.

79 ASTRONOMY AND ASTROPHYSICS↗