Researcher profile

Andrew M. Hopkins

Andrew M. Hopkins contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
60works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2026arXiv

The Secret Lives of Open Clusters: a Multiwavelength Examination of Three Open Clusters

Star clusters are well known for their dynamical interactions, an outcome of their high stellar densities; in this paper we use multiwavelength observations to search for the unique outcomes of these interactions in three nearby Galactic open clusters: IC 2602 (30 Myr), NGC 2632 (750 Myr) and M67 (4 Gyr). We compared X-ray observations from all-sky surveys like eROSITA, plus archival observations from Chandra X-ray Observatory, survey radio observations from ASKAP's Evolutionary Map of the Universe survey plus archival VLA observations, in conjunction with new cluster catalogs with Gaia. From X-ray, we found 77 X-ray sources likely associated with IC 2602, 31 X-ray sources in NGC 2632, and 31 near M67's central regions. We were further able to classify these X-ray sources based on their optical variability and any radio emission. Three IC 2602 X-ray sources had radio counterparts, which are likely all chromospherically active binary stars. We also identified luminous radio and X-ray variability from a spectroscopic triple system in M67, WOCS 3012/S1077, which is either consistent with a quiescent black hole binary, or due to an active binary stellar system. A recent population study of optical variables by Anderson & Hunt 2025 shows that the population of optical variables in open clusters clearly changes over cluster age; this pilot study gives evidence that the X-ray population also changes with time, and demonstrates the need for a broader multiwavelength study of Galactic open clusters.

preprint2024arXiv

Radio continuum from the most massive early-type galaxies detected with ASKAP RACS

All very massive early-type galaxies contain supermassive blackholes but are these blackholes all sufficiently active to produce detectable radio continuum sources? We have used the 887.5~MHz Rapid ASKAP Continuum Survey DR1 to measure the radio emission from morphological early-type galaxies brighter than $K_S=9.5$ selected from the 2MASS Redshift Survey, HyperLEDA and RC3. In line with previous studies, we find median radio power increases with infrared luminosity, with $P_{1.4} \propto L_K^{2.2}$, although the scatter about this relation spans several orders of magnitude. All 40 of the $M_K<-25.7$ early-type galaxies in our sample have measured radio flux densities that are more than $2σ$ above the background noise, with $1.4~{\rm GHz}$ radio powers spanning $\sim 3 \times 10^{20}$ to $\sim 3\times 10^{25}~{\rm W~Hz^{-1}}$. Cross matching our sample with integral field spectroscopy of early-type galaxies reveals that the most powerful radio sources preferentially reside in galaxies with relatively low angular momentum (i.e. slow rotators). While the infrared colours of most galaxies in our early-type sample are consistent with passive galaxies with negligible star formation and the radio emission produced by active galactic nuclei or AGN remnants, very low levels of star formation could power the weakest radio sources with little effect on many other star formation rate tracers.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Bulge-disk decomposition of KiDS data in the nearby universe

We derive single Sérsic fits and bulge-disk decompositions for 13096 galaxies at redshifts z < 0.08 in the GAMA II equatorial survey regions in the Kilo-Degree Survey (KiDS) g, r and i bands. The surface brightness fitting is performed using the Bayesian two-dimensional profile fitting code ProFit. We fit three models to each galaxy in each band independently with a fully automated Markov-chain Monte Carlo analysis: a single Sérsic model, a Sérsic plus exponential and a point source plus exponential. After fitting the galaxies, we perform model selection and flag galaxies for which none of our models are appropriate (mainly mergers/Irregular galaxies). The fit quality is assessed by visual inspections, comparison to previous works, comparison of independent fits of galaxies in the overlap regions between KiDS tiles and bespoke simulations. The latter two are also used for a detailed investigation of systematic error sources. We find that our fit results are robust across various galaxy types and image qualities with minimal biases. Errors given by the MCMC underestimate the true errors typically by factors 2-3. Automated model selection criteria are accurate to > 90 % as calibrated by visual inspection of a subsample of galaxies. We also present g-r component colours and the corresponding colour-magnitude diagram, consistent with previous works despite our increased fit flexibility. Such reliable structural parameters for the components of a diverse sample of galaxies across multiple bands will be integral to various studies of galaxy properties and evolution. All results are integrated into the GAMA database.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Self-Organizing Map Application on Nearby Galaxies

Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and Sérsic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a more complex population structure, represented by between 2 and 6 clusters. Here we use Self Organizing Maps (SOM), an unsupervised learning technique which can be used to visualize similarity in a higher dimensional space using a 2D representation, to map these five-dimensional clusters in the feature space onto two-dimensional projections. To further analyze these clusters, using the SOM information, we agree with previous results that the sub-populations found in the feature space can be reasonably mapped onto three or five clusters. We explore where the "green valley" galaxies are mapped onto the SOM, indicating multiple interstitial populations within the green valley population. Finally, we use the projection of the SOM to verify whether morphological information provided by GalaxyZoo users, for example, if features are visible, can be mapped onto the SOM-generated map. Voting on whether galaxies are smooth, likely ellipticals, or "featured" can reasonably be separated but smaller morphological features (bar, spiral arms) can not. SOMs promise to be a useful tool to map and identify instructive sub-populations in multidimensional galaxy survey feature space, provided they are large enough.

preprint2022arXiv

The eROSITA view of the Abell 3391/95 field: The Northern Clump. The largest infalling structure in the longest known gas filament observed with eROSITA, XMM-Newton, and Chandra

SRG/eROSITA PV observations revealed the A3391/95 cluster system and the Northern Clump (MCXC J0621.7-5242 galaxy cluster) are aligning along a cosmic filament in soft X-rays, similarly to what has been seen in simulations before. We aim to understand the dynamical state of the Northern Clump as it enters the atmosphere ($3\times R_{200}$) of A3391. We analyzed joint eROSITA, XMM-Newton, and Chandra observations to probe the morphological, thermal, and chemical properties of the Northern Clump from its center out to a radius of 988 kpc ($R_{200}$). We utilized the ASKAP/EMU radio data, DECam optical image, and Planck y-map to study the influence of the WAT radio source on the Northern Clump central ICM. From the Magneticum simulation, we identified an analog of the A3391/95 system along with an infalling group resembling the Northern Clump. The Northern Clump is a WCC cluster centered on a WAT radio galaxy. The gas temperature over $0.2-0.5R_{500}$ is $k_BT_{500}=1.99\pm0.04$ keV. We employed the $M-T$ scaling relation and obtained a mass estimate of $M_{500}=(7.68\pm0.43)\times10^{13}M_{\odot}$ and $R_{500}=(636\pm12)$ kpc. Its atmosphere has a boxy shape and deviates from spherical symmetry. We identify a southern surface brightness edge, likely caused by subsonic motion relative to the filament gas. At $\sim\! R_{500}$, the southern atmosphere appears to be 42% hotter than its northern atmosphere. We detect a downstream tail pointing toward the north with a projected length of $\sim318$ kpc, plausibly the result of ram pressure stripping. The analog group in the Magneticum simulation is experiencing changes in its gas properties and a shift between the position of the halo center and that of the bound gas while approaching the main cluster pair.

preprint2022arXiv

The Subaru HSC weak lensing mass-observable scaling relations of spectroscopic galaxy groups from the GAMA survey

We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) to study the dark matter content of galaxy groups in the Universe using weak lensing. We use galaxy groups from the Galaxy Mass and Assembly galaxy survey in approximately $100$ sq. degrees of the sky that overlap with the HSC survey as lenses. We restrict our analysis to the $1587$ groups with at least five members. We divide these groups into six bins each of group luminosity and group member velocity dispersion and measure the lensing signal with a signal-to-noise ratio of $55$ and $51$ for these two different selections, respectively. We use a Bayesian halo model framework to infer the halo mass distribution of our groups binned in the two different observable properties and constrain the power-law scaling relation, and the scatter between mean halo masses and the two group observable properties. We obtain a 5 percent constraint on the amplitude of the scaling relation between halo mass and group luminosity with $\avg{M} = (0.81\pm 0.04)\times10^{14}\hinvMsun$ for $L_{\rm grp}=10^{11.5}\hinvsqLsun$, and a power-law index of $α=1.01\pm 0.07$. We constrain the amplitude of the scaling relation between halo mass and velocity dispersion to be $\avg{M}=(0.93\pm 0.05)\times10^{14}\hinvMsun$ for $σ=500 \kms$ and a power-law index to be $α=1.52\pm0.10$. However, these scaling relations are sensitive to the exact cuts applied to the number of group members. Comparisons with similar scaling relations from the literature show that our results are consistent and have significantly reduced errors.

preprint2022arXiv

The variation of the gas content of galaxy groups and pairs compared to isolated galaxies

We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (provides environmental and galaxy properties), and (ii) DINGO pilot survey data of one of those fields. The mean $f_{HI}$ for different units (groups/pairs/isolated galaxies) are measured in regions of the log($\langle {\rm M_*} \rangle$) -- log($\langle {\rm SFR} \rangle$) plane, relative to the z $\sim 0$ star-forming main sequence (SFMS) of individual galaxies, by stacking $f_{HI}$ spectra of individual units. For ALFALFA, $f_{HI}$ spectra of units are measured by extracting HI spectra over the full groups/pair areas and dividing by the total stellar mass of member galaxies. For DINGO, $f_{HI}$ spectra of units are measured by co-adding HI spectra of individual member galaxies, followed by division by their total stellar mass. For all units the mean $f_{HI}$ decreases as we move to higher $\langle {\rm M_*} \rangle$ along the SFMS, and as we move from above the SFMS to below it at any $\langle {\rm M_*} \rangle$. From the DINGO-based study, mean $f_{HI}$ in groups appears to be lower compared to isolated galaxies for all $\langle {\rm M_*} \rangle$ along the SFMS. From the ALFALFA-based study we find substantially higher mean $f_{HI}$ in groups compared to isolated galaxies (values for pairs being intermediate) for ${\langle{\rm M_*}\rangle}\lesssim10^{9.5}~{\rm M_{\odot}}$, indicating the presence of substantial amounts of HI not associated with cataloged member galaxies in low mass groups.

preprint2022arXiv

The VMC Survey -- XLIX. Discovery of a population of quasars dominated by nuclear dust emission behind the Magellanic Clouds

Following the discovery of SAGE0536AGN ($z \sim$ 0.14), with the strongest 10-$μ$m silicate emission ever observed for an Active Galactic Nucleus (AGN), we discovered SAGE0534AGN ($z \sim$ 1.01), a similar AGN but with less extreme silicate emission. Both were originally mistaken as evolved stars in the Magellanic Clouds. Lack of far-infrared emission, and therefore star-formation, implies we are seeing the central engine of the AGN without contribution from the host galaxy. They could be a key link in galaxy evolution. We used a dimensionality reduction algorithm, t-SNE (t-distributed Stochastic Neighbourhood Embedding) with multi-wavelength data from Gaia EDR3, VISTA survey of the Magellanic Clouds, AllWISE and the Australian SKA Pathfinder to find these two unusual AGN are grouped with 16 other objects separated from the rest, suggesting a rare class. Our spectroscopy at SAAO/SALT and literature data confirm at least 14 of these objects are extragalactic ($0.13 < z < 1.23$), all hosting AGN. Using spectral energy distribution fitter CIGALE we find that the majority of dust emission ($> 70 \%$) in these sources is due to the AGN. Host galaxies appear to be either in or transitioning into the green valley. There is a trend of a thinning torus, increasing X-ray luminosity and decreasing Eddington ratio as the AGN transition through the green valley, implying that as the accretion supply depletes, the torus depletes and the column density reduces. Also, the near-infrared variability amplitude of these sources correlates with attenuation by the torus, implying the torus plays a role in the variability.

preprint2021arXiv

GAMA/DEVILS: Constraining the cosmic star-formation history from improved measurements of the 0.3-2.2 micron Extragalactic Background Light

We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star-formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey datasets, in 9 multi-wavelength filters to measure the COB in the range from 0.35 micron to 2.2 micron. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and Gamma-ray measurements are now fully consistent to within ~10%, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the cosmic star-formation history, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. We conclude that the peak of star-formation rate lies in the range 0.066-0.076 Msol/yr/Mpc^3 at a lookback time of 9.1 to 10.9 Gyrs.

preprint2020arXiv

Discovery of a Radio Relic in the Massive Merging Cluster SPT-CL 2023-5535 from the ASKAP-EMU PILOT SURVEY

The ASKAP-EMU survey is a deep wide-field radio continuum survey designed to cover the entire southern sky and a significant fraction of the northern sky up to $+30^{\circ}$. Here, we report a discovery of a radio relic in the merging cluster SPT-CL 2023-5535 at z=0.23 from the ASKAP-EMU pilot 300 sq. deg survey (800-1088 MHz). The deep high-resolution data reveal a $\sim2$ Mpc-scale radio halo elongated in the east-west direction, coincident with the intracluster gas. The radio relic is located at the western edge of this radio halo stretched $\sim0.5$ Mpc in the north-south orientation. The integrated spectral index of the radio relic within the narrow bandwidth is $α^{\scriptstyle \rm 1088~MHz}_{\scriptstyle \rm 800~MHz}=-0.76 \pm 0.06$. Our weak-lensing analysis shows that the system is massive ($M_{200}=1.04\pm0.36\times 10^{15} M_{\odot}$) and composed of at least three subclusters. We suggest a scenario, wherein the radio features arise from the collision between the eastern and middle subclusters. Our discovery illustrates the effectiveness of the ASKAP-EMU survey in detecting diffuse emissions in galaxy clusters and when completed, the survey will greatly increase the number of merging cluster detections with diffuse radio emissions.

preprint2020arXiv

GAMA+KiDS: Empirical correlations between halo mass and other galaxy properties near the knee of the stellar-to-halo mass relation

We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely: stellar colour, specific star formation rate, Sersic index, and effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range ($M_* \sim 2$--$5 \times 10^{10}$ Msol). This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sersic index are better predictors of halo mass than either colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more tightly correlated with galaxy structure than either past star formation history or current star formation rate. Our results lead to an approximate lower bound on the dispersion in halo masses among $\log M_* \approx {10.5}$ galaxies: we find that the dispersion is $\gtrsim 0.3$ dex. This would imply either that offsets from the mean SHMR are closely coupled to size/structure, or that the dispersion in the SHMR is larger than past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range.

preprint2020arXiv

KiDS+VIKING+GAMA: Testing semi-analytic models of galaxy evolution with galaxy-galaxy-galaxy lensing

Several semi-analytic models (SAMs) try to explain how galaxies form, evolve and interact inside the dark matter large-scale structure. These SAMs can be tested by comparing their predictions for galaxy-galaxy-galaxy-lensing (G3L), which is weak gravitational lensing around galaxy pairs, with observations. We evaluate the SAMs by Henriques et al. (2015; H15) and by Lagos et al. (2012; L12), implemented in the Millennium Run, by comparing their predictions for G3L to observations at smaller scales than previous studies and also for pairs of lens galaxies from different populations. We compare the G3L signal predicted by the SAMs to measurements in the overlap of the Galaxy And Mass Assembly survey (GAMA), the Kilo-Degree Survey (KiDS), and the VISTA Kilo-degree Infrared Galaxy survey (VIKING), splitting lens galaxies into two colour and five stellar-mass samples. Using an improved G3L estimator, we measure the three-point correlation of the matter distribution for mixed lens pairs with galaxies from different samples, and unmixed lens pairs with galaxies from the same sample. Predictions by the H15 SAM agree with the observations for all colour-selected and all but one stellar-mass-selected sample with 95% confidence. Deviations occur for lenses with stellar masses below $9.5h^{-2}\mathrm{M}_\odot$ at scales below $0.2h^{-1}\mathrm{Mpc}$. Predictions by the L12 SAM for stellar-mass selected samples and red galaxies are significantly higher than observed, while the predicted signal for blue galaxy pairs is too low. The L12 SAM predicts more pairs of small stellar-mass and red galaxies than the H15 SAM and the observations, as well as fewer pairs of blue galaxies. Likely explanations are different treatments of environmental effects by the SAMs and different models of the initial mass function. We conclude that G3L provides a stringent test for models of galaxy formation and evolution.

preprint2020arXiv

Unexpected Circular Radio Objects at High Galactic Latitude

We have found a class of circular radio objects in the Evolutionary Map of the Universe Pilot Survey, using the Australian Square Kilometre Array Pathfinder telescope. The objects appear in radio images as circular edge-brightened discs, about one arcmin diameter, that are unlike other objects previously reported in the literature. We explore several possible mechanisms that might cause these objects, but none seems to be a compelling explanation.

preprint2019arXiv

Assembly bias evidence in close galaxy pairs

The growth channel of massive galaxies involving mergers can be studied via close pairs as putative merger progenitors, where the stellar populations of the satellite galaxies will be eventually incorporated into the massive primaries. We extend our recent analysis of the GAMA-based sample of close pairs presented in Ferreras et al. to the general spectroscopic dataset of SDSS galaxies (DR14), for which the high S/N of the data enables a detailed analysis of the differences between satellite galaxies with respect to the mass of the primary galaxy. A sample of approximately two thousand satellites of massive galaxies is carefully selected within a relatively narrow redshift range (0.07<z<0.14). Two main parameters are considered as major drivers of the star formation history of these galaxies, namely: the stellar velocity dispersion of the satellite ($σ$), as a proxy of "local" drivers, and the ratio between the stellar mass of the satellite and the primary, $μ=M_{\rm SAT}/M_{\rm PRI}$, meant to serve as an indicator of environment. Consistently with the independent, GAMA-based work, we find that satellites around the most massive primaries appear older, at fixed velocity dispersion, with respect to satellites of lower mass primaries. This trend is more marked in lower mass satellites ($σ$~100 km/s), with SSP-equivalent age differences up to ~0.5 Gyr, and can be interpreted as a one-halo assembly bias, so that satellites corresponding to smaller values of the mass ratio $μ$ represent older structures, akin to fossil groups.

preprint2016arXiv

Dependence of GAMA galaxy halo masses on the cosmic web environment from 100 square degrees of KiDS weak lensing data

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): $\mathcal{M_\star}-R_{\rm e}$ relations of $z=0$ bulges, discs and spheroids

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

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): Active Galactic Nuclei in Pairs of Galaxies

There exist conflicting observations on whether or not the environment of broad and narrow line AGN differ and this consequently questions the validity of the AGN unification model. The high spectroscopic completeness of the GAMA survey makes it ideal for a comprehensive analysis of the close environment of galaxies. To exploit this, and conduct a comparative analysis of the environment of broad and narrow line AGN within GAMA, we use a double-Gaussian emission line fitting method to model the more complex line profiles associated with broad line AGN. We select 209 type 1 (i.e., unobscured), 464 type 1.5-1.9 (partially obscured), and 281 type 2 (obscured) AGN within the GAMA II database. Comparing the fractions of these with neighbouring galaxies out to a pair separation of $350\,\text{kpc }h^{-1}$ and $Δz < 0.012$ shows no difference between AGN of different type, except at separations less than $20\,\text{kpc }h^{-1}$ where our observations suggest an excess of type 2 AGN in close pairs. We analyse the properties of the galaxies neighbouring our AGN and find no significant differences in colour or the star formation activity of these galaxies. Further to this we find that $Σ_5$ is also consistent between broad and narrow line AGN. We conclude that the observations presented here are consistent with AGN unification.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): Galaxy colour gradients versus colour, structure and luminosity

Using single-component fits to SDSS/UKIDSS images of galaxies in the G09 region of the GAMA survey we study radial colour gradients across the galaxy population. We use the multiwavelength information provided by MegaMorph analysis of galaxy light profiles to calculate intrinsic colour gradients, and divide into six subsamples split by overall Sérsic index ($n$) and galaxy colour. We find a bimodality in the colour gradients of high- and low-$n$ galaxies in all wavebands, which varies with overall galaxy luminosity. Global trends in colour gradients therefore result from combining the contrasting behaviour of a number of different galaxy populations. The ubiquity of strong negative colour gradients supports the picture of inside-out growth through gas accretion for blue, low-$n$ galaxies, and through dry minor mergers for red, high-$n$ galaxies. An exception is the blue high-n population, with properties indicative of dissipative major mergers.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): the Stellar Mass Budget of Galaxy Spheroids and Disks

We build on a recent photometric decomposition analysis of 7506 Galaxy and Mass Assembly (GAMA) survey galaxies to derive stellar mass function fits to individual spheroid and disk component populations down to a lower mass limit of log(M_*/M_sun)= 8. We find that the spheroid/disk mass distributions for individual galaxy morphological types are well described by single Schechter function forms. We derive estimates of the total stellar mass densities in spheroids (rho_spheroid = 1.24+/-0.49 * 10^8 M_sun Mpc^-3 h_0.7) and disks (rho_disk = 1.20+/-0.45 * 10^8 M_sun Mpc^-3 h_0.7), which translates to approximately 50% of the local stellar mass density in spheroids and 48% in disks. The remaining stellar mass is found in the dwarf "little blue spheroid" class, which is not obviously similar in structure to either classical spheroid or disk populations. We also examine the variation of component mass ratios across galaxy mass and group halo mass regimes, finding the transition from spheroid to disk mass dominance occurs near galaxy stellar mass ~10^11 M_sun and group halo mass ~10^12.5 M_sun/h. We further quantify the variation in spheroid-to-total mass ratio with group halo mass for central and satellite populations as well as the radial variation of this ratio within groups.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): Understanding the wavelength dependence of galaxy structure with bulge-disc decompositions

With a large sample of bright, low-redshift galaxies with optical$-$near-IR imaging from the GAMA survey we use bulge-disc decompositions to understand the wavelength-dependent behavior of single-Sérsic structural measurements. We denote the variation in single-Sérsic index with wavelength as $\mathcal{N}$, likewise for effective radius we use $\mathcal{R}$. We find that most galaxies with a substantial disc, even those with no discernable bulge, display a high value of $\mathcal{N}$. The increase in Sérsic index to longer wavelengths is therefore intrinsic to discs, apparently resulting from radial variations in stellar population and/or dust reddening. Similarly, low values of $\mathcal{R}$ ($<$ 1) are found to be ubiquitous, implying an element of universality in galaxy colour gradients. We also study how bulge and disc colour distributions vary with galaxy type. We find that, rather than all bulges being red and all discs being blue in absolute terms, both components become redder for galaxies with redder total colours. We even observe that bulges in bluer galaxies are typically bluer than discs in red galaxies, and that bulges and discs are closer in colour for fainter galaxies. Trends in total colour are therefore not solely due to the colour or flux dominance of the bulge or disc.

preprint2016arXiv

ZFOURGE catalogue of AGN candidates: an enhancement of 160μm-derived star-formation rates in active galaxies to $z$ = 3.2

We investigate active galactic nuclei (AGN) candidates within the FourStar Galaxy Evolution Survey (ZFOURGE) to determine the impact they have on star-formation in their host galaxies. We first identify a population of radio, X-ray, and infrared-selected AGN by cross-matching the deep $K_{s}$-band imaging of ZFOURGE with overlapping multi-wavelength data. From this, we construct a mass-complete (log(M$_{*}$/M$_{\odot}$) $\ge$ 9.75), AGN luminosity limited sample of 235 AGN hosts over z = 0.2 - 3.2. We compare the rest-frame U - V versus V - J (UVJ) colours and specific star-formation rates (sSFRs) of the AGN hosts to a mass-matched control sample of inactive (non-AGN) galaxies. UVJ diagnostics reveal AGN tend to be hosted in a lower fraction of quiescent galaxies and a higher fraction of dusty galaxies than the control sample. Using 160μm Herschel PACS data, we find the mean specific star-formation rate of AGN hosts to be elevated by 0.34$\pm$0.07 dex with respect to the control sample across all redshifts. This offset is primarily driven by infrared-selected AGN, where the mean sSFR is found to be elevated by as much as a factor of ~5. The remaining population, comprised predominantly of X-ray AGN hosts, is found mostly consistent with inactive galaxies, exhibiting only a marginal elevation. We discuss scenarios that may explain these findings and postulate that AGN are less likely to be a dominant mechanism for moderating galaxy growth via quenching than has previously been suggested.

preprint2015arXiv

Direct Shear Mapping - a new weak lensing tool

We have developed a new technique called Direct Shear Mapping (DSM) to measure gravitational lensing shear directly from observations of a single background source. The technique assumes the velocity map of an un-lensed, stably-rotating galaxy will be rotationally symmetric. Lensing distorts the velocity map making it asymmetric. The degree of lensing can be inferred by determining the transformation required to restore axisymmetry. This technique is in contrast to traditional weak lensing methods, which require averaging an ensemble of background galaxy ellipticity measurements, to obtain a single shear measurement. We have tested the efficacy of our fitting algorithm with a suite of systematic tests on simulated data. We demonstrate that we are in principle able to measure shears as small as 0.01. In practice, we have fitted for the shear in very low redshift (and hence un-lensed) velocity maps, and have obtained null result with an error of $\pm 0.01$. This high sensitivity results from analysing spatially resolved spectroscopic images (i.e. 3D data cubes), including not just shape information (as in traditional weak lensing measurements) but velocity information as well. Spirals and rotating ellipticals are ideal targets for this new technique. Data from any large IFU or radio telescope is suitable, or indeed any instrument with spatially resolved spectroscopy such as SAMI, ALMA, HETDEX and SKA.

preprint2015arXiv

Direct Shear Mapping: Prospects for weak lensing studies of individual galaxy-galaxy lensing systems

We have investigated, using both a theoretical and an empirical approach, the frequency of low redshift galaxy-galaxy lensing systems in which the signature of weak lensing might be directly detectable. We find good agreement between these two approaches. In order to make a theoretical estimate of the weak lensing shear, $γ$, for each galaxy in a catalogue, we have made an estimate of the asymptotic circular velocity from the stellar mass using three different approaches: from a simulation based relation, from an empirically-derived relation, and using the baryonic Tully-Fisher relation. Using data from the Galaxy and Mass Assembly redshift survey we estimate the frequency of detectable weak lensing at low redshift. We find that to a redshift of $z\sim 0.6$, the probability of a galaxy being weakly lensed by at least $γ= 0.02$ is $\sim 0.01$. A scatter in the $M_*-M_h$ relation results in a shift towards higher measured shears for a given population of galaxies. Given this, and the good probability of weak lensing at low redshifts, we have investigated the feasibility of measuring the scatter in the $M_*-M_h$ relation using shear statistics. This is a novel measurement, and is made possible because DSM is able to make individual \itshape direct~\upshape shear measurements, in contrast to traditional weak lensing techniques which can only make statistical measurements. We estimate that for a shear measurement error of $Δγ= 0.02$ (consistent with the sensitivity of DSM), a sample of $\sim$50,000 spatially and spectrally resolved galaxies would allow a measurement of the scatter in the $M_*-M_h$ relation to be made. While there are no currently existing IFU surveys of this size, there are upcoming surveys which will provide this data (e.g The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), surveys with Hector, and the Square Kilometre Array (SKA)).

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): Panchromatic Data Release (far-UV --- far-IR) and the low-z energy budget

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$μ$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): the Stellar Mass Budget by Galaxy Type

We report an expanded sample of visual morphological classifications from the Galaxy and Mass Assembly (GAMA) survey phase two, which now includes 7,556 objects (previously 3,727 in phase one). We define a local (z <0.06) sample and classify galaxies into E, S0-Sa, SB0-SBa, Sab-Scd, SBab-SBcd, Sd-Irr, and "little blue spheroid" types. Using these updated classifications, we derive stellar mass function fits to individual galaxy populations divided both by morphological class and more general spheroid- or disk-dominated categories with a lower mass limit of log(Mstar/Msun) = 8 (one dex below earlier morphological mass function determinations). We find that all individual morphological classes and the combined spheroid-/bulge-dominated classes are well described by single Schechter stellar mass function forms. We find that the total stellar mass densities for individual galaxy populations and for the entire galaxy population are bounded within our stellar mass limits and derive an estimated total stellar mass density of rho_star = 2.5 x 10^8 Msun Mpc^-3 h_0.7, which corresponds to an approximately 4% fraction of baryons found in stars. The mass contributions to this total stellar mass density by galaxies that are dominated by spheroidal components (E and S0-Sa classes) and by disk components (Sab-Scd and Sd-Irr classes) are approximately 70% and 30%, respectively.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): the wavelength dependence of galaxy structure versus redshift and luminosity

We study how the sizes and radial profiles of galaxies vary with wavelength, by fitting Sersic functions simultaneously to imaging in nine optical and near-infrared bands. To quantify the wavelength dependence of effective radius we use the ratio, $\mathcal{R}$, of measurements in two restframe bands. The dependence of Sersic index on wavelength, $\mathcal{N}$, is computed correspondingly. Vulcani et al. (2014) have demonstrated that different galaxy populations present sharply contrasting behaviour in terms of $\mathcal{R}$ and $\mathcal{N}$. Here we study the luminosity dependence of this result. We find that at higher luminosities, early-type galaxies display a more substantial decrease in effective radius with wavelength, whereas late-types present a more pronounced increase in Sersic index. The structural contrast between types thus increases with luminosity. By considering samples at different redshifts, we demonstrate that lower data quality reduces the apparent difference between the main galaxy populations. However, our conclusions remain robust to this effect. We show that accounting for different redshift and luminosity selections partly reconciles the size variation measured by Vulcani et al. with the weaker trends found by other recent studies. Dividing galaxies by visual morphology confirms the behaviour inferred using morphological proxies, although the sample size is greatly reduced. Finally, we demonstrate that varying dust opacity and disc inclination can account for features of the joint distribution of $\mathcal{R}$ and $\mathcal{N}$ for late-type galaxies. However, dust does not appear to explain the highest values of $\mathcal{R}$ and $\mathcal{N}$. The bulge-disc nature of galaxies must also contribute to the wavelength-dependence of their structure.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): Trends in galaxy colours, morphology, and stellar populations with large scale structure, group, and pair environments

We explore trends in galaxy properties with Mpc-scale structures using catalogues of environment and large scale structure from the Galaxy And Mass Assembly (GAMA) survey. Existing GAMA catalogues of large scale structure, group and pair membership allow us to construct galaxy stellar mass functions for different environmental types. To avoid simply extracting the known underlying correlations between galaxy properties and stellar mass, we create a mass matched sample of galaxies with stellar masses between $9.5 \leq \log{M_*/h^{-2} M_{\odot}} \leq 11$ for each environmental population. Using these samples, we show that mass normalised galaxies in different large scale environments have similar energy outputs, $u-r$ colours, luminosities, and morphologies. Extending our analysis to group and pair environments, we show galaxies that are not in groups or pairs exhibit similar characteristics to each other regardless of broader environment. For our mass controlled sample, we fail to see a strong dependence of Sérsic index or galaxy luminosity on halo mass, but do find that it correlates very strongly with colour. Repeating our analysis for galaxies that have not been mass controlled introduces and amplifies trends in the properties of galaxies in pairs, groups, and large scale structure, indicating that stellar mass is the most important predictor of the galaxy properties we examine, as opposed to environmental classifications.

preprint2015arXiv

Matching Radio Catalogs with Realistic Geometry: Application to SWIRE and ATLAS

Crossmatching catalogs at different wavelengths is a difficult problem in astronomy, especially when the objects are not point-like. At radio wavelengths an object can have several components corresponding, for example, to a core and lobes. {Considering not all radio detections correspond to visible or infrared sources, matching these catalogs can be challenging.} Traditionally this is done by eye for better quality, which does not scale to the large data volumes expected from the next-generation of radio telescopes. We present a novel automated procedure, using Bayesian hypothesis testing, to achieve reliable associations by explicit modelling of a particular class of radio-source morphology. {The new algorithm not only assesses the likelihood of an association between data at two different wavelengths, but also tries to assess whether different radio sources are physically associated, are double-lobed radio galaxies, or just distinct nearby objects.} Application to the SWIRE and ATLAS CDF-S catalogs shows that this method performs well without human intervention.

preprint2014arXiv

Galaxy and Mass Assembly (GAMA): Fine filaments of galaxies detected within voids

Based on data from the Galaxy and Mass Assembly (GAMA) survey, we report on the discovery of structures that we refer to as `tendrils' of galaxies: coherent, thin chains of galaxies that are rooted in filaments and terminate in neighbouring filaments or voids. On average, tendrils contain 6 galaxies and span 10 $h^{-1}$ Mpc. We use the so-called line correlation function to prove that tendrils represent real structures rather than accidental alignments. We show that voids found in the SDSS-DR7 survey that overlap with GAMA regions contain a large number of galaxies, primarily belonging to tendrils. This implies that void sizes are strongly dependent on the number density and sensitivity limits of a survey. We caution that galaxies in low density regions, that may be defined as `void galaxies,' will have local galaxy number densities that depend on such observational limits and are likely higher than can be directly measured.

preprint2014arXiv

Galaxy and Mass Assembly (GAMA): merging galaxies and their properties

We derive the close pair fractions and volume merger rates as a function of luminosity and morphology for galaxies in the GAMA survey with -23 < M(r) < -17 at 0.01 < z < 0.22. The merger fraction is about 0.015 at all luminosities (assuming 1/2 of pairs merge) and the volume merger rate is about 0.00035 per cubic Mpc per Gyr. Dry mergers (between red or spheroidal galaxies) are uncommon and decrease with decreasing luminosity. Fainter mergers are wet, between blue or disky galaxies. Damp mergers (one of each type) follow the average of dry and wet mergers. In the brighter luminosity bin (-23 < M(r) < -20) the merger rate evolution is flat, irrespective of colour or morphology. The makeup of the merging population does not change since z = 0.2. Major mergers and dry mergers appear comparatively unimportant in the buildup of the red sequence over the past 2 Gyr. We compare the colour, morphology, environmental density and degree of activity of galaxies in pairs to those of more isolated objects in the same volume. Galaxies in close pairs tend to be both redder and slightly more spheroid-dominated. This may be due to "harassment" in multiple previous passes prior to the current interaction. Galaxy pairs do not appear to prefer significantly denser environments. There is no evidence of an enhancement in the AGN fraction in pairs, compared to other galaxies in the same volume.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): Refining the Local Galaxy Merger Rate using Morphological Information

We use the Galaxy And Mass Assembly (GAMA) survey to measure the local Universe mass dependent merger fraction and merger rate using galaxy pairs and the CAS structural method, which identifies highly asymmetric merger candidate galaxies. Our goals are to determine which types of mergers produce highly asymmetrical galaxies, and to provide a new measurement of the local galaxy major merger rate. We examine galaxy pairs at stellar mass limits down to $M_{*} = 10^{8}M_{\odot}$ with mass ratios of $<$100:1 and line of sight velocity differences of $ΔV<500$ km s$^{-1}$. We find a significant increase in mean asymmetries for projected separations less than the sum of the individual galaxy's Petrosian 90 radii. For systems in major merger pairs with mass ratios of $<$4:1 both galaxies in the pair show a strong increase in asymmetry, while in minor merger systems (with mass ratios of $>$4:1) the lower mass companion becomes highly asymmetric, while the larger galaxy is much less affected. The fraction of highly asymmetric paired galaxies which have a major merger companion is highest for the most massive galaxies and drops progressively with decreasing mass. We calculate that the mass dependent major merger fraction is fairly constant at $\sim1.3-2\%$ between $10^{9.5}<M_{*}<10^{11.5} M_{\odot}$, and increases to $\sim4\%$ at lower masses. When the observability time scales are taken into consideration, the major merger rate is found to approximately triple over the mass range we consider. The total co-moving volume major merger rate over the range $10^{8.0}<M_{*}<10^{11.5} M_{\odot}$ is $(1.2 \pm 0.5) \times 10^{-3}$ $h^{3}_{70}$ Mpc$^{-3}$ Gyr$^{-1}$.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): Stellar mass functions by Hubble type

We present an estimate of the galaxy stellar mass function and its division by morphological type in the local (0.025 < z < 0.06) Universe. Adopting robust morphological classifications as previously presented (Kelvin et al.) for a sample of 3,727 galaxies taken from the Galaxy And Mass Assembly survey, we define a local volume and stellar mass limited sub-sample of 2,711 galaxies to a lower stellar mass limit of M = 10^9.0 M_sun. We confirm that the galaxy stellar mass function is well described by a double Schechter function given by M* = 10^10.64 M_sun, α1 = -0.43, ϕ*1 = 4.18 dex^-1 Mpc^-3, α2 = -1.50 and ϕ*2 = 0.74 dex^-1 Mpc^-3. The constituent morphological-type stellar mass functions are well sampled above our lower stellar mass limit, excepting the faint little blue spheroid population of galaxies. We find approximately 71+3-4% of the stellar mass in the local Universe is found within spheroid dominated galaxies; ellipticals and S0-Sas. The remaining 29+4-3% falls predominantly within late type disk dominated systems, Sab-Scds and Sd-Irrs. Adopting reasonable bulge-to-total ratios implies that approximately half the stellar mass today resides in spheroidal structures, and half in disk structures. Within this local sample, we find approximate stellar mass proportions for E : S0-Sa : Sab-Scd : Sd-Irr of 34 : 37 : 24 : 5.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): ugrizYJHK Sérsic luminosity functions and the cosmic spectral energy distribution by Hubble type

We report the morphological classification of 3727 galaxies from the Galaxy and Mass Assembly survey with M_r < -17.4 mag and in the redshift range 0.025 < z < 0.06 (2.1 x 10^5 Mpc^3 ) into E, S0-Sa, SB0-SBa, Sab-Scd, SBab-SBcd, Sd-Irr and little blue spheroid classes. Approximately 70% of galaxies in our sample are disk dominated systems, with the remaining ~30% spheroid dominated. We establish the robustness of our classifications, and use them to derive morphological-type luminosity functions and luminosity densities in the ugrizYJHK passbands, improving on prior studies that split by global colour or light profile shape alone. We find that the total galaxy luminosity function is best described by a double-Schechter function while the constituent morphological-type luminosity functions are well described by a single-Schechter function. These data are also used to derive the star-formation rate densities for each Hubble class, and the attenuated and unattenuated (corrected for dust) cosmic spectral energy distributions, i.e., the instantaneous energy production budget. While the observed optical/near-IR energy budget is dominated 58:42 by galaxies with a significant spheroidal component, the actual energy production rate is reversed, i.e., the combined disk dominated populations generate ~1.3x as much energy as the spheroid dominated populations. On the grandest scale, this implies that chemical evolution in the local Universe is currently confined to mid-type spiral classes like our Milky Way.

preprint2014arXiv

Galaxy And Mass Assembly: Deconstructing Bimodality - I. Red ones and blue ones

We measure the mass functions for generically red and blue galaxies, using a z < 0.12 sample of log M* > 8.7 field galaxies from the Galaxy And Mass Assembly (GAMA) survey. Our motivation is that, as we show, the dominant uncertainty in existing measurements stems from how 'red' and 'blue' galaxies have been selected/defined. Accordingly, we model our data as two naturally overlapping populations, each with their own mass function and colour-mass relation, which enables us characterise the two populations without having to specify a priori which galaxies are 'red' and 'blue'. Our results then provide the means to derive objective operational definitions for the terms 'red' and 'blue', which are based on the phenomenology of the colour-mass diagrams. Informed by this descriptive modelling, we show that: 1.) after accounting for dust, the stellar colours of 'blue' galaxies do not depend strongly on mass; 2.) the tight, flat 'dead sequence' does not extend much below log M* ~ 10.5; instead, 3.) the stellar colours of 'red' galaxies vary rather strongly with mass, such that lower mass 'red' galaxies have bluer stellar populations; 4.) below log M* ~ 9.3, the 'red' population dissolves into obscurity, and it becomes problematic to talk about two distinct populations; as a consequence, 5.) it is hard to meaningfully constrain the shape, including the possibility of an upturn, of the 'red' galaxy mass function below log M* ~ 9. Points 1-4 provide meaningful targets for models of galaxy formation and evolution to aim for.

preprint2014arXiv

Inferring the Redshift Distribution of the Cosmic Infrared Background

Cross-correlating the Planck High Frequency Instrument (HFI) maps against quasars from the Sloan Digital Sky Survey (SDSS) DR7, we estimate the intensity distribution of the Cosmic Infrared Background (CIB) over the redshift range 0 < z < 5.We detect redshift-dependent spatial cross-correlations between the two datasets using the 857, 545 and 353 GHz channels and we obtain upper limits at 217 GHz consistent with expectations. At all frequencies with detectable signal we infer a redshift distribution peaking around z ~ 1.2 and find the recovered spectrum to be consistent with emission arising from star forming galaxies. By assuming simple modified blackbody and Kennicutt relations, we estimate dust and star formation rate density as a function of redshift, finding results consistent with earlier multiwavelength measurements over a large portion of cosmic history. However, we note that, lacking mid-infrared coverage, we are not able to make an accurate determination of the mean temperature for the dust responsible for the CIB. Our results demonstrate that clustering-based redshift inference is a valuable tool for measuring the entire evolution history of the cosmic star formation rate from a single and homogeneous dataset.

preprint2013arXiv

A Merger Shock in Abell 2034

We present a $250\,$ks Chandra observation of the cluster merger A2034 with the aim of understanding the nature of a sharp edge previously characterized as a cold front. The new data reveal that the edge is coherent over a larger opening angle and is significantly more bow-shock-shaped than previously thought. Within $\sim 27\,$degrees about the axis of symmetry of the edge the density, temperature and pressure drop abruptly by factors of $1.83^{+0.09}_{-0.08}$, $1.85^{+0.41}_{-0.41}$ and $3.4^{+0.8}_{-0.7}$, respectively. This is inconsistent with the pressure equilibrium expected of a cold front and we conclude that the edge is a shock front. We measure a Mach number $M = 1.59^{+0.06}_{-0.07}$ and corresponding shock velocity $v_{\rm shock}\simeq 2057\,$km/s. Using spectra collected at the MMT with the Hectospec multi-object spectrograph we identify 328 spectroscopically confirmed cluster members. Significantly, we find a local peak in the projected galaxy density associated with a bright cluster galaxy which is located just ahead of the nose of the shock. The data are consistent with a merger viewed within $\sim 23\,$degrees of the plane of the sky. The merging subclusters are now moving apart along a north-south axis approximately $0.3\,$Gyr after a small impact parameter core passage. The gas core of the secondary subcluster, which was driving the shock, appears to have been disrupted by the merger. Without a driving 'piston' we speculate that the shock is dying. Finally, we propose that the diffuse radio emission near the shock is due to the revival of pre-existing radio plasma which has been overrun by the shock.

preprint2013arXiv

Galaxy And Mass Assembly (GAMA): improved cosmic growth measurements using multiple tracers of large-scale structure

We present the first application of a "multiple-tracer" redshift-space distortion (RSD) analysis to an observational galaxy sample, using data from the Galaxy and Mass Assembly survey (GAMA). Our dataset is an r < 19.8 magnitude-limited sample of 178,579 galaxies covering redshift interval z < 0.5 and area 180 deg^2. We obtain improvements of 10-20% in measurements of the gravitational growth rate compared to a single-tracer analysis, deriving from the correlated sample variance imprinted in the distributions of the overlapping galaxy populations. We present new expressions for the covariances between the auto-power and cross-power spectra of galaxy samples that are valid for a general survey selection function and weighting scheme. We find no evidence for a systematic dependence of the measured growth rate on the galaxy tracer used, justifying the RSD modelling assumptions, and validate our results using mock catalogues from N-body simulations. For multiple tracers selected by galaxy colour, we measure normalized growth rates in two independent redshift bins f*sigma_8(z=0.18) = 0.36 +/- 0.09 and f*sigma_8(z=0.38) = 0.44 +/- 0.06, in agreement with standard GR gravity and other galaxy surveys at similar redshifts.

preprint2013arXiv

On the 3D structure of the mass, metallicity, and SFR space for SF galaxies

We demonstrate that the space formed by the star-formation rate (SFR), gas-phase metallicity (Z), and stellar mass (M), can be reduced to a plane, as first proposed by Lara-Lopez et al. We study three different approaches to find the best representation of this 3D space, using a principal component analysis, a regression fit, and binning of the data. The PCA shows that this 3D space can be adequately represented in only 2 dimensions, i.e., a plane. We find that the plane that minimises the chi^2 for all variables, and hence provides the best representation of the data, corresponds to a regression fit to the stellar mass as a function of SFR and $Z$, M=f(Z,SFR). We find that the distribution resulting from the median values in bins for our data gives the highest chi^2. We also show that the empirical calibrations to the oxygen abundance used to derive the Fundamental Metallicity Relation (Nagao et al.) have important limitations, which contribute to the apparent inconsistencies. The main problem is that these empirical calibrations do not consider the ionization degree of the gas. Furthermore, the use of the N2 index to estimate oxygen abundances cannot be applied for ~8.8 because of the saturation of the [NII]6584 line in the high-metallicity regime. Finally we provide an update of the Fundamental Plane derived by Lara-Lopez et al.

preprint2013arXiv

The Cosmic Star Formation Rate from the Faintest Galaxies in the Unobservable Universe

Observations of high-z galaxies and gamma-ray bursts now allow for empirical studies during reionization. However, even deep surveys see only the brightest galaxies at any epoch and must extrapolate to arbitrary lower limits to estimate the total rate of star formation. We first argue that the galaxy populations seen in LBG surveys yield a GRB rate at z > 8 that is an order of magnitude lower than observed. We find that integrating the inferred UV luminosity functions down to M_UV ~ -10 brings LBG- and GRB-inferred SFR density values into agreement up to z ~ 8. GRBs, however, favor a far larger amount of as yet unseen star formation at z > 9. We suggest that the SFR density may only slowly decline out to z ~ 11, in accord with WMAP and Planck reionization results, and that GRBs may be useful in measuring the scale of this multitude of dwarf galaxies.

preprint2013arXiv

The dependency of AGN infrared color-selection on source luminosity and obscuration: An observational perspective in CDFS and COSMOS

AIMS. This work addresses the AGN IR-selection dependency on intrinsic source luminosity and obscuration, in order to identify and charaterise biases which could affect conclusions in studies. METHODS. We study IR-selected AGN in the Chandra Deep Field South (CDFS) survey and in the Cosmological Survey (COSMOS). The AGN sample is divided into low and high X-ray luminosity classes, and into unobscured (type-1) and obscured (type-2) classes by means of X-ray and optical spectroscopy data. Specifically in the X-rays regime, we adopt the intrinsic luminosity taking into account the estimated column density (N_H). We also take the chance to highlight important differences resulting from adopting different methods to assess AGN obscuration. RESULTS. In agreement with previous studies, we also find that AGN IR-selection efficiency shows a decrease with decreasing source AGN X-ray luminosity. For the intermediate-luminosity AGN population (43.3<~log(L_X[erg/s])<~44$), the efficiency also worsens with increasing obscuration. The same sample also shows an evolution with cosmic time of the obscured fraction at the highest X-ray luminosities, independently of adopted type-1/type-2 classification method. CONCLUSIONS. We confirm that AGN IR-selection is genuinely biased toward unobscured AGN, but only at intermediate luminosities. At the highest luminosities, where AGN IR-selection is more efficient, there is no obscuration bias. We show that type-1 AGN are intrinsically more luminous than type-2 AGN only at z<~1.6, thus resulting in more type-1 AGN being selected the shallower the IR survey is. [ABRIDGED: Conclusion on survey strategy. See manuscript.] Finally, when the James Webb Space Telescope comes online, the broad-band filters 2.0um, 4.4um, 7.7um, and 18um will be key to disentangle AGN from non-AGN dominated SEDs at depths where spectroscopy becomes impractical.

preprint2012arXiv

A new infrared color criterion for the selection of 0<z<7 AGN: application to deep fields and implications for JWST surveys

[Abridged] It is widely accepted that observations at mid-infrared (mid-IR) wavelengths enable the selection of galaxies with nuclear activity, which may not be revealed even in the deepest X-ray surveys. In this work new near- and mid-IR color diagnostics are explored, aiming for improved efficiency - better completeness and less contamination - in selecting AGN out to very high redshifts. We restrict our study to the James Webb Space Telescope wavelength range (0.6-27um). The criteria are created based on the predictions by state-of-the-art galaxy and AGN templates covering a wide variety of galaxy properties, and tested against control samples with deep multi-wavelength coverage (ranging from the X-rays to radio frequencies). We show that the colors Ks-[4.5], [4.5]-[8.0], and [8.0]-[24] are ideal as AGN/non-AGN diagnostics at, respectively, z<~1, 1<~z<~2.5, and z>~2.5-3. However, when the source redshift is unknown, these colors should be combined. We thus develop an improved IR criterion (using Ks and IRAC bands, KI) as a new alternative at z<~2.5. KI does not show improved completeness (50-60% overall) in comparison to commonly used IRAC-based AGN criteria, but is less affected by non-AGN contamination (revealing a >50-90% level of successful AGN selection). We also propose KIM (using Ks, IRAC, and MIPS-24um bands, KIM), which aims to select AGN hosts from local distances to as far back as the end of reionization (0<z<~7) with reduced non-AGN contamination. However, the necessary testing-constraints and the small control-sample sizes prevent the confirmation of its improved efficiency at z<~2.5. Overall, KIM shows a ~30-40% completeness and a >70-90% level of successful AGN selection. KI and KIM are built to be reliable against a ~10-20% error in flux, are based on existing filters, and are suitable for immediate use.

preprint2012arXiv

First Science with SAMI: A Serendipitously Discovered Galactic Wind in ESO 185-G031

We present the first scientific results from the Sydney-AAO Multi-Object IFS (SAMI) at the Anglo-Australian Telescope. This unique instrument deploys 13 fused fibre bundles (hexabundles) across a one-degree field of view allowing simultaneous spatially-resolved spectroscopy of 13 galaxies. During the first SAMI commissioning run, targeting a single galaxy field, one object (ESO 185-G031) was found to have extended minor axis emission with ionisation and kinematic properties consistent with a large-scale galactic wind. The importance of this result is two-fold: (i) fibre bundle spectrographs are able to identify low-surface brightness emission arising from extranuclear activity; (ii) such activity may be more common than presently assumed because conventional multi-object spectrographs use single-aperture fibres and spectra from these are nearly always dominated by nuclear emission. These early results demonstrate the extraordinary potential of multi-object hexabundle spectroscopy in future galaxy surveys.

preprint2012arXiv

Galaxy And Mass Assembly: Estimating galaxy group masses via caustic analysis

We have generated complementary halo mass estimates for all groups in the Galaxy And Mass Assembly Galaxy Group Catalogue (GAMA G3Cv1) using a modified caustic mass estimation algorithm, originally developed by Diaferio & Geller (1997). We calibrate the algorithm by applying it on a series of 9 GAMA mock galaxy light cones and investigate the effects of using different definitions for group centre and size. We select the set of parameters that provide median-unbiased mass estimates when tested on mocks, and generate mass estimates for the real group catalogue. We find that on average, the caustic mass estimates agree with dynamical mass estimates within a factor of 2 in 90.8 +/- 6.1% groups and compares equally well to velocity dispersion based mass estimates for both high and low multiplicity groups over the full range of masses probed by the G3Cv1.

preprint2012arXiv

MegaMorph - multi-wavelength measurement of galaxy structure: complete Sérsic profile information from modern surveys

In this paper, we demonstrate a new method for fitting galaxy profiles which makes use of the full multi-wavelength data provided by modern large optical-near-infrared imaging surveys. We present a new version of GALAPAGOS, which utilises a recently-developed multi-wavelength version of GALFIT, and enables the automated measurement of wavelength dependent Sérsic profile parameters for very large samples of galaxies. Our new technique is extensively tested to assess the reliability of both pieces of software, GALFIT and GALAPAGOS on both real ugrizY JHK imaging data from the GAMA survey and simulated data made to the same specifications. We find that fitting galaxy light profiles with multi-wavelength data increases the stability and accuracy of the measured parameters, and hence produces more complete and meaningful multi-wavelength photometry than has been available previously. The improvement is particularly significant for magnitudes in low S/N bands and for structural parameters like half-light radius re and Sérsic index n for which a prior is used by constraining these parameters to a polynomial as a function of wavelength. This allows the fitting routines to push the magnitude of galaxies for which sensible values can be derived to fainter limits. The technique utilises a smooth transition of galaxy parameters with wavelength, creating more physically meaningful transitions than single-band fitting and allows accurate interpolation between passbands, perfect for derivation of rest-frame values.

preprint2012arXiv

Multiscale probability mapping: groups, clusters and an algorithmic search for filaments in SDSS

We have developed a multiscale structure identification algorithm for the detection of overdensities in galaxy data that identifies structures having radii within a user-defined range. Our "multiscale probability mapping" technique combines density estimation with a shape statistic to identify local peaks in the density field. This technique takes advantage of a user-defined range of scale sizes, which are used in constructing a coarse-grained map of the underlying fine-grained galaxy distribution, from which overdense structures are then identified. In this study we have compiled a catalogue of groups and clusters at 0.025 < z < 0.24 based on the Sloan Digital Sky Survey, Data Release 7, quantifying their significance and comparing with other catalogues. Most measured velocity dispersions for these structures lie between 50 and 400 km/s. A clear trend of increasing velocity dispersion with radius from 0.2 to 1 Mpc/h is detected, confirming the lack of a sharp division between groups and clusters. A method for quantifying elongation is also developed to measure the elongation of group and cluster environments. By using our group and cluster catalogue as a coarse-grained representation of the galaxy distribution for structure sizes of <~ 1 Mpc/h, we identify 53 filaments (from an algorithmically-derived set of 100 candidates) as elongated unions of groups and clusters at 0.025 < z < 0.13. These filaments have morphologies that are consistent with previous samples studied.

preprint2012arXiv

The Australia Telescope Large Area Survey: Spectroscopic Catalogue and Radio Luminosity Functions

The Australia Telescope Large Area Survey (ATLAS) has surveyed seven square degrees of sky around the Chandra Deep Field South (CDFS) and the European Large Area ISO Survey - South 1 (ELAIS-S1) fields at 1.4 GHz. ATLAS aims to reach a uniform sensitivity of $10 μ$Jy beam$^{-1}$ rms over the entire region with data release 1 currently reaching $\sim30 μ$Jy beam$^{-1}$ rms. Here we present 466 new spectroscopic redshifts for radio sources in ATLAS as part of our optical follow-up program. Of the 466 radio sources with new spectroscopic redshifts, 142 have star-forming optical spectra, 282 show evidence for AGN in their optical spectra, 10 have stellar spectra and 32 have spectra revealing redshifts, but with insufficient features to classify. We compare our spectroscopic classifications with two mid-infrared diagnostics and find them to be in broad agreement. We also construct the radio luminosity function for star-forming galaxies to z $= 0.5$ and for AGN to z $= 0.8$. The radio luminosity function for star-forming galaxies appears to be in good agreement with previous studies. The radio luminosity function for AGN appears higher than previous studies of the local AGN radio luminosity function. We explore the possibility of evolution, cosmic variance and classification techniques affecting the AGN radio luminosity function. ATLAS is a pathfinder for the forthcoming EMU survey and the data presented in this paper will be used to guide EMU's survey design and early science papers.

preprint2011arXiv

An unbiased sample of bright southern Compact Steep Spectrum and Gigahertz Peaked Spectrum Sources

Compact Steep Spectrum (CSS) and Gigahertz Peaked Spectrum (GPS) sources are classes of compact, powerful, extragalactic objects. These sources are thought to be the earliest stages in the evolution of radio galaxies, capturing the ignition (or, in some cases, re-ignition) of the AGN. As well as serving as probes of the early stages of large-scale radio sources, these sources are good, stable, amplitude calibrators for radio telescopes. We present an unbiased flux density limited (>1.5 Jy at 2.7 GHz) catalogue of these objects in the Southern Hemisphere, including tabulated data, radio spectra, and where available, optical images and measurements. The catalogue contains 26 sources, consisting of 2 new candidate and 15 known CSS sources, and 9 known GPS sources. We present new Australia Telescope Compact Array (ATCA) data on ten of these 26 sources, and data on a further 42 sources which were excluded from our final sample. This bright sample will serve as a reference sample for comparison with subsequent faint (mJy level) samples of CSS and GPS candidates currently being compiled.

preprint2011arXiv

Galaxy And Mass Assembly (GAMA): Structural Investigation of Galaxies via Model Analysis (SIGMA)

We present single-Sérsic two-dimensional model fits to 167,600 galaxies modelled independently in the ugrizYJHK bandpasses using reprocessed Sloan Digital Sky Survey Data Release Seven (SDSS DR7) and UKIRT Infrared Deep Sky Survey Large Area Survey (UKIDSS-LAS) imaging data available from the GAMA database. In order to facilitate this study we developed SIGMA, an R wrapper around several contemporary astronomy software packages including Source Extractor, PSF Extractor and GALFIT 3. SIGMA produces realistic 2D model fits to galaxies, employing automatic adaptive background subtraction and empirical PSF measurements on the fly for each galaxy in GAMA. Using these results, we define a common coverage area across the three GAMA regions containing 138,269 galaxies. We provide Sérsic magnitudes truncated at 10 re which show good agreement with SDSS Petrosian and GAMA photometry for low Sérsic index systems (n < 4), and much improved photometry for high Sérsic index systems (n > 4), recovering as much as Δm = 0.5 magnitudes in the r band. We employ a K band Sérsic index/u - r colour relation to delineate the massive (n > ~2) early-type galaxies (ETGs) from the late-type galaxies (LTGs). The mean Sérsic index of these ETGs shows a smooth variation with wavelength, increasing by 30% from g through K. LTGs exhibit a more extreme change in Sérsic index, increasing by 52% across the same range. In addition, ETGs and LTGs exhibit a 38% and 25% decrease respectively in half-light radius from g through K. These trends are shown to arise due to the effects of dust attenuation and stellar population/metallicity gradients within galaxy populations.

preprint2011arXiv

Pixel-z: Studying Substructure and Stellar Populations in Galaxies out to z~3 using Pixel Colors I. Systematics

We perform a pixel-by-pixel analysis of 467 galaxies in the GOODS-VIMOS survey to study systematic effects in extracting properties of stellar populations (age, dust, metallicity and SFR) from pixel colors using the pixel-z method. The systematics studied include the effect of the input stellar population synthesis model, passband limitations and differences between individual SED fits to pixels and global SED-fitting to a galaxy's colors. We find that with optical-only colors, the systematic errors due to differences among the models are well constrained. The largest impact on the age and SFR e-folding time estimates in the pixels arises from differences between the Maraston models and the Bruzual&Charlot models, when optical colors are used. This results in systematic differences larger than the 2σ uncertainties in over 10 percent of all pixels in the galaxy sample. The effect of restricting the available passbands is more severe. In 26 percent of pixels in the full sample, passband limitations result in systematic biases in the age estimates which are larger than the 2σ uncertainties. Systematic effects from model differences are reexamined using Near-IR colors for a subsample of 46 galaxies in the GOODS-NICMOS survey. For z > 1, the observed optical/NIR colors span the rest frame UV-optical SED, and the use of different models does not significantly bias the estimates of the stellar population parameters compared to using optical-only colors. We then illustrate how pixel-z can be applied robustly to make detailed studies of substructure in high redshift galaxies such as (a) radial gradients of age, SFR, sSFR and dust and (b) the distribution of these properties within subcomponents such as spiral arms and clumps. Finally, we show preliminary results from the CANDELS survey illustrating how the new HST/WFC3 data can be exploited to probe substructure in z~1-3 galaxies.

preprint2011arXiv

The Sydney-AAO Multi-object Integral field spectrograph (SAMI)

We demonstrate a novel technology that combines the power of the multi-object spectrograph with the spatial multiplex advantage of an integral field spectrograph (IFS). The Sydney-AAO Multi-object IFS (SAMI) is a prototype wide-field system at the Anglo-Australian Telescope (AAT) that allows 13 imaging fibre bundles ("hexabundles") to be deployed over a 1-degree diameter field of view. Each hexabundle comprises 61 lightly-fused multimode fibres with reduced cladding and yields a 75 percent filling factor. Each fibre core diameter subtends 1.6 arcseconds on the sky and each hexabundle has a field of view of 15 arcseconds diameter. The fibres are fed to the flexible AAOmega double-beam spectrograph, which can be used at a range of spectral resolutions (R=lambda/delta(lambda) ~ 1700-13000) over the optical spectrum (3700-9500A). We present the first spectroscopic results obtained with SAMI for a sample of galaxies at z~0.05. We discuss the prospects of implementing hexabundles at a much higher multiplex over wider fields of view in order to carry out spatially--resolved spectroscopic surveys of 10^4 to 10^5 galaxies.

preprint2011arXiv

The Telescopes and Processes of the Australian Astronomical Observatory

The Australian Astronomical Observatory operates the Anglo-Australian Telescope and the United Kingdom Schmidt Telescope in Australia, as well as coordinating access for the Australian community to the Gemini, Magellan, and other international telescope facilities. We review here the processes involved within the AAO related to allocating observing time on these facilities, as well as the impact on telescope use of both the Large Program projects and the AAO's instrument program.

preprint2010arXiv

Galaxy and Mass Assembly (GAMA): FUV, NUV, ugrizYJHK Petrosian, Kron and Sèrsic photometry

In order to generate credible 0.1-2 μm SEDs, the GAMA project requires many Gigabytes of imaging data from a number of instruments to be re-processed into a standard format. In this paper we discuss the software infrastructure we use, and create self-consistent ugrizYJHK photometry for all sources within the GAMA sample. Using UKIDSS and SDSS archive data, we outline the pre-processing necessary to standardise all images to a common zeropoint, the steps taken to correct for seeing bias across the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA regions in each filter. From these mosaics, we extract source catalogues for the GAMA regions using elliptical Kron and Petrosian matched apertures. We also calculate Sérsic magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component modelling wrapper for GALFIT 3. We compare the resultant photometry directly, and also calculate the r band galaxy LF for all photometric datasets to highlight the uncertainty introduced by the photometric method. We find that (1) Changing the object detection threshold has a minor effect on the best-fitting Schechter parameters of the overall population (M* +/- 0.055mag, α +/- 0.014, Φ* +/- 0.0005 h^3 Mpc^{-3}). (2) An offset between datasets that use Kron or Petrosian photometry regardless of the filter. (3) The decision to use circular or elliptical apertures causes an offset in M* of 0.20mag. (4) The best-fitting Schechter parameters from total-magnitude photometric systems (such as SDSS modelmag or Sérsic magnitudes) have a steeper faint-end slope than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's total luminosity density, when calculated using Kron or Petrosian r-band photometry, is underestimated by at least 15%.

preprint2009arXiv

Determining the escape fraction of ionizing photons during reionization with the GRB derived star-formation rate

The fraction of ionizing photons that escape their host galaxies and so are able to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era and early galaxy formation. Studies of the reionization history normally suffer from a degeneracy between the unknown values for the efficiency with which high redshift galaxies turn mass into stars and the escape fraction of ionizing photons. Recent gamma-ray burst (GRB) measurements of the star formation rate density during reionization provide the first opportunity to break this degeneracy. We confront a semi-analytic model for reionization with the GRB-derived star formation rate, as well as observations of the Ly-alpha forest and the CMB. Assuming that UV photons produced in star-forming galaxies dominate the reionization process, we show that the escape fraction of ionizing photons from high redshift galaxies is ~5% [log(f_esc)=-1.35+/-0.15 (68%)] for our fiducial model. This value is reasonably stable against uncertainties in the modeling, including the implementation of radiative feedback, the possibility of an evolving escape fraction, and the unknown shape of the IMF, which in sum contribute ~0.2 dex of additional systematic uncertainty on the value of escape fraction.

preprint2009arXiv

Discovery of Large-Scale Gravitational Infall in a Massive Protostellar Cluster

We report Mopra (ATNF), Anglo-Australian Telescope, and Atacama Submillimeter Telescope Experiment observations of a molecular clump in Carina, BYF73 = G286.21+0.17, which give evidence of large-scale gravitational infall in the dense gas. From the millimetre and far-infrared data, the clump has mass ~ 2 x 10^4 Msun, luminosity ~ 2-3 x 10^4 Lsun, and diameter ~ 0.9 pc. From radiative transfer modelling, we derive a mass infall rate ~ 3.4 x 10^-2 Msun yr-1. If confirmed, this rate for gravitational infall in a molecular core or clump may be the highest yet seen. The near-infrared K-band imaging shows an adjacent compact HII region and IR cluster surrounded by a shell-like photodissociation region showing H2 emission. At the molecular infall peak, the K imaging also reveals a deeply embedded group of stars with associated H2 emission. The combination of these features is very unusual and we suggest they indicate the ongoing formation of a massive star cluster. We discuss the implications of these data for competing theories of massive star formation.

preprint2009arXiv

GAMA: towards a physical understanding of galaxy formation

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

preprint2009arXiv

The Star Formation Rate in the Reionization Era as Indicated by Gamma-ray Bursts

High-redshift gamma-ray bursts (GRBs) offer an extraordinary opportunity to study aspects of the early Universe, including the cosmic star formation rate (SFR). Motivated by the two recent highest-z GRBs, GRB 080913 at z = 6.7 and GRB 090423 at z = 8.1, and more than four years of Swift observations, we first confirm that the GRB rate does not trace the SFR in an unbiased way. Correcting for this, we find that the implied SFR to beyond z = 8 is consistent with LBG-based measurements after accounting for unseen galaxies at the faint end of the UV luminosity function. We show that this provides support for the integrated star formation in the range 6 < z < 8 to have been alone sufficient to reionize the Universe.

preprint2008arXiv

Matching the Local and Cosmic Star Formation Histories

Given the many recent advances in our understanding of the star formation history (SFH) of the Local Group and other nearby galaxies, and in the evolution of star formation with redshift, we present a new comparison of the comoving space density of the star formation rate as a function of look-back time for the Local and Distant Universe. We update the Local SFH derived from the analysis of resolved stellar populations (``fossil records'') in individual nearby galaxies, based on our own estimations as well as available in the literature. While the preliminary comparison of SFHs is found to be broadly consistent, some discrepancies still remain, including an excess of the Local SFR density in the most recent epoch.

preprint2008arXiv

Revealing the High-Redshift Star Formation Rate with Gamma-Ray Bursts

While the high-z frontier of star formation rate (SFR) studies has advanced rapidly, direct measurements beyond z ~ 4 remain difficult, as shown by significant disagreements among different results. Gamma-ray bursts, owing to their brightness and association with massive stars, offer hope of clarifying this situation, provided that the GRB rate can be properly related to the SFR. The Swift GRB data reveal an increasing evolution in the GRB rate relative to the SFR at intermediate z; taking this into account, we use the highest-z GRB data to make a new determination of the SFR at z = 4-7. Our results exceed the lowest direct SFR measurements, and imply that no steep drop exists in the SFR up to at least z ~ 6.5. We discuss the implications of our result for cosmic reionization, the efficiency of the universe in producing stellar-mass black holes, and ``GRB feedback'' in star-forming hosts.

preprint2006arXiv

On the normalisation of the cosmic star formation history

Strong constraints on the cosmic star formation history (SFH) have recently been established using ultraviolet and far-infrared measurements, refining the results of numerous measurements over the past decade. Taken together, the most recent and robust data indicate a compellingly consistent picture of the SFH out to redshift z~6, with especially tight constraints for z<~1. We fit these data with simple analytical forms, and derive conservative bands to indicate possible variations from the best fits. Since the z<~1 SFH data are quite precise, we investigate the sequence of assumptions and corrections that together affect the SFH normalisation, to test their accuracy, both in this redshift range and beyond. As lower limits on this normalisation, we consider the evolution in stellar mass density, metal mass density, and supernova rate density, finding it unlikely that the SFH normalisation is much lower than indicated by our direct fit. Additionally, predictions from the SFH for supernova type Ia rate densities tentatively suggests delay times of ~3 Gyr. As a corresponding upper limit on the SFH normalisation, we consider the Super-Kamiokande (SK) limit on the electron antineutrino flux from past core-collapse supernovae, which applies primarily to z<~1. We find consistency with the SFH only if the neutrino temperatures from SN events are relatively modest. Constraints on the assumed initial mass function (IMF) also become apparent. The traditional Salpeter IMF, assumed for convenience by many authors, is known to be a poor representation at low stellar masses (<~ 1 solar mass), and we show that recently favoured IMFs are also constrained. In particular, somewhat shallow, or top-heavy, IMFs may be preferred, although they cannot be too top-heavy. (Abridged)