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

Comparison of the Stellar Populations of Bulges and Discs using the MaNGA Survey

We use the MaNGA integral-field spectroscopic survey of low-redshift galaxies to compare the stellar populations of the bulge and disc components, identified from their Sersic profiles, for various samples of galaxies. Bulge dominated regions tend to be more metal-rich and have slightly older stellar ages than their associated disc dominated regions. The metallicity difference is consistent with the deeper gravitational potential in bulges relative to discs, which allows bulges to retain more of the metals produced by stars. The age difference is due to star formation persisting longer in discs relative to bulges. Relative to galaxies with lower stellar masses, galaxies with higher stellar masses tend to have bulge dominated regions that are more metal-rich and older (in light-weighted measurements) than their disc dominated regions. This suggests high-mass galaxies quench from the inside out, while lower-mass galaxies quench across the whole galaxy simultaneously. Early-type galaxies tend to have bulge dominated regions the same age as their disc dominated regions, while late-type galaxies tend to have disc dominated regions significantly younger than their bulge dominated regions. Central galaxies tend to have a greater metallicity difference between their bulge dominated regions and disc dominated regions than satellite galaxies at similar stellar mass. This difference may be explained by central galaxies being subject to mergers or extended gas accretion bringing new, lower-metallicity gas to the disc, thereby reducing the average metallicity and age of the stars; quenching of satellite discs may also play a role.Comment: Accepted by PAS

Galaxy And Mass Assembly (GAMA): the dependence of star formation on surface brightness in low-redshift galaxies

The star-formation rate in galaxies is well known to correlate with stellar mass (the ‘star-forming main sequence’). Here, we extend this further to explore any additional dependence on galaxy surface brightness, a proxy for stellar mass surface density. We use a large sample of low-redshift (z ≤ 0.08) galaxies from the Galaxy And Mass Assembly survey which have both spectral energy distribution (SED) derived star-formation rates and photometric bulge-disc decompositions, the latter providing measures of disc surface brightness and disc masses. Using two samples, one of galaxies fitted by a single component with Sérsic index below 2 and one of the discs from two-component fits, we find that once the overall mass dependence of star-formation rate is accounted for, there is no evidence in either sample for a further dependence on stellar surface density

ProPane: Image Warping with Fire

In this paper we introduce the software package ProPane, written for the R data analysis language. ProPane combines the full range of wcslib projections with the C++ image manipulation routines provided by the CImg library. ProPane offers routines for image warping and combining (including stacking), and various related tasks such as image alignment tweaking and pixel masking. It can stack an effectively unlimited number of target frames using multiple parallel cores, and offers threading for many lower level routines. It has been used for a number of current and upcoming large surveys, and we present a range of its capabilities and features. ProPane is already available under a permissive open-source LGPL-3 license at github.com/asgr/ProPane (DOI: 10.5281/zenodo.10057053).Comment: 19 pages, 17 figures, 5 tables, accepted to MNRA

Galapagos-2/Galfitm/Gama – Multi-wavelength measurement of galaxy structure: Separating the properties of spheroid and disk components in modern surveys

Aims. We present the capabilities of GALAPAGOS-2 and GALFITM in the context of fitting two-component profiles – bulge–disk decompositions – to galaxies, with the ultimate goal of providing complete multi-band, multi-component fitting of large samples of galaxies in future surveys. We also release both the code and the fit results to 234 239 objects from the DR3 of the GAMA survey, a sample significantly deeper than in previous works. Methods. We use stringent tests on both simulated and real data, as well as comparison to public catalogues to evaluate the advantages of using multi-band over single-band data. Results. We show that multi-band fitting using GALFITM provides significant advantages when trying to decompose galaxies into their individual constituents, as more data are being used, by effectively being able to use the colour information buried in the individual exposures to its advantage. Using simulated data, we find that multi-band fitting significantly reduces deviations from the real parameter values, allows component sizes and Sérsic indices to be recovered more accurately, and – by design – constrains the band-to-band variations of these parameters to more physical values. On both simulated and real data, we confirm that the spectral energy distributions (SEDs) of the two main components can be recovered to fainter magnitudes compared to using single-band fitting, which tends to recover ‘disks’ and ‘bulges’ with – on average – identical SEDs when the galaxies become too faint, instead of the different SEDs they truly have. By comparing our results to those provided by other fitting codes, we confirm that they agree in general, but measurement errors can be significantly reduced by using the multi-band tools developed by the MEGAMORPH project. Conclusions. We conclude that the multi-band fitting employed by GALAPAGOS-2 and GALFITM significantly improves the accuracy of structural galaxy parameters and enables much larger samples to be be used in a scientific analysis

Galaxy and Mass Assembly (GAMA): Demonstrating the Power of WISE in the Study of Galaxy Groups to z \u3c 0.1

Combining high-fidelity group characterization from the Galaxy and Mass Assembly survey and source-tailored z \u3c 0.1 photometry from the Wide-Field Infrared Survey Explorer (WISE) survey, we present a comprehensive study of the properties of ungrouped galaxies, compared to 497 galaxy groups (4 ≤ N FoF ≤ 20) as a function of stellar and halo mass. Ungrouped galaxies are largely unimodal in WISE color, the result of being dominated by star-forming, late-type galaxies. Grouped galaxies, however, show a clear bimodality in WISE color, which correlates strongly with stellar mass and morphology. We find evidence for an increasing early-type fraction, in stellar mass bins between 1010 M o˙ ≲ M stellar ≲ 1011 M o˙, with increasing halo mass. Using ungrouped, late-type galaxies with star-forming colors (W2-W3 \u3e 3), we define a star-forming main sequence (SFMS), which we use to delineate systems that have moved below the sequence ( quenched for the purposes of this work). We find that with increasing halo mass, the relative number of late-type systems on the SFMS decreases, with a corresponding increase in early-type, quenched systems at high stellar mass (M stellar \u3e 1010.5 M o˙), consistent with mass quenching. Group galaxies with masses M stellar \u3c 1010.5 M o˙ show evidence of quenching consistent with environmentally driven processes. The stellar mass distribution of late-type, quenched galaxies suggests that it may be an intermediate population as systems transition from being star-forming and late-type to the red sequence. Finally, we use the projected area of groups on the sky to extract groups that are (relatively) compact for their halo mass. Although these show a marginal increase in their proportion of high-mass and early-type galaxies compared to nominal groups, a clear increase in quenched fraction is not evident

Galaxy And Mass Assembly (GAMA): bulge-disc decomposition of KiDS data in the nearby Universe

We derive single Sérsic fits and bulge-disc decompositions for 13 096 galaxies at redshifts z \u3c 0.08 in the GAMA II equatorial survey regions in the Kilo-Degree Survey (KiDS) g, r, and i bands. The surface brightness fitting is performed using the Bayesian two-dimensional profile fitting code PROFIT. We fit three models to each galaxy in each band independently with a fully automated Markov chain Monte Carlo analysis: a single Sérsic model, a Sérsic plus exponential and a point source plus exponential. After fitting the galaxies, we perform model selection and flag galaxies for which none of our models are appropriate (mainly mergers/Irregular galaxies). The fit quality is assessed by visual inspections, comparison to previous works, comparison of independent fits of galaxies in the overlap regions between KiDS tiles and bespoke simulations. The latter two are also used for a detailed investigation of systematic error sources. We find that our fit results are robust across various galaxy types and image qualities with minimal biases. Errors given by the MCMC underestimate the true errors typically by factors 2–3. Automated model selection criteria are accurate to \u3e90 per cent role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3e\u3e90 per cent\u3e90 per cent as calibrated by visual inspection of a subsample of galaxies. We also present g−r component colours and the corresponding colour–magnitude diagram, consistent with previous works despite our increased fit flexibility. Such reliable structural parameters for the components of a diverse sample of galaxies across multiple bands will be integral to various studies of galaxy properties and evolution. All results are integrated into the GAMA database

Galaxy And Mass Assembly (GAMA): the dependence of star formation on surface brightness in low-redshift galaxies

The star-formation rate in galaxies is well known to correlate with stellar mass (the 'star-forming main sequence'). Here, we extend this further to explore any additional dependence on galaxy surface brightness, a proxy for stellar mass surface density. We use a large sample of low-redshift (z ≤ 0.08) galaxies from the Galaxy And Mass Assembly survey which have both spectral energy distribution (SED) derived star-formation rates and photometric bulge-disc decompositions, the latter providing measures of disc surface brightness and disc masses. Using two samples, one of galaxies fitted by a single component with Sersic index below 2 and one of the discs from two-component fits, we find that once the overall mass dependence of star-formation rate is accounted for, there is no evidence in either sample for a further dependence on stellar surface density

The SAMI Galaxy Survey: Decomposed stellar kinematics of galaxy bulges and disks

We investigate the stellar kinematics of the bulge and disk components in 826 galaxies with a wide range of morphology from the Sydney-AAO Multi-object Integral-field spectroscopy Galaxy Survey. The spatially resolved rotation velocity (V) and velocity dispersion (σ ) of bulge and disk components have been simultaneously estimated using the penalized pixel fitting (PPXF) method with photometrically defined weights for the two components. We introduce a new subroutine of PPXF for dealing with degeneracy in the solutions. We show that the V and σ distributions in each galaxy can be reconstructed using the kinematics and weights of the bulge and disk components. The combination of two distinct components provides a consistent description of the major kinematic features of galaxies over a wide range of morphological types. We present Tully-Fisher and Faber-Jackson relations showing that the galaxy stellar mass scales with both V and σ for both components of all galaxy types. We find a tight Faber-Jackson relation even for the disk component. We show that the bulge and disk components are kinematically distinct: (1) the two components show scaling relations with similar slopes, but different intercepts; (2) the spin parameter λR indicates bulges are pressure-dominated systems and disks are supported by rotation; and (3) the bulge and disk components have, respectively, low and high values in intrinsic ellipticity. Our findings suggest that the relative contributions of the two components explain, at least to first order, the complex kinematic behaviour of galaxies

The SAMI Galaxy Survey: Decomposed Stellar Kinematics of Galaxy Bulges and Disks

We investigate the stellar kinematics of the bulge and disk components in 826 galaxies with a wide range of morphology from the Sydney-AAO Multi-object Integral-field spectroscopy (SAMI) Galaxy Survey. The spatially-resolved rotation velocity (V) and velocity dispersion ($\sigma$) of bulge and disk components have been simultaneously estimated using the penalized pixel fitting (pPXF) method with photometrically defined weights for the two components. We introduce a new subroutine of pPXF for dealing with degeneracy in the solutions. We show that the V and $\sigma$ distributions in each galaxy can be reconstructed using the kinematics and weights of the bulge and disk components. The combination of two distinct components provides a consistent description of the major kinematic features of galaxies over a wide range of morphological types. We present Tully-Fisher and Faber-Jackson relations showing that the galaxy stellar mass scales with both V and $\sigma$ for both components of all galaxy types. We find a tight Faber-Jackson relation even for the disk component. We show that the bulge and disk components are kinematically distinct: (1) the two components show scaling relations with similar slopes, but different intercepts; (2) the spin parameter $\lambda_R$ indicates bulges are pressure-dominated systems and disks are supported by rotation; (3) the bulge and disk components have, respectively, low and high values in intrinsic ellipticity. Our findings suggest that the relative contributions of the two components explain, at least to first order, the complex kinematic behaviour of galaxies.Comment: 22 pages, 21 figures; Accepted for publication in MNRA