We present a two-component Machine Learning (ML) based approach for classifying astronomical images by data-quality via an examination of sources detected in the images and image pixel values from representative sources within those images. The first component, which uses a clustering algorithm, creates a proper and small fraction of the image pixels to determine the quality of the observation. The representative images (and associated tables) are ~800 times smaller than the original images, significantly reducing the time required to train our algorithm. The useful information in the images is preserved, permitting them to be classified in different categories, but the required storage is reduced. The second component, which is a deep neural network model, classifies the representative images. Using ground-based telescope imaging data, we demonstrate that the method can be used to separate 'usable' images from those that present some problems for scientific projects -- such as images that were taken in sub-optimal conditions. This method uses two different data sets as input to a deep model and provides better performance than if we only used the images' pixel information. The method may be used in cases where large and complex data sets should be examined using deep models. Our automated classification approach achieves 97% agreement when compared to classification generated via manual image inspection. We compare our method with traditional results and show that the method improves the results by about 10%, and also presents more comprehensive outcomes.

Typically, binary classification lens-finding schemes are used to discriminate between lens candidates and non-lenses. However, these models often suffer from substantial false-positive classifications. Such false positives frequently occur due to images containing objects such as crowded sources, galaxies with arms, and also images with a central source and smaller surrounding sources. Therefore, a model might confuse the stated circumstances with an Einstein ring. It has been proposed that by allowing such commonly misclassified image types to constitute their own classes, machine learning models will more easily be able to learn the difference between images that contain real lenses, and images that contain lens imposters. Using Hubble Space Telescope (HST) images, in the F814W filter, we compare the usage of binary and multi-class classification models applied to the lens finding task. From our findings, we conclude there is not a significant benefit to using the multi-class model over a binary model. We will also present the results of a simple lens search using a multi-class machine learning model, and potential new lens candidates.

Using a sample of 11,478 spaxels in 34 galaxies with molecular gas, star formation and stellar maps taken from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we investigate the parameters that correlate with variations in star formation rates on kpc scales. We use a combination of correlation statistics and an artificial neural network to quantify the parameters that drive both the absolute star formation rate surface density (Sigma_SFR), as well as its scatter around the resolved star forming main sequence (Delta Sigma_SFR). We find that Sigma_SFR is primarily regulated by molecular gas surface density (Sigma_H2) with a secondary dependence on stellar mass surface density (Sigma_*), as expected from an `extended Kennicutt-Schmidt relation'. However, Delta Sigma_SFR is driven primarily by changes in star formation efficiency (SFE), with variations in gas fraction playing a secondary role. Taken together, our results demonstrate that whilst the absolute rate of star formation is primarily set by the amount of molecular gas, the variation of star formation rate above and below the resolved star forming main sequence (on kpc scales) is primarily due to changes in SFE.

We present an analysis of star formation and quenching in the SDSS-IV MaNGA-DR15, utilising over 5 million spaxels from $\sim$3500 local galaxies. We estimate star formation rate surface densities ($\Sigma_{\rm SFR}$) via dust corrected $H\alpha$ flux where possible, and via an empirical relationship between specific star formation rate (sSFR) and the strength of the 4000 Angstrom break (D4000) in all other cases. We train a multi-layered artificial neural network (ANN) and a random forest (RF) to classify spaxels into `star forming' and `quenched' categories given various individual (and groups of) parameters. We find that global parameters (pertaining to the galaxy as a whole) perform collectively the best at predicting when spaxels will be quenched, and are substantially superior to local/ spatially resolved and environmental parameters. Central velocity dispersion is the best single parameter for predicting quenching in central galaxies. We interpret this observational fact as a probable consequence of the total integrated energy from AGN feedback being traced by the mass of the black hole, which is well known to correlate strongly with central velocity dispersion. Additionally, we train both an ANN and RF to estimate $\Sigma_{\rm SFR}$ values directly via regression in star forming regions. Local/ spatially resolved parameters are collectively the most predictive at estimating $\Sigma_{\rm SFR}$ in these analyses, with stellar mass surface density at the spaxel location ($\Sigma_*$) being by far the best single parameter. Thus, quenching is fundamentally a global process but star formation is governed locally by processes within each spaxel.

The application of machine learning (ML) methods to the analysis of astrophysical datasets is on the rise, particularly as the computing power and complex algorithms become more powerful and accessible. As the field of ML enjoys a continuous stream of breakthroughs, its applications demonstrate the great potential of ML, ranging from achieving tens of millions of times increase in analysis speed (e.g., modeling of gravitational lenses or analysing spectroscopic surveys) to solutions of previously unsolved problems (e.g., foreground subtraction or efficient telescope operations). The number of astronomical publications that include ML has been steadily increasing since 2010. With the advent of extremely large datasets from a new generation of surveys in the 2020s, ML methods will become an indispensable tool in astrophysics. Canada is an unambiguous world leader in the development of the field of machine learning, attracting large investments and skilled researchers to its prestigious AI Research Institutions. This provides a unique opportunity for Canada to also be a world leader in the application of machine learning in the field of astrophysics, and foster the training of a new generation of highly skilled researchers.

Machine learning is becoming a popular tool to quantify galaxy morphologies and identify mergers. However, this technique relies on using an appropriate set of training data to be successful. By combining hydrodynamical simulations, synthetic observations and convolutional neural networks (CNNs), we quantitatively assess how realistic simulated galaxy images must be in order to reliably classify mergers. Specifically, we compare the performance of CNNs trained with two types of galaxy images, stellar maps and dust-inclusive radiatively transferred images, each with three levels of observational realism: (1) no observational effects (idealized images), (2) realistic sky and point spread function (semi-realistic images), (3) insertion into a real sky image (fully realistic images). We find that networks trained on either idealized or semi-real images have poor performance when applied to survey-realistic images. In contrast, networks trained on fully realistic images achieve 87.1% classification performance. Importantly, the level of realism in the training images is much more important than whether the images included radiative transfer, or simply used the stellar maps (87.1% compared to 79.6% accuracy, respectively). Therefore, one can avoid the large computational and storage cost of running radiative transfer with a relatively modest compromise in classification performance. Making photometry-based networks insensitive to colour incurs a very mild penalty to performance with survey-realistic data (86.0% with r-only compared to 87.1% with gri). This result demonstrates that while colour can be exploited by colour-sensitive networks, it is not necessary to achieve high accuracy and so can be avoided if desired. We provide the public release of our statistical observational realism suite, RealSim, as a companion to this paper.

We present a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]$<-1.2$, the algorithm succeeds in identifying more than 70% of the giants in the training/test set, with a dwarf contamination fraction below 30% (with respect to the SDSS/SEGUE dataset). The photometric metallicity estimates have uncertainties better than 0.2 dex when compared with the spectroscopic measurements. The distances estimated by the algorithm are valid out to a distance of at least $\sim 80$ kpc without requiring any prior on the stellar distribution, and have fully independent uncertainities that take into account both random and systematic errors. These advances allow us to estimate these stellar parameters for approximately 12 million stars in the photometric dataset. This will enable studies involving the chemical mapping of the distant outer disc and the stellar halo, including their kinematics using the Gaia proper motions. This type of algorithm can be applied in the Southern hemisphere to the first release of LSST data, thus providing an almost complete view of the external components of our Galaxy out to at least $\sim 80$ kpc. Critical to the success of these efforts will be ensuring well-defined spectroscopic training sets that sample a broad range of stellar parameters with minimal biases. A catalogue containing the training/test set and all relevant parameters within the public footprint of CFIS is available online.

When multiple star-forming gas structures overlap along the line-of-sight and emit optically thin emission at significantly different radial velocities, the emission can become non-Gaussian and often exhibits two distinct peaks. Traditional line-fitting techniques can fail to account adequately for these double-peaked profiles, providing inaccurate cloud kinematics measurements. We present a new method called Convnet Line-fitting Of Velocities in Emission-line Regions (CLOVER) for distinguishing between one-component, two-component, and noise-only emission lines using 1D convolutional neural networks trained with synthetic spectral cubes. CLOVER utilizes spatial information in spectral cubes by predicting on $3\times3$ pixel sub-cubes, using both the central pixel's spectrum and the average spectrum over the $3\times3$ grid as input. On an unseen set of 10,000 synthetic spectral cubes in each predicted class, CLOVER has classification accuracies of $\sim99\%$ for the one-component class and $\sim97\%$ for the two-component class. For the noise-only class, which is analogous to a signal-to-noise cutoff of four for traditional line-fitting methods, CLOVER has classification accuracy of $100\%$. CLOVER also has exceptional performance on real observations, correctly distinguishing between the three classes across a variety of star-forming regions. In addition, CLOVER quickly and accurately extracts kinematics directly from spectra identified as two-component class members. Moreover, we show that CLOVER is easily scalable to emission lines with hyperfine splitting, making it an attractive tool in the new era of large-scale NH$_3$ and N$_2$H$^+$ mapping surveys.

We investigate the dependence of galaxy structure on a variety of galactic and environmental parameters for ~500,000 galaxies at z<0.2, taken from the Sloan Digital Sky Survey data release 7 (SDSS-DR7). We utilise bulge-to-total stellar mass ratio, (B/T)_*, as the primary indicator of galactic structure, which circumvents issues of morphological dependence on waveband. We rank galaxy and environmental parameters in terms of how predictive they are of galaxy structure, using an artificial neural network approach. We find that distance from the star forming main sequence (Delta_SFR), followed by stellar mass (M_*), are the most closely connected parameters to (B/T)_*, and are significantly more predictive of galaxy structure than global star formation rate (SFR), or any environmental metric considered (for both central and satellite galaxies). Additionally, we make a detailed comparison to the Illustris hydrodynamical simulation and the LGalaxies semi-analytic model. In both simulations, we find a significant lack of bulge-dominated galaxies at a fixed stellar mass, compared to the SDSS. This result highlights a potentially serious problem in contemporary models of galaxy evolution.

The spectrum of a galaxy contains information about its physical properties. Classifying spectra using templates helps elucidate the nature of a galaxy's energy sources. In this paper, we investigate the use of self-organizing maps in classifying galaxy spectra against templates. We trained semi-supervised self-organizing map networks using a set of templates covering the wavelength range from far ultraviolet to near infrared. The trained networks were used to classify the spectra of a sample of 142 galaxies with 0.5 < z < 1 and the results compared to classifications performed using K-means clustering, a supervised neural network, and chi-squared minimization. Spectra corresponding to quiescent galaxies were more likely to be classified similarly by all methods while starburst spectra showed more variability. Compared to classification using chi-squared minimization or the supervised neural network, the galaxies classed together by the self-organizing map had more similar spectra. The class ordering provided by the one-dimensional self-organizing maps corresponds to an ordering in physical properties, a potentially important feature for the exploration of large datasets.

In the absence of the two emission lines H$\alpha$ and [NII] (6584\AA) in a BPT diagram, we show that other spectral information is sufficiently informative to distinguish AGN galaxies from star-forming galaxies. We use pattern recognition methods and a sample of galaxy spectra from the Sloan Digital Sky Survey (SDSS) to show that, in this survey, the flux and equivalent width of [OIII] (5007\AA) and H$\beta$, along with the 4000\AA-break, can be used to classify galaxies in a BPT diagram. This method provides a higher accuracy of predictions than those which use stellar mass and [OIII]/H$\beta$. First, we use BPT diagrams and various physical parameters to re-classify the galaxies. Next, using confusion matrices, we determine the `correctly' predicted classes as well as confused cases. In this way, we investigate the effect of each parameter in the confusion matrices and rank the physical parameters used in the discrimination of the different classes. We show that in this survey, for example, $\rm{g - r}$ colour can provide the same accuracy as galaxy stellar mass to predict whether or not a galaxy hosts an AGN. Finally, with the same information, we also rank the parameters involved in the discrimination of Seyfert and LINER galaxies.

The application of artificial neural networks (ANNs) for the estimation of HI gas mass fraction (\fgas) is investigated, based on a sample of 13,674 galaxies in the Sloan Digital Sky Survey (SDSS) with HI detections or upper limits from the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA). We show that, for an example set of fixed input parameters ($g-r$ colour and $i$-band surface brightness), a multidimensional quadratic model yields \fgas\ scaling relations with a smaller scatter (0.22 dex) than traditional linear fits (0.32 dex), demonstrating that non-linear methods can lead to an improved performance over traditional approaches. A more extensive ANN analysis is performed using 15 galaxy parameters that capture variation in stellar mass, internal structure, environment and star formation. Of the 15 parameters investigated, we find that $g-r$ colour, followed by stellar mass surface density, bulge fraction and specific star formation rate have the best connection with \fgas. By combining two control parameters, that indicate how well a given galaxy in SDSS is represented by the ALFALFA training set (\pr) and the scatter in the training procedure (\sigf), we develop a strategy for quantifying which SDSS galaxies our ANN can be adequately applied to, and the associated errors in the \fgas\ estimation. In contrast to previous works, our \fgas\ estimation has no systematic trend with galactic parameters such as M$_{\star}$, $g-r$ and SFR. We present a catalog of \fgas\ estimates for more than half a million galaxies in the SDSS, of which $\sim$ 150,000 galaxies have a secure selection parameter with average scatter in the \fgas\ estimation of 0.22 dex.

We quantify the impact that a variety of galactic and environmental properties have on the quenching of star formation. We collate a sample of $\sim$ 400,000 central and $\sim$ 100,000 satellite galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Specifically, we consider central velocity dispersion ($\sigma_{c}$), stellar, halo, bulge and disk mass, local density, bulge-to-total ratio, group-centric distance and galaxy-halo mass ratio. We develop and apply a new statistical technique to quantify the impact on the quenched fraction ($f_{\rm Quench}$) of varying one parameter, while keeping the remaining parameters fixed. For centrals, we find that the $f_{\rm Quench} - \sigma_{c}$ relationship is tighter and steeper than for any other variable considered. We compare to the Illustris hydrodynamical simulation and the Munich semi-analytic model (L-Galaxies), finding that our results for centrals are qualitatively consistent with their predictions for quenching via radio-mode AGN feedback, hinting at the viability of this process in explaining our observational trends. However, we also find evidence that quenching in L-Galaxies is too efficient and quenching in Illustris is not efficient enough, compared to observations. For satellites, we find strong evidence that environment affects their quenched fraction at fixed central velocity dispersion, particularly at lower masses. At higher masses, satellites behave identically to centrals in their quenching. Of the environmental parameters considered, local density affects the quenched fraction of satellites the most at fixed central velocity dispersion.

Using artificial neural network (ANN) predictions of total infra-red luminosities (LIR), we compare the host galaxy star formation rates (SFRs) of ~21,000 optically selected active galactic nuclei (AGN), 466 low excitation radio galaxies (LERGs) and 721 mid-IR selected AGN. SFR offsets (Delta SFR) relative to a sample of star-forming `main sequence' galaxies (matched in M*, z and local environment) are computed for the AGN hosts. Optically selected AGN exhibit a wide range of Delta SFR, with a distribution skewed to low SFRs and a median Delta SFR = -0.06 dex. The LERGs have SFRs that are shifted to even lower values with a median Delta SFR = -0.5 dex. In contrast, mid-IR selected AGN have, on average, SFRs enhanced by a factor ~1.5. We interpret the different distributions of Delta SFR amongst the different AGN classes in the context of the relative contribution of triggering by galaxy mergers. Whereas the LERGs are predominantly fuelled through low accretion rate secular processes which are not accompanied by enhancements in SFR, mergers, which can simultaneously boost SFRs, most frequently lead to powerful, obscured AGN.

We present a novel technique for ranking the relative importance of galaxy properties in the process of quenching star formation. Specifically, we develop an artificial neural network (ANN) approach for pattern recognition and apply it to a population of over 400,000 central galaxies taken from the Sloan Digital Sky Survey Data Release 7. We utilise a variety of physical galaxy properties for training the pattern recognition algorithm to recognise star forming and passive systems, for a `training set' of $\sim$100,000 galaxies. We then apply the ANN model to a `verification set' of $\sim$100,000 different galaxies, randomly chosen from the remaining sample. The success rate of each parameter singly, and in conjunction with other parameters, is taken as an indication of how important the parameters are to the process(es) of central galaxy quenching. We find that central velocity dispersion, bulge mass and B/T are excellent predictors of the passive state of the system, indicating that properties related to the central mass of the galaxy are most closely linked to the cessation of star formation. Larger scale galaxy properties (total or disk stellar masses), or those linked to environment (halo masses or $\delta_5$) perform significantly less well. Our results are plausibly explained by AGN feedback driving the quenching of central galaxies, although we discuss other possibilities as well.

The total infra-red (IR) luminosity (L_IR) can be used as a robust measure of a galaxy's star formation rate (SFR), even in the presence of an active galactic nucleus (AGN), or when optical emission lines are weak. Unfortunately, existing all sky far-IR surveys, such as the Infra-red Astronomical Satellite (IRAS) and AKARI, are relatively shallow and are biased towards the highest SFR galaxies and lowest redshifts. More sensitive surveys with the Herschel Space Observatory are limited to much smaller areas. In order to construct a large sample of L_IR measurements for galaxies in the nearby universe, we employ artificial neural networks (ANNs), using 1136 galaxies in the Herschel Stripe 82 sample as the training set. The networks are validated using two independent datasets (IRAS and AKARI) and demonstrated to predict the L_IR with a scatter sigma ~ 0.23 dex, and with no systematic offset. Importantly, the ANN performs well for both star-forming galaxies and those with an AGN. A public catalog is presented with our L_IR predictions which can be used to determine SFRs for 331,926 galaxies in the Sloan Digital Sky Survey (SDSS), including ~ 129,000 SFRs for AGN-dominated galaxies for which SDSS SFRs have large uncertainties.

We derive the dependence of the fraction of passive central galaxies on the mass of their supermassive black holes for a sample of over 400,000 SDSS galaxies at z < 0.2. Our large sample contains galaxies in a wide range of environments, with stellar masses 8 < log(M*/Msun) < 12, spanning the entire morphological spectrum from pure disks to spheroids. We derive estimates for the black hole masses from measured central velocity dispersions and bulge masses, using a variety of published empirical relationships. We find a very strong dependence of the passive fraction on black hole mass, which is largely unaffected by the details of the black hole mass estimate. Moreover, the passive fraction relationship with black hole mass remains strong and tight even at fixed values of galaxy stellar mass (M*), dark matter halo mass (Mhalo), and bulge-to-total stellar mass ratio (B/T). Whereas, the passive fraction dependence on M*, Mhalo and B/T is weak at fixed MBH. These observations show that, for central galaxies, MBH is the strongest correlator with the passive fraction, consistent with quenching from AGN feedback.

The artificial neural network (ANN) is a well-established mathematical technique for data prediction, based on the identification of correlations and pattern recognition in input training sets. We present the application of ANNs to predict the emission line luminosities of H$\alpha$ and [NII] $\lambda$6584 in galaxies. These important spectral diagnostics are used for metallicities, active galactic nuclei (AGN) classification and star formation rates, yet are shifted into the infra-red for galaxies above $z \sim 0.5$, or may not be covered in spectra with limited wavelength coverage. The ANN is trained with a large sample of emission line galaxies selected from the Sloan Digital Sky Survey using various combinations of emission lines and stellar mass. The ANN is tested for galaxies dominated by both star formation and AGN; in both cases the H$\alpha$ and [NII] $\lambda$6584 line luminosities can be predicted with a scatter $\sigma <$ 0.1 dex. We also show that the performance of the ANN does not depend significantly on the covering fraction, mass or metallicity of the data. Polynomial functions are derived that allow easy application of the ANN predictions to determine H$\alpha$ and [NII] $\lambda$6584 line luminosities. An ANN calibration for the Balmer decrement (H$\alpha$/H$\beta$) based on line equivalent widths and colours is also presented. The effectiveness of the ANN calibration is demonstrated with an independent dataset (the Galaxy Mass and Assembly Survey). We demonstrate the application of our line luminosities to the determination of gas-phase metallicities and AGN classification. The ANN technique yields a significant improvement in the measurement of metallicities that require [NII] and H$\alpha$ when compared with the function based conversions of Kewley & Ellison. The AGN classification is successful for 86 per cent of SDSS galaxies.

We will introduce a new method to derive the orbital elements of the spectroscopic binary stars. Fitting the adopted radial velocity curves on the measured experimental data, the corresponding radial acceleration data are obtained. The orbital parameters are estimated by the nonlinear least squares of the radial velocity-acceleration curves. This method is implemented for the three double line spectroscopic binary systems RZ CAS, CC Cas and V1130 Tau. Numerical calculations show that the result are in good agreement with the those obtained by others using the method of Lehmann-Filh\'{e}s.

We present the full data set of the VIMOS spectroscopic campaign of the ESO/GOODS program in the CDFS, which complements the FORS2 ESO/GOODS spectroscopic campaign. The GOODS/VIMOS spectroscopic campaign is structured in two separate surveys using two different VIMOS grisms. The VIMOS Low Resolution Blue (LR-Blue) and Medium Resolution (MR) orange grisms have been used to cover different redshift ranges. The LR-Blue campaign is aimed at observing galaxies mainly at 1.8<z<3.5, while the MR campaign mainly aims at galaxies at z<1 and Lyman Break Galaxies (LBGs) at z>3.5. The full GOODS/VIMOS spectroscopic campaign consists of 20 VIMOS masks. This release adds 8 new masks to the previous release (12 masks, Popesso et al. 2009). In total we obtained 5052 spectra, 3634 from the 10 LR-Blue masks and 1418 from the 10 MR masks. A significant fraction of the extracted spectra comes from serendipitously observed sources: ~21% in the LR-Blue and ~16% in the MR masks. We obtained 2242 redshifts in the LR-Blue campaign and 976 in the MR campaign for a total success rate of 62% and 69% respectively, which increases to 66% and 73% if only primary targets are considered. The typical redshift uncertainty is estimated to be ~0.0012 (~255 km/s) for the LR-Blue grism and ~0.00040 (~120 km/s) for the MR grism. By complementing our VIMOS spectroscopic catalog with all existing spectroscopic redshifts publicly available in the CDFS, we compiled a redshift master catalog with 7332 entries, which we used to investigate large scale structures out to z~3.7. We produced stacked spectra of LBGs in a few bins of equivalent width (EW) of the Ly-alpha and found evidence for a lack of bright LBGs with high EW of the Ly-alpha. Finally, we obtained new redshifts for 12 X-ray sources of the CDFS and extended-CDFS.

[abridged] From a sample of Ly-alpha emitters in the GOODS-S field with uncontaminated photometry and optical (red) spectroscopy, we select a spatially compact object at a redshift of 5.563 (Ly-alpha) that shows a second emission line, identified as N IV] 1486 A. The SED is modelled in a way that accounts for both the N IV] line emission and the photometry in a self-consistent way. The photoionization code CLOUDY is used to calculate a range of nebular models as a function of stellar ionizing source temperature, ionization parameter, density and nebular metallicity. We compare the theoretical and observed magnitudes and search for the model parameters that also reproduce the observed N IV] luminosity and equivalent width. A nebular model with a hot blackbody ionizing source of around 100 kK and a nebular metallicity of ~5% of solar is able to fit the observed SED and, in particular, explain the large apparent Balmer break which is inferred from the pure stellar population model fitting conventionally applied to multi-band photometric observations. In our model, an apparent spectral break is produced by strong [O III] 4959, 5007 A emission falling in one of the IR bands (IRAC1 in this case). A lower limit on the total baryonic mass of a model of this type is 3.2 x 10^8 Msun . It is argued that objects with Ly-alpha emission at high redshift that show an apparent Balmer break may have their SED dominated by nebular emission and so could possibly be identified with very young starbursting galaxies rather than massive evolved stellar populations. Detailed studies of these emission nebulae with large telescopes will provide a unique insight into very early chemical evolution.