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.

This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to stellar population synthesis models. This requires spectra with high signal-to-noise which we achieve by stacking in bins of luminosity (L$_r$) and central velocity dispersion ($\sigma_0$). We find that: 1) There is a bimodality in the S0 population: S0s more massive than $3\times 10^{10}M_\odot$ show stronger velocity dispersion and age gradients (age and $\sigma_r$ decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above $2\times10^{11}M_\odot$ the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. 2) S0s have steeper velocity dispersion profiles than fast rotating elliptical galaxies (E-FRs) of the same luminosity and velocity dispersion. The kinematic profiles and stellar population gradients of E-FRs are both more similar to those of slow rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. 3) At fixed $\sigma_0$, more luminous S0s and E-FRs are younger, more metal rich and less $\alpha$-enhanced. Evidently for these galaxies, the usual statement that 'massive galaxies are older' is not true if $\sigma_0$ is held fixed.

Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $\sigma_e$. Since $I_e$ and $\sigma_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $\sigma_e$ can be used to estimate sizes $R_e$. We show that: 1) The thickness of the FP depends strongly on morphology. If the sample only includes E-SRs, then the observed scatter in $R_e$ is $\sim 16\%$, of which only $\sim 9\%$ is intrinsic. Removing galaxies with $M_*<10^{11}M_\odot$ further reduces the observed scatter to $\sim 13\%$ ($\sim 4\%$ intrinsic). The observed scatter increases to the $\sim 25\%$ usually quoted in the literature if E-FRs and S0s are added. If the FP is defined using the eigenvectors of the covariance matrix of the observables, then the E-SRs again define an exceptionally thin FP, with intrinsic scatter of only $5\%$ orthogonal to the plane. 2) The structure within the FP is most easily understood as arising from the fact that $I_e$ and $\sigma_e$ are nearly independent, whereas the $R_e-I_e$ and $R_e-\sigma_e$ correlations are nearly equal and opposite. 3) If the coefficients of the FP differ from those associated with the virial theorem the plane is said to be `tilted'. If we multiply $I_e$ by the global stellar mass-to-light ratio $M_*/L$ and we account for non-homology across the population by using S\'ersic photometry, then the resulting stellar mass FP is less tilted. Accounting self-consistently for $M_*/L$ gradients will change the tilt. The tilt we currently see suggests that the efficiency of turning baryons into stars increases and/or the dark matter fraction decreases as stellar surface brightness increases.

With the advent of future big-data surveys, automated tools for unsupervised discovery are becoming ever more necessary. In this work, we explore the ability of deep generative networks for detecting outliers in astronomical imaging datasets. The main advantage of such generative models is that they are able to learn complex representations directly from the pixel space. Therefore, these methods enable us to look for subtle morphological deviations which are typically missed by more traditional moment-based approaches. We use a generative model to learn a representation of expected data defined by the training set and then look for deviations from the learned representation by looking for the best reconstruction of a given object. In this first proof-of-concept work, we apply our method to two different test cases. We first show that from a set of simulated galaxies, we are able to detect $\sim90\%$ of merging galaxies if we train our network only with a sample of isolated ones. We then explore how the presented approach can be used to compare observations and hydrodynamic simulations by identifying observed galaxies not well represented in the models.

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

We present new measurements of the evolution of the galaxy stellar mass functions (GSMF) and UV luminosity functions (UV LF) for galaxies from $z=6-9$ within the Frontier Field cluster MACSJ0416.1-2403 and its parallel field. To obtain these results, we derive the stellar masses of our sample by fitting synthetic stellar population models to their observed spectral energy distribution with the inclusion of nebular emission lines. This is the deepest and farthest in distance mass function measured to date and probes down to a level of M$_{*} = 10^{6.8}M_{\odot}$. The main result of this study is that the low-mass end of our GSMF to these limits and redshifts appears to become steeper from $-1.98_{-0.07}^{+0.07}$ at $z=6$ to $-2.38_{-0.88}^{+0.72}$ at $z=9$, steeper than previously observed mass functions at slightly lower redshifts, and we find no evidence of turnover in the mass range probed. We furthermore demonstrate that the UV LF for these system also appears to show a steepening at the highest redshifts, without any evidence of turnover in the luminosity range probed. Our $M_{\mathrm{UV}}-M_{*}$ relation exhibit shallower slopes than previously observed and are in accordance with a constant mass-to-light ratio. Integrating our GSMF, we find that the stellar mass density increases by a factor of $\sim15_{-6}^{+21}$ from $z=9$ to $z=6$. We estimate the dust-corrected star formation rates (SFRs) to calculate the specific star formation rates ($\mathrm{sSFR}=\mathrm{SFR/M_{*}}$) of our sample, and find that for a fixed stellar mass of $5\times10^{9}M_{\odot}$, sSFR $\propto(1+z)^{2.01\pm0.16}$. Finally, from our new measurements, we estimate the UV luminosity density ($\rho_{\textrm{UV}}$) and find that our results support a smooth decline of $\rho_{\textrm{UV}}$ towards high redshifts.

We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

We present a multi-component structural analysis of the internal structure of $1074$ high redshift massive galaxies at $1<z<3$ from the CANDELS HST Survey. In particular we examine galaxies best-fit by two structural components, and thus likely forming discs and bulges. We examine the stellar mass, star formation rates, and colours of both the inner `bulge' and outer `disc' components for these systems using SED information from the resolved ACS+WFC3 HST imaging. We find that the majority of both inner and outer components lie in the star-forming region of UVJ space ($68$ and $90$ per cent respectively). However, the inner portions, or the likely forming bulges, are dominated by dusty star formation. Furthermore, we show that the outer components of these systems have a higher star formation rate than their inner regions, and the ratio of star formation rate between `disc' and `bulge' increases at lower redshifts. Despite the higher star formation rate of the outer component, the stellar mass ratio of inner to outer component remains constant through this epoch. This suggests that there is mass transfer from the outer to inner components for typical two component forming systems, thus building bulges from disks. Finally, using Chandra data we find that the presence of an AGN is more common in both $1$-component spheroid-like galaxies and $2$-component systems ($13\pm3$ and $11\pm2$ per cent) than in $1$-component disc-like galaxies ($3\pm1$ per cent), demonstrating that the formation of a central inner-component likely triggers the formation of central massive black holes in these galaxies.

We perform automated bulge + disc decomposition on a sample of $\sim$7500 galaxies from the Galaxy And Mass Assembly (GAMA) survey in the redshift range of 0.002$<$z$<$0.06 using SIGMA, a wrapper around GALFIT3. To achieve robust profile measurements we use a novel approach of repeatedly fitting the galaxies, varying the input parameters to sample a large fraction of the input parameter space. Using this method we reduce the catastrophic failure rate significantly and verify the confidence in the fit independently of $\chi^2$. Additionally, using the median of the final fitting values and the 16$^{th}$ and 84$^{th}$ percentile produces more realistic error estimates than those provided by GALFIT, which are known to be underestimated. We use the results of our decompositions to analyse the stellar mass - half-light radius relations of bulges, discs and spheroids. We further investigate the association of components with a parent disc or elliptical relation to provide definite $z=0$ disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ relations. We conclude by comparing our local disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ to simulated data from EAGLE and high redshift data from CANDELS-UDS. We show the potential of using the mass-size relation to study galaxy evolution in both cases but caution that for a fair comparison all data sets need to be processed and analysed in the same manner.

We examine a sample of 1495 galaxies in the CANDELS fields to determine the evolution of two component galaxies, including bulges and discs, within massive galaxies at the epoch 1 < z < 3 when the Hubble sequence forms. We fit all of our galaxies' light profiles with a single S\'ersic fit, as well as with a combination of exponential and S\'ersic profiles. The latter is done in order to describe a galaxy with an inner and an outer component, or bulge and disc component. We develop and use three classification methods (visual, F-test and the RFF) to separate our sample into 1-component galaxies (disc/spheroids-like galaxies) and 2-component galaxies (galaxies formed by an 'inner part' or bulge and an 'outer part' or disc). We then compare the results from using these three different ways to classify our galaxies. We find that the fraction of galaxies selected as 2-component galaxies increases on average 50 per cent from the lowest mass bin to the most massive galaxies, and decreases with redshift by a factor of four from z = 1 to z = 3. We find that single S\'ersic 'disc-like' galaxies have the highest relative number densities at all redshifts, and that 2-component galaxies have the greatest increase and become at par with S\'ersic discs by z = 1. We also find that the systems we classify as 2-component galaxies have an increase in the sizes of their outer components, or 'discs', by about a factor of three from z = 3 to z = 1.5, while the inner components or 'bulges' stay roughly the same size. This suggests that these systems are growing from the inside out, whilst the bulges or protobulges are in place early in the history of these galaxies. This is also seen to a lesser degree in the growth of single 'disc-like' galaxies vs. 'spheroid-like' galaxies over the same epoch.

The particle and event horizons are widely known and studied concepts, but the study of their properties, in particular their evolution, have only been done so far considering a single state equation in a decelerating universe. This paper is the first of two where we study this problem from a general point of view. Specifically, this paper is devoted to the study of the evolution of these cosmological horizons in an accelerated universe with two state equations, cosmological constant and dust. We have obtained closed-form expressions for the horizons, which have allowed us to compute their velocities in terms of their respective recession velocities that generalize the previous results for one state equation only. With the equations of state considered, it is proved that both velocities remain always positive.

This paper is the second of two papers devoted to the study of the evolution of the cosmological horizons (particle and event horizons). Specifically, in this paper we consider the extremely general case of an accelerated universe with countably infinitely many constant state equations, and we obtain simple expressions in terms of their respective recession velocities that generalize the previous results for one and two state equations. We also provide a qualitative study of the values of the horizons and their velocities at the origin of the universe and at the far future, and we prove that these values only depend on one dominant state equation. Finally, we compare both horizons and determine when one is larger that the other.