No Arabic abstract
Numerous ongoing and future large area surveys (e.g. DES, EUCLID, LSST, WFIRST), will increase by several orders of magnitude the volume of data that can be exploited for galaxy morphology studies. The full potential of these surveys can only be unlocked with the development of automated, fast and reliable analysis methods. In this paper we present DeepLeGATo, a new method for two-dimensional photometric galaxy profile modeling, based on convolutional neural networks. Our code is trained and validated on analytic profiles (HST/CANDELS F160W filter) and it is able to retrieve the full set of parameters of one- component Sersic models: total magnitude, effective radius, Sersic index, axis ratio. We show detailed comparisons between our code and GALFIT. On simulated data, our method is more accurate than GALFIT and 3000 time faster on GPU (50 times when run on the same CPU). On real data, DeepLeGATo trained on simulations behaves similarly to GALFIT on isolated galaxies. With a fast domain adaptation step made with the 0.1 - 0.8 per cent the size of the training set, our code is easily capable to reproduce the results obtained with GALFIT even on crowded regions. DeepLeGATo does not require any human intervention beyond the training step, rendering it much automated than traditional profiling methods. The development of this method for more complex models (two-component galaxies, variable PSF, dense sky regions) could constitute a fundamental tool in the era of big data in astronomy.
Searches for low-surface-brightness galaxies (LSBGs) in galaxy surveys are plagued by the presence of a large number of artifacts (e.g., objects blended in the diffuse light from stars and galaxies, Galactic cirrus, star-forming regions in the arms of spiral galaxies, etc.) that have to be rejected through time consuming visual inspection. In future surveys, which are expected to collect hundreds of petabytes of data and detect billions of objects, such an approach will not be feasible. We investigate the use of convolutional neural networks (CNNs) for the problem of separating LSBGs from artifacts in survey images. We take advantage of the fact that, for the first time, we have available a large number of labeled LSBGs and artifacts from the Dark Energy Survey, that we use to train, validate, and test a CNN model. That model, which we call DeepShadows, achieves a test accuracy of $92.0 %$, a significant improvement relative to feature-based machine learning models. We also study the ability to use transfer learning to adapt this model to classify objects from the deeper Hyper-Suprime-Cam survey, and we show that after the model is retrained on a very small sample from the new survey, it can reach an accuracy of $87.6%$. These results demonstrate that CNNs offer a very promising path in the quest to study the low-surface-brightness universe.
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 Sersic, Core-Sersic, 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 Serrsic profile for the most common values of the Serrsic 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-fit parameters from fitting images from different sources (SDSS vs KiDS) than from using different codes (ProFit vs 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.
A tool for representation of the one-dimensional astrometric signal of Gaia is described and investigated in terms of fit discrepancy and astrometric performance with respect to number of parameters required. The proposed basis function is based on the aberration free response of the ideal telescope and its derivatives, weighted by the source spectral distribution. The influence of relative position of the detector pixel array with respect to the optical image is analysed, as well as the variation induced by the source spectral emission. The number of parameters required for micro-arcsec level consistency of the reconstructed function with the detected signal is found to be 11. Some considerations are devoted to the issue of calibration of the instrument response representation, taking into account the relevant aspects of source spectrum and focal plane sampling. Additional investigations and other applications are also suggested.
In an effort to probe the origin of surface brightness profile (SBP) breaks widely observed in nearby disk galaxies, we carry out a comparative study of stellar population profiles of 635 disk galaxies selected from the MaNGA spectroscopic survey. We classify our galaxies into single exponential (TI), down-bending (TII) and up-bending (TIII) SBP types, and derive their spin parameters and radial profiles of age/metallicity-sensitive spectral features. Most TII (TIII) galaxies have down-bending (up-bending) star formation rate (SFR) radial profiles, implying that abrupt radial changes of SFR intensities contribute to the formation of both TII and TIII breaks. Nevertheless, a comparison between our galaxies and simulations suggests that stellar migration plays a significant role in weakening down-bending $Sigma_{star}$ profile breaks. While there is a correlation between the break strengths of SBPs and age/metallicity-sensitive spectral features for TII galaxies, no such correlation is found for TIII galaxies, indicating that stellar migration may not play a major role in shaping TIII breaks, as is evidenced by a good correspondence between break strengths of $Sigma_{star}$ and surface brightness profiles of TIII galaxies. We do not find evidence for galaxy spin being a relevant parameter for forming different SBP types, nor do we find significant differences between the asymmetries of galaxies with different SBP types, suggesting that environmental disturbances or satellite accretion in the recent past do not significantly influence the break formation. By dividing our sample into early and late morphological types, we find that galaxies with different SBP types follow nearly the same tight stellar mass-$R_{25}$ relation, which makes the hypothesis that stellar migration alone can transform SBP types from TII to TI and then to TIII highly unlikely.
We present zELDA(redshift Estimator for Line profiles of Distant Lyman-Alpha emitters), an open source code to fit Lyman-Alpha (Lya) line profiles. The main motivation is to provide the community with an easy to use and fast tool to analyze Lya line profiles uniformly to improve the understating of Lya emitting galaxies. zELDA is based on line profiles of the commonly used shell-model pre-computed with the full Monte Carlo radiative transfer code LyaRT. Via interpolation between these spectra and the addition of noise, we assemble a suite of realistic Lya spectra which we use to train a deep neural network. We show that the neural network can predict the model parameters to high accuracy (e.g.,.0.34 dex HI column density for R=12000) and thus allows for a significant speedup over existing fitting methods. As a proof of concept, we demonstrate the potential of zELDA by fitting 97 observed Lya line profiles from the LASD data base. Comparing the fitted value with the measured systemic redshift of these sources, we find that Lya determines their rest frame Lya wavelength with a remarkable good accuracy of 0.3A (75 km/s). Comparing the predicted outflow properties and the observed Lya luminosity and equivalent width, we find several possible trends. For example, we find an anticorrelation between the Lya luminosity and the outflow neutral hydrogen column density, which might be explained by the radiative transfer process within galaxies