Search

We present in this paper the discovery, properties, and a catalog of 1165 high redshift $6.5 < z < 18$ galaxies found in deep JWST NIRCam imaging from the GTO PEARLS survey combined with data from JWST public fields. We describe our bespoke homogeneous reduction process and our analysis of these areas including the NEP, CEERS, GLASS, NGDEEP, JADES, and ERO SMACS-0723 fields with over 214 arcmin$^{2}$ imaged to depths of $\sim 30$ mag. We describe our rigorous methods for identifying these galaxies, involving the use of Lyman-break strength, detection significance criteria, visual inspection, and integrated photometric redshifts probability distributions predominately at high redshift. Our sample is a robust and highly pure collection of distant galaxies from which we also remove brown dwarf stars, and calculate completeness and contamination from simulations. We include a summary of the basic properties of these $z > 6.5$ galaxies, including their redshift distributions, UV absolute magnitudes, and star formation rates. Our study of these young galaxies reveals a wide range of stellar population properties as seen in their colors and SED fits which we compare to stellar population models, indicating a range of star formation histories, dust, AGN and/or nebular emission. We find a strong trend exists between stellar mass and $(U-V)$ color, as well as the existence of the `main-sequence' of star formation for galaxies as early as $z \sim 12$. This indicates that stellar mass, or an underlying variable correlating with stellar mass, is driving galaxy formation, in agreement with simulation predictions. We also discover ultra-high redshift candidates at $z > 12$ in our sample and describe their properties. Finally, we note a significant observed excess of galaxies compared to models at $z > 12$, revealing a tension between predictions and our observations., Comment: Submitted to AAS Journals, 37 pages

We measure the broad impact of galaxy structure on galaxy formation by examining the ongoing star formation and integrated star formation history as revealed through the stellar masses of galaxies at $z < 7$ based on JWST CEERS data from the Extended Groth Strip (EGS). Using the morphological catalog of 3965 visually classified JWST galaxies from Ferreira et al. (2023), we investigate the evolution of stars, and when they form, as a function of morphological type as well as galaxies classified as passive and starburst through spectral energy distributions. Although disk galaxies dominate the structures of galaxies at $z < 7$, we find that these disks are in general either `passive', or on the main-sequence of star formation, and do not contain a large population of starburst galaxies. We also find no significant correlation between morphological type and the star formation rate or colours of galaxies at $z < 7$. In fact, we find that the morphologically classified `spheroids' tend to be blue and are not found to be predominately passive systems at $z > 1.5$. We also find that the stellar mass function for disk galaxies does not evolve significantly during this time, whereas other galaxy types, such as the peculiar population, evolve dramatically, declining at lower redshifts. This indicates that massive peculiars are more common at higher redshifts. We further find that up to $z \sim 7$, the specific star formation rate (sSFR) does not vary with visual morphology, but strongly depends on stellar mass and internal galaxy mass density. This demonstrates that at early epochs galaxy assembly is a mass-driven, rather than a morphologically-driven, process. Quenching of star formation is therefore a mass-dominated process throughout the universe's history, likely due to the presence of supermassive black holes., Comment: MNRAS in press, 20 pages

Star formation quenching is one of the key processes that shape the evolution of galaxies. In this study, we investigate the changes in molecular gas and star formation properties as galaxies transit from the star-forming main sequence to the passive regime. Our analysis reveals that as galaxies move away from the main sequence towards the green valley the radial profile of specific star formation rate surface density ($\Sigma_\mathrm{sSFR}$) is suppressed compared with main sequence galaxies out to a galactocentric radius of 1.5 $R_{e}$ ($\sim$ 7 kpc for our sample). By combining radial profiles of gas fraction ($f_\mathrm{gas}$) and star formation efficiency (SFE), we can discern the underlying mechanism that determines $\Sigma_\mathrm{sSFR}$ at different galactocentric radii. Analysis of relative contributions of $f_\mathrm{gas}$ and SFE to $\Sigma_\mathrm{sSFR}$ uncovers a diverse range of quenching modes. Star formation in approximately half of our quenching galaxies is primarily driven by a single mode (i.e. either $f_\mathrm{gas}$ or SFE), or a combination of both. A collective analysis of all galaxies reveals that the reduction in star formation within the central regions ($R$ $<$ 0.5 $R_{e}$) is primarily attributable to a decrease in SFE. Conversely, in the disk regions ($R$ $>$ 0.5 $R_{e}$), both $f_\mathrm{gas}$ and SFE contribute to the suppression of star formation. Our findings suggest that multiple quenching mechanisms may be at play in our sample galaxies, and even within a single galaxy. We also compare our observational outcomes with those from galaxy simulations and discuss the implications of our data., Comment: 27 pages, 8 figures, 1 table. Accepted for publication in The Astrophysical Journal

We present an analysis of the quenching of local observed and simulated galaxies, including an investigation of the dependence of quiescence on both intrinsic and environmental parameters. We apply an advanced machine learning technique utilizing random forest classification to predict when galaxies are star forming or quenched. We perform separate classification analyses for three groups of galaxies: (a) central galaxies; (b) high-mass satellites ($M_{*} > 10^{10.5}{\rm M_{\odot}}$); and (c) low-mass satellites ($M_{*} < 10^{10}{\rm M_{\odot}}$) for three cosmological hydrodynamical simulations (EAGLE, Illustris, and IllustrisTNG), and observational data from the SDSS. The simulation results are unanimous and unambiguous: quiescence in centrals and high-mass satellites is best predicted by intrinsic parameters (specifically central black hole mass), whilst it is best predicted by environmental parameters (specifically halo mass) for low-mass satellites. In observations, we find black hole mass to best predict quiescence for centrals and high mass satellites, exactly as predicted by the simulations. However, local galaxy over-density is found to be most predictive parameter for low-mass satellites. Nonetheless, both simulations and observations do agree that it is environment which quenches low mass satellites. We provide evidence which suggests that the dominance of local over-density in classifying low mass systems may be due to the high uncertainty in halo mass estimation from abundance matching, rather than it being fundamentally a more predictive parameter. Finally, we establish that the qualitative trends with environment predicted in simulations are recoverable in the observation space. This has important implications for future wide-field galaxy surveys., Comment: Accepted to MNRAS; 32 pages; 22 figures

We present the extended ALMA MaNGA QUEnching and STar formation survey, a combination of the original 46 ALMaQUEST galaxies plus new ALMA observations for a further 20 interacting galaxies. Three well-studied scaling relations are fit to the 19,999 star-forming spaxels in the extended sample, namely the resolved Schmidt-Kennicutt (rSK) relation, the resolved star forming main sequence (rSFMS) and the resolved molecular gas main sequence (rMGMS). We additionally investigate the relationship between the dynamical equilibrium pressure (PDE) and star formation rate surface density (Sigma_SFR), which we refer to as the resolved PDE (rPDE) relation. Contrary to previous studies that have focussed on normal star-forming galaxies and found an approximately linear rPDE relation, the presence of more vigourously star-forming galaxies in the extended ALMaQUEST sample reveals a marked turnover in the relation at high pressures. Although the scatter around the linear fit to the rPDE relation is similar to the other three relations, a random forest analysis, which can extract non-linear dependences, finds that PDE is unambiguously more important than either Sigma_H2 or Sigma_star for predicting Sigma_SFR. We compare the observed rPDE relation to the prediction of the pressure-regulated feedback-modulated (PRFM) model of star formation, finding that galaxies residing on the global SFMS do indeed closely follow the rPDE relation predicted by the PRFM theory. However, galaxies above and below the global SFMS show significant deviations from the model. Galaxies with high SFR are instead consistent with models that include other contributions to turbulence in addition to the local star formation feedback., Comment: Accepted for publication in MNRAS

We present an analysis of the quenching of star formation in massive galaxies ($M_* > 10^{9.5} M_\odot$) within the first 0.5 - 3 Gyr of the Universe's history utilizing JWST-CEERS data. We utilize a combination of advanced statistical methods to accurately constrain the intrinsic dependence of quenching in a multi-dimensional and inter-correlated parameter space. Specifically, we apply Random Forest (RF) classification, area statistics, and a partial correlation analysis to the JWST-CEERS data. First, we identify the key testable predictions from two state-of-the-art cosmological simulations (IllustrisTNG & EAGLE). Both simulations predict that quenching should be regulated by supermassive black hole mass in the early Universe. Furthermore, both simulations identify the stellar potential ($\phi_*$) as the optimal proxy for black hole mass in photometric data. In photometric observations, where we have no direct constraints on black hole masses, we find that the stellar potential is the most predictive parameter of massive galaxy quenching at all epochs from $z = 0 - 8$, exactly as predicted by simulations for this sample. The stellar potential outperforms stellar mass, galaxy size, galaxy density, and S\'ersic index as a predictor of quiescence at all epochs probed in JWST-CEERS. Collectively, these results strongly imply a stable quenching mechanism operating throughout cosmic history, which is closely connected to the central gravitational potential in galaxies. This connection is explained in cosmological models via massive black holes forming and growing in deep potential wells, and subsequently quenching galaxies through a mix of ejective and preventative active galactic nucleus (AGN) feedback., Comment: Accepted to ApJ. 32 pages; 10 figures

One of the most important questions in astrophysics is what causes galaxies to stop forming stars. Previous studies have shown a tight link between quiescence and black hole mass. Other studies have revealed that quiescence is also associated with 'starvation', the halting of gas inflows, which results in the remaining gas being used up by star formation and in rapid chemical enrichment. In this work, we find the missing link between these two findings. Using a large sample of galaxies, we uncover the intrinsic dependencies of the stellar metallicity on galaxy properties. In the case of star-forming galaxies, stellar metallicity is primarily driven by stellar mass. However, for passive galaxies, the stellar metallicity is primarily driven by the stellar velocity dispersion. The latter is known to be tightly correlated with black hole mass. This result can be seen as connecting previous studies, where the integrated effect of black hole feedback (i.e. black hole mass, traced by the velocity dispersion) prevents gas inflows, starving the galaxy, which is seen by the rapid increase in the stellar metallicity, and leading to the galaxy becoming passive., Comment: 10 pages, 6 figures, Accepted for publication in MNRAS

We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass ($M_\star = 10^{11}\,{\rm M}_\odot$), and matching the stellar-mass distributions of our samples, we find 7 $\sigma$ evidence of $h_4$ increasing with cosmic time, from a median value of 0.019$\pm$0.002 at z=0.8 to 0.059$\pm$0.004 at z=0.06. Alternatively, we use a physically motivated sample selection, based on the mass distribution of the progenitors of local quiescent galaxies as inferred from numerical simulations; in this case, we find 10 $\sigma$ evidence. This evolution suggests that, over the last 7 Gyr, there has been a gradual decrease in the rotation-to-dispersion ratio and an increase in the radial anisotropy of the stellar velocity distribution, qualitatively consistent with accretion of gas-poor satellites. These findings demonstrate that massive galaxies continue to accrete mass and increase their dispersion support after becoming quiescent., Comment: 19 pages, 9 figures Accepted for publication in MNRAS

We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star > 10^{10.5} {\rm M}_\odot$) at redshifts z = 0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys, to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. We find that at all redshifts where we have large enough samples, $h_4$ anti-correlates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anti-correlation with rotation-to-dispersion ratio, we also find a correlation between $h_4$ and $M_\star$, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between $h_4$ and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, $h_4$ also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower $h_4$ compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity., Comment: 26 pages, 15 figures Accepted for publication in MNRAS

We identify the intrinsic dependence of star formation quenching on a variety of galactic and environmental parameters, utilizing a machine learning approach with Random Forest classification. We have previously demonstrated the power of this technique to isolate causality, not mere correlation, in complex astronomical data. First, we analyze three cosmological hydrodynamical simulations (Eagle, Illustris, and IllustrisTNG), selecting snapshots spanning the bulk of cosmic history from comic noon ($z \sim 2$) to the present epoch, with stellar masses in the range $9 < \log(M_*/M_{\odot}) < 12$. In the simulations, black hole mass is unanimously found to be the most predictive parameter of central galaxy quenching at all epochs. Perhaps surprisingly, black hole accretion rate (and hence the bolometric luminosity of active galactic nuclei, AGN) is found to be of little predictive power over quenching. This theoretical result is important for observational studies of galaxy quenching as it cautions against using the current AGN state of a galaxy as a useful proxy for the cumulative impact of AGN feedback on a galactic system. The latter is traced by black hole mass not AGN luminosity. Additionally, we explore a sub-set of 'observable' parameters, which can be readily measured in extant wide-field galaxy surveys targeting $z = 0 - 2$, at $9 < \log(M_*/M_{\odot}) < 12$. All three simulations predict that in lieu of black hole mass, the stellar gravitational potential will outperform the other parameters in predicting quenching. We confirm this theoretical prediction observationally in the SDSS (at low redshifts) and in CANDELS (at intermediate and high redshifts)., Comment: Accepted to ApJ. 26 pages; 12 figures

We investigate the relationship between the star formation rate (SFR), stellar mass ($M_*$) and molecular gas mass ($M_{H_2}$) for local star-forming galaxies. We further investigate these relationships for high-z (z=1-3) galaxies and for the hosts of a local sample of Active Galactic Nuclei (AGN). We explore which of these dependencies are intrinsic and which are an indirect by-product by employing partial correlation coefficients and random forest regression. We find that for local star-forming galaxies, high-z galaxies, and AGN host galaxies, the Schmidt-Kennicutt relation (SK, between $M_{H_2}$ and SFR), and the Molecular Gas Main Sequence (MGMS, between $M_{H_2}$ and $M_*$) are intrinsic primary relations, while the relationship between $M_*$ and SFR, i.e. the Star-Forming Main Sequence (SFMS), is an indirect by-product of the former two. Hence the Star-Forming Main Sequence is not a fundamental scaling relation for local or high-redshift galaxies. We find evidence for both the evolution of the MGMS and SK relation over cosmic time, where, at a given stellar mass, the higher the redshift, the greater the molecular gas mass and the star formation efficiency. We offer a parameterisation of both the MGMS and SK relation's evolution with redshift, showing how they combine to form the observed evolution of the SFMS. In addition, we find that the local AGN host galaxies follow an AGN-MGMS relation (as well as a AGN-SK relation), where the MGMS is offset to lower $M_{H_2}$ for a given $M_*$ compared to local SF galaxies., Comment: 16 pages, 14 figures, Accepted for publication in MNRAS

Post-starburst galaxies (PSBs) have recently and rapidly quenched their star-formation, thus they are an important way to understand how galaxies transition from star-forming late-types to quiescent early-types. The recent discovery of large cold gas reservoirs in PSBs calls into question the theory that galaxies must lose their gas to become quiescent. Optical Integral Field Spectroscopy (IFS) surveys have revealed two classes of PSBs: central PSBs with central quenching regions and ring PSBs with quenching in their outskirts. We analyze a sample of 13 nearby (z < 0.1) PSBs with spatially resolved optical IFS data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey and matched resolution Atacama Large (sub-)Millimeter Array (ALMA) observations of $^{12}$CO(1-0). Disturbed stellar kinematics in 7/13 of our PSBs and centrally concentrated molecular gas is consistent with a recent merger for most of our sample. In galaxies without merger evidence, alternate processes may funnel gas inwards and suppress star-formation, which may include outflows, stellar bars, and minor mergers or interactions. The star-formation efficiencies of the post-starburst regions in nearly half our galaxies are suppressed while the gas fractions are consistent with star-forming galaxies. AGN feedback may drive this stabilization, and we observe AGN-consistent emission in the centers of 5/13 galaxies. Finally, our central and ring PSBs have similar properties except the ionized and molecular gas in central PSBs is more disturbed. Overall, the molecular gas in our PSBs tends to be compact and highly disturbed, resulting in concentrated gas reservoirs unable to form stars efficiently., Comment: Accepted for publication in The Astrophysical Journal

We study the scaling relations between gas-phase metallicity, stellar mass surface density ($\Sigma _*$), star formation rate surface density ($\Sigma _{SFR}$), and molecular gas surface density ($\Sigma_{H_2}$) in local star-forming galaxies on scales of a kpc. We employ optical integral field spectroscopy from the MaNGA survey, and ALMA data for a subset of MaNGA galaxies. We use Partial Correlation Coefficients and Random Forest regression to determine the relative importance of local and global galactic properties in setting the gas-phase metallicity. We find that the local metallicity depends primarily on $\Sigma _*$ (the resolved mass-metallicity relation, rMZR), and has a secondary anti-correlation with $\Sigma _{SFR}$ (i.e. a spatially-resolved version of the 'Fundamental Metallicity Relation', rFMR). We find that $\Sigma_{H_2}$ is less important than $\Sigma_{SFR}$ in determining the local metallicity. This result indicates that gas accretion, resulting in local metallicity dilution and local boosting of star formation, is unlikely to be the primary origin of the rFMR. The local metallicity depends also on the global properties of galaxies. We find a strong dependence on the total stellar mass ($M_*$) and a weaker (inverse) dependence on the total SFR. The global metallicity scaling relations, therefore, do not simply stem out of their resolved counterparts; global properties and processes, such as the global gravitational potential well, galaxy-scale winds and global redistribution/mixing of metals, likely contribute to the local metallicity, in addition to local production and retention., Comment: 22 pages, 24 figures, Accepted for publication in MNRAS

Galaxy mergers are known to trigger both extended and central star formation. However, what remains to be understood is whether this triggered star formation is facilitated by enhanced star formation efficiencies, or an abundance of molecular gas fuel. This work presents spatially resolved measurements of CO emission collected with the Atacama Large Millimetre Array (ALMA) for 20 merging galaxies (either pairs or post-mergers) selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. Eleven additional merging galaxies are selected from the ALMA MaNGA QUEnching and STar formation (ALMaQUEST) survey, resulting in a set of 31 mergers at various stages of interaction and covering a broad range of star formation rates (SFR). We investigate galaxy-to-galaxy variations in the resolved Kennicutt-Schmidt relation (rKS: $\Sigma_{H_2}$ vs. $\Sigma_{SFR}$), the resolved molecular gas main sequence (rMGMS: $\Sigma_{\star}$ vs. $\Sigma_{H_2}$), and the resolved star-forming main sequence (rSFMS: $\Sigma_{\star}$ vs. $\Sigma_{SFR}$). We quantify offsets from these resolved relations to determine if star formation rate, molecular gas fraction, and/or star formation efficiency (SFE) is enhanced in different regions of an individual galaxy. By comparing offsets in all three parameters we can discern whether gas fraction or SFE powers an enhanced $\Sigma_{SFR}$. We find that merger-induced star formation can be driven by a variety of mechanisms, both within a galaxy and between different mergers, regardless of interaction stage., Comment: 28 pages, 28 figures, accepted for publication in MNRAS

We present an analysis of the quenching of star formation in galaxies, bulges, and disks throughout the bulk of cosmic history, from $z=2-0$. We utilise observations from the SDSS and MaNGA at low redshifts. We complement these data with observations from CANDELS at high redshifts. Additionally, we compare the observations to detailed predictions from the LGalaxies semi-analytic model. To analyse the data, we developed a machine learning approach utilising a Random Forest classifier. We first demonstrate that this technique is extremely effective at extracting causal insight from highly complex and inter-correlated model data, before applying it to various observational surveys. Our primary observational results are as follows: At all redshifts studied in this work, we find bulge mass to be the most predictive parameter of quenching, out of the photometric parameter set (incorporating bulge mass, disk mass, total stellar mass, and $B/T$ structure). Moreover, we also find bulge mass to be the most predictive parameter of quenching in both bulge and disk structures, treated separately. Hence, intrinsic galaxy quenching must be due to a stable mechanism operating over cosmic time, and the same quenching mechanism must be effective in both bulge and disk regions. Despite the success of bulge mass in predicting quenching, we find that central velocity dispersion is even more predictive (when available in spectroscopic data sets). In comparison to the LGalaxies model, we find that all of these observational results may be consistently explained through quenching via preventative `radio-mode' active galactic nucleus (AGN) feedback. Furthermore, many alternative quenching mechanisms (including virial shocks, supernova feedback, and morphological stabilisation) are found to be inconsistent with our observational results and those from the literature., Comment: Accepted to A&A; 40 pages; 20 figures

We investigate the nature of the scaling relations between the surface density of star formation rate ($\Sigma _{\rm SFR}$), stellar mass ($\Sigma _*$), and molecular gas mass ($\Sigma _{\rm H_2}$), aiming at distinguishing between the relations that are primary, i.e. more fundamental, and those which are instead an indirect by-product of the other relations. We use the ALMaQUEST survey and analyse the data by using both partial correlations and Random Forest regression techniques. We unambiguously find that the strongest intrinsic correlation is between $\Sigma _{\rm SFR}$ and $\Sigma_{\rm H_2}$ (i.e. the resolved Schmidt-Kennicutt relation), followed by the correlation between $\Sigma _{\rm H_2}$ and $\Sigma _*$ (resolved Molecular Gas Main Sequence, rMGMS). Once these two correlations are taken into account, we find that there is no evidence for any intrinsic correlation between $\Sigma _{\rm SFR}$ and $\Sigma _*$, implying that SFR is entirely driven by the amount of molecular gas, while its dependence on stellar mass (i.e. the resolved Star Forming Main Sequence, rSFMS) simply emerges as a consequence of the relationship between molecular gas and stellar mass., Comment: Accepted for publication in MNRAS, 7 pages, 4 figures

We develop a 2D inclined rotating disc model, which we apply to the stellar velocity maps of 1862 galaxies taken from the MaNGA survey (SDSS public Data Release 15). We use a random forest classifier to identify the kinematic parameters that are most connected to galaxy quenching. We find that kinematic parameters that relate predominantly to the disc (such as the mean rotational velocity) and parameters that characterise whether a galaxy is rotation- or dispersion-dominated (such as the ratio of rotational velocity to velocity dispersion) are not fundamentally linked to the quenching of star formation. Instead, we find overwhelmingly that it is the absolute level of velocity dispersion (a property that relates primarily to a galaxy's bulge/spheroidal component) that is most important for separating star forming and quenched galaxies. Furthermore, a partial correlation analysis shows that many commonly discussed correlations between galaxy properties and quenching are spurious, and that the fundamental correlation is between quenching and velocity dispersion. In particular, we find that at fixed velocity dispersion, there is only a very weak dependence of quenching on the disc properties, whereby more discy galaxies are slightly more likely to be forming stars. By invoking the tight relationship between black hole mass and velocity dispersion, and noting that black hole mass traces the total energy released by AGN, we argue that these data support a scenario in which quenching occurs by preventive feedback from AGN. The kinematic measurements from this work are publicly available., Comment: Accepted to MNRAS; 29 pages; 17 figures

In this paper we investigate how massive central galaxies cease their star formation by comparing theoretical predictions from cosmological simulations: EAGLE, Illustris and IllustrisTNG with observations of the local Universe from the Sloan Digital Sky Survey (SDSS). Our machine learning (ML) classification reveals supermassive black hole mass ($M_{\rm BH}$) as the most predictive parameter in determining whether a galaxy is star forming or quenched at redshift $z=0$ in all three simulations. This predicted consequence of active galactic nucleus (AGN) quenching is reflected in the observations, where it is true for a range of indirect estimates of $M_{\rm BH}$ via proxies as well as its dynamical measurements. Our partial correlation analysis shows that other galactic parameters lose their strong association with quiescence, once their correlations with $M_{\rm BH}$ are accounted for. In simulations we demonstrate that it is the integrated power output of the AGN, rather than its instantaneous activity, which causes galaxies to quench. Finally, we analyse the change in molecular gas content of galaxies from star forming to passive populations. We find that both gas fractions ($f_{\rm gas}$) and star formation efficiencies (SFEs) decrease upon transition to quiescence in the observations but SFE is more predictive than $f_{\rm gas}$ in the ML passive/star-forming classification. These trends in the SDSS are most closely recovered in IllustrisTNG and are in direct contrast with the predictions made by Illustris. We conclude that a viable AGN feedback prescription can be achieved by a combination of preventative feedback and turbulence injection which together quench star formation in central galaxies., Comment: 29 pages, 16 figures, accepted for publication in MNRAS

Modern spectroscopic surveys of galaxies such as MaNGA consist of millions of diverse spectra covering different regions of thousands of galaxies. We propose and implement a deep unsupervised machine learning method to summarize the entire diversity of MaNGA spectra onto a 15x15 map (DESOM-1), where neighbouring points on the map represent similar spectra. We demonstrate our method as an alternative to conventional full spectral fitting for deriving physical quantities, as well as their full probability distributions, much more efficiently than traditional resource-intensive Bayesian methods. Since spectra are grouped by similarity, the distribution of spectra onto the map for a single galaxy, i.e., its "fingerprint", reveals the presence of distinct stellar populations within the galaxy indicating smoother or episodic star-formation histories. We further map the diversity of galaxy fingerprints onto a second map (DESOM-2). Using galaxy images and independent measures of galaxy morphology, we confirm that galaxies with similar fingerprints have similar morphologies and inclination angles. Since morphological information was not used in the mapping algorithm, relating galaxy morphology to the star-formation histories encoded in the fingerprints is one example of how the DESOM maps can be used to make scientific inferences., Comment: 19 pages, 14 figures. Accepted in The Astronomical Journal

Quantitative morphologies, such as asymmetry and concentration, have long been used as an effective way to assess the distribution of galaxy starlight in large samples. Application of such quantitative indicators to other data products could provide a tool capable of capturing the 2-dimensional distribution of a range of galactic properties, such as stellar mass or star-formation rate maps. In this work, we utilize galaxies from the Illustris and IllustrisTNG simulations to assess the applicability of concentration and asymmetry indicators to the stellar mass distribution in galaxies. Specifically, we test whether the intrinsic values of concentration and asymmetry (measured directly from the simulation stellar mass particle maps) are recovered after the application of measurement uncertainty and a point spread function (PSF). We find that random noise has a non-negligible systematic effect on asymmetry that scales inversely with signal-to-noise, particularly at signal-to-noise less than 100. We evaluate different methods to correct for the noise contribution to asymmetry at very low signal-to-noise, where previous studies have been unable to explore due to systematics. We present algebraic corrections for noise and resolution to recover the intrinsic morphology parameters. Using Illustris as a comparison dataset, we evaluate the robustness of these fits in the presence of a different physics model, and confirm these correction methods can be applied to other datasets. Lastly, we provide estimations for the uncertainty on different correction methods at varying signal-to-noise and resolution regimes., Comment: 19 pages, 18 figures, accepted for publication in MNRAS

In order to investigate the role of gas in the demise of star formation on kpc-scales, we compare the resolved molecular gas main sequence (rMGMS: Sigma_* vs Sigma_H2) of star-forming regions to the sequence of `retired' regions that have ceased to form new stars. Using data from the ALMaQUEST survey, we find that retired spaxels form a rMGMS that is distinct from that of star-forming spaxels, offset to lower Sigma_H2 at fixed Sigma_* by a factor of ~5. We study the rMGMS of star-forming and retired spaxels on a galaxy-by-galaxy basis for eight individual ALMaQUEST galaxies. Six of these galaxies have their retired spaxels concentrated within the central few kpc. Molecular gas is detected in 40-100% of retired spaxels in the eight galaxies in our sample. Both the star-forming and retired rMGMS show a diversity in normalization from galaxy-to-galaxy. However, in any given galaxy, the rMGMS for retired regions is found to be distinct from the star-forming sequence and gas fractions of retired spaxels are up to an order of magnitude lower than the star-forming spaxels. We conclude that quenching is associated with a depletion (but not absence) of molecular gas via a mechanism that typically begins in the centre of the galaxy., Comment: MNRAS, in press

Using a sample of ~15,000 kpc-scale star-forming spaxels in 28 galaxies drawn from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we investigate the galaxy-to-galaxy variation of the `resolved' Schmidt-Kennicutt relation (rSK; Sigma_H2 - Sigma_SFR), the `resolved' star forming main sequence (rSFMS; Sigma_* - Sigma_SFR) and the `resolved' molecular gas main sequence (rMGMS; Sigma_* - Sigma_H2). The rSK relation, rSFMS and rMGMS all show significant galaxy-to-galaxy variation in both shape and normalization, indicating that none of these relations is universal between galaxies. The rSFMS shows the largest galaxy-to-galaxy variation and the rMGMS the least. By defining an `offset' from the average relations, we compute a Delta_rSK, Delta_rSFMS, Delta_rMGMS for each galaxy, to investigate correlations with global properties. We find the following correlations with at least 2 sigma significance: the rSK is lower (i.e. lower star formation efficiency) in galaxies with higher M_*, larger Sersic index and lower specific SFR (sSFR); the rSFMS is lower (i.e. lower sSFR) in galaxies with higher M_* and larger Sersic index; the rMGMS is lower (i.e. lower gas fraction) in galaxies with lower sSFR. In the ensemble of all 15,000 data points, the rSK relation and rMGMS show equally tight scatters and strong correlation coefficients, compared with a larger scatter and weaker correlation in the rSFMS. Moreover, whilst there is no correlation between Delta_rSK and Delta_rMGMS in the sample, the offset of a galaxy's rSFMS does correlate with both of the other two offsets. Our results therefore indicate that the rSK and rMGMS are independent relations, whereas the rSFMS is a result of their combination., Comment: MNRAS, in press

We investigate the spatial structure and evolution of star formation and the interstellar medium (ISM) in interacting galaxies. We use an extensive suite of parsec-scale galaxy merger simulations (stellar mass ratio = 2.5:1), which employs the "Feedback In Realistic Environments-" model (fire-2). This framework resolves star formation, feedback processes, and the multi-phase structure of the ISM. We focus on the galaxy-pair stages of interaction. We find that close encounters substantially augment cool (HI) and cold-dense (H2) gas budgets, elevating the formation of new stars as a result. This enhancement is centrally-concentrated for the secondary galaxy, and more radially extended for the primary. This behaviour is weakly dependent on orbital geometry. We also find that galaxies with elevated global star formation rate (SFR) experience intense nuclear SFR enhancement, driven by high levels of either star formation efficiency (SFE) or available cold-dense gas fuel. Galaxies with suppressed global SFR also contain a nuclear cold-dense gas reservoir, but low SFE levels diminish SFR in the central region. Concretely, in the majority of cases, SFR-enhancement in the central kiloparsec is fuel-driven (55% for the secondary, 71% for the primary) -- whilst central SFR-suppression is efficiency-driven (91% for the secondary, 97% for the primary). Our numerical predictions underscore the need of substantially larger, and/or merger-dedicated, spatially-resolved galaxy surveys -- capable of examining vast and diverse samples of interacting systems -- coupled with multi-wavelength campaigns aimed to capture their internal ISM structure., Comment: 21 pages, 9 figures, MNRAS accepted. For videos: http://research.pomona.edu/galaxymergers/videos/

We investigate how star formation quenching proceeds within central and satellite galaxies using spatially resolved spectroscopy from the SDSS-IV MaNGA DR15. We adopt a complete sample of star formation rate surface densities ($\Sigma_{\rm SFR}$), derived in Bluck et al. (2020), to compute the distance at which each spaxel resides from the resolved star forming main sequence ($\Sigma_{\rm SFR} - \Sigma_*$ relation): $\Delta \Sigma_{\rm SFR}$. We study galaxy radial profiles in $\Delta \Sigma_{\rm SFR}$, and luminosity weighted stellar age (${\rm Age_L}$), split by a variety of intrinsic and environmental parameters. Via several statistical analyses, we establish that the quenching of central galaxies is governed by intrinsic parameters, with central velocity dispersion ($\sigma_c$) being the most important single parameter. High mass satellites quench in a very similar manner to centrals. Conversely, low mass satellite quenching is governed primarily by environmental parameters, with local galaxy over-density ($\delta_5$) being the most important single parameter. Utilising the empirical $M_{BH}$ - $\sigma_c$ relation, we estimate that quenching via AGN feedback must occur at $M_{BH} \geq 10^{6.5-7.5} M_{\odot}$, and is marked by steeply rising $\Delta \Sigma_{\rm SFR}$ radial profiles in the green valley, indicating `inside-out' quenching. On the other hand, environmental quenching occurs at over-densities of 10 - 30 times the average galaxy density at z$\sim$0.1, and is marked by steeply declining $\Delta \Sigma_{\rm SFR}$ profiles, indicating `outside-in' quenching. Finally, through an analysis of stellar metallicities, we conclude that both intrinsic and environmental quenching must incorporate significant starvation of gas supply., Comment: Accepted for publication in MNRAS. 40 pages; 24 figures

Starburst galaxies have elevated star formation rates (SFRs) for their stellar mass. In Ellison et al. (2018) we used integral field unit (IFU) maps of star formation rate surface density (Sigma_SFR) and stellar mass surface density (Sigma_*) to show that starburst galaxies in the local universe are driven by SFRs that are preferentially boosted in their central regions. Here, we present molecular gas maps obtained with the Atacama Large Millimeter Array (ALMA) observatory for 12 central starburst galaxies at z~0 drawn from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The ALMA and MaNGA data are well matched in spatial resolution, such that the ALMA maps of molecular gas surface density (Sigma_H2) can be directly compared with MaNGA maps at kpc-scale resolution. The combination of Sigma_H2, Sigma_* and Sigma_SFR at the same resolution allow us to investigate whether central starbursts are driven primarily by enhancements in star formation efficiency (SFE) or by increased gas fractions. By computing offsets from the resolved Kennicutt-Schmidt relation (Sigma_H2 vs. Sigma_SFR) and the molecular gas main sequence (Sigma_* vs. Sigma_H2), we conclude that the primary driver of the central starburst is an elevated SFE. We also show that the enhancement in Sigma_SFR is accompanied by a dilution in O/H, consistent with a triggering that is induced by metal poor gas inflow. These observational signatures are found in both undisturbed (9/12 galaxies in our sample) and recently merged galaxies, indicating that both interactions and secular mechanisms contribute to central starbursts., Comment: Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. 45 pages; 20 figures

In order to investigate the importance of different proposed quenching mechanisms, we use an indirect method to estimate gas masses for ~62,000 SDSS DR7 galaxies. We infer gas surface densities from dust column densities as traced by extinction within the fibre, applying a metallicity correction to account for varying dust-to-gas ratios. We find that both gas fraction and star formation efficiency (SFE) decrease moving away from the star forming main sequence (MS) towards quiescence for all galaxy masses. We further show that both quantities correlate similarly strongly with the departure from the MS, implying the need for any physical model of quenching to invoke a change in $\textit{both}$ gas fraction and SFE. Our results call for a better understanding of the physical processes driving the decrease in star formation efficiency, which has received relatively little attention in the theory of quenching until now., Comment: 7 pages, 4 figures, Accepted for publication in MNRAS Letters. In order to ensure reproducibility, all analysis is available at https://hub.docker.com/u/jpiotrowska

There are several supervised machine learning methods used for the application of automated morphological classification of galaxies; however, there has not yet been a clear comparison of these different methods using imaging data, or a investigation for maximising their effectiveness. We carry out a comparison between several common machine learning methods for galaxy classification (Convolutional Neural Network (CNN), K-nearest neighbour, Logistic Regression, Support Vector Machine, Random Forest, and Neural Networks) by using Dark Energy Survey (DES) data combined with visual classifications from the Galaxy Zoo 1 project (GZ1). Our goal is to determine the optimal machine learning methods when using imaging data for galaxy classification. We show that CNN is the most successful method of these ten methods in our study. Using a sample of $\sim$2,800 galaxies with visual classification from GZ1, we reach an accuracy of $\sim$0.99 for the morphological classification of Ellipticals and Spirals. The further investigation of the galaxies that have a different ML and visual classification but with high predicted probabilities in our CNN usually reveals an the incorrect classification provided by GZ1. We further find the galaxies having a low probability of being either spirals or ellipticals are visually Lenticulars (S0), demonstrating that supervised learning is able to rediscover that this class of galaxy is distinct from both Es and Spirals. We confirm that $\sim$2.5\% galaxies are misclassified by GZ1 in our study. After correcting these galaxies' labels, we improve our CNN performance to an average accuracy of over 0.99 (accuracy of 0.994 is our best result)., Comment: 20 pages, 19 figures. Accepted to MNRAS

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., Comment: Accepted to MNRAS. 31 pages, 15 figures

We study the correlation of galaxy structural properties with their location relative to the SFR-M* correlation, also known as the star formation "main sequence" (SFMS), in the CANDELS and GAMA surveys and in a semi-analytic model (SAM) of galaxy formation. We first study the distribution of median Sersic index, effective radius, star formation rate (SFR) density and stellar mass density in the SFR-M* plane. We then define a redshift dependent main sequence and examine the medians of these quantities as a function of distance from this main sequence, both above (higher SFRs) and below (lower SFRs). Finally, we examine the distributions of distance from the main sequence in bins of these quantities. We find strong correlations between all of these galaxy structural properties and the distance from the SFMS, such that as we move from galaxies above the SFMS to those below it, we see a nearly monotonic trend towards higher median Sersic index, smaller radius, lower SFR density, and higher stellar density. In the semi-analytic model, bulge growth is driven by mergers and disk instabilities, and is accompanied by the growth of a supermassive black hole which can regulate or quench star formation via Active Galactic Nucleus (AGN) feedback. We find that our model qualitatively reproduces the trends described above, supporting a picture in which black holes and bulges co-evolve, and AGN feedback plays a critical role in moving galaxies off of the SFMS., Comment: 25 pages, 17 figures. Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. 28 pages, 18 figures, 2 tables. This is the final accepted version of an earlier publication: arXiv:1412.3862

We describe a statistical approach for measuring the influence that a galaxy's closest companion has on the galaxy's properties out to arbitrarily wide separations. We begin by identifying the closest companion for every galaxy in a large spectroscopic sample of Sloan Digital Sky Survey galaxies. We then characterize the local environment of each galaxy by using the number of galaxies within 2 Mpc and by determining the isolation of the galaxy pair from other neighbouring galaxies. We introduce a sophisticated algorithm for creating a statistical control sample for each galaxy, matching on stellar mass, redshift, local density and isolation. Unlike traditional studies of close galaxy pairs, this approach is effective in a wide range of environments, regardless of how far away the closest companion is (although a very distant closest companion is unlikely to have a measurable influence on the galaxy in question). We apply this methodology to measurements of galaxy asymmetry, and find that the presence of nearby companions drives a clear enhancement in galaxy asymmetries. The asymmetry excess peaks at the smallest projected separations (< 10 kpc), where the mean asymmetry is enhanced by a factor of 2.0 $\pm$ 0.2. Enhancements in mean asymmetry decline as pair separation increases, but remain statistically significant (1-2$\sigma$) out to projected separations of at least 50 kpc., Comment: Accepted for publication in MNRAS. 17 pages, 14 figures

We present the stellar mass functions (SMF) and mass densities of galaxies, and their spheroid and disk components in the local (z~0.1) universe over the range 8.9 <= log(M/M_solar) <= 12 from spheroid+disk decompositions and corresponding stellar masses of a sample of over 600,000 galaxies in the SDSS-DR7 spectroscopic sample. The galaxy SMF is well represented by a single Schechter function (M* = 11.116+/-0.011, alpha = -1.145+/-0.008), though with a hint of a steeper faint end slope. The corresponding stellar mass densities are (2.670+/-0.110), (1.687+/-0.063) and (0.910+/-0.029)x10^8 M_solar Mpc^-3 for galaxies, spheroids and disks respectively. We identify a crossover stellar mass of log(M/M_solar) = 10.3+/-0.030 at which the spheroid and disk SMFs are equal. Relative contributions of four distinct spheroid/disk dominated sub-populations to the overall galaxy SMF are also presented. The mean disk-to-spheroid stellar mass ratio shows a five fold disk dominance at the low mass end, decreasing monotonically with a corresponding increase in the spheroidal fraction till the two are equal at a galaxy stellar mass, log(M/M_solar)=10.479+/-0.013, the dominance of spheroids then grows with increasing stellar mass. The relative numbers of composite disk and spheroid dominated galaxies show peaks in their distributions, perhaps indicative of a preferred galaxy mass. Our characterization of the low redshift galaxy population provides stringent constraints for numerical simulations to reproduce., Comment: 30 pages, 18 figures, 5 tables (2 online), Accepted for publication in MNRAS

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., Comment: 22 pages, 22 figures, Accepted for publication in MNRAS

We employ a suite of 75 simulations of galaxies in idealised major mergers (stellar mass ratio ~2.5:1), with a wide range of orbital parameters, to investigate the spatial extent of interaction-induced star formation. Although the total star formation in galaxy encounters is generally elevated relative to isolated galaxies, we find that this elevation is a combination of intense enhancements within the central kpc and moderately suppressed activity at large galacto-centric radii. The radial dependence of the star formation enhancement is stronger in the less massive galaxy than in the primary, and is also more pronounced in mergers of more closely aligned disc spin orientations. Conversely, these trends are almost entirely independent of the encounter's impact parameter and orbital eccentricity. Our predictions of the radial dependence of triggered star formation, and specifically the suppression of star formation beyond kph-scales, will be testable with the next generation of integral-field spectroscopic surveys., Comment: 12 pages, 8 figures, accepted by MNRAS

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., Comment: Submitted to MNRAS Letters. 6 pages, 2 figures. Comments welcome

As galaxy formation and evolution over long cosmic time-scales depends to a large degree on the structure of the universe, the assembly history of galaxies is potentially a powerful approach for learning about the universe itself. In this paper we examine the merger history of dark matter halos based on the Extended Press-Schechter formalism as a function of cosmological parameters, redshift and halo mass. We calculate how major halo mergers are influenced by changes in the cosmological values of $\Omega_{\rm m}$, $\Omega_{\Lambda}$, $\sigma_{8}$, the dark matter particle temperature (warm vs. cold dark matter), and the value of a constant and evolving equation of state parameter $w(z)$. We find that the merger fraction at a given halo mass varies by up to a factor of three for halos forming under the assumption of Cold Dark Matter, within different underling cosmological parameters. We find that the current measurements of the merger history, as measured through observed galaxy pairs as well as through structure, are in agreement with the concordance cosmology with the current best fit giving $1 - \Omega_{\rm m} = \Omega_{\rm \Lambda} = 0.84^{+0.16}_{-0.17}$. To obtain a more accurate constraint competitive with recently measured cosmological parameters from Planck and WMAP requires a measured merger accuracy of $\delta f_{\rm m} \sim 0.01$, implying surveys with an accurately measured merger history over 2 - 20 deg$^{2}$, which will be feasible with the next generation of imaging and spectroscopic surveys such as Euclid and LSST., Comment: MNRAS, accepted for publication, 22 pages

We investigate the origin of galaxy bimodality by quantifying the relative role of intrinsic and environmental drivers to the cessation (or `quenching') of star formation in over half a million local Sloan Digital Sky Survey (SDSS) galaxies. Our sample contains a wide variety of galaxies at z=0.02-0.2, with stellar masses of 8 < log(M*/M_sun) < 12, spanning the entire morphological range from pure disks to spheroids, and over four orders of magnitude in local galaxy density and halo mass. We utilise published star formation rates and add to this recent GIM2D photometric and stellar mass bulge + disk decompositions from our group. We find that the passive fraction of galaxies increases steeply with stellar mass, halo mass, and bulge mass, with a less steep dependence on local galaxy density and bulge-to-total stellar mass ratio (B/T). At fixed internal properties, we find that central and satellite galaxies have different passive fraction relationships. For centrals, we conclude that there is less variation in the passive fraction at a fixed bulge mass, than for any other variable, including total stellar mass, halo mass, and B/T. This implies that the quenching mechanism must be most tightly coupled to the bulge. We argue that radio-mode AGN feedback offers the most plausible explanation of the observed trends., Comment: Accepted to MNRAS. 32 pages, 27 figures. [This version is virtually identical to v1]

We use the Millennium Simulation, and an abundance-matching framework, to investigate the dynamical behaviour of galaxy pairs embedded in a cosmological context. Our main galaxy-pair sample, selected to have separations under 250 kpc/h, consists of over 1.3 million pairs at redshift z = 0, with stellar masses greater than 10^9 Msun, probing mass ratios down to 1:1000. We use dark matter halo membership and energy to classify our galaxy pairs. In terms of halo membership, central-satellite pairs tend to be in isolation (in relation to external more massive galaxies), are energetically- bound to each other, and are also weakly-bound to a neighbouring massive galaxy. Satellite-satellite pairs, instead, inhabit regions in close proximity to a more massive galaxy, are energetically-unbound, and are often bound to that neighbour. We find that 60% of our paired galaxies are bound to both their companion and to a third external object. Moreover, only 9% of our pairs resemble the kind of systems described by idealised binary merger simulations in complete isolation. In sum, we demonstrate the importance of properly connecting galaxy pairs to the rest of the Universe., Comment: 25 pages, 14 figures, accepted by MNRAS

In this paper we present a detailed study of the structures and morphologies of a sample of 1188 massive galaxies with Mstar>10^10Msun between redshifts z=1-3 within the Ultra Deep Survey (UDS) region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. Using this sample we determine how galaxy structure and morphology evolve with time. We visually classify our sample into disks, ellipticals and peculiar systems and correct for redshift effects on our classifications through simulations. We find evolution in the fractions of galaxies at a given visual classification as a function of redshift. The peculiar population is dominant at z>2 with a substantial spheroid population, and a negligible disk population. We compute the transition redshift, ztrans, where the combined fraction of spheroids and disks is equal to that of peculiar galaxies, as ztrans=1.86+/-0.62 for galaxies in our stellar mass range. We find that this changes as a function of stellar mass, with Hubble-type systems becoming dominant at higher redshifts for higher mass galaxies (ztrans=2.22+/-0.82), than for the lower mass galaxies (ztrans=1.73+/-0.57). Higher mass galaxies become morphologically settled before their lower mass counterparts, a form of morphological downsizing. We furthermore compare our visual classifications with Sersic index, the concentration, asymmetry and clumpiness (CAS) parameters, star formation rate and rest frame U-B colour. We find links between the colour of a galaxy, its star formation rate and how extended or peculiar it appears. Finally, we discuss the negligible z>2 disk fraction based on visual morphologies and speculate that this is an effect of forming disks appearing peculiar through processes such as violent disk instabilities or mergers. We conclude that to properly define high redshift morphology a new and more exact classification scheme is needed., Comment: 18 pages, 14 figures, Accepted for publication in MNRAS

The ability to resolve all processes which drive galaxy formation is one of the most fundamental goals in extragalactic astronomy. While star formation rates and the merger history are now measured with increasingly high certainty, the role of gas accretion from the intergalactic medium in supplying gas for star formation still remains largely unknown. We present in this paper indirect evidence for the accretion of gas into massive galaxies with initial stellar masses M_*>10^{11} M_sol and following the same merger adjusted co-moving number density at lower redshifts during the epoch 1.5 < z < 3, using results from the GOODS NICMOS Survey (GNS). We show that the measured gas mass fractions of these massive galaxies are inconsistent with the observed star formation history for the same galaxy population. We further demonstrate that this additional gas mass cannot be accounted for by cold gas delivered through minor and major mergers. We also consider the effects of gas outflows and gas recycling due to stellar evolution in these calculations. We argue that to sustain star formation at the observed rates there must be additional methods for increasing the cold gas mass, and that the likeliest method for establishing this supply of gas is by accretion from the intergalactic medium. We calculate that the average gas mass accretion rate into these massive galaxies between 1.5 < z < 3.0, is \dot{M} = 96+/-19 M_sol/yr after accounting for outflowing gas. We show that during this epoch, and for these very massive galaxies, 49+/-20% of baryonic mass assembly is a result of gas accretion and unresolved mergers. However, 66+/-20% of all star formation in this epoch is the result of gas accretion. This reveals that for the most massive galaxies at 1.5< z< 3 gas accretion is the dominant method for instigating new stellar mass assembly., Comment: MNRAS in press, 11 pages, 5 figures

We investigate the total major (> 1:4 by stellar mass) and minor (> 1:100 by stellar mass) merger history of a population of 80 massive (M_* > 10^11 M_sol) galaxies at high redshifts (z = 1.7 - 3). We utilize extremely deep and high resolution HST H-band imaging from the GOODS NICMOS Survey (GNS), which corresponds to rest-frame optical wavelengths at the redshifts probed. We find that massive galaxies at high redshifts are often morphologically disturbed, with a CAS deduced merger fraction f_m = 0.23 +/- 0.05 at z = 1.7 - 3. We find close accord between close pair methods (within 30 kpc apertures) and CAS methods for deducing major merger fractions at all redshifts. We deduce the total (minor + major) merger history of massive galaxies with M_* > 10^9 M_sol galaxies, and find that this scales roughly linearly with log-stellar-mass and magnitude range. We test our close pair methods by utilizing mock galaxy catalogs from the Millennium Simulation. We compute the total number of mergers to be (4.5 +/- 2.9) / <\tau_m> from z = 3 to the present, to a stellar mass sensitivity threshold of ~ 1:100 (where \tau_m is the merger timescale in Gyr which varies as a function of mass). This corresponds to an average mass increase of (3.4 +/- 2.2) x 10^11 M_sol over the past 11.5 Gyrs due to merging. We show that the size evolution observed for these galaxies may be mostly explained by this merging., Comment: 19 pages, 10 figures, re-submitted to ApJ after a positive referee report, originally submitted on Sept 20 2011

We present a study of the star-forming properties of a stellar mass-selected sample of galaxies in the GOODS NICMOS Survey (GNS), based on deep Hubble Space Telescope imaging of the GOODS North and South fields. Using a stellar mass selected sample, combined with HST/ACS and Spitzer data to measure both UV and infrared derived star formation rates (SFR), we investigate the star forming properties of a complete sample of ~1300 galaxies down to log M*=9.5 at redshifts 1.510^11 Msun. We derive optical colours, dust extinctions, and ultraviolet and infrared SFR to determine how the star formation rate changes as a function of both stellar mass and time. Our results show that SFR increases at higher stellar mass such that massive galaxies nearly double their stellar mass from star formation alone over the redshift range studied, but the average value of SFR for a given stellar mass remains constant over this 2 Gyr period. Furthermore, we find no strong evolution in the SFR for our sample as a function of mass over our redshift range of interest, in particular we do not find a decline in the SFR among massive galaxies, as is seen at z < 1. The most massive galaxies in our sample (log M*>11) have high average SFRs with values, SFR(UV,corr) = 103+/-75 Msun/yr, yet exhibit red rest-frame (U-B) colours at all redshifts. We conclude that the majority of these red high-redshift massive galaxies are red due to dust extinction. We find that A(2800) increases with stellar mass, and show that between 45% and 85% of massive galaxies harbour dusty star formation. These results show that even just a few Gyr after the first galaxies appear, there are strong relations between the global physical properties of galaxies, driven by stellar mass or another underlying feature of galaxies strongly related to the stellar mass., Comment: 18 pages, 10 figures, accepted for publication in MNRAS

We use a sample of 8298 galaxies observed in the HST GOODS NICMOS Survey (GNS) to construct the galaxy stellar mass function as a function of both redshift and stellar mass up to z=3.5 and down to masses of Mstar=10^8.5 Msun at z~1. We discover that a significant fraction of all massive Mstar>10^11 Msun galaxies are in place up to the highest redshifts we probe, with a decreasing fraction of lower mass galaxies present at all redshifts. This is an example of `galaxy mass downsizing', and is the result of massive galaxies forming before lower mass ones, and not just simply ending their star formation earlier as in traditional downsizing scenarios. We find that the faint end slope is significantly steeper than what is found in previous investigations. We demonstrate that this steeper mass function better matches the stellar mass added due to star formation, thereby alleviating some of the mismatch between these two measures of the evolution of galaxy mass. We furthermore examine the stellar mass function divided into blue/red systems, as well as for star forming and non-star forming galaxies. We find a similar mass downsizing present for both blue/red and star-forming/non-star forming galaxies, and that the low mass galaxies are mostly all blue, and are therefore creating the steep mass functions. We furthermore show that, although there is a downsizing such that high mass galaxies are nearer their z=0 values at high redshift, this turns over at masses Mstar~10^10 Msun, such that the lowest mass galaxies are more common than galaxies at slight higher masses, creating a `dip' in the observed galaxy mass function. We argue that the galaxy assembly process may be driven by different mechanisms at low and high masses, and that the efficiency of the galaxy formation process is lowest at masses Mstar~10^10 Msun at 1

We present a study of the relationship between galaxy colour, stellar mass, and local galaxy density in a deep near-infrared imaging survey up to a redshift of z~3 using the GOODS NICMOS Survey (GNS). The GNS is a very deep, near-infrared Hubble Space Telescope survey imaging a total of 45 arcmin^2 in the GOODS fields, reaching a stellar mass completeness limit of M* = 10^9.5 M_sun at z=3. Using this data we measure galaxy local densities based on galaxy counts within a fixed aperture, as well as the distance to the 3rd, 5th and 7th nearest neighbour. We find a strong correlation between colour and stellar mass at all redshifts up to z~3. We do not find a strong correlation between colour and local density, however, the highest overdensities might be populated by a higher fraction of blue galaxies than average or underdense areas, indicating a possible reversal of the colour-density relation at high redshift. Our data suggests that the possible higher blue fraction at extreme overdensities might be due to a lack of massive red galaxies at the highest local densities., Comment: 5 pages, 3 figures, to appear in the proceedings of JENAM 2010, Symposium 2 "Environment and the Formation of Galaxies: 30 years later"