-
The Farmer: A reproducible profile-fitting photometry package for deep galaxy surveys
Authors:
J. R. Weaver,
L. Zalesky,
V. Kokorev,
C. J. R. McPartland,
N. Chartab,
K. M. L. Gould,
M. Shuntov,
I. Davidzon,
A. Faisst,
N. Stickley,
P. L. Capak,
S. Toft,
D. Masters,
B. Mobasher,
D. B. Sanders,
O. B. Kauffmann,
H. J. McCracken,
O. Ilbert,
G. Brammer,
A. Moneti
Abstract:
While space-borne optical and near-infrared facilities have succeeded in delivering a precise and spatially resolved picture of our Universe, their small survey area is known to under-represent the true diversity of galaxy populations. Ground-based surveys have reached comparable depths but at lower spatial resolution, resulting in source confusion that hampers accurate photometry extractions. Wha…
▽ More
While space-borne optical and near-infrared facilities have succeeded in delivering a precise and spatially resolved picture of our Universe, their small survey area is known to under-represent the true diversity of galaxy populations. Ground-based surveys have reached comparable depths but at lower spatial resolution, resulting in source confusion that hampers accurate photometry extractions. What once was limited to the infrared regime has now begun to challenge ground-based ultra-deep surveys, affecting detection and photometry alike. Failing to address these challenges will mean forfeiting a representative view into the distant Universe. We introduce The Farmer: an automated, reproducible profile-fitting photometry package that pairs a library of smooth parametric models from The Tractor (Lang et al. 2016) with a decision tree that determines the best-fit model in concert with neighboring sources. Photometry is measured by fitting the models on other bands leaving brightness free to vary. The resulting photometric measurements are naturally total, and no aperture corrections are required. Supporting diagnostics (e.g. $χ^2$) enable measurement validation. As fitting models is relatively time intensive, The Farmer is built with high-performance computing routines. We benchmark The Farmer on a set of realistic COSMOS-like images and find accurate photometry, number counts, and galaxy shapes. The Farmer is already being utilized to produce catalogs for several large-area deep extragalactic surveys where it has been shown to tackle some of the most challenging optical and near-infrared data available, with the promise of extending to other ultra-deep surveys expected in the near future. The Farmer is available to download from GitHub and Zenodo.
△ Less
Submitted 11 October, 2023;
originally announced October 2023.
-
Joint Survey Processing II: Stellar Proper Motions in the COSMOS Field, from Hubble Space Telescope ACS and Subaru Telescope HSC Observations
Authors:
Sergio B. Fajardo-Acosta,
Andreas Faisst,
Carl J. Grillmair,
Ranga-Ram Chary,
Roberta Paladini,
Ben Rusholme,
Nathaniel Stickley
Abstract:
We analyze stellar proper motions in the COSMOS field to assess the presence of bulk motions. At bright magnitudes (G-band 18.5--20.76 AB), we use the proper motions of 1,010 stars in the Gaia DR2 catalog. At the faint end, we computed proper motions of 11,519 point-like objects at i-band magnitudes 19--25 AB using Hubble ACS and Subaru HSC, which span two epochs about 11 years apart. In order to…
▽ More
We analyze stellar proper motions in the COSMOS field to assess the presence of bulk motions. At bright magnitudes (G-band 18.5--20.76 AB), we use the proper motions of 1,010 stars in the Gaia DR2 catalog. At the faint end, we computed proper motions of 11,519 point-like objects at i-band magnitudes 19--25 AB using Hubble ACS and Subaru HSC, which span two epochs about 11 years apart. In order to measure these proper motions with unprecedented accuracy at faint magnitudes, we developed a foundational set of astrometric tools which will be required for Joint Survey Processing (JSP) of data from the next generation of optical/infrared surveys. The astrometric grids of Hubble ACS and Subaru HSC mosaics were corrected at the catalog level, using proper motion-propagated and parallax-corrected Gaia DR2 sources. These astrometric corrections were verified using compact extragalactic sources. Upon comparison of our measured proper motions with Gaia DR2, we estimate the uncertainties in our measurements to be ~2--3 mas/yr per axis, down to 25.5 AB mag. We corrected proper motions for the mean motion of the Sun, and we find that late-type main-sequence stars predominantly in the thin disk in the COSMOS field have space velocities mainly towards the Galactic center. We detect candidate high-velocity (> 220 km/s) stars, 6 of them at ~0.4-6 kpc from the Gaia sample, and 5 of them at ~20 kpc from the faint star HSC and ACS sample. The sources from the faint star sample may be candidate halo members of the Sangarius stream.
△ Less
Submitted 10 March, 2022;
originally announced March 2022.
-
Joint Survey Processing I: Compact oddballs in the COSMOS field -- low-luminosity Quasars at z > 6?
Authors:
Andreas L. Faisst,
Ranga Ram Chary,
Sergio Fajardo-Acosta,
Roberta Paladini,
Benjamin Rusholme,
Nathaniel Stickley,
George Helou,
John R. Weaver,
Gabriel Brammer,
Anton M. Koekemoer,
Hironao Miyatake
Abstract:
The faint-end slope of the quasar luminosity function at z~6 and its implication on the role of quasars in reionizing the intergalactic medium at early times has been an outstanding problem for some time. The identification of faint high-redshift quasars with luminosities of <1e44.5 erg/s is challenging. They are rare (few per square-degree) and the separation of these unresolved quasars from late…
▽ More
The faint-end slope of the quasar luminosity function at z~6 and its implication on the role of quasars in reionizing the intergalactic medium at early times has been an outstanding problem for some time. The identification of faint high-redshift quasars with luminosities of <1e44.5 erg/s is challenging. They are rare (few per square-degree) and the separation of these unresolved quasars from late-type stars and compact star-forming galaxies is difficult from ground-based observations alone. In addition, source confusion becomes significant at >25mag, with 30% of sources having their flux contaminated by foreground objects when the seeing resolution is ~0.7". We mitigate these issues by performing a pixel-level joint processing of ground and space-based data from Subaru/HSC and HST/ACS. We create a deconfused catalog over the 1.64 square-degrees of the COSMOS field, after accounting for spatial varying PSFs and astrometric differences between the two datasets. We identify twelve low-luminosity (M_UV ~ -21 mag) z>6 quasar candidates through (i) their red color measured between ACS/F814W and HSC/i-band and (ii) their compactness in the space-based data. Non-detections of our candidates in Hubble DASH data argues against contamination from late-type stars. Our constraints on the faint end of the quasar luminosity function at z~6.4 suggests a negligibly small contribution to reionization compared to the star-forming galaxy population. The confirmation of our candidates and the evolution of number density with redshift could provide better insights into how supermassive galaxies grew in the first billion years of cosmic time.
△ Less
Submitted 26 March, 2022; v1 submitted 17 March, 2021;
originally announced March 2021.
-
Joint Survey Processing of Euclid, Rubin and Roman: Final Report
Authors:
R. Chary,
G. Helou,
G. Brammer,
P. Capak,
A. Faisst,
D. Flynn,
S. Groom,
H. C. Ferguson,
C. Grillmair,
S. Hemmati,
A. Koekemoer,
B. Lee,
S. Malhotra,
H. Miyatake,
P. Melchior,
I. Momcheva,
J. Newman,
J. Masiero,
R. Paladini,
A. Prakash,
B. Rusholme,
N. R. Stickley,
A. Smith,
W. M. Wood-Vasey,
H. I. Teplitz
Abstract:
The Euclid, Rubin/LSST and Roman (WFIRST) projects will undertake flagship optical/near-infrared surveys in the next decade. By mapping thousands of square degrees of sky and covering the electromagnetic spectrum between 0.3 and 2 microns with sub-arcsec resolution, these projects will detect several tens of billions of sources, enable a wide range of astrophysical investigations by the astronomic…
▽ More
The Euclid, Rubin/LSST and Roman (WFIRST) projects will undertake flagship optical/near-infrared surveys in the next decade. By mapping thousands of square degrees of sky and covering the electromagnetic spectrum between 0.3 and 2 microns with sub-arcsec resolution, these projects will detect several tens of billions of sources, enable a wide range of astrophysical investigations by the astronomical community and provide unprecedented constraints on the nature of dark energy and dark matter. The ultimate cosmological, astrophysical and time-domain science yield from these missions will require joint survey processing (JSP) functionality at the pixel level that is outside the scope of the individual survey projects. The JSP effort scoped here serves two high-level objectives: 1) provide precise concordance multi-wavelength images and catalogs over the entire sky area where these surveys overlap, which accounts for source confusion and mismatched isophotes, and 2) provide a science platform to analyze concordance images and catalogs to enable a wide range of astrophysical science goals to be formulated and addressed by the research community. For the cost of about 200WY, JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also allow for the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.
△ Less
Submitted 24 August, 2020;
originally announced August 2020.
-
Joint Survey Processing of LSST, Euclid and WFIRST: Enabling a broad array of astrophysics and cosmology through pixel level combinations of datasets
Authors:
R. Chary,
G. Brammer,
P. Capak,
W. Dawson,
A. Faisst,
S. Fajardo-Acosta,
H. C. Ferguson,
C. J. Grillmair,
S. Hemmati,
A. Koekemoer,
B. Lee,
R. Lupton,
S. Malhotra,
P. Melchior,
I. Momcheva,
J. Newman,
J. Masiero,
R. Paladini,
A. Prakash,
J. Rhodes,
B. Rusholme,
M. Schneider,
N. Stickley,
A. Smith,
W. M. Wood-Vasey
, et al. (1 additional authors not shown)
Abstract:
Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper…
▽ More
Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.
△ Less
Submitted 2 October, 2019;
originally announced October 2019.
-
Stellar Velocity Dispersion in Mergers: The Effects of Dust and Star Formation
Authors:
Nathaniel R. Stickley,
Gabriela Canalizo
Abstract:
We investigate the effects of stellar evolution and dust on measurements of stellar velocity dispersion in mergers of disk galaxies. $N$-body simulations and radiative transfer analysis software are used to obtain mass-weighted and flux-weighted measurements of stellar velocity dispersion. We find that the distribution of dust with respect to the distribution of young stars in such systems is more…
▽ More
We investigate the effects of stellar evolution and dust on measurements of stellar velocity dispersion in mergers of disk galaxies. $N$-body simulations and radiative transfer analysis software are used to obtain mass-weighted and flux-weighted measurements of stellar velocity dispersion. We find that the distribution of dust with respect to the distribution of young stars in such systems is more important than the total degree of attenuation. The presence of dust typically causes flux-weighted measurements of stellar velocity dispersion to be elevated with respect to mass-weighted measurements because dust preferentially obscures young stars, which tend to be dynamically cooler than older stellar populations in such systems. In exceptional situations, in which young stars are not preferentially obscured by dust, flux-weighted velocity dispersion measurements tend to be negatively offset with respect to mass-weighted measurements because the dynamically cool young stellar populations are more luminous, per unit mass, than older stellar populations. Our findings provide a context for comparing observationally-obtained measurements of velocity dispersion with measurements of velocity dispersion obtained from galaxy merger simulations.
△ Less
Submitted 8 October, 2016;
originally announced October 2016.
-
An Empirical Approach to Cosmological Galaxy Survey Simulation: Application to SPHEREx Low-Resolution Spectroscopy
Authors:
Nathaniel R. Stickley,
Peter Capak,
Daniel Masters,
Roland de Putter,
Olivier Doré,
Jamie Bock
Abstract:
Highly accurate models of the galaxy population over cosmological volumes are necessary in order to predict the performance of upcoming cosmological missions. We present a data-driven model of the galaxy population constrained by deep 0.1-8 $\rm μm$ imaging and spectroscopic data in the COSMOS survey, with the immediate goal of simulating the spectroscopic redshift performance of the proposed SPHE…
▽ More
Highly accurate models of the galaxy population over cosmological volumes are necessary in order to predict the performance of upcoming cosmological missions. We present a data-driven model of the galaxy population constrained by deep 0.1-8 $\rm μm$ imaging and spectroscopic data in the COSMOS survey, with the immediate goal of simulating the spectroscopic redshift performance of the proposed SPHEREx mission. SPHEREx will obtain over the full-sky $R\sim41$ spectrophotometry at moderate spatial resolution ($\sim6"$) over the wavelength range 0.75-4.18 $\rm μm$ and $R\sim135$ over the wavelength range 4.18-5 $\rm μm$. We show that our simulation accurately reproduces a range of known galaxy properties, encapsulating the full complexity of the galaxy population and enables realistic, full end-to-end simulations to predict mission performance. Finally, we discuss potential applications of the simulation framework to future cosmology missions and give a description of released data products.
△ Less
Submitted 20 June, 2016;
originally announced June 2016.
-
NebulOS: A Big Data Framework for Astrophysics
Authors:
Nathaniel R. Stickley,
Miguel A. Aragon-Calvo
Abstract:
We introduce NebulOS, a Big Data platform that allows a cluster of Linux machines to be treated as a single computer. With NebulOS, the process of writing a massively parallel program for a datacenter is no more complicated than writing a Python script for a desktop computer. The platform enables most pre-existing data analysis software to be used, as scale, in a datacenter without modification. T…
▽ More
We introduce NebulOS, a Big Data platform that allows a cluster of Linux machines to be treated as a single computer. With NebulOS, the process of writing a massively parallel program for a datacenter is no more complicated than writing a Python script for a desktop computer. The platform enables most pre-existing data analysis software to be used, as scale, in a datacenter without modification. The shallow learning curve and compatibility with existing software greatly reduces the time required to develop distributed data analysis pipelines. The platform is built upon industry-standard, open-source Big Data technologies, from which it inherits several fault tolerance features. NebulOS enhances these technologies by adding an intuitive user interface, automated task monitoring, and other usability features. We present a summary of the architecture, provide usage examples, and discuss the system's performance scaling.
△ Less
Submitted 14 September, 2016; v1 submitted 7 March, 2015;
originally announced March 2015.
-
Stellar Velocity Dispersion in Dissipative Galaxy Mergers with Star Formation
Authors:
Nathaniel R. Stickley,
Gabriela Canalizo
Abstract:
In order to better understand stellar dynamics in merging systems, such as NGC 6240, we examine the evolution of central stellar velocity dispersion (\sig) in dissipative galaxy mergers using a suite of binary disk merger simulations that include feedback from stellar formation and active galactic nuclei (AGNs). We find that $σ_*$ undergoes the same general stages of evolution that were observed i…
▽ More
In order to better understand stellar dynamics in merging systems, such as NGC 6240, we examine the evolution of central stellar velocity dispersion (\sig) in dissipative galaxy mergers using a suite of binary disk merger simulations that include feedback from stellar formation and active galactic nuclei (AGNs). We find that $σ_*$ undergoes the same general stages of evolution that were observed in our previous dissipationless simulations: coherent oscillation, then phase mixing, followed by dynamical equilibrium. We also find that measurements of $σ_*$ that are based only upon the youngest stars in simulations consistently yield lower values than measurements based upon the total stellar population. This finding appears to be consistent with the so-called "$σ_*$ discrepancy," observed in real galaxies. We note that quasar-level AGN activity is much more likely to occur when $σ_*$ is near its equilibrium value rather than during periods of extreme $σ_*$. Finally, we provide estimates of the scatter inherent in measuring $σ_*$ in ongoing mergers.
△ Less
Submitted 19 March, 2014; v1 submitted 25 October, 2013;
originally announced October 2013.
-
The Evolution of Stellar Velocity Dispersion During Dissipationless Galaxy Mergers
Authors:
Nathaniel R. Stickley,
Gabriela Canalizo
Abstract:
Using N-body simulations, we studied the detailed evolution of central stellar velocity dispersion, σ, during dissipationless binary mergers of galaxies. Stellar velocity dispersion was measured using the common mass-weighting method as well as a flux-weighting method designed to simulate the technique used by observers. A toy model for dust attenuation was introduced in order to study the effect…
▽ More
Using N-body simulations, we studied the detailed evolution of central stellar velocity dispersion, σ, during dissipationless binary mergers of galaxies. Stellar velocity dispersion was measured using the common mass-weighting method as well as a flux-weighting method designed to simulate the technique used by observers. A toy model for dust attenuation was introduced in order to study the effect of dust attenuation on measurements of σ. We found that there are three principal stages in the evolution of σ in such mergers: oscillation, phase mixing, and dynamical equilibrium. During the oscillation stage, σ undergoes damped oscillations of increasing frequency. The oscillation stage is followed by a phase mixing stage during which the amplitude of the variations in σ is smaller and more chaotic than in the oscillation stage. Upon reaching dynamical equilibrium, σ assumes a stable value. We used our data regarding the evolution of σ during mergers to characterize the scatter inherent in making measurements of σ in non-quiescent systems. In particular, we found that σ does not fall below 70% nor exceed 200% of its final, quiescent value during a merger and that a random measurement of σ in such a system is much more likely to fall near the equilibrium value than near an extremum. Our toy model of dust attenuation suggested that dust can systematically reduce observational measurements of σ and increase the scatter in σ measurements.
△ Less
Submitted 22 January, 2012;
originally announced January 2012.
