Search

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 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

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., 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

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)., 20 pages, 19 figures. Accepted to MNRAS

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/

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., MNRAS, in press

Cheng, Ting-Yun, et al., We present in this paper one of the largest galaxy morphological classification catalogues to date, including over 20 million galaxies, using the Dark Energy Survey (DES) Year 3 data based on convolutional neural networks (CNNs). Monochromatic i-band DES images with linear, logarithmic, and gradient scales, matched with debiased visual classifications from the Galaxy Zoo 1 (GZ1) catalogue, are used to train our CNN models. With a training set including bright galaxies (16 ≤ i < 18) at low redshift (z < 0.25), we furthermore investigate the limit of the accuracy of our predictions applied to galaxies at fainter magnitude and at higher redshifts. Our final catalogue covers magnitudes 16 ≤ i < 21, and redshifts z < 1.0, and provides predicted probabilities to two galaxy types – ellipticals and spirals (disc galaxies). Our CNN classifications reveal an accuracy of over 99 per cent for bright galaxies when comparing with the GZ1 classifications (i < 18). For fainter galaxies, the visual classification carried out by three of the co-authors shows that the CNN classifier correctly categorizes discy galaxies with rounder and blurred features, which humans often incorrectly visually classify as ellipticals. As a part of the validation, we carry out one of the largest examinations of non-parametric methods, including ∼100 ,000 galaxies with the same coverage of magnitude and redshift as the training set from our catalogue. We find that the Gini coefficient is the best single parameter discriminator between ellipticals and spirals for this data set., The DES data management system is supported by the National Science Foundation under grant numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015- 71825, ESP2015-66861, FPA2015-68048, SEV-2016-0588, SEV2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA programme of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020, and the Brazilian Instituto Nacional de Ciencia e Tecnologia (INCT) e-Universe (CNPq grant 465376/2014-2).

We demonstrate that highly accurate joint redshift-stellar mass probability distribution functions (PDFs) can be obtained using the Random Forest (RF) machine learning (ML) algorithm, even with few photometric bands available. As an example, we use the Dark Energy Survey (DES), combined with the COSMOS2015 catalogue for redshifts and stellar masses. We build two ML models: one containing deep photometry in the $griz$ bands, and the second reflecting the photometric scatter present in the main DES survey, with carefully constructed representative training data in each case. We validate our joint PDFs for $10,699$ test galaxies by utilizing the copula probability integral transform and the Kendall distribution function, and their univariate counterparts to validate the marginals. Benchmarked against a basic set-up of the template-fitting code BAGPIPES, our ML-based method outperforms template fitting on all of our predefined performance metrics. In addition to accuracy, the RF is extremely fast, able to compute joint PDFs for a million galaxies in just under $6$ min with consumer computer hardware. Such speed enables PDFs to be derived in real time within analysis codes, solving potential storage issues. As part of this work we have developed GALPRO, a highly intuitive and efficient Python package to rapidly generate multivariate PDFs on-the-fly. GALPRO is documented and available for researchers to use in their cosmology and galaxy evolution studies., Comment: 18 pages, 8 figures, Accepted by MNRAS

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., MNRAS, in press

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., 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., 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., ERC Advanced Grant: 695671 "Quench"

Monthly Notices of the Royal Astronomical Society, 465 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We present a comprehensive suite of high-resolution (parsec-scale), idealised (non-cosmological) galaxy merger simulations (24 runs, stellar mass ratio ~2.5:1) to investigate the connection between interaction-induced star formation and the evolution of the interstellar medium (ISM) in various temperature-density regimes. We use the GIZMO code and the second version of the 'Feedback in Realistic Environments' model (FIRE-2), which captures the multi-phase structure of the ISM. Our simulations are designed to represent galaxy mergers in the local Universe. In this work, we focus on the 'galaxy-pair period' between first and second pericentric passage. We split the ISM into four regimes: hot, warm, cool and cold-dense, motivated by the hot, ionised, atomic and molecular gas phases observed in real galaxies. We find that, on average, interactions enhance the star formation rate of the pair (~30%, merger-suite sample average) and elevate their cold-dense gas content (~18%). This is accompanied by a decrease in warm gas (~11%), a negligible change in cool gas (~4% increase), and a substantial increase in hot gas (~400%). The amount of cold-dense gas with densities above 1000 cm^3 (the cold ultra-dense regime) is elevated significantly (~240%), but only accounts for 0.15% (on average) of the cold-dense gas budget., Comment: 21 pages, 17 figures, accepted by MNRAS

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

We investigate the dependence of galaxy structure on a variety of galactic and environmental parameters for ~500,000 galaxies at z, Accepted to MNRAS. 31 pages, 15 figures

We use the Henriques et al. (2015) version of the Munich galaxy formation model (L-GALAXIES) to investigate why the halo and stellar mass scales above which galaxies are quenched are constant with redshift and coincide with the scale where baryons are most efficiently converted into stars. This model assumes that central galaxies are quenched by AGN feedback when hot halo gas accretes onto a supermassive black hole. Nevertheless, we find that supernova (SN) feedback sets both mass scales. As haloes grow above a threshold mass, SNe can no longer eject material so their hot gas content increases, enhancing the cooling rate onto the central galaxy, its cold gas content, its star formation rate and the growth rate of its central black hole. Strong AGN feedback terminates this short-lived phase by suppressing the fuel supply for star formation. Despite strong evolution of the halo mass - temperature relation, quenching occurs at a redshift-independent halo and stellar mass which coincides with the mass where baryons have been converted into stars with maximal efficiency. These regularities and coincidences are a result of the specific parameters selected by MCMC tuning of the model to fit the observed abundance and passive fraction of galaxies over the redshift range 0, Comment: 12 pages, 8 figures, submitted to 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 present a structural and morphological catalogue for 45 million objects selected from the first year of data from the Dark Energy Survey (DES). Single Sersic fits and non-parametric measurements are produced for g, r and i filters. The parameters from the best-fitting Sersic model (total magnitude, half-light radius, Sersic index, axis ratio and position angle) are measured with Galfit; the non-parametric coefficients (concentration, asymmetry, clumpiness, Gini, M20) are provided using the Zurich Estimator of Structural Types (ZEST+). To study the statistical uncertainties, we consider a sample of state-of-the-art image simulations with a realistic distribution in the input parameter space and then process and analyse them as we do with real data: this enables us to quantify the observational biases due to PSF blurring and magnitude effects and correct the measurements as a function of magnitude, galaxy size, Sersic index (concentration for the analysis of the non-parametric measurements) and ellipticity. We present the largest structural catalogue to date: we find that accurate and complete measurements for all the structural parameters are typically obtained for galaxies with SExtractor MAG AUTO I < 21. Indeed, the parameters in the filters i and r can be overall well recovered up to MAG AUTO < 21.5, corresponding to a fitting completeness of ~90% below this threshold, for a total of 25 million galaxies. The combination of parametric and non-parametric structural measurements makes this catalogue an important instrument to explore and understand how galaxies form and evolve. The catalogue described in this paper will be publicly released alongside the Dark Energy Survey collaboration Y1 cosmology data products at the following URL: https://des.ncsa.illinois.edu/releases/y1a1/gold/morphology., Accepted for publication in MNRAS. 26 pages, 16 figures. Catalogue data are available at https://des.ncsa.illinois.edu/releases/y1a1/gold/morphology

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

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 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 investigate the relation between star formation rates and local galaxy environment for a stellar-mass-selected galaxy sample in the redshift range 1.5 ≤z≤ 3. We use near-infrared imaging from an extremely deep Hubble Space Telescope survey, the GOODS NICMOS Survey (GNS), to measure local galaxy densities based on the nearest neighbour approach, while star formation rates are estimated from rest-frame ultraviolet fluxes. Due to our imaging depth, we can examine galaxies down to a colour-independent stellar mass completeness limit of log M*= 9.5 M⊙ at z∼ 3. We find a strong dependence of star formation activity on galaxy stellar mass over the whole redshift range, which does not depend on local environment. The average star formation rates are largely independent of local environment, apart from those in the highest relative overdensities. Galaxies in overdensities of a factor of >5 have on average lower star formation rates by a factor of 2–3, but only up to redshifts of z∼ 2. We do not see any evidence for active galactic nucleus activity influencing these relations. We also investigate the influence of the very local environment on star formation activity by counting neighbours within 30 kpc radius. This shows that galaxies with two or more close neighbours have on average significantly lower star formation rates as well as lower specific star formation rates up to z∼ 2.5. We suggest that this might be due to star formation quenching induced by galaxy merging processes.

We examine in this paper a stellar mass selected sample of galaxies at 1 1. By using a locally defined definition of disk and elliptical galaxies based on structure and spectral-type, we find no true ellipticals at z > 2, and a fraction of 3.2+/-2.3% at 1.5 10^{10} M_sol, consistent with previous results from ACS and NICMOS. We compare our results to semi-analytical model results and find a relatively good agreement between our morphological break-down and the predictions. Finally, we argue that the peculiars, ellipticals and peculiar ellipticals have similar properties, suggesting similar formation modes, likely driven by major mergers.

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 the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

[Abridged] To investigate the evolution in the relation between galaxy stellar and central black hole mass we construct a volume limited complete sample of 85 AGN with host galaxy stellar masses M_{*} > 10^{10.5} M_{sol}, and specific X-ray luminosities L_{X} > 2.35 x 10^{43} erg s^{-1} at 0.4 2.35 x 10^{43} erg s^{-1}) rises with redshift from 1.2 +/- 0.2 % at z = 0.7 to 7.4 +/- 2.0 % at z = 2.5. We use our results to calculate the maximum timescales for which our sample of AGN can continue to accrete at their observed rates before surpassing the local galaxy-black hole mass relation. We use these timescales to calculate the total fraction of massive galaxies which will be active (with L_{X} > 2.35 x 10^{43} erg s^{-1}) since z = 3, finding that at least ~ 40% of all massive galaxies will be Seyfert luminosity AGN or brighter during this epoch. Further, we calculate the energy density due to AGN activity in the Universe as 1.0 (+/- 0.3) x 10^{57} erg Mpc^{-3} Gyr^{-1}, potentially providing a significant source of energy for AGN feedback on star formation. We also use this method to compute the evolution in the X-ray luminosity density of AGN with redshift, finding that massive galaxy Seyfert luminosity AGN are the dominant source of X-ray emission in the Universe at z < 3.

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, 30 pages, 18 figures, 5 tables (2 online), 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., 12 pages, 8 figures, accepted by MNRAS

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., Accepted to MNRAS. 32 pages, 27 figures. [This version is virtually identical to v1]

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 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

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., 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 * > 1011 M ⊙) galaxies at high redshifts (z = 1.7-3). We utilize extremely deep and high-resolution Hubble Space Telescope H-band imaging from the GOODS NICMOS Survey, 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 (concentration, C; asymmetry, A; clumpiness, S) deduced merger fraction fm = 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 * > 109 M ⊙ 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)/langτ m rang from z = 3 to the present, to a stellar mass sensitivity threshold of ~1:100 (where τ 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) × 1011 M ⊙ over the past 11.5 Gyr due to merging. We show that the size evolution observed for these galaxies may be mostly explained by this merging., The Astrophysical Journal, 747 (1), ISSN:0004-637X, ISSN:2041-8213

We present our results on the structure and activity of massive galaxies at z=1-3 using one of the largest (166 with M_star>=5e10 M_sun) and most diverse samples of massive galaxies derived from the GOODS-NICMOS survey: (1) Sersic fits to deep NIC3/F160W images indicate that the rest-frame optical structures of massive galaxies are very different at z=2-3 compared to z~0. Approximately 40% of massive galaxies are ultra-compact (r_e2. The large fraction of n2. In order for massive galaxies at z=2-3 to evolve into present-day massive E/S0s, they need to significantly increase (n, r_e). Dry minor and major mergers may play an important role in this process., 49 pages, 22 figures; accepted by ApJ, minor changes from version v1

Utilizing deep Hubble Space Telescope imaging from the two largest field galaxy surveys, the Extended Groth Strip (EGS) and the COSMOS survey, we examine the structural properties, and derive the merger history for 21,902 galaxies with M_*>10^{10} M_0 at z10^{10} M_0, the merger fraction can be parameterised by f_m = f_0*(1+z)^m with the power-law slope m=2.3+/-0.4. By using the best available z = 0 prior the slope increases to m=3.8+/-0.2, showing how critical the measurement of local merger properties are for deriving the evolution of the merger fraction. We furthermore show that the merger fraction derived through structure is roughly a factor of 3-6 higher than pair fractions. Based on the latest cosmological simulations of mergers we show that this ratio is predicted, and that both methods are likely tracing the merger fraction and rate properly. We calculate, utilising merger time scales from simulations, and previously published merger fractions that the merger rate of galaxies with M_*>10^{10} M_0 increases linearly between z = 0.7 and z = 3, and that a typical M_*>10^{10} M_0 galaxy undergoes between 1-2 major mergers at z, Comment: MNRAS, in press, 20 pages

We calculate the major pair fraction and derive the major merger fraction and rate for 82 massive ($M_{*}>10^{11}M_{\odot}$) galaxies at $1.7 < z < 3.0$ utilising deep HST NICMOS data taken in the GOODS North and South fields. For the first time, our NICMOS data provides imaging with sufficient angular resolution and depth to collate a sufficiently large sample of massive galaxies at z $>$ 1.5 to reliably measure their pair fraction history. We find strong evidence that the pair fraction of massive galaxies evolves with redshift. We calculate a pair fraction of $f_{m}$ = 0.29 +/- 0.06 for our whole sample at $1.7 < z < 3.0$. Specifically, we fit a power law function of the form $f_{m}=f_{0}(1+z)^{m}$ to a combined sample of low redshift data from Conselice et al. (2007) and recently acquired high redshift data from the GOODS NICMOS Survey. We find a best fit to the free parameters of $f_{0}$ = 0.008 +/- 0.003 and $m$ = 3.0 +/- 0.4. We go on to fit a theoretically motivated Press-Schechter curve to this data. This Press-Schechter fit, and the data, show no sign of levelling off or turning over, implying that the merger fraction of massive galaxies continues to rise with redshift out to z $\sim$ 3. Since previous work has established that the merger fraction for lower mass galaxies turns over at z $\sim$ 1.5 - 2.0, this is evidence that higher mass galaxies experience more mergers earlier than their lower mass counterparts, i.e. a galaxy assembly downsizing. Finally, we calculate a merger rate at z = 2.6 of $\Re$ $10^{11}M_{\odot}$ galaxy experiencing 1.7 +/- 0.5 mergers between z = 3 and z = 0., Comment: 5 pages, 3 figures, accepted to MNRAS

Understanding the origin and diversity of emission processes responsible for gamma-ray bursts (GRBs) remains a pressing challenge. While prompt and contemporaneous panchromatic observations have the potential to test predictions of the internal-external shock model, extensive multiband imaging has been conducted for only a few GRBs. We present rich, early-time, multiband data sets for two Swift events, GRB 110205A and GRB 110213A. The former shows optical emission since the early stages of the prompt phase, followed by the steep rising in flux up to 1000s after the burst (t -α with α = -6.13 ± 0.75). We discuss this feature in the context of the reverse-shock scenario and interpret the following single power-law decay as being forward-shock dominated. Polarization measurements, obtained with the RINGO2 instrument mounted on the Liverpool Telescope, also provide hints on the nature of the emitting ejecta. The latter event, instead, displays a very peculiar optical to near-infrared light curve, with two achromatic peaks. In this case, while the first peak is probably due to the onset of the afterglow, we interpret the second peak to be produced by newly injected material, signifying a late-time activity of the central engine. © 2011. The American Astronomical Society. All rights reserved.