A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.

Maximizing the public impact of astronomy projects in the next decade requires NSF-funded centers to support the development of online, mobile-friendly outreach and education activities. EPO teams with astronomy, education, and web development expertise should be in place to build accessible programs at scale and support astronomers doing outreach.

We present a study of the dust, stars and atomic gas (HI) in an HI-selected sample of local galaxies (z<0.035) in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) fields. This HI-selected sample reveals a population of very high gas fraction (>80 per cent), low stellar mass sources that appear to be in the earliest stages of their evolution. We compare this sample with dust and stellar mass selected samples to study the dust and gas scaling relations over a wide range of gas fraction (proxy for evolutionary state of a galaxy). The most robust scaling relations for gas and dust are those linked to NUV-r (SSFR) and gas fraction, these do not depend on sample selection or environment. At the highest gas fractions, our additional sample shows the dust content is well below expectations from extrapolating scaling relations for more evolved sources, and dust is not a good tracer of the gas content. The specific dust mass for local galaxies peaks at a gas fraction of ~75 per cent. The atomic gas depletion time is also longer for high gas fraction galaxies, opposite to the trend found for molecular gas depletion timescale. We link this trend to the changing efficiency of conversion of HI to H2 as galaxies increase in stellar mass surface density as they evolve. Finally, we show that galaxies start out barely obscured and increase in obscuration as they evolve, yet there is no clear and simple link between obscuration and global galaxy properties.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

Commodity cloud computing, as provided by commercial vendors such as Amazon, Google, and Microsoft, has revolutionized computing in many sectors. With the advent of a new class of big data, public access astronomical facility such as LSST, DKIST, and WFIRST, there exists a real opportunity to combine these missions with cloud computing platforms and fundamentally change the way astronomical data is collected, processed, archived, and curated. Making these changes in a cross-mission, coordinated way can provide unprecedented economies of scale in personnel, data collection and management, archiving, algorithm and software development and, most importantly, science.

We present Hubble Space Telescope WFC3-IR imaging in the fields of six apparently bright dusty star-forming galaxies (DSFGs) at $z$ = 2-4 identified by their rest-frame far-infrared colors using the Planck and Herschel space facilities. We detect near-infrared counterparts for all six submillimeter sources, allowing us to undertake strong-lensing analyses. One field in particular stands out for its prominent giant arcs, PLCK G165.7+67.0 (G165). After combining the color and morphological information, we identify 11 sets of image multiplicities in this one field. We construct a strong-lensing model constrained by this lensing evidence, which uncovers a bimodal spatial mass distribution, and from which we measure a mass of $(2.6 \pm 0.11)$ $\times$ $10^{14}$ $M_{\odot}$ within $\sim$250 kpc. The bright ($S_{350}$ $\approx$ 750 mJy) DSFG appears as two images: a giant arc with a spatial extent of 4.5" that is merging with the critical curve, and a lower-magnification counterimage that is detected in our new longer-wavelength ground- and space-based imaging data. Using our ground-based spectroscopy, we calculate a dynamical mass of $1.3^{+0.04}_{-0.70} \times 10^{15}$ $M_{\odot}$ to the same fixed radius, although this value may be inflated relative to the true value if the velocity distribution is enhanced in the line-of-sight direction. We suggest that the bimodal mass taken in combination with the weak X-ray flux and low SZ decrement may be explained as a pre-merger for which the intracluster gas is diluted along the line of sight, while the integrated surface mass density is supercritical to strong-lensing effects.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We present the first major release of data from the SAMI Galaxy Survey. This data release focuses on the emission-line physics of galaxies. Data Release One includes data for 772 galaxies, about 20% of the full survey. Galaxies included have the redshift range 0.004 < z < 0.092, a large mass range (7.6 < log(Mstellar/M$_\odot$) < 11.6), and star-formation rates of 10^-4 to 10^1\ M$_\odot$/yr. For each galaxy, we include two spectral cubes and a set of spatially resolved 2D maps: single- and multi-component emission-line fits (with dust extinction corrections for strong lines), local dust extinction and star-formation rate. Calibration of the fibre throughputs, fluxes and differential-atmospheric-refraction has been improved over the Early Data Release. The data have average spatial resolution of 2.16 arcsec (FWHM) over the 15~arcsec diameter field of view and spectral (kinematic) resolution R=4263 (sigma=30km/s) around Halpha. The relative flux calibration is better than 5\% and absolute flux calibration better than $\pm0.22$~mag, with the latter estimate limited by galaxy photometry. The data are presented online through the Australian Astronomical Observatory's Data Central.

We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.

We describe the selection of galaxies targeted in eight low redshift clusters (APMCC0917, A168, A4038, EDCC442, A3880, A2399, A119 and A85; $0.029 < z < 0.058$) as part of the Sydney-AAO Multi-Object integral field Spectrograph Galaxy Survey (SAMI-GS). We have conducted a redshift survey of these clusters using the AAOmega multi-object spectrograph on the 3.9m Anglo-Australian Telescope. The redshift survey is used to determine cluster membership and to characterise the dynamical properties of the clusters. In combination with existing data, the survey resulted in 21,257 reliable redshift measurements and 2899 confirmed cluster member galaxies. Our redshift catalogue has a high spectroscopic completeness ($\sim 94\%$) for $r_{\rm petro} \leq 19.4$ and clustercentric distances $R< 2\rm{R}_{200}$. We use the confirmed cluster member positions and redshifts to determine cluster velocity dispersion, $\rm{R}_{200}$, virial and caustic masses, as well as cluster structure. The clusters have virial masses $14.25 \leq {\rm log }({\rm M}_{200}/\rm{M}_{\odot}) \leq 15.19$. The cluster sample exhibits a range of dynamical states, from relatively relaxed-appearing systems, to clusters with strong indications of merger-related substructure. Aperture- and PSF-matched photometry are derived from SDSS and VST/ATLAS imaging and used to estimate stellar masses. These estimates, in combination with the redshifts, are used to define the input target catalogue for the cluster portion of the SAMI-GS. The primary SAMI-GS cluster targets have $R< \rm{R}_{200}$, velocities $|v_{\rm pec}| < 3.5\sigma_{200}$ and stellar masses $9.5 \leq {\rm log(M}^*_{approx}/\rm{M}_{\odot}) \leq 12$. Finally, we give an update on the SAMI-GS progress for the cluster regions.

We use data from the Sydney-AAO Multi-Object Integral Field Spectrograph (SAMI) Galaxy Survey and the Galaxy And Mass Assembly (GAMA) survey to investigate the spatially-resolved signatures of the environmental quenching of star formation in galaxies. Using dust-corrected measurements of the distribution of H$\alpha$ emission we measure the radial profiles of star formation in a sample of 201 star-forming galaxies covering three orders of magnitude in stellar mass (M$_{*}$; $10^{8.1}$-$10^{10.95}\, $M$_{\odot}$) and in $5^{th}$ nearest neighbour local environment density ($\Sigma_{5}$; $10^{-1.3}$-$10^{2.1}\,$Mpc$^{-2}$). We show that star formation rate gradients in galaxies are steeper in dense ($\log_{10}(\Sigma_{5}/$Mpc$^{2})>0.5$) environments by $0.58\pm 0.29\, dex\, $r$_{e}^{-1}$ in galaxies with stellar masses in the range $10^{10}<$M$_{*}/$M$_{\odot}<10^{11}$ and that this steepening is accompanied by a reduction in the integrated star formation rate. However, for any given stellar mass or environment density the star-formation morphology of galaxies shows large scatter. We also measure the degree to which the star formation is centrally concentrated using the unitless scale-radius ratio ($r_{50,H\alpha}/r_{50,cont}$), which compares the extent of ongoing star formation to previous star formation. With this metric we find that the fraction of galaxies with centrally concentrated star formation increases with environment density, from $\sim 5\pm 4\%$ in low-density environments ($\log_{10}(\Sigma_{5}/$Mpc$^{2})<0.0$) to $30\pm 15\%$ in the highest density environments ($\log_{10}(\Sigma_{5}/$Mpc$^{2})>1.0$). These lines of evidence strongly suggest that with increasing local environment density the star formation in galaxies is suppressed, and that this starts in their outskirts such that quenching occurs in an outside-in fashion in dense environments and is not instantaneous.

The Dark Energy Survey (DES) is currently undertaking an observational program imaging $1/4$ of the southern hemisphere sky with unprecedented photometric accuracy. In the process of observing millions of faint stars and galaxies to constrain the parameters of the dark energy equation of state, the DES will obtain pre-discovery images of the regions surrounding an estimated 100 gamma-ray bursts (GRBs) over five years. Once GRBs are detected by, e.g., the Swift satellite, the DES data will be extremely useful for follow-up observations by the transient astronomy community. We describe a recently-commissioned suite of software that listens continuously for automated notices of GRB activity, collates useful information from archival DES data, and disseminates relevant data products back to the community in near-real-time. Of particular importance are the opportunities that non-public DES data provide for relative photometry of the optical counterparts of GRBs, as well as for identifying key characteristics (e.g., photometric redshifts) of potential GRB host galaxies. We provide the functional details of the DESAlert software as it presently operates, as well as the data products that it produces, and we show sample results from the application of DESAlert to several previously-detected GRBs.

We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7sigma. Cosmological constraints from the measurements in this work are presented in a companion paper (DES et al. 2015).

We present the Lambda Adaptive Multi-Band Deblending Algorithm in R (LAMBDAR), a novel code for calculating matched aperture photometry across images that are neither pixel- nor PSF-matched, using prior aperture definitions derived from high resolution optical imaging. The development of this program is motivated by the desire for consistent photometry and uncertainties across large ranges of photometric imaging, for use in calculating spectral energy distributions. We describe the program, specifically key features required for robust determination of panchromatic photometry: propagation of apertures to images with arbitrary resolution, local background estimation, aperture normalisation, uncertainty determination and propagation, and object deblending. Using simulated images, we demonstrate that the program is able to recover accurate photometric measurements in both high-resolution, low-confusion, and low-resolution, high-confusion, regimes. We apply the program to the 21-band photometric dataset from the Galaxy And Mass Assembly (GAMA) Panchromatic Data Release (PDR; Driver et al. 2016), which contains imaging spanning the far-UV to the far-IR. We compare photometry derived from LAMBDAR with that presented in Driver et al. (2016), finding broad agreement between the datasets. Nonetheless, we demonstrate that the photometry from LAMBDAR is superior to that from the GAMA PDR, as determined by a reduction in the outlier rate and intrinsic scatter of colours in the LAMBDAR dataset. We similarly find a decrease in the outlier rate of stellar masses and star formation rates using LAMBDAR photometry. Finally, we note an exceptional increase in the number of UV and mid-IR sources able to be constrained, which is accompanied by a significant increase in the mid-IR colour-colour parameter-space able to be explored.

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.

The nature of dark matter, dark energy and large-scale gravity pose some of the most pressing questions in cosmology today. These fundamental questions require highly precise measurements, and a number of wide-field spectroscopic survey instruments are being designed to meet this requirement. A key component in these experiments is the development of a simulation tool to forecast science performance, define requirement flow-downs, optimize implementation, demonstrate feasibility, and prepare for exploitation. We present SPOKES (SPectrOscopic KEn Simulation), an end-to-end simulation facility for spectroscopic cosmological surveys designed to address this challenge. SPOKES is based on an integrated infrastructure, modular function organization, coherent data handling and fast data access. These key features allow reproducibility of pipeline runs, enable ease of use and provide flexibility to update functions within the pipeline. The cyclic nature of the pipeline offers the possibility to make the science output an efficient measure for design optimization and feasibility testing. We present the architecture, first science, and computational performance results of the simulation pipeline. The framework is general, but for the benchmark tests, we use the Dark Energy Spectrometer (DESpec), one of the early concepts for the upcoming project, the Dark Energy Spectroscopic Instrument (DESI). We discuss how the SPOKES framework enables a rigorous process to optimize and exploit spectroscopic survey experiments in order to derive high-precision cosmological measurements optimally.

We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range $0.2 < z < 1.2$, a cross-correlation signal is detected at $6 \sigma$ and $4\sigma$ with SPT and Planck respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant ($>$$2 \sigma$) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the $z<1.2$ universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck cosmology, a $1.7\sigma$ deviation.

Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.

We study the clustering of galaxies detected at $i<22.5$ in the Science Verification observations of the Dark Energy Survey (DES). Two-point correlation functions are measured using $2.3\times 10^6$ galaxies over a contiguous 116 deg$^2$ region in five bins of photometric redshift width $\Delta z = 0.2$ in the range $0.2 < z < 1.2.$ The impact of photometric redshift errors are assessed by comparing results using a template-based photo-$z$ algorithm (BPZ) to a machine-learning algorithm (TPZ). A companion paper (Leistedt et al 2015) presents maps of several observational variables (e.g. seeing, sky brightness) which could modulate the galaxy density. Here we characterize and mitigate systematic errors on the measured clustering which arise from these observational variables, in addition to others such as Galactic dust and stellar contamination. After correcting for systematic effects we measure galaxy bias over a broad range of linear scales relative to mass clustering predicted from the Planck $\Lambda$CDM model, finding agreement with CFHTLS measurements with $\chi^2$ of 4.0 (8.7) with 5 degrees of freedom for the TPZ (BPZ) redshifts. We test a "linear bias" model, in which the galaxy clustering is a fixed multiple of the predicted non-linear dark-matter clustering. The precision of the data allow us to determine that the linear bias model describes the observed galaxy clustering to $2.5\%$ accuracy down to scales at least $4$ to $10$ times smaller than those on which linear theory is expected to be sufficient.

We measure the weak lensing shear around galaxy troughs, i.e. the radial alignment of background galaxies relative to underdensities in projections of the foreground galaxy field over a wide range of redshift in Science Verification data from the Dark Energy Survey. Our detection of the shear signal is highly significant (10 to 15$\sigma$ for the smallest angular scales) for troughs with the redshift range z in [0.2,0.5] of the projected galaxy field and angular diameters of 10 arcmin...1{\deg}. These measurements probe the connection between the galaxy, matter density, and convergence fields. By assuming galaxies are biased tracers of the matter density with Poissonian noise, we find agreement of our measurements with predictions in a fiducial Lambda cold dark matter model. The prediction for the lensing signal on large trough scales is virtually independent of the details of the underlying model for the connection of galaxies and matter. Our comparison of the shear around troughs with that around cylinders with large galaxy counts is consistent with a symmetry between galaxy and matter over- and underdensities. In addition, we measure the two-point angular correlation of troughs with galaxies which, in contrast to the lensing signal, is sensitive to galaxy bias on all scales. The lensing signal of troughs and their clustering with galaxies is therefore a promising probe of the statistical properties of matter underdensities and their connection to the galaxy field.

Using the science verification data of the Dark Energy Survey (DES) for a new sample of 106 X-Ray selected clusters and groups, we study the stellar mass growth of Bright Central Galaxies (BCGs) since redshift 1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncertainties associated with cluster mass, redshift, and BCG stellar mass measurements into analysis of a redshift-dependent BCG-cluster mass relation, $m_{*}\propto(\frac{M_{200}}{1.5\times 10^{14}M_{\odot}})^{0.24\pm 0.08}(1+z)^{-0.19\pm0.34}$, and compare the observed relation to the model prediction. We estimate the average growth rate since $z = 1.0$ for BCGs hosted by clusters of $M_{200, z}=10^{13.8}M_{\odot}$, at $z=1.0$: $m_{*, BCG}$ appears to have grown by $0.13\pm0.11$ dex, in tension at $\sim 2.5 \sigma$ significance level with the $0.40$ dex growth rate expected from the semi-analytic model. We show that the buildup of extended intra-cluster light after $z=1.0$ may alleviate this tension in BCG growth rates.

We combine Herschel/SPIRE sub-millimeter (submm) observations with existing multi-wavelength data to investigate the characteristics of low redshift, optically red galaxies detected in submm bands. We select a sample of galaxies in the redshift range 0.01$\leq$z$\leq$0.2, having >5$\sigma$ detections in the SPIRE 250 micron submm waveband. Sources are then divided into two sub-samples of $red$ and $blue$ galaxies, based on their UV-optical colours. Galaxies in the $red$ sample account for $\approx$4.2 per cent of the total number of sources with stellar masses M$_{*}\gtrsim$10$^{10}$ Solar-mass. Following visual classification of the $red$ galaxies, we find that $\gtrsim$30 per cent of them are early-type galaxies and $\gtrsim$40 per cent are spirals. The colour of the $red$-spiral galaxies could be the result of their highly inclined orientation and/or a strong contribution of the old stellar population. It is found that irrespective of their morphological types, $red$ and $blue$ sources occupy environments with more or less similar densities (i.e., the $\Sigma_5$ parameter). From the analysis of the spectral energy distributions (SEDs) of galaxies in our samples based on MAGPHYS, we find that galaxies in the $red$ sample (of any morphological type) have dust masses similar to those in the $blue$ sample (i.e. normal spiral/star-forming systems). However, in comparison to the $red$-spirals and in particular $blue$ systems, $red$-ellipticals have lower mean dust-to-stellar mass ratios. Besides galaxies in the $red$-elliptical sample have much lower mean star-formation/specific-star-formation rates in contrast to their counterparts in the $blue$ sample. Our results support a scenario where dust in early-type systems is likely to be of an external origin.

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$\mu$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/

Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high order discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed second order finite volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than a FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall-back to first order can be avoided. Finally, DG requires no wide stencils at high order, and offers an improved data locality and a focus on local computations, which is favourable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.

Hydrodynamical simulations of galaxy formation such as the Illustris simulations have progressed to a state where they approximately reproduce the observed stellar mass function from high to low redshift. This in principle allows self-consistent models of reionization that exploit the accurate representation of the diffuse gas distribution together with the realistic growth of galaxies provided by these simulations, within a representative cosmological volume. In this work, we apply and compare two radiative transfer algorithms implemented in a GPU-accelerated code to the $106.5\,{\rm Mpc}$ wide volume of Illustris in postprocessing in order to investigate the reionization transition predicted by this model. We find that the first generation of galaxies formed by Illustris is just about able to reionize the universe by redshift $z\sim 7$, provided quite optimistic assumptions about the escape fraction and the resolution limitations are made. Our most optimistic model finds an optical depth of $\tau\simeq 0.065$, which is in very good agreement with recent Planck 2015 determinations. Furthermore, we show that moment-based approaches for radiative transfer with the M1 closure give broadly consistent results with our angular-resolved radiative transfer scheme. In our favoured fiducial model, 20\% of the hydrogen is reionized by redshift $z=9.20$, and this rapidly climbs to 80\% by redshift $z=6.92$. It then takes until $z=6.24$ before 99\% of the hydrogen is ionized. On average, reionization proceeds `inside-out' in our models, with a size distribution of reionized bubbles that progressively features regions of ever larger size while the abundance of small bubbles stays fairly constant.

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

We present photometric redshift estimates for galaxies used in the weak lensing analysis of the Dark Energy Survey Science Verification (DES SV) data. Four model- or machine learning-based photometric redshift methods -- ANNZ2, BPZ calibrated against BCC-Ufig simulations, SkyNet, and TPZ -- are analysed. For training, calibration, and testing of these methods, we construct a catalogue of spectroscopically confirmed galaxies matched against DES SV data. The performance of the methods is evaluated against the matched spectroscopic catalogue, focusing on metrics relevant for weak lensing analyses, with additional validation against COSMOS photo-zs. From the galaxies in the DES SV shear catalogue, which have mean redshift $0.72\pm0.01$ over the range $0.3<z<1.3$, we construct three tomographic bins with means of $z=\{0.45, 0.67, 1.00\}$. These bins each have systematic uncertainties $\delta z \lesssim 0.05$ in the mean of the fiducial SkyNet photo-z $n(z)$. We propagate the errors in the redshift distributions through to their impact on cosmological parameters estimated with cosmic shear, and find that they cause shifts in the value of $\sigma_8$ of approx. 3%. This shift is within the one sigma statistical errors on $\sigma_8$ for the DES SV shear catalog. We further study the potential impact of systematic differences on the critical surface density, $\Sigma_{\mathrm{crit}}$, finding levels of bias safely less than the statistical power of DES SV data. We recommend a final Gaussian prior for the photo-z bias in the mean of $n(z)$ of width $0.05$ for each of the three tomographic bins, and show that this is a sufficient bias model for the corresponding cosmology analysis.

We present a mass map reconstructed from weak gravitational lensing shear measurements over 139 sq. deg from the Dark Energy Survey (DES) Science Verification data. The mass map probes both luminous and dark matter, thus providing a tool for studying cosmology. We find good agreement between the mass map and the distribution of massive galaxy clusters identified using a red-sequence cluster finder. Potential candidates for super-clusters and voids are identified using these maps. We measure the cross-correlation between the mass map and a magnitude-limited foreground galaxy sample and find a detection at the 5-7 sigma level on a large range of scales. These measurements are consistent with simulated galaxy catalogs based on LCDM N-body simulations, suggesting low systematics uncertainties in the map. We summarize our key findings in this letter; the detailed methodology and tests for systematics are presented in a companion paper.

We introduce redMaGiC, an automated algorithm for selecting Luminous Red Galaxies (LRGs). The algorithm was specifically developed to minimize photometric redshift uncertainties in photometric large-scale structure studies. redMaGiC achieves this by self-training the color-cuts necessary to produce a luminosity-thresholded LRG sample of constant comoving density. We demonstrate that redMaGiC photozs are very nearly as accurate as the best machine-learning based methods, yet they require minimal spectroscopic training, do not suffer from extrapolation biases, and are very nearly Gaussian. We apply our algorithm to Dark Energy Survey (DES) Science Verification (SV) data to produce a redMaGiC catalog sampling the redshift range $z\in[0.2,0.8]$. Our fiducial sample has a comoving space density of $10^{-3}\ (h^{-1} Mpc)^{-3}$, and a median photoz bias ($z_{spec}-z_{photo}$) and scatter $(\sigma_z/(1+z))$ of 0.005 and 0.017 respectively. The corresponding $5\sigma$ outlier fraction is 1.4%. We also test our algorithm with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) and Stripe 82 data, and discuss how spectroscopic training can be used to control photoz biases at the 0.1% level.

Weak gravitational lensing allows one to reconstruct the spatial distribution of the projected mass density across the sky. These "mass maps" provide a powerful tool for studying cosmology as they probe both luminous and dark matter. In this paper, we present a weak lensing mass map reconstructed from shear measurements in a 139 sq. deg area from the Dark Energy Survey (DES) Science Verification (SV) data. We compare the distribution of mass with that of the foreground distribution of galaxies and clusters. The overdensities in the reconstructed map correlate well with the distribution of optically detected clusters. We demonstrate that candidate superclusters and voids along the line of sight can be identified, exploiting the tight scatter of the cluster photometric redshifts. We cross-correlate the mass map with a foreground magnitude-limited galaxy sample from the same data. Our measurement gives results consistent with mock catalogs from N-body simulations that include the primary sources of statistical uncertainties in the galaxy, lensing, and photo-z catalogs. The statistical significance of the cross-correlation is at the 6.8-sigma level with 20 arcminute smoothing. A major goal of this study is to investigate systematic effects arising from a variety of sources, including PSF and photo-z uncertainties. We make maps derived from twenty variables that may characterize systematics and find the principal components. We find that the contribution of systematics to the lensing mass maps is generally within measurement uncertainties. In this work, we analyze less than 3% of the final area that will be mapped by the DES; the tools and analysis techniques developed in this paper can be applied to forthcoming larger datasets from the survey.

Spatially-varying depth and characteristics of observing conditions, such as seeing, airmass, or sky background, are major sources of systematic uncertainties in modern galaxy survey analyses, in particular in deep multi-epoch surveys. We present a framework to extract and project these sources of systematics onto the sky, and apply it to the Dark Energy Survey (DES) to map the observing conditions of the Science Verification (SV) data. The resulting distributions and maps of sources of systematics are used in several analyses of DES SV to perform detailed null tests with the data, and also to incorporate systematics in survey simulations. We illustrate the complementarity of these two approaches by comparing the SV data with the BCC-UFig, a synthetic sky catalogue generated by forward-modelling of the DES SV images. We analyse the BCC-UFig simulation to construct galaxy samples mimicking those used in SV galaxy clustering studies. We show that the spatially-varying survey depth imprinted in the observed galaxy densities and the redshift distributions of the SV data are successfully reproduced by the simulation and well-captured by the maps of observing conditions. The combined use of the maps, the SV data and the BCC-UFig simulation allows us to quantify the impact of spatial systematics on $N(z)$, the redshift distributions inferred using photometric redshifts. We conclude that spatial systematics in the SV data are mainly due to seeing fluctuations and are under control in current clustering and weak lensing analyses. The framework presented here is relevant to all multi-epoch surveys, and will be essential for exploiting future surveys such as the Large Synoptic Survey Telescope (LSST), which will require detailed null-tests and realistic end-to-end image simulations to correctly interpret the deep, high-cadence observations of the sky.

The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low-redshift galaxies. Covering an area of ~286 deg^2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238,000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm - 1 m. Here we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component S\'ersic fits, stellar masses, H$\alpha$-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72,225 objects in total). The database serving these data is available at http://www.gama-survey.org/.

We cross-match galaxy cluster candidates selected via their Sunyaev-Zel'dovich effect (SZE) signatures in 129.1 deg$^2$ of the South Pole Telescope 2500d SPT-SZ survey with optically identified clusters selected from the Dark Energy Survey (DES) science verification data. We identify 25 clusters between $0.1\lesssim z\lesssim 0.8$ in the union of the SPT-SZ and redMaPPer (RM) samples. RM is an optical cluster finding algorithm that also returns a richness estimate for each cluster. We model the richness $\lambda$-mass relation with the following function $\langle\ln\lambda|M_{500}\rangle\propto B_\lambda\ln M_{500}+C_\lambda\ln E(z)$ and use SPT-SZ cluster masses and RM richnesses $\lambda$ to constrain the parameters. We find $B_\lambda= 1.14^{+0.21}_{-0.18}$ and $C_\lambda=0.73^{+0.77}_{-0.75}$. The associated scatter in mass at fixed richness is $\sigma_{\ln M|\lambda} = 0.18^{+0.08}_{-0.05}$ at a characteristic richness $\lambda=70$. We demonstrate that our model provides an adequate description of the matched sample, showing that the fraction of SPT-SZ selected clusters with RM counterparts is consistent with expectations and that the fraction of RM selected clusters with SPT-SZ counterparts is in mild tension with expectation. We model the optical-SZE cluster positional offset distribution with the sum of two Gaussians, showing that it is consistent with a dominant, centrally peaked population and a sub-dominant population characterized by larger offsets. We also cross-match the RM catalog with SPT-SZ candidates below the official catalog threshold significance $\xi=4.5$, using the RM catalog to provide optical confirmation and redshifts for additional low-$\xi$ SPT-SZ candidates.In this way, we identify 15 additional clusters with $\xi\in [4,4.5]$ over the redshift regime explored by RM in the overlapping region between DES science verification data and the SPT-SZ survey.

We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.

We present the catalogue of blended galaxy spectra from the Galaxy And Mass Assembly (GAMA) survey. These are cases where light from two galaxies are significantly detected in a single GAMA fibre. Galaxy pairs identified from their blended spectrum fall into two principal classes: they are either strong lenses, a passive galaxy lensing an emission-line galaxy; or occulting galaxies, serendipitous overlaps of two galaxies, of any type. Blended spectra can thus be used to reliably identify strong lenses for follow-up observations (high resolution imaging) and occulting pairs, especially those that are a late-type partly obscuring an early-type galaxy which are of interest for the study of dust content of spiral and irregular galaxies. The GAMA survey setup and its autoz automated redshift determination were used to identify candidate blended galaxy spectra from the cross-correlation peaks. We identify 280 blended spectra with a minimum velocity separation of 600 km/s, of which 104 are lens pair candidates, 71 emission-line-passive pairs, 78 are pairs of emission-line galaxies and and 27 are pairs of galaxies with passive spectra. We have visually inspected the candidates in the Sloan Digital Sky Survey (SDSS) and Kilo Degree Survey (KiDS) images. Many blended objects are ellipticals with blue fuzz (Ef in our classification). These latter "Ef" classifications are candidates for possible strong lenses, massive ellipticals with an emission-line galaxy in one or more lensed images. The GAMA lens and occulting galaxy candidate samples are similar in size to those identified in the entire SDSS. This blended spectrum sample stands as a testament of the power of this highly complete, second-largest spectroscopic survey in existence and offers the possibility to expand e.g., strong gravitational lens surveys.

We report the discovery of eight new Milky Way companions in ~1,800 deg^2 of optical imaging data collected during the first year of the Dark Energy Survey (DES). Each system is identified as a statistically significant over-density of individual stars consistent with the expected isochrone and luminosity function of an old and metal-poor stellar population. The objects span a wide range of absolute magnitudes (M_V from -2.2 mag to -7.4 mag), physical sizes (10 pc to 170 pc), and heliocentric distances (30 kpc to 330 kpc). Based on the low surface brightnesses, large physical sizes, and/or large Galactocentric distances of these objects, several are likely to be new ultra-faint satellite galaxies of the Milky Way and/or Magellanic Clouds. We introduce a likelihood-based algorithm to search for and characterize stellar over-densities, as well as identify stars with high satellite membership probabilities. We also present completeness estimates for detecting ultra-faint galaxies of varying luminosities, sizes, and heliocentric distances in the first-year DES data.

OzDES is a five-year, 100-night, spectroscopic survey on the Anglo-Australian Telescope, whose primary aim is to measure redshifts of approximately 2,500 Type Ia supernovae host galaxies over the redshift range 0.1 < z < 1.2, and derive reverberation-mapped black hole masses for approximately 500 active galactic nuclei and quasars over 0.3 < z < 4.5. This treasure trove of data forms a major part of the spectroscopic follow-up for the Dark Energy Survey for which we are also targeting cluster galaxies, radio galaxies, strong lenses, and unidentified transients, as well as measuring luminous red galaxies and emission line galaxies to help calibrate photometric redshifts. Here we present an overview of the OzDES program and our first-year results. Between Dec 2012 and Dec 2013, we observed over 10,000 objects and measured more than 6,000 redshifts. Our strategy of retargeting faint objects across many observing runs has allowed us to measure redshifts for galaxies as faint as m_r=25 mag. We outline our target selection and observing strategy, quantify the redshift success rate for different types of targets, and discuss the implications for our main science goals. Finally, we highlight a few interesting objects as examples of the fortuitous yet not totally unexpected discoveries that can come from such a large spectroscopic survey.

We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.02<z<0.2. Far-IR spectra centred on this redshifted line were taken with the PACS instrument on-board the Herschel Space Observatory. The galaxies span 10<log(L_IR/Lo)<12 (where L_IR=L_IR[8-1000um]) and 7.3<log(L_[CII]/Lo)<9.3, covering a variety of optical galaxy morphologies. The sample exhibits the so-called [CII] deficit at high IR luminosities, i.e. L_[CII]/L_IR (hereafter [CII]/IR) decreases at high L_IR. We find significant differences between those galaxies presenting [CII]/IR>2.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5<S_12um/S_22um<1.0; (5) low surface brightness Sigma_IR~10^8-9 Lo kpc^-2, (6) and specific star-formation rates of sSFR~0.05-3 Gyr^-1. We suggest that the strength of the far-UV radiation fields (<G_O>) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high <G_O> creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11<log(L_IR/Lo)<12, the low [CII]/IR ratio is unlikely to be modified by [CII] self absorption or controlled by the presence of a moderately luminous AGN (identified via the BPT diagram).

The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes r$_{pet}$ < 19.4, stellar masses $10^{7} - 10^{12}$ M$_{sol}$, and environments from isolated field galaxies through groups to clusters of $10^{15}$ M$_{sol}$.

We present a methodology for the regularisation and combination of sparse sampled and irregularly gridded observations from fibre-optic multi-object integral-field spectroscopy. The approach minimises interpolation and retains image resolution on combining sub-pixel dithered data. We discuss the methodology in the context of the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey underway at the Anglo-Australian Telescope. The SAMI instrument uses 13 fibre bundles to perform high-multiplex integral-field spectroscopy across a one degree diameter field of view. The SAMI Galaxy Survey is targeting 3000 galaxies drawn from the full range of galaxy environments. We demonstrate the subcritical sampling of the seeing and incomplete fill factor for the integral-field bundles results in only a 10% degradation in the final image resolution recovered. We also implement a new methodology for tracking covariance between elements of the resulting datacubes which retains 90% of the covariance information while incurring only a modest increase in the survey data volume.

We present the Early Data Release of the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The SAMI Galaxy Survey is an ongoing integral field spectroscopic survey of ~3400 low-redshift (z<0.12) galaxies, covering galaxies in the field and in groups within the Galaxy And Mass Assembly (GAMA) survey regions, and a sample of galaxies in clusters. In the Early Data Release, we publicly release the fully calibrated datacubes for a representative selection of 107 galaxies drawn from the GAMA regions, along with information about these galaxies from the GAMA catalogues. All datacubes for the Early Data Release galaxies can be downloaded individually or as a set from the SAMI Galaxy Survey website. In this paper we also assess the quality of the pipeline used to reduce the SAMI data, giving metrics that quantify its performance at all stages in processing the raw data into calibrated datacubes. The pipeline gives excellent results throughout, with typical sky subtraction residuals in the continuum of 0.9-1.2 per cent, a relative flux calibration uncertainty of 4.1 per cent (systematic) plus 4.3 per cent (statistical), and atmospheric dispersion removed with an accuracy of 0."09, less than a fifth of a spaxel.

We use a highly complete subset of the GAMA-II redshift sample to fully describe the stellar mass dependence of close-pairs and mergers between 10^8 Msun and 10^12 Msun. Using the analytic form of this fit we investigate the total stellar mass accreting onto more massive galaxies across all mass ratios. Depending on how conservatively we select our robust merging systems, the fraction of mass merging onto more massive companions is 2.0%-5.6%. Using the GAMA-II data we see no significant evidence for a change in the close-pair fraction between redshift $z = 0.05-0.2$. However, we find a systematically higher fraction of galaxies in similar mass close-pairs compared to published results over a similar redshift baseline. Using a compendium of data and the function $\gamma_M =A(1+z)m$ to predict the major close-pair fraction, we find fitting parameters of $A = 0.021 \pm 0.001$ and $m = 1.53 \pm 0.08$, which represents a higher low-redshift normalisation and shallower power-law slope than recent literature values. We find that the relative importance of in-situ star-formation versus galaxy merging is inversely correlated, with star-formation dominating the addition of stellar material below Mstar and merger accretion events dominating beyond Mstar. We find mergers have a measurable impact on the whole extent of the GSMF, manifest as a deepening of the dip in the GSMF over the next Gyr and an increase in Mstar by as much as 0.01-0.05 dex.

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, {\alpha}1 = -0.43, {\phi}*1 = 4.18 dex^-1 Mpc^-3, {\alpha}2 = -1.50 and {\phi}*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.

We investigate the relationship between colour and structure within galaxies using a large, volume-limited sample of bright, low-redshift galaxies with optical to near-infrared imaging from the GAMA survey. We fit single-component, wavelength-dependent, elliptical Sersic models to all passbands simultaneously, using software developed by the MegaMorph project. Dividing our sample by Sersic index and colour, the recovered wavelength variations in effective radius (R_e) and Sersic index (n) reveal the internal structure, and hence formation history, of different types of galaxies. All these trends depend on n; some have an additional dependence on galaxy colour. Late-type galaxies (n_r < 2.5) show a dramatic increase in Sersic index with wavelength. This might be a result of their two-component (bulge-disk) nature, though stellar population gradients within each component and dust attenuation are likely to play a role. All galaxies show a substantial decrease in R_e with wavelength. This is strongest for early-types (n_r > 2.5), even though they maintain constant n with wavelength, revealing that ellipticals are a superimposition of different stellar populations associated with multiple collapse and merging events. Processes leading to structures with larger R_e must be associated with lower metallicity or younger stellar populations. This appears to rule out the formation of young cores through dissipative gas accretion as an important mechanism in the recent lives of luminous elliptical galaxies.

The Galaxy And Mass Assembly (GAMA) survey has obtained spectra of over 230000 targets using the Anglo-Australian Telescope. To homogenise the redshift measurements and improve the reliability, a fully automatic redshift code was developed (autoz). The measurements were made using a cross-correlation method for both absorption-line and emission-line spectra. Large deviations in the high-pass filtered spectra are partially clipped in order to be robust against uncorrected artefacts and to reduce the weight given to single-line matches. A single figure of merit (FOM) was developed that puts all template matches onto a similar confidence scale. The redshift confidence as a function of the FOM was fitted with a tanh function using a maximum likelihood method applied to repeat observations of targets. The method could be adapted to provide robust automatic redshifts for other large galaxy redshift surveys. For the GAMA survey, there was a substantial improvement in the reliability of assigned redshifts and in the lowering of redshift uncertainties with a median velocity uncertainty of 33 km/s.

The magnification effect of gravitational lensing is a powerful probe of the distribution of matter in the universe, yet it is frequently overlooked due to the fact that its signal to noise is smaller than that of lensing shear. Because its systematic errors are quite different from those of shear, magnification is nevertheless an important approach with which to study the distribution of large scale structure. We present lensing mass profiles of spectroscopic luminous red galaxies (LRGs) and galaxy clusters determined through measurements of the weak lensing magnification of photometric LRGs in their background. We measure the change in detected galaxy counts as well as the increased average galaxy flux behind the lenses. In addition, we examine the average change in source color due to extinction by dust in the lenses. By simultaneously fitting these three probes we constrain the mass profiles and dust-to-mass ratios of the lenses in six bins of lens richness. For each richness bin we fit an NFW halo mass, brightest cluster galaxy (BCG) mass, second halo term, and dust-to-mass ratio. The resulting mass-richness relation is consistent with previous analyses of the catalogs, and limits on the dust-to-mass ratio in the lenses are in agreement with expectations. We explore the effects of including the (low signal-to-noise) flux magnification and reddening measurements in the analysis compared to using only the counts magnification data; the additional probes significantly improve the agreement between our measured mass-richness relation and previous results.

We have analysed the growth of Brightest Group Galaxies and Brightest Cluster Galaxies (BGGs/BCGs) over the last 3 billion years using a large sample of 883 galaxies from the Galaxy And Mass Assembly Survey. By comparing the stellar mass of BGGs and BCGs in groups and clusters of similar dynamical masses, we find no significant growth between redshift $z=0.27$ and $z=0.09$. We also examine the number of BGGs/BCGs that have line emission, finding that approximately 65 per cent of BGGs/BCGs show H$\alpha$ in emission. From the galaxies where the necessary spectroscopic lines were accurately recovered (54 per cent of the sample), we find that half of this (i.e. 27 per cent of the sample) harbour on-going star formation with rates up to $10\,$M$_{\odot}$yr$^{-1}$, and the other half (i.e. 27 per cent of the sample) have an active nucleus (AGN) at the centre. BGGs are more likely to have ongoing star formation, while BCGs show a higher fraction of AGN activity. By examining the position of the BGGs/BCGs with respect to their host dark matter halo we find that around 13 per cent of them do not lie at the centre of the dark matter halo. This could be an indicator of recent cluster-cluster mergers. We conclude that BGGs and BCGs acquired their stellar mass rapidly at higher redshifts as predicted by semi-analytic models, mildly slowing down at low redshifts.

The Galaxy And Mass Assembly (GAMA) survey furnishes a deep redshift catalog that, when combined with the Wide-field Infrared Explorer ($WISE$), allows us to explore for the first time the mid-infrared properties of $> 110, 000$ galaxies over 120 deg$^2$ to $z\simeq 0.5$. In this paper we detail the procedure for producing the matched GAMA-$WISE$ catalog for the G12 and G15 fields, in particular characterising and measuring resolved sources; the complete catalogs for all three GAMA equatorial fields will be made available through the GAMA public releases. The wealth of multiwavelength photometry and optical spectroscopy allows us to explore empirical relations between optically determined stellar mass (derived from synthetic stellar population models) and 3.4micron and 4.6micron WISE measurements. Similarly dust-corrected Halpha-derived star formation rates can be compared to 12micron and 22micron luminosities to quantify correlations that can be applied to large samples to $z<0.5$. To illustrate the applications of these relations, we use the 12micron star formation prescription to investigate the behavior of specific star formation within the GAMA-WISE sample and underscore the ability of WISE to detect star-forming systems at $z\sim0.5$. Within galaxy groups (determined by a sophisticated friends-of-friends scheme), results suggest that galaxies with a neighbor within 100$\,h^{-1} $kpc have, on average, lower specific star formation rates than typical GAMA galaxies with the same stellar mass.

Photometric galaxy surveys are an essential tool to further our understanding of the large-scale structure of the universe, its matter and energy content and its evolution. These surveys necessitate the determination of the galaxy redshifts using photometric techniques (photo-z). Oftentimes, it is advantageous to remove from the galaxy sample those for which one suspects that the photo-z estimation might be unreliable. In this paper, we show that applying these photo-z quality cuts blindly can grossly bias the measured galaxy correlations within and across photometric redshift bins. We then extend the work of Ho et al. (2012) and Ross et al. (2011) to develop a simple and effective method to correct for this using the data themselves. Finally, we apply the method to the Mega-Z catalog, containing about a million luminous red galaxies in the redshift range 0.45 < z < 0.65. After splitting the sample into four \Delta z = 0.05 photo-z bins using the BPZ algorithm, we see how our corrections bring the measured galaxy auto- and cross-correlations into agreement with expectations. We then look for the BAO feature in the four bins, with and without applying the photo-z quality cuts, and find a broad agreement between the BAO scales extracted in both cases. Intriguingly, we observe a correlation between galaxy density and photo-z quality even before any photo-z quality cuts are applied. This may be due to uncorrected observational effects that result in correlated gradients across the sky of the galaxy density and the galaxy photo-z precision. Our correction procedure could also help to mitigate some of these systematic effects.

From a volume limited sample of 45,542 galaxies and 6,000 groups with $z \leq 0.213$ we use an adapted minimal spanning tree algorithm to identify and classify large scale structures within the Galaxy and Mass Assembly (GAMA) survey. Using galaxy groups, we identify 643 filaments across the three equatorial GAMA fields that span up to 200 $h^{-1}$ Mpc in length, each with an average of 8 groups within them. By analysing galaxies not belonging to groups we identify a secondary population of smaller coherent structures composed entirely of galaxies, dubbed `tendrils' that appear to link filaments together, or penetrate into voids, generally measuring around 10 $h^{-1}$ Mpc in length and containing on average 6 galaxies. Finally we are also able to identify a population of isolated void galaxies. By running this algorithm on GAMA mock galaxy catalogues we compare the characteristics of large scale structure between observed and mock data; finding that mock filaments reproduce observed ones extremely well. This provides a probe of higher order distribution statistics not captured by the popularly used two-point correlation function.

We present observations of 18 galaxies from the Galaxy And Mass Assembly (GAMA) survey made with the SPIRAL optical integral field unit (IFU) on the Anglo-Australian Telescope. The galaxies are selected to have a narrow range in stellar mass (6x10^9Msolar < M* <2x10^10 Msolar) in order to focus on the effects of environment. Local galaxy environments are measured quantitatively using 5th nearest neighbour surface densities. We find that the total star formation rates (SFR) measured from the IFU data are consistent with total SFRs measured from aperture correcting either GAMA or Sloan Digital Sky Survey single-fibre observations. The mean differences are SFR_GAMA/SFR_IFU = 1.26+/-0.23, sigma=0.90 and for the Sloan Digital Sky Survey we similarly find SFR_Brinchmann/SFR_IFU = 1.34+/-0.17, sigma=0.67. Examining the relationships with environment, we find off-centre and clumpy Halpha emission is not significantly dependent on environment, being present in 2/7 (29^+20_-11 per cent) galaxies in high-density environments (>0.77 Mpc^-2), and 5/11 (45^+15_-13 per cent) galaxies in low-density environments (<0.77 Mpc^-2). We find a weak but not significant relationship of the total star formation rates of star-forming galaxies with environment. Due to the size of our sample and the scatter observed we do not draw a definitive conclusion about a possible SFR dependence on environment. Examining the spatial distribution of the Halpha emission, we find no evidence for a change in shape or amplitude of the radial profile of star-forming galaxies with environment. If these observations are borne out in larger samples this would infer that any environment-driven star-formation suppression must either act very rapidly (the `infall-and-quench' model) or that galaxies must evolve in a density-dependent manner (an `in-situ evolution' model).

We present a combined optical and X-ray analysis of the rich cluster Abell 1882 with the aim of identifying merging substructure and understanding the recent assembly history of this system. Our optical data consist of spectra drawn from the Galaxy and Mass Assembly (GAMA) survey, which lends itself to this kind of detailed study thanks to its depth and high spectroscopic completeness. We use 283 spectroscopically confirmed cluster members to detect and characterize substructure. We complement the optical data with X-ray data taken with both Chandra and XMM. Our analysis reveals that A1882 harbors two main components, A1882A and A1882B, which have a projected separation of 2Mpc and a line of sight velocity difference of v_{los}=-428km/s. The primary system, A1882A, has velocity dispersion sigma_v=500km/s and Chandra (XMM) temperature kT=3.57keV (3.31keV) while the secondary, A1882B, has sigma_v=457km/s and Chandra (XMM) temperature kT=2.39keV (2.12keV). The optical and X-ray estimates for the masses of the two systems are consistent within the uncertainties and indicate that there is twice as much mass in A1882A (M_{500}=1.5-1.9x10^{14}Msolar) when compared with A1882B (M_{500}=0.8-1.0x10^{14}Msolar). We interpret the A1882A/A1882B system as being observed prior to a core passage. Supporting this interpretation is the large projected separation of A1882A and A1882B and the dearth of evidence for a recent (<2Gyr) major interaction in the X-ray data. Two-body analyses indicate that A1882A and A1882B form a bound system with bound incoming solutions strongly favored. We compute blue fractions of f_b=0.28 and 0.18 for the spectroscopically confirmed member galaxies within r_{500} of the centers of A1882A and A1882B, respectively. These blue fractions do not differ significantly from the blue fraction measured from an ensemble of 20 clusters with similar mass and redshift.

We present evidence for stochastic star formation histories in low-mass (M* < 10^10 Msun) galaxies from observations within the Galaxy And Mass Assembly (GAMA) survey. For ~73,000 galaxies between 0.05<z<0.32, we calculate star formation rates (SFR) and specific star formation rates (SSFR = SFR/M*) from spectroscopic Halpha measurements and apply dust corrections derived from Balmer decrements. We find a dependence of SSFR on stellar mass, such that SSFRs decrease with increasing stellar mass for star-forming galaxies, and for the full sample, SSFRs decrease as a stronger function of stellar mass. We use simple parametrizations of exponentially declining star formation histories to investigate the dependence on stellar mass of the star formation timescale and the formation redshift. We find that parametrizations previously fit to samples of z~1 galaxies cannot recover the distributions of SSFRs and stellar masses observed in the GAMA sample between 0.05<z<0.32. In particular, a large number of low-mass (M* < 10^10 Msun) galaxies are observed to have much higher SSFRs than can be explained by these simple models over the redshift range of 0.05<z<0.32, even when invoking mass-dependent staged evolution. For such a large number of galaxies to maintain low stellar masses, yet harbour such high SSFRs, requires the late onset of a weak underlying exponentially declining SFH with stochastic bursts of star formation superimposed.

A full appreciation of the role played by gas metallicity (Z), star-formation rate (SFR), and stellar mass is fundamental to understanding how galaxies form and evolve. The connections between these three parameters at different redshifts significantly affect galaxy evolution, and thus provide important constraints for galaxy evolution models. Using data from the Sloan Digital Sky Survey-Data Release 7 (SDSS-DR7) and the Galaxy and Mass Assembly (GAMA) surveys we study the relationships and dependencies between SFR, Z, and stellar mass, as well as the Fundamental Plane for star-forming galaxies. We combine both surveys using volume-limited samples up to a redshift of z ~ 0.36. The GAMA and SDSS surveys complement each other when analyzing the relationships between SFR, Mass and Z. We present evidence for SFR and metallicity evolution to z ~ 0.2. We study the dependencies between SFR, Mass, Z, and specific star-formation rate (SSFR) on the M-Z, M-SFR, M-SSFR, Z-SFR, and Z-SSFR relations, finding strong correlations between all. Based on those dependencies, we propose a simple model that allows us to explain the different behaviour observed between low and high mass galaxies. Finally, our analysis allows us to confirm the existence of a Fundamental Plane, for which Mass=f(Z, SFR) in star-forming galaxies.

Measurements of the low-z Halpha luminosity function have a large dispersion in the local number density of sources, and correspondingly in the SFR density. The possible causes for these discrepancies include limited volume sampling, biases arising from survey sample selection, different methods of correcting for dust obscuration and AGN contamination. The Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey (SDSS) provide deep spectroscopic observations over a wide sky area enabling detection of a large sample of star-forming galaxies spanning 0.001<SFR(Halpha)<100 with which to robustly measure the evolution of the SFR density in the low-z universe. The large number of high SFR galaxies present in our sample allow an improved measurement of the bright end of the luminosity function, indicating that the decrease in number density of sources at bright luminosities is best described by a Saunders functional form rather than the traditional Schechter function. This result is consistent with other published luminosity functions in the FIR and radio. For GAMA and SDSS we find the r-band apparent magnitude limit, combined with the subsequent requirement for Halpha detection leads to an incompleteness due to missing bright Halpha sources with faint r-band magnitudes.

Recent independent results from numerical simulations and observations have shown that brightest cluster galaxies (BCGs) have increased their stellar mass by a factor of almost two between z~0.9 and z~0.2. The numerical simulations further suggest that more than half this mass is accreted through major mergers. Using a sample of 18 distant galaxy clusters with over 600 spectroscopically confirmed cluster members between them, we search for observational evidence that major mergers do play a significant role. We find a major merger rate of 0.38 +/- 0.14 mergers per Gyr at z~1. While the uncertainties, which stem from the small size of our sample, are relatively large, our rate is consistent with the results that are derived from numerical simulations. If we assume that this rate continues to the present day and that half of the mass of the companion is accreted onto the BCG during these mergers, then we find that this rate can explain the growth in the stellar mass of the BCGs that is observed and predicted by simulations. Major mergers therefore appear to be playing an important role, perhaps even the dominant one, in the build up of stellar mass in these extraordinary galaxies.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

In this work we investigate in detail the effects local environment (groups and pairs) has on galaxies with stellar mass similar to the Milky-Way (L* galaxies). A volume limited sample of 6,150 galaxies is classified to determine emission features, morphological type and presence of a disk. This sample allows for characteristics of galaxies to be isolated (e.g. stellar mass and group halo mass), and their codependencies determined. We observe that galaxy-galaxy interactions play the most important role in shaping the evolution within a group halo, the main role of halo mass is in gathering the galaxies together to encourage such interactions. Dominant pair galaxies find their overall star formation enhanced when the pair's mass ratio is close to 1, otherwise we observe the same galaxies as we would in an unpaired system. The minor galaxy in a pair is greatly affected by its companion galaxy, and whilst the star forming fraction is always suppressed relative to equivalent stellar mass unpaired galaxies, it becomes lower still when the mass ratio of a pair system increases. We find that, in general, the close galaxy-galaxy interaction rate drops as a function of halo mass for a given amount of stellar mass. We find evidence of a local peak of interactions for Milky-Way stellar mass galaxies in Milky-Way halo mass groups. Low mass halos, and in particular Local Group mass halos, are an important environment for understanding the typical evolutionary path of a unit of stellar mass. We find compelling evidence for galaxy conformity in both groups and pairs, where morphological type conformity is dominant in groups, and emission class conformity is dominant in pairs. This suggests that group scale conformity is the result of many galaxy encounters over an extended period of time, whilst pair conformity is a fairly instantaneous response to a transitory interaction.

The goal of the .Astronomy conference series is to bring together astronomers, educators, developers and others interested in using the Internet as a medium for astronomy. Attendance at the event is limited to approximately 50 participants, and days are split into mornings of scheduled talks, followed by 'unconference' afternoons, where sessions are defined by participants during the course of the event. Participants in unconference sessions are discouraged from formal presentations, with discussion, workshop-style formats or informal practical tutorials encouraged. The conference also designates one day as a 'hack day', in which attendees collaborate in groups on day-long projects for presentation the following morning. These hacks are often a way of concentrating effort, learning new skills, and exploring ideas in a practical fashion. The emphasis on informal, focused interaction makes recording proceedings more difficult than for a normal meeting. While the first .Astronomy conference is preserved formally in a book, more recent iterations are not documented. We therefore, in the spirit of .Astronomy, report 'unproceedings' from .Astronomy 4, which was held in Heidelberg in July 2012.

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.

We present the results of an extended narrow-band H{\alpha} study of the massive galaxy cluster XMMU J2235.3-2557 at z = 1.39. This paper represents a follow up study to our previous investigation of star-formation in the cluster centre, extending our analysis out to a projected cluster radius of 1.5 Mpc. Using the Near InfraRed Imager and Spectrograph (NIRI) on Gemini North we obtained deep H narrow-band imaging corresponding to the rest-frame wavelength of H{\alpha} at the cluster's redshift. We identify a total of 163 potential cluster members in both pointings, excluding stars based on their near-IR colours derived from VLT/HAWK-I imaging. Of these 163 objects 14 are spectroscopically confirmed cluster members, and 20% are excess line-emitters. We find no evidence of star formation activity within a radius of 200 kpc of the brightest cluster galaxy in the cluster core. Dust-corrected star formation rates (SFR) of excess emitters outside this cluster quenching radius, RQ \sim 200 kpc, are on average <SFR> = 2.7 \pm 1.0 M\odot yr-1, but do not show evidence of increasing star-formation rates toward the extreme 1.5 Mpc radius of the cluster. No individual cluster galaxy exceeds an SFR of 6 M\odot yr-1 . Massive galaxies (log M\ast /M\odot > 10.75) all have low specific SFRs (SSFRs, i.e. SFR per unit stellar mass). At fixed stellar mass, galaxies in the cluster centre have lower SSFRs than the rest of the cluster galaxies, which in turn have lower SSFRs than field galaxies at the same redshift by a factor of a few to 10. For the first time we can demonstrate through measurements of individual SFRs that already at very early epochs (at an age of the Universe of \sim4.5 Gyr) the suppression of star-formation is an effect of the cluster environment which persists at fixed galaxy stellar mass. [Erratum added after the original paper]

We study the SLACS strong lensing system SDSSJ1430+4105 at z=0.285. The lensed source (z=0.575) of this system has a complex morphology with several subcomponents. Its subcomponents span a radial range from 4 to 10kpc in the lens plane. Therefore we can constrain the slope of the total projected mass profile around the Einstein radius (R_E) from lensing alone. We measure a density profile that is slightly but not significantly shallower than isothermal at R_E. We decompose the mass of the lensing galaxy into a de Vaucouleurs (deV) component to trace the stars and an additional dark component. The spread of multiple image components over a large radial range also allows us to determine the amplitude of the deV and dark matter components separately. We get a mass to light ratio of M_deV/L_B~5.5\pm1.5M/L_sun,B and a dark matter fraction within R_E of ~20 to 40%. Modelling the star formation history assuming composite stellar populations at solar metallicity to the galaxy's photometry yields a mass to light ratio of M_star,salp/L_B~4.0_{-1.3}^{+0.6}M/L_sun,B and M_star,chab/L_B~2.3_{-0.8}^{+0.3}M/L_sun,B for Salpeter and Chabrier IMFs, respectively. Hence, the mass to light ratio derived from lensing is more Salpeter-like, in agreement with results for massive Coma galaxies and other nearby massive early type galaxies. We examine the consequences of the galaxy group in which the lensing galaxy is embedded, showing that it has little influence on the mass to light ratio obtained for the deV component of the lensing galaxy. Finally, we decompose the projected, azimuthally averaged 2D density distribution of the deV and dark matter component of the lensing signal into spherically averaged 3D density profiles. We can show that the 3D dark and luminous matter density within R_E~0.6R_eff of this SLACS galaxy is similar to the values of Coma galaxies with the same velocity dispersions.

Context: The mass-metallicity relationship (MMR) of star-forming galaxies is well-established, however there is still some disagreement with respect to its exact shape and its possible dependence on other observables. Aims: We measure the MMR in the Galaxy And Mass Assembly (GAMA) survey. We compare our measured MMR to that measured in the Sloan Digital Sky Survey (SDSS) and study the dependence of the MMR on various selection criteria to identify potential causes for disparities seen in the literature. Methods: We use strong emission line ratio diagnostics to derive oxygen abundances. We then apply a range of selection criteria for the minimum signal-to-noise in various emission lines, as well as the apparent and absolute magnitude to study variations in the inferred MMR. Results: The shape and position of the MMR can differ significantly depending on the metallicity calibration and selection used. After selecting a robust metallicity calibration amongst those tested, we find that the mass-metallicity relation for redshifts 0.061< z<0.35 in GAMA is in reasonable agreement with that found in the SDSS despite the difference in the luminosity range probed. Conclusions: In view of the significant variations of the MMR brought about by reasonable changes in the sample selection criteria and method, we recommend that care be taken when comparing the MMR from different surveys and studies directly. We also conclude that there could be a modest level of evolution over 0.06<z<0.35 within the GAMA sample.

Analysing all Galaxy and Mass Assembly (GAMA) galaxies within a factor two (+/- 0.3 dex) of the stellar mass of the Milky Way (MW), there is a 11.9% chance that one of these galaxies will have a close companion (within a projected separation of 70 kpc and radial separation of 400 km/s) that is at least as massive as the Large Magellanic Cloud (LMC). Two close companions at least as massive as the Small Magellanic Cloud (SMC) are rare at the 3.4% level. Two full analogues to the MW- LMC-SMC system were found in GAMA (all galaxies late-type and star forming), suggesting such a combination of close together, late-type, star-forming galaxies is rare: only 0.4% of MW mass galaxies (in the range where we could observe both the LMC and SMC) have such a system. In summary, the MW-LMC-SMC system is a 2.7? event (when recast into Gaussian statistics). Using cross-correlation comparisons we find that there is a preference for SMC- LMC binary pair analogues to be located within 2 Mpc of a range of different lumi- nosity groups. There is a particular preference is for such binaries to be located near LG luminosity systems. When these groups are subdivided into small magnitude gap and large magnitude gap subsets, the binaries prefer to be spatially associated with the small magnitude gap systems. These systems will be dynamically less evolved, but still offer the same amount of gravitational dark matter. This suggests that binaries such as the SMC-LMC might be transient systems, usually destroyed during vigorous merger events. Details of a particularly striking analogue to the MW-SMC-LMC and M31 complex are included.

We present a detailed investigation into the effects of galaxy environment on their star formation rates (SFR) using galaxies observed in the Galaxy and Mass Assembly Survey (GAMA). We use three independent volume-limited samples of galaxies within z < 0.2 and Mr < -17.8. We investigate the known SFR-density relationship and explore in detail the dependence of SFR on stellar mass and density. We show that the SFR-density trend is only visible when we include the passive galaxy population along with the star-forming population. This SFR-density relation is absent when we consider only the star-forming population of galaxies, consistent with previous work. While there is a strong dependence of the EWH?a on density we find, as in previous studies, that these trends are largely due to the passive galaxy population and this relationship is absent when considering a "star-forming" sample of galaxies. We find that stellar mass has the strongest influence on SFR and EWH?a with the environment having no significant effect on the star-formation properties of the star forming population. We also show that the SFR-density relationship is absent for both early and late-type star-forming galaxies. We conclude that the stellar mass has the largest impact on the current SFR of a galaxy, and any environmental effect is not detectable. The observation that the trends with density are due to the changing morphology fraction with density implies that the timescales must be very short for any quenching of the SFR in infalling galaxies. Alternatively galaxies may in fact undergo predominantly in-situ evolution where the infall and quenching of galaxies from the field into dense environments is not the dominant evolutionary mode.

Highly supersonic, compressible turbulence is thought to be of tantamount importance for star formation processes in the interstellar medium. Likewise, cosmic structure formation is expected to give rise to subsonic turbulence in the intergalactic medium, which may substantially modify the thermodynamic structure of gas in virialized dark matter halos and affect small-scale mixing processes in the gas. Numerical simulations have played a key role in characterizing the properties of astrophysical turbulence, but thus far systematic code comparisons have been restricted to the supersonic regime, leaving it unclear whether subsonic turbulence is faithfully represented by the numerical techniques commonly employed in astrophysics. Here we focus on comparing the accuracy of smoothed particle hydrodynamics (SPH) and our new moving-mesh technique AREPO in simulations of driven subsonic turbulence. To make contact with previous results, we also analyze simulations of transsonic and highly supersonic turbulence. We find that the widely employed standard formulation of SPH yields problematic results in the subsonic regime. Instead of building up a Kolmogorov-like turbulent cascade, large-scale eddies are quickly damped close to the driving scale and decay into small-scale velocity noise. Reduced viscosity settings improve the situation, but the shape of the dissipation range differs compared with expectations for a Kolmogorov cascade. In contrast, our moving-mesh technique does yield power-law scaling laws for the power spectra of velocity, vorticity and density, consistent with expectations for fully developed isotropic turbulence. We show that large errors in SPH's gradient estimate and the associated subsonic velocity noise are ultimately responsible for producing inaccurate results in the subsonic regime. In contrast, SPH's performance is much better for supersonic turbulence. [Abridged]

Accurate measurement of galaxy cluster masses is an essential component not only in studies of cluster physics, but also for probes of cosmology. However, different mass measurement techniques frequently yield discrepant results. The SDSS MaxBCG catalog's mass-richness relation has previously been constrained using weak lensing shear, Sunyaev-Zeldovich (SZ), and X-ray measurements. The mass normalization of the clusters as measured by weak lensing shear is >~25% higher than that measured using SZ and X-ray methods, a difference much larger than the stated measurement errors in the analyses. We constrain the mass-richness relation of the MaxBCG galaxy cluster catalog by measuring the gravitational lensing magnification of type I quasars in the background of the clusters. The magnification is determined using the quasars' variability and the correlation between quasars' variability amplitude and intrinsic luminosity. The mass-richness relation determined through magnification is in agreement with that measured using shear, confirming that the lensing strength of the clusters implies a high mass normalization, and that the discrepancy with other methods is not due to a shear-related systematic measurement error. We study the dependence of the measured mass normalization on the cluster halo orientation. As expected, line-of-sight clusters yield a higher normalization; however, this minority of haloes does not significantly bias the average mass-richness relation of the catalog.

We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (M\ast \approx 2 \times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalog on maps: at 24 {\mu}m (Spitzer/MIPS); 70, 100, and 160{\mu}m (Herschel/PACS); 250, 350, 500{\mu}m (BLAST); and 870{\mu}m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n \leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

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.

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

We present a near-infrared spectroscopic study of a stellar mass selected sample of galaxies at z~1 utilising the LIRIS multi-object spectrograph on the WHT. We detect continuum, and the H alpha line for our sample, which is one of the better direct tracers of star formation in external galaxies. We spectroscopically measure the H alpha emission from 41 massive (M_{*}>10^{10.5} Msol) galaxies taken from the POWIR Survey with spectroscopic redshifts 0.4<z_{spec}<1.4. We correct our H alpha fluxes for dust extinction by using multi-wavelength data, and investigate SFR trends with mass and colour. We find a drop in the fraction of massive galaxies with M_{*}>10^{11} Msol which are detected in H alpha emission at z<0.9. We furthermore find that the fraction of galaxies with H alpha emission drops steadily and significantly with redder (U-B) colours at z~1, and that the SSFR drops with increasing (U-B) colour for galaxies at all masses. By investigating the SFR-mass relation we find that the SFR is roughly constant with mass, in possible contrast to previous work, and that the specific star formation rate (SSFR) is lower in the most massive galaxies. The scatter in the SFR vs. mass relationship is very small for those systems with ongoing star formation which suggests that star formation in the most massive galaxies at z~1 shuts off rather abruptly over <1 Gyr, without an obvious gradual decline. We furthermore investigate the SFR as a function of (U-B) colour divided into different mass bins, revealing a tracer of the epoch of transition from star forming to passive, as a form of star formation "downsizing". This suggests that the shut off of star formation occurs before the change in a galaxy's colour. We find that galaxy stellar mass is the primary driving mechanisms behind the star formation history for these galaxies and discuss several possible mechanisms for regulating this process.

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_e<=2 kpc), compared to less than 1% at z~0. Furthermore, most (~65%) systems at z=2-3 have a low Sersic index n<=2, compared to ~13% at z~0. We present evidence that the n<=2 systems at z=2-3 likely contain prominent disks, unlike most massive z~0 systems. (2) There is a correlation between structure and star formation rates (SFR). The majority (~85%) of non-AGN massive galaxies at z=2-3, with SFR high enough to yield a 5 sigma (30 micro Jy) 24 micron Spitzer detection have low n<=2. Such n<=2 systems host the highest SFR. (3) The frequency of AGN is ~40% at z=2-3. Most (~65%) AGN hosts have disky (n<=2) morphologies. Ultra-compact galaxies appear quiescent in terms of both AGN activity and star formation. (4) Large stellar surface densities imply massive galaxies at z=2-3 formed via rapid, highly dissipative events at z>2. The large fraction of n<=2 disky systems suggests cold mode accretion complements gas-rich major mergers at z>2. 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.

Magnetic fields play an important role in astrophysics on a wide variety of scales, ranging from the Sun and compact objects to galaxies and galaxy clusters. Here we discuss a novel implementation of ideal magnetohydrodynamics (MHD) in the moving mesh code AREPO which combines many of the advantages of Eulerian and Lagrangian methods in a single computational technique. The employed grid is defined as the Voronoi tessellation of a set of mesh-generating points which can move along with the flow, yielding an automatic adaptivity of the mesh and a substantial reduction of advection errors. Our scheme solves the MHD Riemann problem in the rest frame of the Voronoi interfaces using the HLLD Riemann solver. To satisfy the divergence constraint of the magnetic field in multiple dimensions, the Dedner divergence cleaning method is applied. In a set of standard test problems we show that the new code produces accurate results, and that the divergence of the magnetic field is kept sufficiently small to closely preserve the correct physical solution. We also apply the code to two first application problems, namely supersonic MHD turbulence and the spherical collapse of a magnetized cloud. We verify that the code is able to handle both problems well, demonstrating the applicability of this MHD version of AREPO to a wide range of problems in astrophysics.

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-infra-red 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 UV-fluxes. Due to our imaging depth we can examine galaxies down to a colour-independent stellar mass completeness limit of log M\ast = 9.5 M\odot 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 in the highest relative over-densities. Galaxies in over-densities 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 AGN 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.

Here we report on the Women in Astronomy Workshop (http://asawomeninastronomy.org/meetings/wia2011/), which was held on 13 May 2011 in Sydney, Australia. The workshop was organised by the Astronomical Society of Australia's Chapter on Women in Astronomy, to discuss some of the issues that face women in astronomy and make recommendations to help support the success of women in Australian astronomy but came to broader conclusions that have value for the whole astronomical community. The workshop consisted of four sessions, with presentations by invited speakers on demographics, leadership, varied career paths, and how institutions & individuals can help. The workshop ended with a discussion panel that summarised the day's debate and presented a list of recommendations for the Australian astronomical community (both individuals and institutions) that are provided in this report.

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

We introduce a technique to measure gravitational lensing magnification using the variability of type I quasars. Quasars' variability amplitudes and luminosities are tightly correlated, on average. Magnification due to gravitational lensing increases the quasars' apparent luminosity, while leaving the variability amplitude unchanged. Therefore, the mean magnification of an ensemble of quasars can be measured through the mean shift in the variability-luminosity relation. As a proof of principle, we use this technique to measure the magnification of quasars spectroscopically identified in the Sloan Digital Sky Survey, due to gravitational lensing by galaxy clusters in the SDSS MaxBCG catalog. The Palomar-QUEST Variability Survey, reduced using the DeepSky pipeline, provides variability data for the sources. We measure the average quasar magnification as a function of scaled distance (r/R200) from the nearest cluster; our measurements are consistent with expectations assuming NFW cluster profiles, particularly after accounting for the known uncertainty in the clusters' centers. Variability-based lensing measurements are a valuable complement to shape-based techniques because their systematic errors are very different, and also because the variability measurements are amenable to photometric errors of a few percent and to depths seen in current wide-field surveys. Given the data volume expected from current and upcoming surveys, this new technique has the potential to be competitive with weak lensing shear measurements of large scale structure.

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<z<3. (Abridged)

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

We study the relationship between galaxy colour, stellar mass, and local galaxy density in a deep near-infrared imaging survey up to a redshift of z~3 using the GOODS NICMOS Survey (GNS). The GNS is a deep near-infrared Hubble Space Telescope survey imaging a total of 45 square arcminutes of the GOODS fields, reaching a stellar mass completeness limit of log M* = 9.5 M_sun at z=3. Using this data we measure galaxy local densities based on galaxy counts within a fixed aperture, as well as the distance to the 3rd, 5th and 7th nearest neighbour. We compare the average rest-frame (U-B) colour and fraction of blue galaxies in different local densities and at different stellar masses. We find a strong correlation between colour and stellar mass at all redshifts up to z~3. Massive red galaxies are already in place at z~3 at the expected location of the red-sequence in the colour-magnitude diagram, although they are star forming. We do not find a strong correlation between colour and local density, however, there may be evidence that the highest overdensities are populated by a higher fraction of blue galaxies than average or underdense areas. This could indicating that the colour-density relation at high redshift is reversed with respect to lower redshifts (z<1), where higher densities are found to have lower blue fractions. Our data suggests that the possible higher blue fraction at extreme overdensities might be due to a lack of massive red galaxies at the highest local densities.

We present the first results of a narrow-band photometric study of the massive galaxy cluster XMMU J2235.3-2557 at z=1.39. We obtained deep $H$ narrow-band imaging with NIRI on Gemini North, corresponding to H-alpha emission at the cluster's redshift. Our sample consists of 82 galaxies within a radius of ~500 kpc, ten of which are spectroscopically confirmed cluster members. Sixteen galaxies are identified as excess line-emitters. Among just the excess line-emitting galaxies we find an average SFR of 3.6 +/- 1.3 Msun/yr. For spectroscopically confirmed cluster members we find a correlation between H broad-band magnitude and SFR such that brighter galaxies have lower SFRs. The probability that SFR and magnitude of confirmed members are uncorrelated is 0.7%. We also find a correlation between SFR and distance from the cluster centre for both confirmed and excess line-emitting candidate members, with a probability of 5% for there to be no correlation among confirmed members. All excess line-emitting candidate cluster members are located outside a radius of 200 kpc. We conclude that star formation is effectively shut off within the central 200 kpc radius (R_QUENCH ~ 200 kpc) of this massive galaxy cluster at z=1.39, when the universe was only 4.5 Gyr old.

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 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 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.

We present a lensing and photometric study of the exceptional system SDSS J1538+5817, identified by the SLACS survey. The lens is a luminous elliptical at redshift z=0.143. Using HST public images in two different filters, the presence of two background sources lensed into an Einstein ring and a double system is ascertained. Our new spectroscopic observations, performed at the NOT, reveal that the two sources are located at the same redshift z=0.531. We investigate the total mass distribution of the lens between 1 and 4 kpc from the galaxy center by means of parametric and non-parametric lensing codes that describe the multiple images as point-like objects. Several disparate lensing models agree on: (1) reproducing accurately the observed image positions; (2) predicting a nearly axisymmetric total mass distribution, centered and oriented as the light distribution; (3) measuring a value of 8.11 x 10^{10} M_{Sun} for the total mass projected within the Einstein radius of 2.5 kpc; (4) estimating a total mass density profile slightly steeper than an isothermal one. A fit of the SDSS multicolor photometry with CSP models provides a value of 20 x 10^{10} M_{Sun} for the total stellar mass of the galaxy and of 0.9 for the fraction of projected luminous over total mass enclosed inside the Einstein radius. By combining lensing and photometric mass measurements, we differentiate the lens mass content in terms of luminous and dark matter components. This two-component modeling, which is viable only in extraordinary systems like SDSS J1538+5817, leads to a description of the global properties of the galaxy dark matter halo. Extending these results to a larger number of lenses would improve considerably our understanding of galaxy formation and evolution processes in the LCDM scenario.

We identify 3,113 highly variable objects in 7,200 square degrees of the Palomar-QUEST Survey, which each varied by more than 0.4 magnitudes simultaneously in two broadband optical filters on timescales from hours to roughly 3.5 years. The primary goal of the selection is to find blazars by their well-known violent optical variability. Because most known blazars have been found in radio and/or X-ray wavelengths, a sample discovered through optical variability may have very different selection effects, elucidating the range of behavior possible in these systems. A set of blazars selected in this unusual manner will improve our understanding of the physics behind this extremely variable and diverse class of AGN. The object positions, variability statistics, and color information are available using the Palomar-QUEST CasJobs server. The time domain is just beginning to be explored over large sky areas; we do not know exactly what a violently variable sample will hold. About 20% of the sample has been classified in the literature; over 70% of those objects are known or likely AGN. The remainder largely consists of a variety of variable stars, including a number of RR Lyrae and cataclysmic variables.

We study the ensemble optical variability of 276 FSRQs and 86 BL Lacs in the Palomar-QUEST Survey with the goal of searching for common fluctuation properties, examining the range of behavior across the sample, and characterizing the appearance of blazars in such a survey so that future work can more easily identify such objects. The survey, which covers 15,000 square degrees multiple times over 3.5 years, allows for the first ensemble blazar study of this scale. Variability amplitude distributions are shown for the FSRQ and BL Lac samples for numerous time lags, and also studied through structure function analyses. Individual blazars show a wide range of variability amplitudes, timescales, and duty cycles. Of the best sampled objects, 35% are seen to vary by more than 0.4 magnitudes; for these, the fraction of measurements contributing to the high amplitude variability ranges constantly from about 5% to 80%. Blazar variability has some similarities to that of type I quasars but includes larger amplitude fluctuations on all timescales. FSRQ variability amplitudes are particularly similar to those of QSOs on timescales of several months, suggesting significant contributions from the accretion disk to the variable flux at these timescales. Optical variability amplitudes are correlated with the maximum apparent velocities of the radio jet for the subset of FSRQs with MOJAVE VLBA measurements, implying that the optically variable flux's strength is typically related to that of the radio emission. We also study CRATES radio-selected FSRQ candidates, which show similar variability characteristics to known FSRQs; this suggests a high purity for the CRATES sample.

The ensemble variability properties of nearly 23,000 quasars are studied using the Palomar-QUEST Survey. The survey has covered 15,000 square degrees multiple times over 3.5 years using 7 optical filters, and has been calibrated specifically for variability work. Palomar-QUEST allows for the study of rare objects using multiple epochs of consistently calibrated, homogeneous data, obviating the common problem of generating comparable measurements from disparate datasets. A power law fit to the quasar structure function versus time yields an index of 0.432 +/- 0.024 for our best measured sample. We see the commonly reported anticorrelation between average optical variability amplitude and optical luminosity, and measure the logarithmic decrease in variability amplitude to scale as the logarithm of the luminosity times 0.205 +/- 0.002. Black hole mass is positively correlated with variability amplitude over three orders of magnitude in mass. Quasar variability amplitude is seen to decrease with Eddington ratio as a step function with a transition around Eddington ratio of 0.5. The higher variability at low Eddington ratios is due to excess power at timescales shorter than roughly 300 days. X-ray and radio measurements exist for subsets of the quasar sample. We observe an anticorrelation between optical variability amplitude and X-ray luminosity. No significant correlation is seen between average optical variability properties and radio luminosity. The timescales of quasar fluctuations are suggestive of accretion disk instabilities. The relationships seen between variability, Eddington ratio, and radio and X-ray emission are discussed in terms of a possible link between the behavior of quasars and black hole X-ray binaries.

We report here on the discovery of a red quasar, J004929.4+351025.7 at a redshift of z = 2.48, situated within a large Lya emission-line halo. The radio spectral energy distribution implies that the radio jets were triggered < 10^4 years prior to the time at which the object is observed, suggesting that the jet triggering of the active galactic nucleus is recent. The loosely biconical structure of the emission-line halo suggests that it is ionised by photons emitted by the central quasar nucleus and that the central nucleus is obscured by a dusty torus with Av ~ 3.0. The large spatial extent of the Lya halo relative to the radio emission means this could only have occurred if the radio jets emerged from an already established highly-accreting black hole. This suggests that the radio-jet triggering is delayed with respect to the onset of accretion activity on to the central supermassive black hole.

We describe briefly the Palomar-Quest (PQ) digital synoptic sky survey, including its parameters, data processing, status, and plans. Exploration of the time domain is now the central scientific and technological focus of the survey. To this end, we have developed a real-time pipeline for detection of transient sources. We describe some of the early results, and lessons learned which may be useful for other, similar projects, and time-domain astronomy in general. Finally, we discuss some issues and challenges posed by the real-time analysis and scientific exploitation of massive data streams from modern synoptic sky surveys.

A program that we call the QUEST Data Processing Software Pipeline has been written to process the large volumes of data produced by the QUEST camera on the Samuel Oschin Schmidt Telescope at the Palomar Observatory. The program carries out both aperture and PSF photometry, combines data from different repeated observations of the same portion of sky, and produces a Master Object Catalog. A rough calibration of the data is carried out. This program, as well as the calibration procedures and quality checks on the output are described.

We have designed, constructed and put into operation a very large area CCD camera that covers the field of view of the 1.2 m Samuel Oschin Schmidt Telescope at the Palomar Observatory. The camera consists of 112 CCDs arranged in a mosaic of four rows with 28 CCDs each. The CCDs are 600 x 2400 pixel Sarnoff thinned, back illuminated devices with 13 um x 13 um pixels. The camera covers an area of 4.6 deg x 3.6 deg on the sky with an active area of 9.6 square degrees. This camera has been installed at the prime focus of the telescope, commissioned, and scientific quality observations on the Palomar-QUEST Variability Sky Survey were started in September of 2003. The design considerations, construction features, and performance parameters of this camera are described in this paper.

We present four spectra of the Type Ia supernova (SN Ia) 2006D extending from -7 to +13 days with respect to B-band maximum. The spectra include the strongest signature of unburned material at photospheric velocities observed in a SN Ia to date. The earliest spectrum exhibits C II absorption features below 14,000 km/s, including a distinctive C II \lambda 6580 absorption feature. The carbon signatures dissipate as the SN approaches peak brightness. In addition to discussing implications of photospheric-velocity carbon for white dwarf explosion models, we outline some factors that may influence the frequency of its detection before and around peak brightness. Two effects are explored in this regard, including depopulation of the C II optical levels by non-LTE effects, and line-of-sight effects resulting from a clumpy distribution of unburned material with low volume-filling factor.

We report Nearby Supernova Factory observations of SN 2005gj, the second confirmed case of a "hybrid" Type Ia/IIn supernova. Our early-phase photometry of SN 2005gj shows that the interaction is much stronger than for the prototype, SN 2002ic. Our first spectrum shows a hot continuum with broad and narrow H-alpha emission. Later spectra, spanning over 4 months from outburst, show clear Type Ia features combined with broad and narrow H-gamma, H-beta, H-alpha and HeI 5876,7065 in emission. At higher resolution, P Cygni profiles are apparent. Surprisingly, we also observe an inverted P Cygni profile for [OIII] 5007. We find that the lightcurve and measured velocity of the unshocked circumstellar material imply mass loss as recently as 8 years ago. The early lightcurve is well-described by a flat radial density profile for the circumstellar material. However, our decomposition of the spectra into Type Ia and shock emission components allows for little obscuration of the supernova, suggesting an aspherical or clumpy distribution for the circumstellar material. We suggest that the emission line velocity profiles arise from electron scattering rather than the kinematics of the shock. This is supported by the inferred high densities, and the lack of evidence for evolution in the line widths. Ground- and space-based photometry, and Keck spectroscopy, of the host galaxy are used to ascertain that the host galaxy has low metallicity Z/Zsun < 0.3; (95% confidence) and that this galaxy is undergoing a significant star formation event that began roughly 200+/-70 Myr ago. We discuss the implications of these observations for progenitor models and cosmology using Type Ia supernovae.

We present results from a study to determine how star formation contributes to galaxy growth since redshift z=1.5. Using galaxies from the MUnich Near-Infrared Cluster Survey (MUNICS) and the FORS Deep Field (FDF), we investigate the specific star formation rate (SSFR, star formation rate [SFR] per unit galaxy stellar mass) as a function of galaxy stellar mass and redshift. We test the compatibility of our results with a sample drawn from a larger volume using the Sloan Digital Sky Survey. We find that the SSFR decreases as galaxy stellar mass increases, suggesting that star formation contributes more to the growth of low-mass galaxies than high-mass galaxies at all redshifts in this study. We also find a ridge in the SSFR that runs parallel to lines of constant SFR and decreases by a factor of 10 from z=1 to today, matching the results of the evolution in SFR density seen in the ``Lilly-Madau'' diagram. The ridge evolves independently of galaxy stellar mass to a particular turnover mass at the high mass end. Galaxies above the turnover mass show a sharp decrease in SSFR compared to the average at that epoch, and the turnover mass increases with redshift.

We present a study to determine how star formation contributes to galaxy growth since z=1.5 over five decades in galaxy stellar mass. We investigate the specific star formation rate (SSFR; star formation rate [SFR] per unit galaxy stellar mass) as a function of galaxy stellar mass and redshift. A sample of 175 K-band selected galaxies from the MUnich Near-Infrared Cluster Survey spectroscopic dataset provide intermediate to high mass galaxies (mostly M* > 10^10 Msun) to z=1. The FORS Deep Field provides 168 low mass galaxies (mostly M* < 10^10 Msun) to z=1.5. We use a Sloan Digital Sky Survey galaxy sample to test the compatibility of our results with data drawn from a larger volume. We find that at all redshifts, the SSFR decreases with increasing galaxy stellar mass suggesting that star formation contributes more to the growth of low mass galaxies than to the growth of high mass galaxies, and that high mass galaxies formed the bulk of their stellar content before z=1. At each epoch we find a ridge in SSFR versus stellar mass that is parallel to lines of constant SFR and evolves independently of galaxy stellar mass to a particular turnover mass. Galaxies above this turnover mass show a sharp decrease in the SFR compared to the average at each epoch and the turnover mass increases with redshift. The SFR along the SSFR ridge decreases by roughly a factor of 10, from 10 Msun/yr at z=1.5 to 1 Msun/yr at z=0. High mass galaxies could sustain the observed rates of star formation over the 10 Gyr observed, but low mass galaxies likely undergo episodic starbursts.

One of the new frontiers of astronomical research is the exploration of time variability on the sky at different wavelengths and flux levels. We have carried out a pilot project using DPOSS data to study strong variables and transients, and are now extending it to the new Palomar-QUEST synoptic sky survey. We report on our early findings and outline the methodology to be implemented in preparation for a real-time transient detection pipeline. In addition to large numbers of known types of highly variable sources (e.g., SNe, CVs, OVV QSOs, etc.), we expect to find numerous transients whose nature may be established by a rapid follow-up. Whereas we will make all detected variables publicly available through the web, we anticipate that email alerts would be issued in the real time for a subset of events deemed to be the most interesting. This real-time process entails many challenges, in an effort to maintain a high completeness while keeping the contamination low. We will utilize distributed Grid services developed by the GRIST project, and implement a variety of advanced statistical and machine learning techniques.

Exploration of the time variability on the sky over a broad range of flux levels and wavelengths is rapidly becoming a new frontier of astronomical research. We describe here briefly the Palomar-QUEST survey being carried out from the Samuel Oschin 48-inch Schmidt telescope at Palomar. The following features make the survey an attractive candidate for studying time variability: anticipated survey area of 12,000 - 15,000 sq. degrees in the drift scan mode, point source depth of 21st mag. in I under good conditions, near simultaneous observations in four filters, and at least four passes per year at each location covered. The survey will yield a large number of transients and highly variable sources in the near future and in that sense is a prototype of LSST and Pan-STARRS. We briefly outline our strategy for searching such objects and the proposed pipeline for detecting transients in real-time.

The Universe became fully reionized, and observable optically, at a time corresponding to redshift z ~ 6.5, so it is only by studying the HI and molecular absorption lines against higher-redshift, radio-loud sources that one can hope to make detailed studies of the earliest stages of galaxy formation. At present no targets for such studies are known. In these proceedings we describe a survey which is underway to find radio-loud quasars at z > 6.5, and present broad-band SEDs of our most promising candidates.

The advent of wide-area multicolour synoptic sky surveys is leading to data sets unprecedented in size, complexity and data throughput. VO technology offers a way to exploit these to the full but requires changes in design philosophy. The Palomar-QUEST survey is a major new survey being undertaken by Caltech, Yale, JPL and Indiana University to repeatedly observe 1/3 of the sky (~15000 sq. deg. between -27 < Dec <27 in seven passbands. Utilising the 48-inch Oschin Schmidt Telescope at the Palomar Observatory with the 112-CCD QUEST camera covering the full 4 x 4 sq. deg. field of view, it will generate \~1TB of data per month. In this paper, we review the design of QUEST as a VO resource, a federated data set and an exemplar of VO standards.

The Universe became fully reionized, and observable optically, at a time corresponding to redshift z~6.5, so it is only by studying the HI and molecular absorption lines against higher-redshift, radio-loud sources that one can hope to make detailed studies of the earliest stages of galaxy formation. At present no targets for such studies are known. In these proceedings we describe a survey which is underway to find radio-loud quasars at z > 6.5.