Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sample”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A pure-sampling quantum Monte Carlo algorithm

The objective of pure-sampling quantum Monte Carlo is to calculate physical properties that are independent of the importance sampling function being employed in the calculation, save for the mismatch of its nodal hypersurface with that of the exact wave function. To achieve this objective, we report a pure-sampling algorithm that combines features of forward walking methods of pure-sampling and reptation quantum Monte Carlo (RQMC). The new algorithm accurately samples properties from the mixed and pure distributions simultaneously in runs performed at a single set of time-steps, over which extrapolation to zero time-step is performed. In a detailed comparison, we found RQMC to be less efficient. It requires different sets of time-steps to accurately determine the energy and other properties, such as the dipole moment. We implement our algorithm by systematically increasing an algorithmic parameter until the properties converge to statistically equivalent values. As a proof in principle, we calculated the fixed-node energy, static α polarizability, and other one-electron expectation values for the ground-states of LiH and water molecules. These quantities are free from importance sampling bias, population control bias, time-step bias, extrapolation-model bias, and the finite-field approximation. We found excellent agreement with the accepted values for the energy and a variety of other properties for those systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Enhanced configurational sampling with hybrid non-equilibrium molecular dynamics–Monte Carlo propagator

Molecular dynamics (MD) trajectories based on classical equations of motion can be used to sample the configurational space of complex molecular systems. However, brute-force MD often converges slowly due to the ruggedness of the underlying potential energy surface. Several schemes have been proposed to address this problem by effectively smoothing the potential energy surface. However, in order to recover the proper Boltzmann equilibrium probability distribution, these approaches must then rely on statistical reweighting techniques or generate the simulations within a Hamiltonian tempering replica-exchange scheme. The present work puts forth a novel hybrid sampling propagator combining Metropolis-Hastings Monte Carlo (MC) with proposed moves generated by non-equilibrium MD (neMD). This hybrid neMD-MC propagator comprises three elementary elements: (i) an atomic system is dynamically propagated for some period of time using standard equilibrium MD on the correct potential energy surface; (ii) the system is then propagated for a brief period of time during what is referred to as a “boosting phase,” via a time-dependent Hamiltonian that is evolved toward the perturbed potential energy surface and then back to the correct potential energy surface; (iii) the resulting configuration at the end of the neMD trajectory is then accepted or rejected according to a Metropolis criterion before returning to step 1. A symmetric two-end momentum reversal prescription is used at the end of the neMD trajectories to guarantee that the hybrid neMD-MC sampling propagator obeys microscopic detailed balance and rigorously yields the equilibrium Boltzmann distribution. The hybrid neMD-MC sampling propagator is designed and implemented to enhance the sampling by relying on the accelerated MD and solute tempering schemes. It is also combined with the adaptive biased force sampling algorithm to examine. Illustrative tests with specific biomolecular systems indicate that the method can yield a significant speedup.

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

Thermal deflection in neutron scattering sample environments at Oak Ridge National Laboratory

The neutron sources at Oak Ridge National Laboratory use a wide suite of sample environment equipment to deliver extreme conditions for a number of experiments. Much of this instrumentation focuses on extremes of temperature, such as cryostats, closed-cycle refrigerators in both low and high temperature configurations, and radiant heating furnaces. When the temperature is controlled across a large range, thermal deflection effects can notably move the sample and affect its alignment in the beam. Here, we combine these sample environments with neutron imaging and machine vision to determine the motion of a representative sample with respect to the neutron beam. We find vertical sample displacement on the order of 1–2 mm and horizontal displacement that varies from near-negligible to 1.2 mm. While these deflections are not relevant for some of the beamlines at the first target station at the spallation neutron source and the high flux isotope reactor, they will become critical for upcoming instrumentation at the second target station, as well as any instruments targeting sub-mm samples, as neutron sources and optics evolve to smaller and more focused beams. We discuss mitigation protocols and potential modifications to the environment to minimize the effect of misalignment due to thermal deflection.

47 OTHER INSTRUMENTATION↗
DOE OSTI2011

THE zCOSMOS-SINFONI PROJECT. I. SAMPLE SELECTION AND NATURAL-SEEING OBSERVATIONS

The zCOSMOS-SINFONI project is aimed at studying the physical and kinematical properties of a sample of massive z {approx} 1.4-2.5 star-forming galaxies, through SINFONI near-infrared integral field spectroscopy (IFS), combined with the multiwavelength information from the zCOSMOS (COSMOS) survey. The project is based on one hour of natural-seeing observations per target, and adaptive optics (AO) follow-up for a major part of the sample, which includes 30 galaxies selected from the zCOSMOS/VIMOS spectroscopic survey. This first paper presents the sample selection, and the global physical characterization of the target galaxies from multicolor photometry, i.e., star formation rate (SFR), stellar mass, age, etc. The H{alpha} integrated properties, such as, flux, velocity dispersion, and size, are derived from the natural-seeing observations, while the follow-up AO observations will be presented in the next paper of this series. Our sample appears to be well representative of star-forming galaxies at z {approx} 2, covering a wide range in mass and SFR. The H{alpha} integrated properties of the 25 H{alpha} detected galaxies are similar to those of other IFS samples at the same redshifts. Good agreement is found among the SFRs derived from H{alpha} luminosity and other diagnostic methods, provided the extinction affecting the H{alpha} luminosity is about twice that affecting the continuum. A preliminary kinematic analysis, based on the maximum observed velocity difference across the source and on the integrated velocity dispersion, indicates that the sample splits nearly 50-50 into rotation-dominated and velocity-dispersion-dominated galaxies, in good agreement with previous surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

REDSHIFTS, SAMPLE PURITY, AND BCG POSITIONS FOR THE GALAXY CLUSTER CATALOG FROM THE FIRST 720 SQUARE DEGREES OF THE SOUTH POLE TELESCOPE SURVEY

We present the results of the ground- and space-based optical and near-infrared (NIR) follow-up of 224 galaxy cluster candidates detected with the Sunyaev-Zel'dovich (SZ) effect in the 720 deg{sup 2} of the South Pole Telescope (SPT) survey completed in the 2008 and 2009 observing seasons. We use the optical/NIR data to establish whether each candidate is associated with an overdensity of galaxies and to estimate the cluster redshift. Most photometric redshifts are derived through a combination of three different cluster redshift estimators using red-sequence galaxies, resulting in an accuracy of {Delta}z/(1 + z) = 0.017, determined through comparison with a subsample of 57 clusters for which we have spectroscopic redshifts. We successfully measure redshifts for 158 systems and present redshift lower limits for the remaining candidates. The redshift distribution of the confirmed clusters extends to z = 1.35 with a median of z{sub med} = 0.57. Approximately 18% of the sample with measured redshifts lies at z > 0.8. We estimate a lower limit to the purity of this SPT SZ-selected sample by assuming that all unconfirmed clusters are noise fluctuations in the SPT data. We show that the cumulative purity at detection significance {xi} > 5({xi} > 4.5) is {>=}95% ({>=}70%). We present the red brightest cluster galaxy (rBCG) positions for the sample and examine the offsets between the SPT candidate position and the rBCG. The radial distribution of offsets is similar to that seen in X-ray-selected cluster samples, providing no evidence that SZ-selected cluster samples include a different fraction of recent mergers from X-ray-selected cluster samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE ATACAMA COSMOLOGY TELESCOPE: DYNAMICAL MASSES AND SCALING RELATIONS FOR A SAMPLE OF MASSIVE SUNYAEV-ZEL'DOVICH EFFECT SELECTED GALAXY CLUSTERS {sup ,}

We present the first dynamical mass estimates and scaling relations for a sample of Sunyaev-Zel'dovich effect (SZE) selected galaxy clusters. The sample consists of 16 massive clusters detected with the Atacama Cosmology Telescope (ACT) over a 455 deg{sup 2} area of the southern sky. Deep multi-object spectroscopic observations were taken to secure intermediate-resolution (R {approx} 700-800) spectra and redshifts for Almost-Equal-To 60 member galaxies on average per cluster. The dynamical masses M{sub 200c} of the clusters have been calculated using simulation-based scaling relations between velocity dispersion and mass. The sample has a median redshift z = 0.50 and a median mass M{sub 200c}{approx_equal}12 Multiplication-Sign 10{sup 14} h{sub 70}{sup -1} M{sub sun} with a lower limit M{sub 200c}{approx_equal}6 Multiplication-Sign 10{sup 14} h{sub 70}{sup -1} M{sub sun}, consistent with the expectations for the ACT southern sky survey. These masses are compared to the ACT SZE properties of the sample, specifically, the match-filtered central SZE amplitude y{sub 0}-tilde, the central Compton parameter y{sub 0}, and the integrated Compton signal Y{sub 200c}, which we use to derive SZE-mass scaling relations. All SZE estimators correlate with dynamical mass with low intrinsic scatter ({approx}< 20%), in agreement with numerical simulations. We explore the effects of various systematic effects on these scaling relations, including the correlation between observables and the influence of dynamically disturbed clusters. Using the three-dimensional information available, we divide the sample into relaxed and disturbed clusters and find that {approx}50% of the clusters are disturbed. There are hints that disturbed systems might bias the scaling relations, but given the current sample sizes, these differences are not significant; further studies including more clusters are required to assess the impact of these clusters on the scaling relations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A topological analysis of large-scale structure, studied using the CMASS sample of SDSS-III

We study the three-dimensional genus topology of large-scale structure using the northern region of the CMASS Data Release 10 (DR10) sample of the SDSS-III Baryon Oscillation Spectroscopic Survey. We select galaxies with redshift 0.452 < z < 0.625 and with a stellar mass M {sub stellar} > 10{sup 11.56} M {sub ☉}. We study the topology at two smoothing lengths: R {sub G} = 21 h {sup –1} Mpc and R {sub G} = 34 h {sup –1} Mpc. The genus topology studied at the R {sub G} = 21 h {sup –1} Mpc scale results in the highest genus amplitude observed to date. The CMASS sample yields a genus curve that is characteristic of one produced by Gaussian random phase initial conditions. The data thus support the standard model of inflation where random quantum fluctuations in the early universe produced Gaussian random phase initial conditions. Modest deviations in the observed genus from random phase are as expected from shot noise effects and the nonlinear evolution of structure. We suggest the use of a fitting formula motivated by perturbation theory to characterize the shift and asymmetries in the observed genus curve with a single parameter. We construct 54 mock SDSS CMASS surveys along the past light cone from the Horizon Run 3 (HR3) N-body simulations, where gravitationally bound dark matter subhalos are identified as the sites of galaxy formation. We study the genus topology of the HR3 mock surveys with the same geometry and sampling density as the observational sample and find the observed genus topology to be consistent with ΛCDM as simulated by the HR3 mock samples. We conclude that the topology of the large-scale structure in the SDSS CMASS sample is consistent with cosmological models having primordial Gaussian density fluctuations growing in accordance with general relativity to form galaxies in massive dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The Lyα reference sample. I. Survey outline and first results for Markarian 259

The Lyα Reference Sample (LARS) is a substantial program with the Hubble Space Telescope (HST) that provides a sample of local universe laboratory galaxies in which to study the detailed astrophysics of the visibility and strength of the Lyαline of neutral hydrogen. Lyα is the dominant spectral line in use for characterizing high-redshift (z) galaxies. This paper presents an overview of the survey, its selection function, and HST imaging observations. The sample was selected from the combined GALEX+Sloan Digital Sky Survey catalog at z = 0.028-0.19, in order to allow Lyα to be captured with combinations of long-pass filters in the Solar Blind Channel (SBC) of the Advanced Camera for Surveys (ACS) onboard HST. In addition, LARS utilizes Hα and Hβ narrowband and u, b, i broadband imaging with ACS and the Wide Field Camera 3 (WFC3). In order to study galaxies in which large numbers of Lyα photons are produced (whether or not they escape), we demanded an Hα equivalent width W(Hα) ≥100 Å. The final sample of 14 galaxies covers far-UV (FUV, λ ∼ 1500 Å) luminosities that overlap with those of high-z Lyα emitters (LAEs) and Lyman break galaxies (LBGs), making LARS a valid comparison sample. We present the reduction steps used to obtain the Lyα images, including our LARS eXtraction software (LaXs), which utilizes pixel-by-pixel spectral synthesis fitting of the energy distribution to determine and subtract the continuum at Lyα. We demonstrate that the use of SBC long-pass-filter combinations increase the signal-to-noise ratio by an order of magnitude compared to the nominal Lyα filter available in SBC. To exemplify the science potential of LARS, we also present some first results for a single galaxy, Mrk 259 (LARS #1). This irregular galaxy shows bright and extended (indicative of resonance scattering) but strongly asymmetric Lyα emission. Spectroscopy from the Cosmic Origins Spectrograph on board HST centered on the brightest UV knot shows a moderate outflow in the neutral interstellar medium (probed by low ionization stage absorption features) and Lyα emission with an asymmetric profile. Radiative transfer modeling is able to reproduce the essential features of the Lyα line profile and confirms the presence of an outflow. From the integrated photometry we measure an Lyα luminosity of L {sub Lyα}=1.3 × 10{sup 42} erg s{sup –1} an equivalent width W(Lyα) = 45 Å and an FUV absolute magnitude M {sub FUV} = –19.2 (AB). Mrk 259 would hence be detectable in high-z Lyα and LBG surveys. The total Lyα escape fraction is 12%. This number is higher than the low-z average, but similar to that at z > 4, demonstrating that LARS provides a valid comparison sample for high-z galaxy studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

THE IMPORTANCE OF THE MAGNETIC FIELD FROM AN SMA-CSO-COMBINED SAMPLE OF STAR-FORMING REGIONS

Submillimeter dust polarization measurements of a sample of 50 star-forming regions, observed with the Submillimeter Array (SMA) and the Caltech Submillimeter Observatory (CSO) covering parsec-scale clouds to milliparsec-scale cores, are analyzed in order to quantify the magnetic field importance. The magnetic field misalignment δ—the local angle between magnetic field and dust emission gradient—is found to be a prime observable, revealing distinct distributions for sources where the magnetic field is preferentially aligned with or perpendicular to the source minor axis. Source-averaged misalignment angles (|δ|) fall into systematically different ranges, reflecting the different source-magnetic field configurations. Possible bimodal (|δ|) distributions are found for the separate SMA and CSO samples. Combining both samples broadens the distribution with a wide maximum peak at small (|δ|) values. Assuming the 50 sources to be representative, the prevailing source-magnetic field configuration is one that statistically prefers small magnetic field misalignments |δ|. When interpreting |δ| together with a magnetohydrodynamics force equation, as developed in the framework of the polarization-intensity gradient method, a sample-based log-linear scaling fits the magnetic field tension-to-gravity force ratio (Σ {sub B}) versus (|δ|) with (Σ {sub B}) = 0.116 · exp (0.047 · (|δ|)) ± 0.20 (mean error), providing a way to estimate the relative importance of the magnetic field, only based on measurable field misalignments |δ|. The force ratio Σ {sub B} discriminates systems that are collapsible on average ((Σ {sub B}) < 1) from other molecular clouds where the magnetic field still provides enough resistance against gravitational collapse ((Σ {sub B}) > 1). The sample-wide trend shows a transition around (|δ|) ≈ 45°. Defining an effective gravitational force ∼1 – (Σ {sub B}), the average magnetic-field-reduced star formation efficiency is at least a factor of two smaller than the free-fall efficiency. For about one fourth of the sources the average efficiency drops to zero. The force ratio Σ {sub B} can further be linked to the normalized mass-to-flux ratio, yielding an estimate for the latter one without the need of field strength measurements. Across the sample, a transition from magnetically supercritical to subcritcal is observed with growing misalignment (|δ|)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Insights on the dusty torus and neutral torus from optical and X-ray obscuration in a complete volume limited hard X-ray AGN sample

We describe a complete volume limited sample of nearby active galaxies selected by their 14–195 keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We also describe a complementary sample of inactive galaxies, selected to match the host galaxy properties. The active sample appears to have no bias in terms of active galactic nucleus (AGN) type, the only difference being the neutral absorbing column, which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, logL{sub 14−195 keV}[erg s{sup −1}]=42.4–43.7, the optically obscured and X-ray absorbed fractions are 50%–65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad-line region (BLR) emission, and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The redshift distribution of cosmological samples: a forward modeling approach

Determining the redshift distribution n ( z ) of galaxy samples is essential for several cosmological probes including weak lensing. For imaging surveys, this is usually done using photometric redshifts estimated on an object-by-object basis. We present a new approach for directly measuring the global n ( z ) of cosmological galaxy samples, including uncertainties, using forward modeling. Our method relies on image simulations produced using \textsc(UFig) (Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation) within the MCCL (Monte-Carlo Control Loops) framework. The galaxy population is modeled using parametric forms for the luminosity functions, spectral energy distributions, sizes and radial profiles of both blue and red galaxies. We apply exactly the same analysis to the real data and to the simulated images, which also include instrumental and observational effects. By adjusting the parameters of the simulations, we derive a set of acceptable models that are statistically consistent with the data. We then apply the same cuts to the simulations that were used to construct the target galaxy sample in the real data. The redshifts of the galaxies in the resulting simulated samples yield a set of n ( z ) distributions for the acceptable models. We demonstrate the method by determining n ( z ) for a cosmic shear like galaxy sample from the 4-band Subaru Suprime-Cam data in the COSMOS field. We also complement this imaging data with a spectroscopic calibration sample from the VVDS survey. We compare our resulting posterior n ( z ) distributions to the one derived from photometric redshifts estimated using 36 photometric bands in COSMOS and find good agreement. This offers good prospects for applying our approach to current and future large imaging surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Forecasting super-sample covariance in future weak lensing surveys with SuperSCRAM

The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constraints for future surveys requires accurately accounting for the super-sample effects. If super-sample covariance is in fact a major component of the survey error budget, it may be necessary to investigate mitigation strategies that constrain the specific realization of the long-wavelength modes. We present a Fisher matrix based formalism for approximating the magnitude of super-sample covariance and the effectiveness of mitigation strategies for realistic survey geometries. We implement our formalism in the public code SuperSCRAM: Super-Sample Covariance Reduction and Mitigation. We illustrate SuperSCRAM with an example application, where the modes contributing to supersample covariance in the WFIRST weak lensing survey are constrained by the low-redshift galaxy number counts in the wider LSST footprint. We find that super-sample covariance increases the volume of the error ellipsoid in 7D cosmological parameter space by a factor of 4.5 relative to Gaussian statistical errors only, but our simple mitigation strategy more than halves the contamination, to a factor of 2.0.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DESI luminous red galaxy samples for cross-correlations

We present two galaxy samples, selected from DESI Legacy Imaging Surveys (LS) DR9, with approximately 20,000 square degrees of coverage and spectroscopic redshift distributions designed for cross-correlations such as with CMB lensing, galaxy lensing, and the Sunyaev-Zel'dovich effect. The first sample is identical to the DESI Luminous Red Galaxy (LRG) sample, and the second sample is an extended LRG sample with 2–3 times the DESI LRG density. We present the improved photometric redshifts, tomographic binning and their spectroscopic redshift distributions and imaging systematics weights, and magnification bias coefficients. The catalogs and related data products will be made publicly available. The cosmological constraints using this sample and Planck lensing maps are presented in a companion paper. Here, we also make public the new set of general-purpose photometric redshifts trained using DESI spectroscopic redshifts, which are used in this work, for all galaxies in LS DR9.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Skip Nav Destination The effect of environment on Type Ia supernovae in the Dark Energy Survey three-year cosmological sample

Analyses of Type Ia supernovae (SNe Ia) have found puzzling correlations between their standardized luminosities and host galaxy properties: SNe Ia in high-mass, passive hosts appear brighter than those in lower mass, star-forming hosts. We examine the host galaxies of SNe Ia in the Dark Energy Survey 3-yr spectroscopically confirmed cosmological sample, obtaining photometry in a series of ‘local’ apertures centred on the SN, and for the global host galaxy. We study the differences in these host galaxy properties, such as stellar mass and rest-frame U − R colours, and their correlations with SN Ia parameters including Hubble residuals. We find all Hubble residual steps to be >3σ in significance, both for splitting at the traditional environmental property sample median and for the step of maximum significance. For stellar mass, we find a maximal local step of 0.098 ± 0.018 mag; ∼0.03 mag greater than the largest global stellar mass step in our sample (0.070 ± 0.017 mag). When splitting at the sample median, differences between local and global U − R steps are small, both ∼0.08 mag, but are more significant than the global stellar mass step (0.057 ± 0.017 mag). We split the data into sub-samples based on SN Ia light-curve parameters: stretch (x_1) and colour (c), finding that redder objects (c > 0) have larger Hubble residual steps, for both stellar mass and U − R, for both local and global measurements, of ∼0.14 mag. Additionally, the bluer (star-forming) local environments host a more homogeneous SN Ia sample, with local U − R rms scatter as low as 0.084 ± 0.017 mag for blue (c < 0) SNe Ia in locally blue U − R environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Probing gravity with the DES-CMASS sample and BOSS spectroscopy

ABSTRACT The DES-CMASS sample (DMASS) is designed to optimally combine the weak lensing measurements from the Dark Energy Survey (DES) and redshift-space distortions (RSD) probed by the CMASS galaxy sample from the Baryonic Oscillation Spectroscopic Survey. In this paper, we demonstrate the feasibility of adopting DMASS as the equivalent of CMASS for a joint analysis of DES and BOSS in the framework of modified gravity. We utilize the angular clustering of the DMASS galaxies, cosmic shear of the DES metacalibration sources, and cross-correlation of the two as data vectors. By jointly fitting the combination of the data with the RSD measurements from the CMASS sample and Planck data, we obtain the constraints on modified gravity parameters $\mu _0=-0.37^{+0.47}_{-0.45}$ and $\Sigma _0=0.078^{+0.078}_{-0.082}$. Our constraints of modified gravity with DMASS are tighter than those with the DES Year 1 redMaGiC sample with the same external data sets by 29 per cent for μ0 and 21 per cent for Σ0, and comparable to the published results of the DES Year 1 modified gravity analysis despite this work using fewer external data sets. This improvement is mainly because the galaxy bias parameter is shared and more tightly constrained by both CMASS and DMASS, effectively breaking the degeneracy between the galaxy bias and other cosmological parameters. Such an approach to optimally combine photometric and spectroscopic surveys using a photometric sample equivalent to a spectroscopic sample can be applied to combining future surveys having a limited overlap such as DESI and LSST.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quantifying the poor purity and completeness of morphological samples selected by galaxy colour

The galaxy population is strongly bimodal in both colour and morphology, and the two measures correlate strongly, with most blue galaxies being late-types (spirals) and most early-types, typically ellipticals, being red. This observation has led to the use of colour as a convenient selection criterion to make samples that are then labelled by morphology. Such use of colour as a proxy for morphology results in necessarily impure and incomplete samples. In this paper, we make use of the morphological labels produced by Galaxy Zoo to measure how incomplete and impure such samples are, considering optical (ugriz), near-ultraviolet (NUV), and near-infrared (NIR; JHK) bands. The best single colour optical selection is found using a threshold of g - r = 0.742, but this still results in a sample where only 56 per cent of red galaxies are smooth and 56 per cent of smooth galaxies are red. Use of the NUV gives some improvement over purely optical bands, particularly for late-types, but still results in low purity/completeness for early-types. No significant improvement is found by adding NIR bands. With any two bands, including NUV, a sample of early-types with greater than two-thirds purity cannot be constructed. In conclusion, advances in quantitative galaxy morphologies have made colour–morphology proxy selections largely unnecessary going forward; where such assumptions are still required, we recommend studies carefully consider the implications of sample incompleteness/impurity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evidence of galaxy assembly bias in SDSS DR7 galaxy samples from count statistics

We present observational constraints on the galaxy–halo connection, focusing particularly on galaxy assembly bias from a novel combination of counts-in-cylinders statistics, P(N CIC ), with the standard measurements of the projected two-point correlation function w p (r p ), and number density n gal of galaxies. We measure n gal , w p (r p ), and P(N CIC ) for volume-limited, luminosity-threshold samples of galaxies selected from SDSS DR7, and use them to constrain halo occupation distribution (HOD) models, including a model in which galaxy occupation depends upon a secondary halo property, namely halo concentration. We detect significant positive central assembly bias for the M r < -20.0 and M r <-19.5 samples. Central galaxies preferentially reside within haloes of high concentration at fixed mass. Positive central assembly bias is also favoured in the M r < -20.5 and M r < -19.0 samples. We find no evidence of central assembly bias in the M r < -21.0 sample. We observe only a marginal preference for negative satellite assembly bias in the M r < -20.0 and M r < -19.0 samples, and non-zero satellite assembly bias is not indicated in other samples. Our findings underscore the necessity of accounting for galaxy assembly bias when interpreting galaxy survey data, and demonstrate the potential of count statistics in extracting information from the spatial distribution of galaxies, which could be applied to both galaxy–halo connection studies and cosmological analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

nautilus : boosting Bayesian importance nested sampling with deep learning

ABSTRACT We introduce a novel approach to boost the efficiency of the importance nested sampling (INS) technique for Bayesian posterior and evidence estimation using deep learning. Unlike rejection-based sampling methods such as vanilla nested sampling (NS) or Markov chain Monte Carlo (MCMC) algorithms, importance sampling techniques can use all likelihood evaluations for posterior and evidence estimation. However, for efficient importance sampling, one needs proposal distributions that closely mimic the posterior distributions. We show how to combine INS with deep learning via neural network regression to accomplish this task. We also introduce nautilus, a reference open-source python implementation of this technique for Bayesian posterior and evidence estimation. We compare nautilus against popular NS and MCMC packages, including emcee, dynesty, ultranest, and pocomc, on a variety of challenging synthetic problems and real-world applications in exoplanet detection, galaxy SED fitting and cosmology. In all applications, the sampling efficiency of nautilus is substantially higher than that of all other samplers, often by more than an order of magnitude. Simultaneously, nautilus delivers highly accurate results and needs fewer likelihood evaluations than all other samplers tested. We also show that nautilus has good scaling with the dimensionality of the likelihood and is easily parallelizable to many CPUs.

97 MATHEMATICS AND COMPUTING↗