ProFit : Bayesian profile fitting of galaxy images

We present ProFit, a new code for Bayesian two-dimensional photometric galaxy profile modelling. ProFit consists of a low-level c++ library (libprofit), accessible via a command-line interface and documented API, along with high-level R (ProFit) and Python (PyProFit) interfaces (available at github.com/ICRAR/libprofit, github.com/ICRAR/ProFit, and github.com/ICRAR/pyprofit, respectively). R ProFit is also available pre-built from cran; however, this version will be slightly behind the latest GitHub version. libprofit offers fast and accurate two-dimensional integration for a useful number of profiles, including Sérsic, Core-Sérsic, broken-exponential, Ferrer, Moffat, empirical King, point-source, and sky, with a simple mechanism for adding new profiles. We show detailed comparisons between libprofit and galfit. libprofit is both faster and more accurate than galfit at integrating the ubiquitous Sérsic profile for the most common values of the Sérsic index n (0.5 < n < 8). The high-level fitting code ProFit is tested on a sample of galaxies with both SDSS and deeper KiDS imaging. We find good agreement in the fit parameters, with larger scatter in best-fitting parameters from fitting images from different sources (SDSS versus KiDS) than from using different codes (ProFit versus galfit). A large suite of Monte Carlo-simulated images are used to assess prospects for automated bulge-disc decomposition with ProFit on SDSS, KiDS, and future LSST imaging. We find that the biggest increases in fit quality come from moving from SDSS- to KiDS-quality data, with less significant gains moving from KiDS to LSST.Publisher PDFPeer reviewe

Architecture of the Andromeda galaxy : a quantitative analysis of clustering in the inner stellar halo

We present a quantitative measurement of the amount of clustering present in the inner ∼30 kpc of the stellar halo of the Andromeda galaxy (M31). For this we analyse the angular positions and radial velocities of the carefully selected planetary nebulae in the M31 stellar halo. We study the cumulative distribution of pairwise distances in angular position and line-of-sight velocity space, and find that the M31 stellar halo contains substantially more stars in the form of close pairs as compared to that of a featureless smooth halo. In comparison to a smoothed/scrambled distribution, we estimate that the clustering excess in the M31 inner halo is roughly 40 per cent at maximum and on average ∼20 per cent. Importantly, comparing against the 11 stellar halo models of Bullock & Johnston, which were simulated within the context of the ΛCDM (Λ cold dark matter) cosmological paradigm, we find that the amount of substructures in the M31 stellar halo closely resembles that of a typical ΛCDM halo.Publisher PDFPeer reviewe

Galaxy And Mass Assembly (GAMA): curation and reanalysis of 16.6k redshifts in the G10/COSMOS region

We discuss the construction of the Galaxy And Mass Assembly (GAMA) 10h region (G10) using publicly available data in the Cosmic Evolution Survey region (COSMOS) in order to extend the GAMA survey to z ∼ 1 in a single deg2 field. In order to obtain the maximum number of high precision spectroscopic redshifts we re-reduce all archival zCOSMOS-bright data and use the GAMA automatic cross-correlation redshift fitting code autoz. We use all available redshift information (autoz, zCOSMOS-bright 10k, PRIMUS, VVDS, SDSS and photometric redshifts) to calculate robust best-fitting redshifts for all galaxies and visually inspect all 1D and 2D spectra to obtain 16 583 robust redshifts in the full COSMOS region. We then define the G10 region to be the central ∼1 deg2 of COSMOS, which has relatively high spectroscopic completeness, and encompasses the CHILES VLA region. We define a combined r < 23.0 mag and i < 22.0 mag G10 sample (selected to have the highest bijective overlap) with which to perform future analysis, containing 9861 sources with reliable high-precision VLT-VIMOS spectra. All tables, spectra and imaging are available at http://ict.icrar.org/cutout/G10

GAMA/G10-COSMOS/3D-HST : evolution of the galaxy stellar mass function over 12.5 Gyr

Using a combined and consistently analysed GAMA, G10-COSMOS, and 3D-HST data set, we explore the evolution of the galaxy stellar mass function over lookback times tL ∈ [0.2, 12.5] h70-1 Gyr. We use a series of volume-limited samples to fit Schechter functions in bins of ~constant lookback time and explore the evolution of the best-fitting parameters in both single and two-component cases. In all cases, we employ a fitting procedure that is robust to the effects of Eddington bias and sample variance. Surprisingly, when fitting a two-component Schechter function, we find essentially no evidence of temporal evolution in M*, the two α slope parameters, or the normalization of the low-mass component. Instead, our fits suggest that the various shape parameters have been exceptionally stable over cosmic time, as has the normalization of the low-mass component, and that the evolution of the stellar mass function is well described by a simple build-up of the high-mass component over time. When fitting a single component Schechter function, there is an observed evolution in both M* and α however, this is interpreted as being an artefact. Finally, we find that the evolution of the stellar mass function, and the observed stellar mass density, can be well described by a simple model of constant growth in the high-mass source density over the last 11 h70-1 Gyr.Publisher PDFPeer reviewe

ProFound : source extraction and application to modern survey data

Parts of this research were conducted by the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020.We introduce ProFound, a source finding and image analysis package. ProFound provides methods to detect sources in noisy images, generate segmentation maps identifying the pixels belonging to each source, and measure statistics like flux, size and ellipticity. These inputs are key requirements of ProFit, our recently released galaxy profiling package, where the design aim is that these two software packages will be used in unison to semi-automatically profile large samples of galaxies. The key novel feature introduced in ProFound is that all photometry is executed on dilated segmentation maps that fully contain the identifiable flux, rather than using more traditional circular or ellipse based photometry. Also, to be less sensitive to pathological segmentation issues, the de-blending is made across saddle points in flux. We apply ProFound in a number of simulated and real world cases, and demonstrate that it behaves reasonably given its stated design goals. In particular, it offers good initial parameter estimation for ProFit, and also segmentation maps that follow the sometimes complex geometry of resolved sources, whilst capturing nearly all of the flux. A number of bulge-disc decomposition projects are already making use of the ProFound and ProFit pipeline, and adoption is being encouraged by publicly releasing the software for the open source R data analysis platform under an LGPL-3 license on GitHub (github.com/asgr/ProFound).PostprintPeer reviewe