Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

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.

The current consensus on the formation and evolution of the brightest cluster galaxies is that their stellar mass forms early ($z \gtrsim 4$) in separate galaxies that then eventually assemble the main structure at late times ($z \lesssim 1$). However, advances in observational techniques have led to the discovery of protoclusters out to $z \sim 7$, suggesting that the late-assembly picture may not be fully complete. If these protoclusters assemble rapidly in the early universe, they should form the brightest cluster galaxies much earlier than suspected by the late-assembly picture. Using a combination of observationally constrained hydrodynamical and dark-matter-only simulations, we show that the stellar assembly time of a sub-set of the brightest cluster galaxies occurs at high redshifts ($z > 3$) rather than at low redshifts ($z < 1$), as is commonly thought. We find, using isolated non-cosmological hydrodynamical simulations, that highly overdense protoclusters assemble their stellar mass into brightest cluster galaxies within $\sim 1$ $\mathrm{Gyr}$ of evolution -- producing massive blue elliptical galaxies at high redshifts ($z \gtrsim 1.5$). We argue that there is a downsizing effect on the cluster scale wherein some of the brightest cluster galaxies in the cores of the most-massive clusters assemble earlier than those in lower-mass clusters. In those clusters with $z = 0$ virial mass $\geqslant 5\times10^{14}$ M$_\mathrm{\odot}$, we find that $9.8$% have their cores assembly early, and a higher fraction of $16.4$% in those clusters above $10^{15}$ M$_\mathrm{\odot}$. The James Webb Space Telescope will be able to detect and confirm our prediction in the near future, and we discuss the implications to constraining the value of $\sigma_\mathrm{8}$.

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.

The question of whether galaxy mergers are linked to the triggering of active galactic nuclei (AGN) continues to be a topic of considerable debate. The issue can be broken down into two distinct questions: 1) Can galaxy mergers trigger AGN? 2) Are galaxy mergers the dominant AGN triggering mechanism? A complete picture of the AGN-merger connection requires that both of these questions are addressed with the same dataset. In previous work, we have shown that galaxy mergers selected from the Sloan Digital Sky Survey (SDSS) show an excess of both optically-selected, and mid-IR colour-selected AGN, demonstrating that the answer to the first of the above questions is affirmative. Here, we use the same optical and mid-IR AGN selection to address the second question, by quantifying the frequency of morphological disturbances in low surface brightness r-band images from the Canada France Imaging Survey (CFIS). Only ~30 per cent of optical AGN host galaxies are morphologically disturbed, indicating that recent interactions are not the dominant trigger. However, almost 60 per cent of mid-IR AGN hosts show signs of visual disturbance, indicating that interactions play a more significant role in nuclear feeding. Both mid-IR and optically selected AGN have interacting fractions that are a factor of two greater than a mass and redshift matched non-AGN control sample, an excess that increases with both AGN luminosity and host galaxy stellar mass.

Quantitative characterization of galaxy morphology is vital in enabling comparison of observations to predictions from galaxy formation theory. However, without significant overlap between the observational footprints of deep and shallow galaxy surveys, the extent to which structural measurements for large galaxy samples are robust to image quality (e.g., depth, spatial resolution) cannot be established. Deep images from the Sloan Digital Sky Survey (SDSS) Stripe 82 co-adds provide a unique solution to this problem - offering $1.6-1.8$ magnitudes improvement in depth with respect to SDSS Legacy images. Having similar spatial resolution to Legacy, the co-adds make it possible to examine the sensitivity of parametric morphologies to depth alone. Using the Gim2D surface-brightness decomposition software, we provide public morphology catalogs for 16,908 galaxies in the Stripe 82 $ugriz$ co-adds. Our methods and selection are completely consistent with the Simard et al. (2011) and Mendel et al. (2014) photometric decompositions. We rigorously compare measurements in the deep and shallow images. We find no systematics in total magnitudes and sizes except for faint galaxies in the $u$-band and the brightest galaxies in each band. However, characterization of bulge-to-total fractions is significantly improved in the deep images. Furthermore, statistics used to determine whether single-S\'ersic or two-component (e.g., bulge+disc) models are required become more bimodal in the deep images. Lastly, we show that asymmetries are enhanced in the deep images and that the enhancement is positively correlated with the asymmetries measured in Legacy images.

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 interpretive power of the newest generation of large-volume hydrodynamical simulations of galaxy formation rests upon their ability to reproduce the observed properties of galaxies. In this second paper in a series, we employ bulge+disc decompositions of realistic dust-free galaxy images from the Illustris simulation in a consistent comparison with galaxies from the Sloan Digital Sky Survey (SDSS). Examining the size-luminosity relations of each sample, we find that galaxies in Illustris are roughly twice as large and $0.7$ magnitudes brighter on average than galaxies in the SDSS. The trend of increasing slope and decreasing normalization of size-luminosity as a function of bulge-fraction is qualitatively similar to observations. However, the size-luminosity relations of Illustris galaxies are quantitatively distinguished by higher normalizations and smaller slopes than for real galaxies. We show that this result is linked to a significant deficit of bulge-dominated galaxies in Illustris relative to the SDSS at stellar masses $\log\mathrm{M}_{\star}/\mathrm{M}_{\odot}\lesssim11$. We investigate this deficit by comparing bulge fraction estimates derived from photometry \emph{and} internal kinematics. We show that photometric bulge fractions are systematically lower than the kinematic fractions at low masses, but with increasingly good agreement as the stellar mass increases.

We present an image-based method for comparing the structural properties of galaxies produced in hydrodynamical simulations to real galaxies in the Sloan Digital Sky Survey. The key feature of our work is the introduction of extensive observational realism, such as object crowding, noise and viewing angle, to the synthetic images of simulated galaxies, so that they can be fairly compared to real galaxy catalogs. We apply our methodology to the dust-free synthetic image catalog of galaxies from the Illustris simulation at $z=0$, which are then fit with bulge+disc models to obtain morphological parameters. In this first paper in a series, we detail our methods, quantify observational biases, and present publicly available bulge+disc decomposition catalogs. We find that our bulge+disc decompositions are largely robust to the observational biases that affect decompositions of real galaxies. However, we identify a significant population of galaxies (roughly 30\% of the full sample) in Illustris that are prone to internal segmentation, leading to systematically reduced flux estimates by up to a factor of 6, smaller half-light radii by up to a factor of $\sim$ 2, and generally erroneous bulge-to-total fractions of (B/T)=0.