-
Extreme Scale Survey Simulation with Python Workflows
Authors:
A. S. Villarreal,
Yadu Babuji,
Tom Uram,
Daniel S. Katz,
Kyle Chard,
Katrin Heitmann
Abstract:
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will soon carry out an unprecedented wide, fast, and deep survey of the sky in multiple optical bands. The data from LSST will open up a new discovery space in astronomy and cosmology, simultaneously providing clues toward addressing burning issues of the day, such as the origin of dark energy and and the nature of dark matter, w…
▽ More
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will soon carry out an unprecedented wide, fast, and deep survey of the sky in multiple optical bands. The data from LSST will open up a new discovery space in astronomy and cosmology, simultaneously providing clues toward addressing burning issues of the day, such as the origin of dark energy and and the nature of dark matter, while at the same time yielding data that will, in turn, pose fresh new questions. To prepare for the imminent arrival of this remarkable data set, it is crucial that the associated scientific communities be able to develop the software needed to analyze it. Computational power now available allows us to generate synthetic data sets that can be used as a realistic training ground for such an effort. This effort raises its own challenges -- the need to generate very large simulations of the night sky, scaling up simulation campaigns to large numbers of compute nodes across multiple computing centers with different architectures, and optimizing the complex workload around memory requirements and widely varying wall clock times. We describe here a large-scale workflow that melds together Python code to steer the workflow, Parsl to manage the large-scale distributed execution of workflow components, and containers to carry out the image simulation campaign across multiple sites. Taking advantage of these tools, we developed an extreme-scale computational framework and used it to simulate five years of observations for 300 square degrees of sky area. We describe our experiences and lessons learned in developing this workflow capability, and highlight how the scalability and portability of our approach enabled us to efficiently execute it on up to 4000 compute nodes on two supercomputers.
△ Less
Submitted 24 September, 2021;
originally announced September 2021.
-
How to Optimally Constrain Galaxy Assembly Bias: Supplement Projected Correlation Functions with Count-in-cells Statistics
Authors:
Kuan Wang,
Yao-Yuan Mao,
Andrew R. Zentner,
Frank C. van den Bosch,
Johannes U. Lange,
Chad M. Schafer,
Antonia S. Villarreal,
Andrew P. Hearin,
Duncan Campbell
Abstract:
Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real U…
▽ More
Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $ΔΣ(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+ΔΣ(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $ΔΣ(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection.
△ Less
Submitted 6 September, 2019; v1 submitted 22 March, 2019;
originally announced March 2019.
-
Updated Results on the Galaxy-Halo Connection from Satellite Kinematics in SDSS
Authors:
Johannes U. Lange,
Frank C. van den Bosch,
Andrew R. Zentner,
Kuan Wang,
Antonia S. Villarreal
Abstract:
We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter…
▽ More
We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter haloes, as well as the radial profile of satellite galaxies. Confirming previous results, we find that red central galaxies live in more massive haloes than blue galaxies at fixed luminosity. Additionally, our results suggest that satellite galaxies have a radial profile less centrally concentrated than dark matter but not as cored as resolved subhaloes in dark matter-only simulations. Compared to previous works using satellite kinematics by More et al., we find much more competitive constraints on the galaxy-halo connection, on par with those derived from a combination of galaxy clustering and galaxy-galaxy lensing. We compare our results on the galaxy-halo connection to other studies using galaxy clustering and group catalogues, showing very good agreement between these different techniques. We discuss future applications of satellite kinematics in the context of constraining cosmology and the relationship between galaxies and dark matter haloes.
△ Less
Submitted 8 November, 2018;
originally announced November 2018.
-
Maturing Satellite Kinematics into a Competitive Probe of the Galaxy-Halo Connection
Authors:
Johannes U. Lange,
Frank C. van den Bosch,
Andrew R. Zentner,
Kuan Wang,
Antonia S. Villarreal
Abstract:
The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on…
▽ More
The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on average. In this work, we test the impact of various modelling assumptions on constraints derived from analysing satellite phase-space distributions in the non-linear, 1-halo regime. We discuss their potential to explain the discrepancy between average halo masses derived from satellite kinematics and gravitational lensing previously reported. Furthermore, we develop an updated, more robust analysis to extract constraints on the galaxy-halo relation from satellite properties in spectroscopic galaxy surveys such as the SDSS. We test the accuracy of this approach using a large number of realistic mock catalogues. Furthermore, we find that constraints derived from such an analysis are complementary and competitive with respect to the commonly used galaxy clustering and galaxy-galaxy lensing observables.
△ Less
Submitted 24 October, 2018;
originally announced October 2018.
-
The Immitigable Nature of Assembly Bias: The Impact of Halo Definition on Assembly Bias
Authors:
Antonia S. Villarreal,
Andrew R. Zentner,
Yao-Yuan Mao,
Chris W. Purcell,
Frank C. van den Bosch,
Benedikt Diemer,
Johannes U. Lange,
Kuan Wang,
Duncan Campbell
Abstract:
Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $Δ$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo m…
▽ More
Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $Δ$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $Δ$. For conventional halo definitions ($Δ\sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $Δ\sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $Δ$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $Δ\gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.
△ Less
Submitted 16 September, 2017; v1 submitted 11 May, 2017;
originally announced May 2017.
-
On the Mass Distribution and Birth Masses of Neutron Stars
Authors:
Feryal Ozel,
Dimitrios Psaltis,
Ramesh Narayan,
Antonia Sierra Villarreal
Abstract:
We investigate the distribution of neutron star masses in different populations of binaries, employing Bayesian statistical techniques. In particular, we explore the differences in neutron star masses between sources that have experienced distinct evolutionary paths and accretion episodes. We find that the distribution of neutron star masses in non-recycled eclipsing high-mass binaries as well as…
▽ More
We investigate the distribution of neutron star masses in different populations of binaries, employing Bayesian statistical techniques. In particular, we explore the differences in neutron star masses between sources that have experienced distinct evolutionary paths and accretion episodes. We find that the distribution of neutron star masses in non-recycled eclipsing high-mass binaries as well as of slow pulsars, which are all believed to be near their birth masses, has a mean of 1.28 M_solar and a dispersion of 0.24 M_solar. These values are consistent with expectations for neutron star formation in core-collapse supernovae. On the other hand, double neutron stars, which are also believed to be near their birth masses, have a much narrower mass distribution, peaking at 1.33 M_solar but with a dispersion of only 0.05 M_solar. Such a small dispersion cannot easily be understood and perhaps points to a particular and rare formation channel. The mass distribution of neutron stars that have been recycled has a mean of 1.48 M_solar and a dispersion of 0.2 M_solar, consistent with the expectation that they have experienced extended mass accretion episodes. The fact that only a very small fraction of recycled neutron stars in the inferred distribution have masses that exceed ~2 M_solar suggests that only a few of these neutron stars cross the mass threshold to form low mass black holes.
△ Less
Submitted 7 September, 2012; v1 submitted 4 January, 2012;
originally announced January 2012.
-
New light curves and ephemeris for the close eclipsing binary V963 PER
Authors:
Andrew P. Odell,
Patrick Wils,
Clarissa Dirks,
Blythe Guvenen,
Jo C. O'Malley,
Antonia S. Villarreal,
Rita M. Weinzettle
Abstract:
We have obtained CCD photometry in 2010-11 of V963 Per (=GSC3355 0394), which is a recently identified close binary star with unequal eclipse depths. The seven new eclipse timings yield an improved ephemeris, but we caution that secondary eclipse can be affected by variation of the light curve. This variation seems to be on a monthly timescale at the few percent level.
We have obtained CCD photometry in 2010-11 of V963 Per (=GSC3355 0394), which is a recently identified close binary star with unequal eclipse depths. The seven new eclipse timings yield an improved ephemeris, but we caution that secondary eclipse can be affected by variation of the light curve. This variation seems to be on a monthly timescale at the few percent level.
△ Less
Submitted 1 August, 2011;
originally announced August 2011.