sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Simpson, C.

Normalized to: Simpson, C.

196 article(s) in total. 1245 co-authors, from 1 to 50 common article(s). Median position in authors list is 6,0.

[1]  [pdf] - 2053407
The Edge of the Galaxy
Comments: 14 pages, 11 figures. Submitted to MNRAS
Submitted: 2020-02-21
We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge'' --- a caustic manifested in a drop in density or radial velocity --- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback'' radius, and a second caustic located at ~0.6r_200m which likely corresponds to the edge of the virialized material which has completed at least two pericentric passages. The splashback radius is ill defined in Local Group type environments where the halos of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and seems a more natural definition for the boundary of the Milky Way halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase-space distribution of dwarf galaxies in the Local Group. Using the current dwarf galaxy population, we predict the edge of the Milky Way halo to be 292 +/- 61 kpc.
[2]  [pdf] - 2050526
Gravitational-wave parameter estimation with autoregressive neural network flows
Comments: 14 pages, 7 figures
Submitted: 2020-02-18
We introduce the use of autoregressive normalizing flows for rapid likelihood-free inference of binary black hole system parameters from gravitational-wave data with deep neural networks. A normalizing flow is an invertible mapping on a sample space that can be used to induce a transformation from a simple probability distribution to a more complex one: if the simple distribution can be rapidly sampled and its density evaluated, then so can the complex distribution. Our first application to gravitational waves uses an autoregressive flow, conditioned on detector strain data, to map a multivariate standard normal distribution into the posterior distribution over system parameters. We train the model on artificial strain data consisting of IMRPhenomPv2 waveforms drawn from a five-parameter $(m_1, m_2, \phi_0, t_c, d_L)$ prior and stationary Gaussian noise realizations with a fixed power spectral density. This gives performance comparable to current best deep-learning approaches to gravitational-wave parameter estimation. We then build a more powerful latent variable model by incorporating autoregressive flows within the variational autoencoder framework. This model has performance comparable to Markov chain Monte Carlo and, in particular, successfully models the multimodal $\phi_0$ posterior. Finally, we train the autoregressive latent variable model on an expanded parameter space, including also aligned spins $(\chi_{1z}, \chi_{2z})$ and binary inclination $\theta_{JN}$, and show that all parameters and degeneracies are well-recovered. In all cases, sampling is extremely fast, requiring less than two seconds to draw $10^4$ posterior samples.
[3]  [pdf] - 2038426
An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: Dust attenuation in high-redshift Lyman break Galaxies
Comments: 17 pages, 12 figures. Resubmitted to MNRAS after referee report
Submitted: 2020-01-02
We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.
[4]  [pdf] - 2005590
Magnetising the circumgalactic medium of disk galaxies
Comments: 12 pages, 12 figures, submitted to MNRAS, comments welcome
Submitted: 2019-11-25
The circumgalactic medium (CGM) is one of the frontiers of galaxy formation and intimately connected to the galaxy via accretion of gas on to the galaxy and gaseous outflows from the galaxy. Here we analyse the magnetic field in the CGM of the Milky Way-like galaxies simulated as part of the \textsc{Auriga} project that constitutes a set of high resolution cosmological magnetohydrodynamical zoom simulations. We show that at high redshift the CGM becomes magnetised via galactic outflows that transport magnetised gas from the disk into the halo. At this time the magnetisation of the CGM closely follows its metal enrichment. We then show that at low redshift an in-situ turbulent dynamo that operates on a timescale of Gigayears further amplifies the magnetic field in the CGM and sets its strength down $z=0$. The magnetic field strength reaches a typical value of $0.1\,\mu G$ at the virial radius at $z=0$ and becomes mostly uniform within the virial radius. Its Faraday rotation signal is in excellent agreement with recent observations. For most of its evolution the magnetic field in the CGM is an unordered small scale field. Only strong coherent outflows at low redshift are able to order the magnetic field in parts of the CGM that are directly displaced by these outflows.
[5]  [pdf] - 2025442
Satellites of Satellites: The Case for Carina and Fornax
Comments: 7 pages, 3 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society
Submitted: 2019-04-01, last modified: 2019-11-22
We use the Auriga cosmological simulations of Milky Way (MW)-mass galaxies and their surroundings to study the satellite populations of dwarf galaxies in $\Lambda$CDM. As expected from prior work, the number of satellites above a fixed stellar mass is a strong function of the mass of the primary dwarf. For galaxies as luminous as the Large Magellanic Cloud (LMC), and for halos as massive as expected for the LMC (determined by its rotation speed), the simulations predict about 3 satellites with stellar masses exceeding $M_*>10^5\, M_\odot$. If the LMC is on its first pericentric passage, then these satellites should be near the LMC and should have orbital angular momenta roughly coincident with that of the LMC. We use 3D positions and velocities from the 2nd data release of the Gaia mission to revisit which of the "classical" MW dwarf spheroidals could plausibly be LMC satellites. The new proper motions of the Fornax and Carina dwarf spheroidals place them on orbits closely aligned with the orbital plane of the Magellanic Clouds, hinting at a potential Magellanic association. Together with the Small Magellanic Cloud (SMC), this result raises to $3$ the number of LMC satellites with $M_*>10^5\, M_\odot$, as expected from simulations. This also fills the 12-mag luminosity gap between the SMC and the ultra-faints Hyi1, Car2, Hor1, and Car3, the few ultra-faint satellites confirmed to have orbits consistent with a Magellanic origin.
[6]  [pdf] - 2026037
Search for Astronomical Neutrinos from Blazar TXS0506+056 in Super-Kamiokande
Hagiwara, K.; Abe, K.; Bronner, C.; Hayato, Y.; Ikeda, M.; Ito, H.; Kameda, J.; Kataoka, Y.; Kato, Y.; Kishimoto, Y.; Marti, Ll.; Miura, M.; Moriyama, S.; Mochizuki, T.; Nakahata, M.; Nakajima, Y.; Nakayama, S.; Okada, T.; Okamoto, K.; Orii, A.; Pronost, G.; Sekiya, H.; Shiozawa, M.; Sonoda, Y.; Takeda, A.; Takenaka, A.; Tanaka, H.; Yano, T.; Akutsu, R.; Kajita, T.; Okumura, K.; Wang, R.; Xia, J.; Bravo-Berguno, D.; Labarga, L.; Fernandez, P.; Blaszczyk, F. d. M.; Kearns, E.; Raaf, J. L.; Stone, J. L.; Wan, L.; Wester, T.; Bian, J.; Griskevich, N. J.; Kropp, W. R.; Locke, S.; Mine, S.; Smy, M. B.; Sobel, H. W.; Takhistov, V.; Weatherly, P.; Ganezer, K. S.; Hill, J.; Kim, J. Y.; Lim, I. T.; Park, R. G.; Bodur, B.; Scholberg, K.; Walter, C. W.; Coffani, A.; Drapier, O.; Gonin, M.; Mueller, Th. A.; Paganini, P.; Ishizuka, T.; Nakamura, T.; Jang, J. S.; Learned, J. G.; Matsuno, S.; Litchfield, R. P.; Sztuc, A. A.; Uchida, Y.; Berardi, V.; Calabria, N. F.; Catanesi, M. G.; Radicioni, E.; De Rosa, G.; Collazuol, G.; Iacob, F.; Ludovici, L.; Nishimura, Y.; Cao, S.; Friend, M.; Hasegawa, T.; Ishida, T.; Kobayashi, T.; Nakadaira, T.; Nakamura, K.; Oyama, Y.; Sakashita, K.; Sekiguchi, T.; Tsukamoto, T.; Hasegawa, M.; Isobe, Y.; Miyabe, H.; Nakano, Y.; Shiozawa, T.; Sugimoto, T.; Suzuki, A. T.; Takeuchi, Y.; Ali, A.; Ashida, Y.; Hirota, S.; Jiang, M.; Kikawa, T.; Mori, M.; Nakamura, KE.; Nakaya, T.; Wendell, R. A.; Anthony, L. H. V.; McCauley, N.; Pritchard, A.; Tsui, K. M.; Fukuda, Y.; Itow, Y.; Niwa, T.; Taani, M.; Tsukada, M.; Mijakowski, P.; Frankiewicz, K.; Jung, C. K.; Vilela, C.; Wilking, M. J.; Yanagisawa, C.; Fukuda, D.; Harada, M.; Horai, T.; Ishino, H.; Ito, S.; Koshio, Y.; Sakuda, M.; Takahira, Y.; Xu, C.; Kuno, Y.; Cook, L.; Simpson, C.; Wark, D.; Di Lodovico, F.; Sedgwick, S. Molina; Richards, B.; Zsoldos, S.; Kim, S. B.; Thiesse, M.; Thompson, L.; Okazawa, H.; Choi, Y.; Nishijima, K.; Koshiba, M.; Yokoyama, M.; Goldsack, A.; Martens, K.; Quilain, B.; Suzuki, Y.; Vagins, M. R.; Kuze, M.; Tanaka, M.; Yoshida, T.; Ishitsuka, M.; Matsumoto, R.; Ohta, K.; Martin, J. F.; Nantais, C. M.; Tanaka, H. A.; Towstego, T.; Hartz, M.; Konaka, A.; de Perio, P.; Chen, S.; Jamieson, B.; Walker, J.; Minamino, A.; Okamoto, K.; Pintaudi, G.
Submitted: 2019-10-16, last modified: 2019-11-18
We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.
[7]  [pdf] - 1996808
The Milky Way total mass profile as inferred from Gaia DR2
Comments: 21 pages, 14 figures, main results: figs 3, 7, 10 and 11; submitted to MNRAS, comments are welcomed
Submitted: 2019-11-11
We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the \textit{Gaia} DR2 Galactic rotation curve and other data. Using various hydrodynamical simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, $r\lesssim 20$ kpc. We provide an analytic expression that relates the baryonic distribution to the change in the DM halo profile. For our galaxy, the contraction increases the enclosed DM halo mass by factors of roughly 1.3, 2 and 4 at radial distances of 20, 8 and 1 kpc, respectively compared to an uncontracted halo. Ignoring this contraction results in systematic biases in the inferred halo mass and concentration. We provide a best-fitting contracted NFW halo model to the MW rotation curve that matches the data very well. The best-fit has a DM halo mass, $M_{200}^{\rm DM}=0.99_{-0.20}^{+0.18}\times10^{12} M_\odot$, and concentration before baryon contraction of $8.2_{-1.5}^{+1.7}$, which lie close to the median halo mass--concentration relation predicted in $\Lambda$CDM. The inferred total mass, $M_{200}^{\rm total}=1.12_{-0.22}^{+0.20} \times 10^{12} M_\odot$, is in good agreement with recent measurements. The model gives a MW stellar mass of $4.99_{-0.50}^{+0.34}\times10^{10} M_\odot$, of which $60\%$ is contained in the thin stellar disc, with a bulge-to-total ratio of 0.2. We infer that the DM density at the Solar position is $\rho_{\odot}^{\rm DM} = 9.0_{-0.4}^{+0.5}\times10^{-3} M_\odot \rm{pc}^{-3} \equiv 0.34_{-0.02}^{+0.02}~\rm{GeV}~\rm{cm}^{-3}$. The rotation curve data can also be fitted with an uncontracted NFW halo model, but with very different DM and stellar parameters. The observations prefer the physically motivated contracted NFW halo, but the measurement uncertainties are too large to rule out the uncontracted NFW halo.
[8]  [pdf] - 1981807
An ALMA survey of the SCUBA-2 CLS UDS field: Physical properties of 707 Sub-millimetre Galaxies
Submitted: 2019-10-16
We analyse the physical properties of a large, homogeneously selected sample of ALMA-located sub-mm galaxies (SMGs) detected in the SCUBA-2 CLS 850-$\mu$m map of the UKIDSS/UDS field. This survey, AS2UDS, identified 707 SMGs across the ~1 sq.deg. field, including ~17 per cent which are undetected in the optical/near-infrared to $K$>~25.7 mag. We interpret the UV-to-radio data using a physically motivated model, MAGPHYS and determine a median photometric redshift of z=2.61+-0.08, with a 68th percentile range of z=1.8-3.4 and just ~6 per cent at z>4. The redshift distribution is well fit by a model combining evolution of the gas fraction in halos with the growth of halo mass past a threshold of ~4x10$^{12}$M$_\odot$, thus SMGs may represent the highly efficient collapse of gas-rich massive halos. Our survey provides a sample of the most massive, dusty galaxies at z>~1, with median dust and stellar masses of $M_d$=(6.8+-0.3)x10$^{8}$M$_\odot$ (thus, gas masses of ~10$^{11}$M$_\odot$) and $M_\ast=$(1.26+-0.05)x10$^{11}$M$_\odot$. These galaxies have gas fractions of $f_{gas}=$0.41+-0.02 with depletion timescales of ~150Myr. The gas mass function evolution at high masses is consistent with constraints at lower masses from blind CO-surveys, with an increase to z~2-3 and then a decline at higher redshifts. The space density and masses of SMGs suggests that almost all galaxies with $M_\ast$>~2x10$^{11}$M$_\odot$ have passed through an SMG-like phase. We find no evolution in dust temperature at a constant far-infrared luminosity across z~1.5-4. We show that SMGs appear to behave as simple homologous systems in the far-infrared, having properties consistent with a centrally illuminated starburst. Our study provides strong support for an evolutionary link between the active, gas-rich SMG population at z>1 and the formation of massive, bulge-dominated galaxies across the history of the Universe.
[9]  [pdf] - 1996732
Sensitivity of Super-Kamiokande with Gadolinium to Low Energy Anti-neutrinos from Pre-supernova Emission
Simpson, C.; Abe, K.; Bronner, C.; Hayato, Y.; Ikeda, M.; Ito, H.; Iyogi, K.; Kameda, J.; Kataoka, Y.; Kato, Y.; Kishimoto, Y.; Marti, Ll.; Miura, M.; Moriyama, S.; Mochizuki, T.; Nakahata, M.; Nakajima, Y.; Nakayama, S.; Okada, T.; Okamoto, K.; Orii, A.; Pronost, G.; Sekiya, H.; Shiozawa, M.; Sonoda, Y.; Takeda, A.; Takenaka, A.; Tanaka, H.; Yano, T.; Akutsu, R.; Kajita, T.; Okumura, K.; Wang, R.; Xia, J.; Bravo-Berguño, D.; Labarga, L.; Fernandez, P.; Blaszczyk, F. d. M.; Kachulis, C.; Kearns, E.; Raaf, J. L.; Stone, J. L.; Wan, L.; Wester, T.; Sussman, S.; Berkman, S.; Bian, J.; Griskevich, N. J.; Kropp, W. R.; Locke, S.; Mine, S.; Smy, M. B.; Sobel, H. W.; Takhistov, V.; Weatherly, P.; Ganezer, K. S.; Hill, J.; Kim, J. Y.; Lim, I. T.; Park, R. G.; Bodur, B.; Scholberg, K.; Walter, C. W.; Coffani, A.; Drapier, O.; Gonin, M.; Imber, J.; Mueller, Th. A.; Paganini, P.; Ishizuka, T.; Nakamura, T.; Jang, J. S.; Choi, K.; Learned, J. G.; Matsuno, S.; Litchfield, R. P.; Sztuc, A. A.; Uchida, Y.; Wascko, M. O.; Berardi, V.; Calabria, N. F.; Catanesi, M. G.; Intonti, R. A.; Radicioni, E.; De Rosa, G.; Collazuol, G.; Iacob, F.; Ludovici, L.; Nishimura, Y.; Cao, S.; Friend, M.; Hasegawa, T.; Ishida, T.; Kobayashi, T.; Nakadaira, T.; Nakamura, K.; Oyama, Y.; Sakashita, K.; Sekiguchi, T.; Tsukamoto, T.; Abe, KE.; Hasegawa, M.; Isobe, Y.; Miyabe, H.; Nakano, Y.; Shiozawa, T.; Sugimoto, T.; Suzuki, A. T.; Takeuchi, Y.; Ali, A.; Ashida, Y.; Hayashino, T.; Hirota, S.; Jiang, M.; Kikawa, T.; Mori, M.; Nakamura, KE.; Nakaya, T.; Wendell, R. A.; Anthony, L. H. V.; McCauley, N.; Pritchard, A.; Tsui, K. M.; Fukuda, Y.; Itow, Y.; Murrase, M.; Niwa, T.; Taani, M.; Tsukada, M.; Mijakowski, P.; Frankiewicz, K.; Jung, C. K.; Li, X.; Palomino, J. L.; Santucci, G.; Vilela, C.; Wilking, M. J.; Yanagisawa, C.; Fukuda, D.; Harada, M.; Hagiwara, K.; Horai, T.; Ishino, H.; Ito, S.; Koshio, Y.; Sakuda, M.; Takahira, Y.; Xu, C.; Kuno, Y.; Cook, L.; Wark, D.; Di Lodovico, F.; Sedgwick, S. Molina; Richards, B.; Zsoldos, S.; Kim, S. B.; Tacik, R.; Thiesse, M.; Thompson, L.; Okazawa, H.; Choi, Y.; Nishijima, K.; Koshiba, M.; Yokoyama, M.; Goldsack, A.; Martens, K.; Murdoch, M.; Quilain, B.; Suzuki, Y.; Vagins, M. R.; Kuze, M.; Okajima, Y.; Tanaka, M.; Yoshida, T.; Ishitsuka, M.; Matsumoto, R.; Ohta, K.; Martin, J. F.; Nantais, C. M.; Tanaka, H. A.; Towstego, T.; Hartz, M.; Konaka, A.; de Perio, P.; Chen, S.; Jamieson, B.; Walker, J.; Minamino, A.; Okamoto, K.; Pintaudi, G.
Comments: 20 pages, 12 figures. Accepted by The Astrophysical Journal
Submitted: 2019-08-20, last modified: 2019-09-26
Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.
[10]  [pdf] - 2025519
Simulating cosmological substructure in the solar neighbourhood
Comments: 6 pages, 4 figures, accepted for publication in MNRAS letters
Submitted: 2019-05-23, last modified: 2019-09-12
We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase space substructure and velocity correlations with the Auriga simulations and Aurigaia mock-Gaia catalogues. We show that at the solar circle the Auriga simulations commonly host phase space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These structures can persist for a few Gyrs, even after coherent streams in position space have been erased. We also explore velocity two-point correlation functions and find this diagnostic is not deterministic for particular clustering patterns in phase space. Finally, we explore these structure diagnostics with the Aurigaia catalogues and show that the current catalogues have the ability to recover some structures in phase space but careful consideration is required to separate physical structures from numerical structures arising from catalogue generation methods.
[11]  [pdf] - 1975247
Radio spectra and sizes of ALMA-identified submillimetre galaxies: evidence of age-related spectral curvature and cosmic ray diffusion?
Comments: 25 pages, 10 colour figures, 2 tables; accepted for publication in ApJ
Submitted: 2019-04-18, last modified: 2019-08-16
We analyse the multi-frequency radio spectral properties of $41$ 6GHz-detected ALMA-identified, submillimetre galaxies (SMGs), observed at 610MHz, 1.4GHz, 6GHz with GMRT and the VLA. Combining high-resolution ($\sim0.5''$) 6GHz radio and ALMA $870\,\mu$m imaging (tracing rest-frame $\sim20$GHz, and $\sim250\,\mu$m dust continuum), we study the far-infrared/radio correlation via the logarithmic flux ratio $q_{\rm IR}$, measuring $\langle q_{\rm IR}\rangle=2.20\pm 0.06$ for our sample. We show that the high-frequency radio sizes of SMGs are $\sim1.9\pm 0.4\times$ ($\sim2$-$3$kpc) larger than those of the cool dust emission, and find evidence for a subset of our sources being extended on $\sim 10$kpc scales at 1.4GHz. By combining radio flux densities measured at three frequencies, we can move beyond simple linear fits to the radio spectra of high-redshift star-forming galaxies, and search for spectral curvature, which has been observed in local starburst galaxies. At least a quarter (10/41) of our sample show evidence of a spectral break, with a median $\langle\alpha^{1.4\,{\rm GHz}}_{610\,{\rm GHz}}\rangle=-0.60\pm 0.06$, but $\langle\alpha^{6\,{\rm GHz}}_{1.4\,{\rm GHz}}\rangle=-1.06\pm 0.04$ -- a high-frequency flux deficit relative to simple extrapolations from the low-frequency data. We explore this result within this subset of sources in the context of age-related synchrotron losses, showing that a combination of weak magnetic fields ($B\sim35\,\mu$G) and young ages ($t_{\rm SB}\sim40$--$80\,$Myr) for the central starburst can reproduce the observed spectral break. Assuming these represent evolved (but ongoing) starbursts and we are observing these systems roughly half-way through their current episode of star formation, this implies starburst durations of $\lesssim100$Myr, in reasonable agreement with estimates derived via gas depletion timescales.
[12]  [pdf] - 1927963
The velocity anisotropy of the Milky Way satellite system
Comments: 16 pages, 8 figures. Key results are summarized in Figures 6 and 7. Accepted to MNRAS
Submitted: 2018-10-24, last modified: 2019-08-05
We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $\beta$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $\beta\sim -2$ at $r\sim20$ kpc to $\beta\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $\beta(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $\beta(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system.
[13]  [pdf] - 1905856
The effects of dynamical substructure on Milky Way mass estimates from the high velocity tail of the local stellar halo
Comments: Accepted to the Monthly Notices of the Royal Astronomical Society, Letters, with minor revision
Submitted: 2019-05-23, last modified: 2019-06-03
We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise and truncation of the high velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by $\sim 20 \%$, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of $1.29 ^{+0.37}_{-0.47} \times 10^{12}$ $\rm M_{\odot}$.
[14]  [pdf] - 1892184
No cores in dark matter-dominated dwarf galaxies with bursty star formation histories
Comments: 15 pages, 9 figures, 2 tables. Updated version accepted for publication in MNRAS
Submitted: 2018-10-08, last modified: 2019-05-30
Measurements of the rotation curves of dwarf galaxies are often interpreted as requiring a constant density core at the centre, at odds with the "cuspy" inner profiles predicted by $N$-body simulations of cold dark matter (CDM) haloes. It has been suggested that this conflict could be resolved by fluctuations in the inner gravitational potential caused by the periodic removal of gas following bursts of star formation. Earlier work has suggested that core formation requires a bursty and extended star formation history (SFH). Here we investigate the structure of CDM haloes of dwarf galaxies ($M_{{\rm DM}} \sim 10^9-5\times10^{10}\,{\rm M}_\odot$) formed in the APOSTLE ('A Project of Simulating the Local Environment') and AURIGA cosmological hydrodynamic simulations. Our simulations have comparable or better resolution than others that make cores ($M_{{\rm gas}} \sim 10^4\,{\rm M}_\odot$, gravitational softening $\sim 150$ pc). Yet, we do not find evidence of core formation at {\it any} mass or any correlation between the inner slope of the DM density profile and temporal variations in the SFH. APOSTLE and AURIGA dwarfs display a similar diversity in their cumulative SFHs to available data for Local Group dwarfs. Dwarfs in both simulations are DM-dominated on all resolved scales at all times, likely limiting the ability of gas outflows to alter significantly the central density profiles of their haloes. We conclude that recurrent bursts of star formation are not sufficient to cause the formation of cores, and that other conditions must also be met for baryons to be able to modify the central DM cusp.
[15]  [pdf] - 1905841
Gas engaged in noncircular motions in LITTLE THINGS dwarf irregular galaxies
Comments: In press in AJ
Submitted: 2019-05-16
We have examined gas engaged in noncircular motions in 22 of the nearby LITTLE THINGS dwarf irregular galaxies. The HI data cubes have been deconvolved into kinematic components - bulk rotation and noncircular motions, to produce maps of integrated gas, velocity field, and velocity dispersion in the different components. We found significant regions of gas engaged in noncircular motions in half of the galaxies, involving 1%-20% of the total HI mass of the galaxy. In one galaxy we found a pattern in the velocity field that is characteristic of streaming motions around the stellar bar potential and star formation at the end of bar. Two galaxies have large-scale filamentary structures found in their outer disks, and these filaments could be transient instabilities in the gas. We found no spatial correlation between noncircular motion gas and enhanced star formation. We found noncircular motion gas in only one galaxy associated with higher HI velocity dispersion.
[16]  [pdf] - 1849199
The star formation histories of dwarf galaxies in Local Group cosmological simulations
Comments: 18 pages, 10 figures, submitted to MNRAS
Submitted: 2018-12-13, last modified: 2019-03-12
We use the APOSTLE and Auriga cosmological simulations to study the star formation histories (SFHs) of field and satellite dwarf galaxies. Despite sizeable galaxy-to-galaxy scatter, the SFHs of APOSTLE and Auriga dwarfs exhibit robust average trends with galaxy stellar mass: faint field dwarfs ($10^5<M_{\rm star}/M_\odot<10^{6.5}$) have, on average, steadily declining SFHs, whereas brighter dwarfs ($10^{7.5}<M_{\rm star}/M_\odot<10^{9}$) show the opposite trend. Intermediate-mass dwarfs have roughly constant SFHs. Satellites exhibit similar average trends, but with substantially suppressed star formation in the most recent $\sim 5$ Gyr, likely as a result of gas loss due to tidal and ram-pressure stripping after entering the haloes of their primaries. These simple mass and environmental trends are in good agreement with the derived SFHs of Local Group (LG) dwarfs whose photometry reaches the oldest main sequence turnoff. SFHs of galaxies with less deep data show deviations from these trends, but this may be explained, at least in part, by the large galaxy-to-galaxy scatter, the limited sample size, and the large uncertainties of the inferred SFHs. Confirming the predicted mass and environmental trends will require deeper photometric data than currently available, especially for isolated dwarfs.
[17]  [pdf] - 1842310
The Auriga Stellar Haloes: Connecting stellar population properties with accretion and merging history
Comments: Accepted to MNRAS. 30 pages, 19 figures
Submitted: 2018-04-20, last modified: 2019-02-21
We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass--metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies.
[18]  [pdf] - 2025354
Efficacy of early stellar feedback in low gas surface density environments
Comments: 17 pages, 18 figures, Submitted to MNRAS
Submitted: 2018-12-04
We present a suite of high resolution radiation hydrodynamic simulations of a small patch ($1 \ {\rm kpc}^2$) of the inter-stellar medium (ISM) performed with Arepo-RT, with the aim to quantify the efficacy of various feedback processes like supernovae explosions (SNe), photoheating and radiation pressure in low gas surface density galaxies ($\Sigma_{\rm gas} \simeq 10 \ {\rm M}_\odot \ {\rm pc}^{-2}$). We show that radiation fields decrease the star formation rate and therefore the total stellar mass formed by a factor of $\sim 2$. This increases the gas depletion timescale and brings the simulated Kennicutt-Schmidt relation closer to the observational estimates. Radiation feedback coupled with SNe is more efficient at driving outflows with the mass and energy loading increasing by a factor of $\sim 10$. This increase is mainly driven by the additional entrainment of medium density ($10^{-2} \leq n< 1 \ {\rm cm}^{-3}$), warm ($300 \ {\rm K}\leq T<8000 \ {\rm K}$) material. Therefore including radiation fields tends to launch colder, denser and higher mass and energy loaded outflows. This is because photoheating of the high density gas around a newly formed star over-pressurises the region, causing it to expand. This reduces the ambient density in which the SNe explode by a factor of $10-100$ which in turn increases their momentum output by a factor of $\sim 1.5-2.5$. Finally, we note that in these low gas surface density environments, radiation fields primarily impact the ISM via photoheating and radiation pressure has only a minimal role in regulating star formation.
[19]  [pdf] - 1790711
Hyper-Kamiokande Design Report
Proto-Collaboration, Hyper-Kamiokande; :; Abe, K.; Abe, Ke.; Aihara, H.; Aimi, A.; Akutsu, R.; Andreopoulos, C.; Anghel, I.; Anthony, L. H. V.; Antonova, M.; Ashida, Y.; Aushev, V.; Barbi, M.; Barker, G. J.; Barr, G.; Beltrame, P.; Berardi, V.; Bergevin, M.; Berkman, S.; Berns, L.; Berry, T.; Bhadra, S.; Bravo-Berguño, D.; Blaszczyk, F. d. M.; Blondel, A.; Bolognesi, S.; Boyd, S. B.; Bravar, A.; Bronner, C.; Avanzini, M. Buizza; Cafagna, F. S.; Cole, A.; Calland, R.; Cao, S.; Cartwright, S. L.; Catanesi, M. G.; Checchia, C.; Chen-Wishart, Z.; Choi, J. H.; Choi, K.; Coleman, J.; Collazuol, G.; Cowan, G.; Cremonesi, L.; Dealtry, T.; De Rosa, G.; Densham, C.; Dewhurst, D.; Drakopoulou, E. L.; Di Lodovico, F.; Drapier, O.; Dumarchez, J.; Dunne, P.; Dziewiecki, M.; Emery, S.; Esmaili, A.; Evangelisti, A.; Fernandez-Martinez, E.; Feusels, T.; Finch, A.; Fiorentini, A.; Fiorillo, G.; Fitton, M.; Frankiewicz, K.; Friend, M.; Fujii, Y.; Fukuda, Y.; Fukuda, D.; Ganezer, K.; Giganti, C.; Gonin, M.; Grant, N.; Gumplinger, P.; Hadley, D. R.; Hartfiel, B.; Hartz, M.; Hayato, Y.; Hayrapetyan, K.; Hill, J.; Hirota, S.; Horiuchi, S.; Ichikawa, A. K.; Iijima, T.; Ikeda, M.; Imber, J.; Inoue, K.; Insler, J.; Intonti, R. A.; Ioannisian, A.; Ishida, T.; Ishino, H.; Ishitsuka, M.; Itow, Y.; Iwamoto, K.; Izmaylov, A.; Jamieson, B.; Jang, H. I.; Jang, J. S.; Jeon, S. H.; Jiang, M.; Jonsson, P.; Joo, K. K.; Kaboth, A.; Kachulis, C.; Kajita, T.; Kameda, J.; Kataoka, Y.; Katori, T.; Kayrapetyan, K.; Kearns, E.; Khabibullin, M.; Khotjantsev, A.; Kim, J. H.; Kim, J. Y.; Kim, S. B.; Kim, S. Y.; King, S.; Kishimoto, Y.; Kobayashi, T.; Koga, M.; Konaka, A.; Kormos, L. L.; Koshio, Y.; Korzenev, A.; Kowalik, K. L.; Kropp, W. R.; Kudenko, Y.; Kurjata, R.; Kutter, T.; Kuze, M.; Labarga, L.; Lagoda, J.; Lasorak, P. J. J.; Laveder, M.; Lawe, M.; Learned, J. G.; Lim, I. T.; Lindner, T.; Litchfield, R. P.; Longhin, A.; Loverre, P.; Lou, T.; Ludovici, L.; Ma, W.; Magaletti, L.; Mahn, K.; Malek, M.; Maret, L.; Mariani, C.; Martens, K.; Marti, Ll.; Martin, J. F.; Marzec, J.; Matsuno, S.; Mazzucato, E.; McCarthy, M.; McCauley, N.; McFarland, K. S.; McGrew, C.; Mefodiev, A.; Mermod, P.; Metelko, C.; Mezzetto, M.; Migenda, J.; Mijakowski, P.; Minakata, H.; Minamino, A.; Mine, S.; Mineev, O.; Mitra, A.; Miura, M.; Mochizuki, T.; Monroe, J.; Moon, D. H.; Moriyama, S.; Mueller, T.; Muheim, F.; Murase, K.; Muto, F.; Nakahata, M.; Nakajima, Y.; Nakamura, K.; Nakaya, T.; Nakayama, S.; Nantais, C.; Needham, M.; Nicholls, T.; Nishimura, Y.; Noah, E.; Nova, F.; Nowak, J.; Nunokawa, H.; Obayashi, Y.; O'Keeffe, H. M.; Okajima, Y.; Okumura, K.; Onishchuk, Yu.; O'Sullivan, E.; O'Sullivan, L.; Ovsiannikova, T.; Owen, R. A.; Oyama, Y.; Pac, M. Y.; Palladino, V.; Palomino, J. L.; Paolone, V.; Parker, W.; Parsa, S.; Payne, D.; Perkin, J. D.; Pidcott, C.; Guerra, E. Pinzon; Playfer, S.; Popov, B.; Posiadala-Zezula, M.; Poutissou, J. M.; Pritchard, A.; Prouse, N. W.; Pronost, G.; Przewlocki, P.; Quilain, B.; Radicioni, E.; Ratoff, P. N.; Retiere, F.; Riccio, C.; Richards, B.; Rondio, E.; Rose, H. J.; Rott, C.; Rountree, S. D.; Ruggeri, A. C.; Rychter, A.; Sacco, R.; Sakuda, M.; Sanchez, M. C.; Scantamburlo, E.; Scott, M.; Sedgwick, S. M.; Seiya, Y.; Sekiguchi, T.; Sekiya, H.; Seo, S. H.; Sgalaberna, D.; Shah, R.; Shaikhiev, A.; Shimizu, I.; Shiozawa, M.; Shitov, Y.; Short, S.; Simpson, C.; Sinnis, G.; Smy, M. B.; Snow, S.; Sobczyk, J.; Sobel, H. W.; Sonoda, Y.; Spina, R.; Stewart, T.; Stone, J. L.; Suda, Y.; Suwa, Y.; Suzuki, Y.; Suzuki, A. T.; Svoboda, R.; Taani, M.; Tacik, R.; Takeda, A.; Takenaka, A.; Taketa, A.; Takeuchi, Y.; Takhistov, V.; Tanaka, H. A.; Tanaka, H. K. M.; Tanaka, H.; Terri, R.; Thiesse, M.; Thompson, L. F.; Thorpe, M.; Tobayama, S.; Touramanis, C.; Towstego, T.; Tsukamoto, T.; Tsui, K. M.; Tzanov, M.; Uchida, Y.; Vagins, M. R.; Vasseur, G.; Vilela, C.; Vogelaar, R. B.; Walding, J.; Walker, J.; Ward, M.; Wark, D.; Wascko, M. O.; Weber, A.; Wendell, R.; Wilkes, R. J.; Wilking, M. J.; Wilson, J. R.; Xin, T.; Yamamoto, K.; Yanagisawa, C.; Yano, T.; Yen, S.; Yershov, N.; Yeum, D. N.; Yokoyama, M.; Yoshida, T.; Yu, I.; Yu, M.; Zalipska, J.; Zaremba, K.; Ziembicki, M.; Zito, M.; Zsoldos, S.
Comments: 325 pages
Submitted: 2018-05-09, last modified: 2018-11-28
On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.
[20]  [pdf] - 1773120
Compact star-forming galaxies preferentially quenched to become PSBs in $z<1$ clusters
Comments: 12 pages, 10 figures, accepted to MNRAS
Submitted: 2018-10-23
We analyse the structure of galaxies with high specific star formation rate (SSFR) in cluster and field environments in the redshift range $0.5<z<1.0$. Recent studies have shown that these galaxies are strongly depleted in dense environments due to rapid environmental quenching, giving rise to post-starburst galaxies (PSBs). We use effective radii and S\'ersic indices as tracers of galaxy structure, determined using imaging from the UKIDSS Ultra Deep Survey (UDS). We find that the high-SSFR galaxies that survive into the cluster environment have, on average, larger effective radii than those in the field. We suggest that this trend is likely to be driven by the most compact star-forming galaxies being preferentially quenched in dense environments. We also show that the PSBs in clusters have stellar masses and effective radii that are similar to the missing compact star-forming population, suggesting that these PSBs are the result of size-dependent quenching. We propose that both strong stellar feedback and the stripping of the extended halo act together to preferentially and rapidly quench the compact and low-mass star-forming systems in clusters to produce PSBs. We test this scenario using the stacked spectra of 124 high-SSFR galaxies, showing that more compact galaxies are more likely to host outflows. We conclude that a combination of environmental and secular processes is the most likely explanation for the appearance of PSBs in galaxy clusters.
[21]  [pdf] - 1767506
Faraday rotation maps of disk galaxies
Comments: 10 pages, 10 figures, accepted by MNRAS
Submitted: 2018-07-05, last modified: 2018-09-19
Faraday rotation is one of the most widely used observables to infer the strength and configuration of the magnetic field in the ionised gas of the Milky Way and nearby spiral galaxies. Here we compute synthetic Faraday rotation maps at $z=0$ for a set of disk galaxies from the Auriga high-resolution cosmological simulations, for different observer positions within and outside the galaxy. We find that the strength of the Faraday rotation of our simulated galaxies for a hypothetic observer at the solar circle is broadly consistent with the Faraday rotation seen for the Milky Way. The same holds for an observer outside the galaxy and the observed signal of the nearby spiral galaxy M51. However, we also find that the structure and angular power spectra of the synthetic all-sky Faraday rotation maps vary strongly with azimuthal position along the solar circle. We argue that this variation is a result of the structure of the magnetic field of the galaxy that is dominated by an azimuthal magnetic field ordered scales of several kpc, but has radial and vertical magnetic field components that are only ordered on scales of 1-2 kpc. Because the magnetic field strength decreases exponentially with height above the disk, the Faraday rotation for an observer at the solar circle is dominated by the local environment. This represents a severe obstacle for attempts to reconstruct the global magnetic field of the Milky Way from Faraday rotation maps alone without including additional observables.
[22]  [pdf] - 1747895
Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations
Comments: Accepted for publication in MNRAS (1st of September, 2018) after improvements to several figures and some discussion. The mock data are available for access and download at and documentation and supplementary material is available at
Submitted: 2018-04-23, last modified: 2018-09-05
We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $\Lambda$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $\Lambda$CDM cosmological model.
[23]  [pdf] - 1717237
The structure of post-starburst galaxies at $0.5 < z < 2$: evidence for two distinct quenching routes at different epochs
Comments: 21 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2018-07-03
We present an analysis of the structure of post-starburst (PSB) galaxies in the redshift range $0.5 < z < 2$, using a photometrically-selected sample identified in the Ultra Deep Survey (UDS) field. We examine the structure of $\sim80$ of these transient galaxies using radial light $\mu(r)$ profiles obtained from CANDELS $\textit{Hubble Space Telescope}$ near-infrared/optical imaging, and compare to a large sample of $\sim2000$ passive and star-forming galaxies. For each population, we determine their typical structural properties (effective radius $r_{\rm e}$, S\'ersic index $n$) and find significant differences in PSB structure at different epochs. At high redshift ($z > 1$), PSBs are typically massive ($M_* > 10^{10}\rm\,M_{\odot}$), very compact and exhibit high S\'ersic indices, with structures that differ significantly from their star-forming progenitors but are similar to massive passive galaxies. In contrast, at lower redshift ($0.5 < z < 1$), PSBs are generally of low mass ($M_* < 10^{10}\rm\,M_{\odot}$) and exhibit compact but less concentrated profiles (i.e. lower S\'ersic indices), with structures similar to low-mass passive discs. Furthermore, for both epochs we find remarkably consistent PSB structure across the optical/near-infrared wavebands (which largely trace different stellar populations), suggesting that any preceding starburst and/or quenching in PSBs was not strongly centralized. Taken together, these results imply that PSBs at $z > 1$ have been recently quenched during a major disruptive event (e.g. merger or protogalactic collapse) which formed a compact remnant, while at $z < 1$ an alternative less disruptive process is primarily responsible. Our results suggest that high-$z$ PSBs are an intrinsically different population to those at lower redshifts, and indicate different quenching routes are active at different epochs.
[24]  [pdf] - 1729063
A machine-learning method for identifying multi-wavelength counterparts of submillimeter galaxies: training and testing using AS2UDS and ALESS
Comments: 25 pages, 10 figures, three tables, accepted for publication in ApJ
Submitted: 2018-06-18
We describe the application of the supervised machine-learning algorithms to identify the likely multi-wavelength counterparts to submillimeter sources detected in panoramic, single-dish submillimeter surveys. As a training set, we employ a sample of 695 ($S_{\rm 870\mu m}$ >1 mJy) submillimeter galaxies (SMGs) with precise identifications from the ALMA follow-up of the SCUBA-2 Cosmology Legacy Survey's UKIDSS-UDS field (AS2UDS). We show that radio emission, near-/mid-infrared colors, photometric redshift, and absolute $H$-band magnitude are effective predictors that can distinguish SMGs from submillimeter-faint field galaxies. Our combined radio+machine-learning method is able to successfully recover $\sim$85 percent of ALMA-identified SMGs which are detected in at least three bands from the ultraviolet to radio. We confirm the robustness of our method by dividing our training set into independent subsets and using these for training and testing respectively, as well as applying our method to an independent sample of $\sim$100 ALMA-identified SMGs from the ALMA/LABOCA ECDF-South Survey (ALESS). To further test our methodology, we stack the 870$\mu$m ALMA maps at the positions of those $K$-band galaxies that are classified as SMG counterparts by the machine-learning but do not have a $>$4.3$ \sigma$ ALMA detection. The median peak flux density of these galaxies is $S_{\rm 870\mu m}=(0.61\pm0.03)$ mJy, demonstrating that our method can recover faint and/or diffuse SMGs even when they are below the detection threshold of our ALMA observations. In future, we will apply this method to samples drawn from panoramic single-dish submillimeter surveys which currently lack interferometric follow-up observations, to address science questions which can only be tackled with large, statistical samples of SMGs.
[25]  [pdf] - 1712958
An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: number counts of submillimeter galaxies
Comments: 11 pages, 4 figures, fixed typo (sorry Rob)
Submitted: 2018-05-14, last modified: 2018-05-16
We report the first results of AS2UDS: an 870 $\mu$m continuum survey with the Atacama Large Millimeter/Submillimeter Array (ALMA) of a total area of $\sim$ 50 arcmin$^2$ comprising a complete sample of 716 submillimeter sources drawn from the SCUBA-2 Cosmology Legacy Survey (S2CLS) map of the UKIDSS/UDS field. The S2CLS parent sample covers a 0.96 degree$^2$ field at $\sigma_{850}=0.90\pm0.05$ mJy beam$^{-1}$. Our deep, high-resolution ALMA observations with $\sigma_{\rm 870}\sim$ 0.25 mJy and a 0.15"--0.30" FWHM synthesized beam, provide precise locations for 695 submillimetre galaxies (SMGs) responsible for the submillimeter emission corresponding to 606 sources in the low resolution, single-dish map. We measure the number counts of SMGs brighter than $S_{\rm 870}\geq$ 4 mJy, free from the effects of blending and show that the normalisation of the counts falls by 28 $\pm$ 2% in comparison to the SCUBA-2 parent sample, but that the shape remains unchanged. We determine that 44$^{+16}_{-14}$% of the brighter single-dish sources with $S_{850}\geq$ 9 mJy consist of a blend of two or more ALMA-detectable SMGs brighter than $S_{870}\sim$ 1 mJy (corresponding to a galaxy with a total-infrared luminosity of $L_{IR}\geq$ 10$^{12}$ L$_\odot$), in comparison to 28 $\pm$ 2% for the single-dish sources at $S_{\rm 850}\geq$ 5 mJy. Using the 46 single-dish submillimeter sources that contain two or more ALMA-detected SMGs with photometric redshifts, we show that there is a significant statistical excess of pairs of SMGs with similar redshifts ($<$1% probability of occurring by chance), suggesting that at least 30% of these blends arise from physically associated pairs of SMGs.
[26]  [pdf] - 1686716
On the relevance of chaos for halo stars in the solar neighbourhood II
Comments: 17 pages, 10 figures, 3 tables. Accepted for publication in MNRAS
Submitted: 2018-01-11, last modified: 2018-05-15
In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of chaotic mixing using multicomponent Galactic potential models and solar neighbourhood-like volumes extracted from fully cosmological hydrodynamic simulations, thus naturally accounting for the gravitational potential associated with baryonic components, such as the bulge and disc. Despite the strong change in the global Galactic potentials with respect to those obtained in dark matter only simulations, our results confirm that a large fraction of halo particles evolving on chaotic orbits exhibit their chaotic behaviour after periods of time significantly larger than a Hubble time. In addition, significant diffusion in phase space is not observed on those particles that do exhibit chaotic behaviour within a Hubble time.
[27]  [pdf] - 1721117
An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: Identifying candidate z~4.5 [CII] emitters
Comments: 16 pages, 10 figures, accepted for publication in ApJ
Submitted: 2018-05-14
We report the results of a search for serendipitous [CII] 157.74$\mu$m emitters at $z\sim4.4$-$4.7$ using the Atacama Large Millimeter/submillimeter Array (ALMA). The search exploits the AS2UDS continuum survey, which covers ~50 arcmin$^2$ of the sky towards 695 luminous ($S_{870}\gtrsim1$mJy) submillimeter galaxies (SMGs), selected from the SCUBA-2 Cosmology Legacy Survey (S2CLS) 0.96deg$^2$ Ultra Deep Survey (UDS) field. We detect ten candidate line emitters, with an expected false detection rate of ten percent. All of these line emitters correspond to 870$\mu$m continuum-detected sources in AS2UDS. The emission lines in two emitters appear to be high-J CO, but the remainder have multi-wavelength properties consistent with [CII] from $z\simeq4.5$ galaxies. Using our sample, we place a lower limit of $>5\times10^{-6}$Mpc$^{-3}$ on the space density of luminous ($L_{\rm IR} \simeq 10^{13}$Lsun) SMGs at $z=4.40$-$4.66$, suggesting $\ge7$percent of SMGs with $S_{870\mu{\rm m}}\gtrsim1$mJy lie at $4<z<5$. From stacking the high-resolution ($\sim0.15"$ full-width half maximum) ALMA $870\mu$m imaging, we show that the [CII] line emission is more extended than the continuum dust emission, with an average effective radius for the [CII] of $r_{\rm e} = 1.7^{+0.1}_{-0.2}$kpc compared to $r_{\rm e} = 1.0\pm0.1$kpc for the continuum (rest-frame $160\mu$m). By fitting the far-infrared photometry for these galaxies from $100$-$870\mu$m, we show that SMGs at $z\sim4.5$ have a median dust temperature of $T_{\rm d}=55\pm4$K. This is systematically warmer than $870\mu$m-selected SMGs at $z\simeq2$, which typically have temperatures around $35$K. These $z\simeq4.5$ SMGs display a steeper trend in the luminosity-temperature plane than $z\le2$ SMGs. We discuss the implications of this result in terms of the selection biases of high redshift starbursts in far-infrared/submillimeter surveys.
[28]  [pdf] - 1670610
Search for Neutrinos in Super-Kamiokande associated with the GW170817 neutron-star merger
Abe, K.; Bronner, C.; Hayato, Y.; Ikeda, M.; Iyogi, K.; Kameda, J.; Kato, Y.; Kishimoto, Y.; Marti, Ll.; Miura, M.; Moriyama, S.; Nakahata, M.; Nakajima, Y.; Nakano, Y.; Nakayama, S.; Orii, A.; Pronost, G.; Sekiya, H.; Shiozawa, M.; Sonoda, Y.; Takeda, A.; Takenaka, A.; Tanaka, H.; Tasaka, S.; Yano, T.; Akutsu, R.; Kajita, T.; Nishimura, Y.; Okumura, K.; Tsui, K. M.; Labarga, L.; Fernandez, P.; Blaszczyk, F. d. M.; Kachulis, C.; Kearns, E.; Raaf, J. L.; Stone, J. L.; Sulak, L. R.; Berkman, S.; Tobayama, S.; Bian, J.; Elnimr, M.; Kropp, W. R.; Locke, S.; Mine, S.; Weatherly, P.; Smy, M. B.; Sobel, H. W.; Takhistov, V.; Ganezer, K. S.; Hill, J.; Kim, J. Y.; Lim, I. T.; Park, R. G.; Li, Z.; O'Sullivan, E.; Scholberg, K.; Walter, C. W.; Gonin, M.; Imber, J.; Mueller, Th. A.; Ishizuka, T.; Nakamura, T.; Jang, J. S.; Choi, K.; Learned, J. G.; Matsuno, S.; Amey, J.; Litchfield, R. P.; Ma, W. Y.; Uchida, Y.; Wascko, M. O.; Catanesi, M. G.; Intonti, R. A.; Radicioni, E.; De Rosa, G.; Ali, A.; Collazuol, G.; Ludovici, L.; Cao, S.; Friend, M.; Hasegawa, T.; Ishida, T.; Ishii, T.; Kobayashi, T.; Nakadaira, T.; Nakamura, K.; Oyama, Y.; Sakashita, K.; Sekiguchi, T.; Tsukamoto, T.; Abe, KE.; Hasegawa, M.; Suzuki, A. T.; Takeuchi, Y.; Hayashino, T.; Hirota, S.; Jiang, M.; Mori, M.; Nakamura, KE.; Nakaya, T.; Wendell, R. A.; Anthony, L. H. V.; McCauley, N.; Pritchard, A.; Fukuda, Y.; Itow, Y.; Murase, M.; Muto, F.; Mijakowski, P.; Frankiewicz, K.; Jung, C. K.; Li, X.; Palomino, J. L.; Santucci, G.; Viela, C.; Wilking, M. J.; Yanagisawa, C.; Fukuda, D.; Ishino, H.; Ito, S.; Kibayashi, A.; Koshio, Y.; Nagata, H.; Sakuda, M.; Xu, C.; Kuno, Y.; Wark, D.; Di Lodovico, F.; Richards, B.; Sedgwick, S. Molina; Tacik, R.; Kim, S. B.; Cole, A.; Thompson, L.; Okazawa, H.; Choi, Y.; Ito, K.; Nishijima, K.; Koshiba, M.; Suda, Y.; Yokoyama, M.; Calland, R. G.; Hartz, M.; Martens, K.; Murdoch, M.; Quilain, B.; Simpson, C.; Suzuki, Y.; Vagins, M. R.; Hamabe, D.; Kuze, M.; Okajima, Y.; Yoshida, T.; Ishitsuka, M.; Martin, J. F.; Nantais, C. M.; Tanaka, H. A.; Towstego, T.; Konaka, A.; Chen, S.; Wan, L.; Minamino, A.
Comments: 8 pages, 4 figures
Submitted: 2018-02-12, last modified: 2018-03-29
We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.
[29]  [pdf] - 1659507
Quenching and ram pressure stripping of simulated Milky Way satellite galaxies
Comments: 22 pages, 17 figures, accepted for publication in MNRAS
Submitted: 2017-05-08, last modified: 2018-03-20
We present predictions for the quenching of star formation in satellite galaxies of the Local Group from a suite of 30 cosmological zoom simulations of Milky Way-like host galaxies. The Auriga simulations resolve satellites down to the luminosity of the classical dwarf spheroidal galaxies of the Milky Way. We find strong mass-dependent and distance-dependent quenching signals, where dwarf systems beyond 600 kpc are only strongly quenched below a stellar mass of $10^7$ M$_\odot$. Ram pressure stripping appears to be the dominant quenching mechanism and 50% of quenched systems cease star formation within 1 Gyr of first infall. We demonstrate that systems within a host galaxy's $R_{200}$ radius are comprised of two populations: (i) a first infall population that has entered the host halo within the past few Gyrs and (ii) a population of returning `backsplash' systems that have had a much more extended interaction with the host. Backsplash galaxies that do not return to the host galaxy by redshift zero exhibit quenching properties similar to galaxies within $R_{200}$ and are distinct from other external systems. The simulated quenching trend with stellar mass has some tension with observations, but our simulations are able reproduce the range of quenching times measured from resolved stellar populations of Local Group dwarf galaxies.
[30]  [pdf] - 1942018
The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test
Comments: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at, The first paper of the AGORA Initiative is at, More information on AGORA is at
Submitted: 2016-10-10, last modified: 2018-03-08
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.
[31]  [pdf] - 1641376
The multiplicity and anisotropy of galactic satellite accretion
Comments: 15 pages, 17 figures. Minor changes. MNRAS in press
Submitted: 2017-12-14, last modified: 2018-02-27
We study the incidence of group and filamentary dwarf galaxy accretion into Milky Way (MW) mass haloes using two types of hydrodynamical simulations: EAGLE, which resolves a large cosmological volume, and the AURIGA suite, which are very high resolution zoom-in simulations of individual MW-sized haloes. The present-day 11 most massive satellites are predominantly (75%) accreted in single events, 14% in pairs and 6% in triplets, with higher group multiplicities being unlikely. Group accretion becomes more common for fainter satellites, with 60% of the top 50 satellites accreted singly, 12% in pairs, and 28% in richer groups. A group similar in stellar mass to the Large Magellanic Cloud (LMC) would bring on average 15 members with stellar mass larger than $10^4{~\rm M_\odot}$. Half of the top 11 satellites are accreted along the two richest filaments. The accretion of dwarf galaxies is highly anisotropic, taking place preferentially perpendicular to the halo minor axis, and, within this plane, preferentially along the halo major axis. The satellite entry points tend to be aligned with the present-day central galaxy disc and satellite plane, but to a lesser extent than with the halo shape. Dwarfs accreted in groups or along the richest filament have entry points that show an even larger degree of alignment with the host halo than the full satellite population. We also find that having most satellites accreted as a single group or along a single filament is unlikely to explain the MW disc of satellites.
[32]  [pdf] - 1630404
The enhancement of rapidly quenched galaxies in distant clusters at 0.5<z<1.0
Comments: 17 pages, 10 figures, Accepted to MNRAS
Submitted: 2018-02-05
We investigate the relationship between environment and galaxy evolution in the redshift range $0.5 < z < 1.0$. Galaxy overdensities are selected using a Friends-of-Friends algorithm, applied to deep photometric data in the Ultra-Deep Survey (UDS) field. A study of the resulting stellar mass functions reveals clear differences between cluster and field environments, with a strong excess of low-mass rapidly quenched galaxies in cluster environments compared to the field. Cluster environments also show a corresponding deficit of young, low-mass star-forming galaxies, which show a sharp radial decline towards cluster centres. By comparing mass functions and radial distributions, we conclude that young star-forming galaxies are rapidly quenched as they enter overdense environments, becoming post-starburst galaxies before joining the red sequence. Our results also point to the existence of two environmental quenching pathways operating in galaxy clusters, operating on different timescales. Fast quenching acts on galaxies with high specific star-formation rates, operating on timescales shorter than the cluster dynamical time ($ < 1$ Gyr). In contrast, slow quenching affects galaxies with moderate specific star-formation rates, regardless of their stellar mass, and acts on longer timescales ($\gtrsim 1$ Gyr). Of the cluster galaxies in the stellar mass range $9.0 < \log(M_{*}/M_{\odot}) < 10.5$ quenched during this epoch, we find that 73% were transformed through fast quenching, while the remaining 27% followed the slow quenching route.
[33]  [pdf] - 1629287
Massive post-starburst galaxies at z > 1 are compact proto-spheroids
Comments: 12 pages, 10 figures, published in MNRAS. References corrected and updated to match journal article
Submitted: 2017-07-31, last modified: 2018-02-05
We investigate the relationship between the quenching of star formation and the structural transformation of massive galaxies, using a large sample of photometrically-selected post-starburst galaxies in the UKIDSS UDS field. We find that post-starburst galaxies at high-redshift ($z>1$) show high S\'ersic indices, significantly higher than those of active star-forming galaxies, but with a distribution that is indistinguishable from the old quiescent population. We conclude that the morphological transformation occurs before (or during) the quenching of star formation. Recently quenched galaxies are also the most compact; we find evidence that massive post-starburst galaxies (M$_{\ast}> 10^{10.5} ~$M$_{\odot}$) at high redshift ($z>1$) are on average smaller than comparable quiescent galaxies at the same epoch. Our findings are consistent with a scenario in which massive passive galaxies are formed from three distinct phases: (1) gas-rich dissipative collapse to very high densities, forming the proto-spheroid; (2) rapid quenching of star formation, to create the "red nugget" with post-starburst features; (3) a gradual growth in size as the population ages, perhaps as a result of minor mergers.
[34]  [pdf] - 1613084
Search for an excess of events in the Super-Kamiokande detector in the directions of the astrophysical neutrinos reported by the IceCube Collaboration
Collaboration, The Super-Kamiokande; :; Abe, K.; Bronner, C.; Pronost, G.; Hayato, Y.; Ikeda, M.; Iyogi, K.; Kameda, J.; Kato, Y.; Kishimoto, Y.; Marti, Ll.; Miura, M.; Moriyama, S.; Nakahata, M.; Nakano, Y.; Nakayama, S.; Okajima, Y.; Orii, A.; Sekiya, H.; Shiozawa, M.; Sonoda, Y.; Takeda, A.; Takenaka, A.; Tanaka, H.; Tasaka, S.; Tomura, T.; Akutsu, R.; Kajita, T.; Kaneyuki, K.; Nishimura, Y.; Okumura, K.; Tsui, K. M.; Labarga, L.; Fernandez, P.; Blaszczyk, F. d. M.; Gustafson, J.; Kachulis, C.; Kearns, E.; Raaf, J. L.; Stone, J. L.; Sulak, L. R.; Berkman, S.; Tobayama, S.; Goldhaber, M.; Elnimr, M.; Kropp, W. R.; Mine, S.; Locke, S.; Weatherly, P.; Smy, M. B.; Sobel, H. W.; Takhistov, V.; Ganezer, K. S.; Hill, J.; Kim, J. Y.; Lim, I. T.; Park, R. G.; Himmel, A.; Li, Z.; O'Sullivan, E.; Scholberg, K.; Walter, C. W.; Ishizuka, T.; Nakamura, T.; Jang, J. S.; Choi, K.; Learned, J. G.; Matsuno, S.; Smith, S. N.; Amey, J.; Litchfield, R. P.; Ma, W. Y.; Uchida, Y.; Wascko, M. O.; Cao, S.; Friend, M.; Hasegawa, T.; Ishida, T.; Ishii, T.; Kobayashi, T.; Nakadaira, T.; Nakamura, K.; Oyama, Y.; Sakashita, K.; Sekiguchi, T.; Tsukamoto, T.; Abe, KE.; Hasegawa, M.; Nakano, Y.; Suzuki, A. T.; Takeuchi, Y.; Yano, T.; Cao, S. V.; Hayashino, T.; Hiraki, T.; Hirota, S.; Huang, K.; Jiang, M.; Minamino, A.; Nakamura, KE.; Nakaya, T.; Quilain, B.; Patel, N. D.; Wendell, R. A.; Anthony, L. H. V.; McCauley, N.; Pritchard, A.; Fukuda, Y.; Itow, Y.; Murase, M.; Muto, F.; Mijakowski, P.; Frankiewicz, K.; Jung, C. K.; Li, X.; Palomino, J. L.; Santucci, G.; Vilela, C.; Wilking, M. J.; Yanagisawa, C.; Ito, S.; Fukuda, D.; Ishino, H.; Kibayashi, A.; Koshio, Y.; Nagata, H.; Sakuda, M.; Xu, C.; Kuno, Y.; Wark, D.; Di Lodovico, F.; Richards, B.; Tacik, R.; Kim, S. B.; Cole, A.; Thompson, L.; Okazawa, H.; Choi, Y.; Ito, K.; Nishijima, K.; Koshiba, M.; Totsuka, Y.; Suda, Y.; Yokoyama, M.; Calland, R. G.; Hartz, M.; Martens, K.; Simpson, C.; Suzuki, Y.; Vagins, M. R.; Hamabe, D.; Kuze, M.; Yoshida, T.; Ishitsuka, M.; Martin, J. F.; Nantais, C. M.; Tanaka, H. A.; Konaka, A.; Chen, S.; Wan, L.; Zhang, Y.; Minamino, A.; Wilkes, R. J.
Comments: 9 pages, 4 figures, Matches version published in Astrophysical Journal
Submitted: 2017-07-26, last modified: 2018-01-05
We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.
[35]  [pdf] - 1705180
A direct calibration of the IRX-{\beta} relation in Lyman-break Galaxies at z=3-5
Comments: 12 pages, 10 figures, submitted to MNRAS
Submitted: 2018-01-02
We use a sample of 4178 Lyman break galaxies (LBGs) at z = 3, 4 and 5 in the UKIRT Infrared Deep Sky Survey (UKIDSS) Ultra Deep Survey (UDS) field to investigate the relationship between the observed slope of the stellar continuum emission in the ultraviolet, {\beta}, and the thermal dust emission, as quantified via the so-called 'infrared excess' (IRX = LIR/LUV). Through a stacking analysis we directly measure the 850-{\mu}m flux density of LBGs in our deep (0.9mJy) James Clerk Maxwell Telescope (JCMT) SCUBA-2 850-{\mu}m map, as well as deep public Herschel/SPIRE 250-, 350- and 500-{\mu}m imaging. We establish functional forms for the IRX-{\beta} relation to z ~ 5, confirming that there is no significant redshift evolution of the relation and that the resulting average IRX-{\beta} curve is consistent with a Calzetti-like attenuation law. We compare our results with recent work in the literature, finding that discrepancies in the slope of the IRX-{\beta} relation are driven by biases in the methodology used to determine the ultraviolet slopes. Consistent results are found when IRX-{\beta} is evaluated by stacking in bins of stellar mass, M, and we argue that the near-linear IRX-M relationship is a better proxy for correcting observed UV luminosities to total star formation rates, provided an accurate handle on M can be had, and also gives clues as to the physical driver of the role of dust-obscured star formation in high-redshift galaxies.
[36]  [pdf] - 1626455
The dependence of cosmic ray driven galactic winds on halo mass
Comments: 16 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2017-12-13
Galactic winds regulate star formation in disk galaxies and help to enrich the circum-galactic medium. They are therefore crucial for galaxy formation, but their driving mechanism is still poorly understood. Recent studies have demonstrated that cosmic rays (CRs) can drive outflows if active CR transport is taken into account. Using hydrodynamical simulations of isolated galaxies with virial masses between $10^{10}$ and $10^{13}\mathrm{~M_\odot}$, we study how the properties of CR-driven winds depend on halo mass. CRs are treated in a two-fluid approximation and their transport is modelled through isotropic or anisotropic diffusion. We find that CRs are only able to drive mass-loaded winds beyond the virial radius in haloes with masses below $10^{12}\mathrm{~M_\odot}$. For our lowest examined halo mass, the wind is roughly spherical and has velocities of $\sim20\mathrm{~km\;s^{-1}}$. With increasing halo mass, the wind becomes biconical and can reach ten times higher velocities. The mass loading factor drops rapidly with virial mass, a dependence that approximately follows a power-law with a slope between $-1$ and $-2$. This scaling is slightly steeper than observational inferences, and also steeper than commonly used prescriptions for wind feedback in cosmological simulations. The slope is quite robust to variations of the CR injection efficiency or the CR diffusion coefficient. In contrast to the mass loading, the energy loading shows no significant dependence on halo mass. While these scalings are close to successful heuristic models of wind feedback, the CR-driven winds in our present models are not yet powerful enough to fully account for the required feedback strength.
[37]  [pdf] - 1608434
Origin of chemically distinct discs in the Auriga cosmological simulations
Comments: 12 pages, MNRAS accepted
Submitted: 2017-08-25, last modified: 2017-11-20
The stellar disk of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a clearly bimodal distribution in the $[\rm \alpha/Fe]$ -- $[\rm Fe/H]$ plane. We find two main pathways to creating this dichotomy which operate in different regions of the galaxies: a) an early ($z>1$) and intense high-$\rm[\alpha/Fe]$ star formation phase in the inner region ($R\lesssim 5$ kpc) induced by gas-rich mergers, followed by more quiescent low-$\rm[\alpha/Fe]$ star formation; and b) an early phase of high-$\rm[\alpha/Fe]$ star formation in the outer disc followed by a shrinking of the gas disc owing to a temporarily lowered gas accretion rate, after which disc growth resumes. In process b), a double-peaked star formation history around the time and radius of disc shrinking accentuates the dichotomy. If the early star formation phase is prolonged (rather than short and intense), chemical evolution proceeds as per process a) in the inner region, but the dichotomy is less clear. In the outer region, the dichotomy is only evident if the first intense phase of star formation covers a large enough radial range before disc shrinking occurs; otherwise, the outer disc consists of only low-$\rm[\alpha/Fe]$ sequence stars. We discuss the implication that both processes occurred in the Milky Way.
[38]  [pdf] - 1588452
Simulating Gamma-ray Emission in Star-forming Galaxies
Comments: Accepted by ApJL, 8 pages, 4 figures
Submitted: 2017-09-15
Star forming galaxies emit GeV- and TeV-gamma rays that are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. To understand the emission, we have used the moving mesh code Arepo to perform magneto-hydrodynamical galaxy formation simulations with self-consistent CR physics. Our galaxy models exhibit a first burst of star formation that injects CRs at supernovae. Once CRs have sufficiently accumulated in our Milky-Way like galaxy, their buoyancy force overcomes the magnetic tension of the toroidal disk field. As field lines open up, they enable anisotropically diffusing CRs to escape into the halo and to accelerate a bubble-like, CR-dominated outflow. However, these bubbles are invisible in our simulated gamma-ray maps of hadronic pion-decay and secondary inverse-Compton emission because of low gas density in the outflows. By adopting a phenomenological relation between star formation rate (SFR) and far-infrared emission and assuming that gamma rays mainly originate from decaying pions, our simulated galaxies can reproduce the observed tight relation between far-infrared and gamma-ray emission, independent of whether we account for anisotropic CR diffusion. This demonstrates that uncertainties in modeling active CR transport processes only play a minor role in predicting gamma-ray emission from galaxies. We find that in starbursts, most of the CR energy is "calorimetrically" lost to hadronic interactions. In contrast, the gamma-ray emission deviates from this calorimetric property at low SFRs due to adiabatic losses, which cannot be identified in traditional one-zone models.
[39]  [pdf] - 1588022
The MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Survey
Comments: Proceedings of Science, "MeerKAT Science: On the Pathway to the SKA", Stellenbosch, 25-27 May 2016
Submitted: 2017-09-06
The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.
[40]  [pdf] - 1585449
A complete distribution of redshifts for sub-millimetre galaxies in the SCUBA-2 Cosmology Legacy Survey UDS field
Comments: 10 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2017-07-03, last modified: 2017-07-05
Sub-milllimetre galaxies (SMGs) are some of the most luminous star-forming galaxies in the Universe, however their properties remain hard to determine due to the difficulty of identifying their optical\slash near-infrared counterparts. One of the key steps to determining the nature of SMGs is measuring a redshift distribution representative of the whole population. We do this by applying statistical techniques to a sample of 761 850$\mu$m sources from the SCUBA-2 Cosmology Legacy Survey observations of the UKIDSS Ultra-Deep Survey (UDS) Field. We detect excess galaxies around $> 98.4$ per cent of the 850$\mu$m positions in the deep UDS catalogue, giving us the first 850$\mu$m selected sample to have virtually complete optical\slash near-infrared redshift information. Under the reasonable assumption that the redshifts of the excess galaxies are representative of the SMGs themselves, we derive a median SMG redshift of $z = 2.05 \pm 0.03$, with 68 per cent of SMGs residing between $1.07 < z < 3.06$. We find an average of $1.52\pm 0.09$ excess $K$-band galaxies within 12 arc sec of an 850$\mu$m position, with an average stellar mass of $2.2\pm 0.1 \times 10^{10}$ M$_\odot$. While the vast majority of excess galaxies are star-forming, $8.0 \pm 2.1$ per cent have passive rest-frame colours, and are therefore unlikely to be detected at sub-millimetre wavelengths even in deep interferometry. We show that brighter SMGs lie at higher redshifts, and use our SMG redshift distribution -- along with the assumption of a universal far-infrared SED -- to estimate that SMGs contribute around 30 per cent of the cosmic star formation rate density between $0.5 < z < 5.0$.
[41]  [pdf] - 1584952
Extragalactic radio surveys in the pre-Square Kilometre Array era
Comments: Invited review for Royal Society Open Science, 22 pages
Submitted: 2017-06-20
The era of the Square Kilometre Array is almost upon us, and pathfinder telescopes are already in operation. This brief review summarizes our current knowledge of extragalactic radio sources, accumulated through six decades of continuum surveys at the low-frequency end of the electromagnetic spectrum and the extensive complementary observations at other wavelengths necessary to gain this understanding. The relationships between radio survey data and surveys at other wavelengths are discussed. Some of the outstanding questions are identified and prospects over the next few years are outlined.
[42]  [pdf] - 1584054
Enhancement of AGN in a protocluster at z=1.6
Comments: 9 pages, 8 figures, Accepted to MNRAS (24th May 2017)
Submitted: 2017-05-30
We investigate the prevalence of AGN in the high-redshift protocluster $\rm{Cl}\,0218.3$-$0510$ at $z=1.62$. Using imaging from the Chandra X-ray Telescope, we find a large overdensity of AGN in the protocluster; a factor of $23\pm9$ times the field density of AGN. Only half of this AGN overdensity is due to the overdensity of massive galaxies in the protocluster (a factor of $11\pm2$), as we find that $17^{+6}_{-5}\%$ of massive galaxies ($M_* > 10^{10}\,\rm{M}_{\odot}$) in the protocluster host an X-ray luminous AGN, compared to $8\pm1\%$ in the field. This corresponds to an enhancement of AGN activity in massive protocluster galaxies by a factor of $2.1\pm0.7$ at $1.6\sigma$ significance. We also find that the AGN overdensity is centrally concentrated, located within 3 arcmin and most pronounced within 1 arcmin of the centre of the protocluster. Our results confirm that there is a reversal in the local anti-correlation between galaxy density and AGN activity, so there is an enhancement of AGN in high-redshift protoclusters. We compare the properties of AGN in the protocluster to the field and find no significant differences in the distributions of their stellar mass, X-ray luminosity, or hardness ratio. We therefore suggest that triggering mechanisms are similar in both environments, and that the mechanisms simply occur more frequently in denser environments.
[43]  [pdf] - 1583851
A deep search for metals near redshift 7: the line-of-sight towards ULAS J1120+0641
Comments: 18 pages, 23 figures, accepted for publication in MNRAS
Submitted: 2017-05-24
We present a search for metal absorption line systems at the highest redshifts to date using a deep (30h) VLT/X-Shooter spectrum of the z = 7.084 quasi-stellar object (QSO) ULAS J1120+0641. We detect seven intervening systems at z > 5.5, with the highest-redshift system being a C IV absorber at z = 6.51. We find tentative evidence that the mass density of C IV remains flat or declines with redshift at z < 6, while the number density of C II systems remains relatively flat over 5 < z < 7. These trends are broadly consistent with models of chemical enrichment by star formation-driven winds that include a softening of the ultraviolet background towards higher redshifts. We find a larger number of weak ( W_rest < 0.3A ) Mg II systems over 5.9 < z < 7.0 than predicted by a power-law fit to the number density of stronger systems. This is consistent with trends in the number density of weak Mg II systems at z = 2.5, and suggests that the mechanisms that create these absorbers are already in place at z = 7. Finally, we investigate the associated narrow Si IV, C IV, and N V absorbers located near the QSO redshift, and find that at least one component shows evidence of partial covering of the continuum source.
[44]  [pdf] - 1732518
Magnetic field formation in the Milky Way-like disk galaxies of the Auriga project
Comments: 16 pages, 16 figures, accepted by MNRAS, section 4.1 (turbulent dynamo) extended substantially, now includes magnetic and kinetic power spectra
Submitted: 2017-01-24, last modified: 2017-05-02
The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{\mu G}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.
[45]  [pdf] - 1582893
Lessons from the Auriga discs: The hunt for the Milky Way's ex-situ disc is not yet over
Comments: 13 pages, 11 figures, submitted to MNRAS
Submitted: 2017-04-26
We characterize the contribution from accreted material to the galactic discs of the Auriga Project, a set of high resolution magnetohydrodynamic cosmological simulations of late-type galaxies performed with the moving-mesh code AREPO. Our goal is to explore whether a significant accreted (or ex-situ) stellar component in the Milky Way disc could be hidden within the near-circular orbit population, which is strongly dominated by stars born in-situ. One third of our models shows a significant ex-situ disc but this fraction would be larger if constraints on orbital circularity were relaxed. Most of the ex-situ material ($\gtrsim 50\%$) comes from single massive satellites ($> 6 \times 10^{10}~M_{\odot}$). These satellites are accreted with a wide range of infall times and inclination angles (up to $85^{\circ}$). Ex-situ discs are thicker, older and more metal-poor than their in-situ counterparts. They show a flat median age profile, which differs from the negative gradient observed in the in-situ component. As a result, the likelihood of identifying an ex-situ disc in samples of old stars on near-circular orbits increases towards the outskirts of the disc. We show three examples that, in addition to ex-situ discs, have a strongly rotating dark matter component. Interestingly, two of these ex-situ stellar discs show an orbital circularity distribution that is consistent with that of the in-situ disc. Thus, they would not be detected in typical kinematic studies.
[46]  [pdf] - 1554062
The Nature of Turbulence in the LITTLE THINGS Dwarf Irregular Galaxies
Comments: AJ, in press
Submitted: 2017-03-01
We present probability density functions and higher order (skewness and kurtosis) analyses of the galaxy-wide and spatially-resolved HI column density distributions in the LITTLE THINGS sample of dwarf irregular galaxies. This analysis follows that of Burkhart et al. (2010) for the Small Magellanic Cloud. About 60% of our sample have galaxy-wide values of kurtosis that are similar to that found for the Small Magellanic Cloud, with a range up to much higher values, and kurtosis increases with integrated star formation rate. Kurtosis and skewness were calculated for radial annuli and for a grid of 32 pixel X 32 pixel kernels across each galaxy. For most galaxies, kurtosis correlates with skewness. For about half of the galaxies, there is a trend of increasing kurtosis with radius. The range of kurtosis and skewness values is modeled by small variations in the Mach number close to the sonic limit and by conversion of HI to molecules at high column density. The maximum HI column densities decrease with increasing radius in a way that suggests molecules are forming in the weak field limit, where H_2 formation balances photodissociation in optically thin gas at the edges of clouds.
[47]  [pdf] - 1550495
The compact, ~1 kpc host galaxy of a quasar at z=7.1
Comments: 10 pages, 4 figures. Accepted for publication in ApJ
Submitted: 2017-02-13
We present ALMA observations of the [CII] fine structure line and the underlying far-infrared (FIR) dust continuum emission in J1120+0641, the most distant quasar currently known (z=7.1). We also present observations targeting the CO(2-1), CO(7-6) and [CI] 369 micron lines in the same source obtained at the VLA and PdBI. We find a [CII] line flux of F_[CII]=1.11+/-0.10 Jy km/s and a continuum flux density of S_227GHz=0.53+/-0.04 mJy/beam, consistent with previous unresolved measurements. No other source is detected in continuum or [CII] emission in the field covered by ALMA (~25"). At the resolution of our ALMA observations (0.23", or 1.2 kpc, a factor ~70 smaller beam area compared to previous measurements), we find that the majority of the emission is very compact: a high fraction (~80%) of the total line and continuum flux is associated with a region 1-1.5 kpc in diameter. The remaining ~20% of the emission is distributed over a larger area with radius <4 kpc. The [CII] emission does not exhibit ordered motion on kpc-scales: applying the virial theorem yields an upper limit on the dynamical mass of the host galaxy of (4.3+/-0.9)x10^10 M_sun, only ~20x higher than the central black hole. The other targeted lines (CO(2-1), CO(7-6) and [CI]) are not detected, but the limits of the line ratios with respect to the [CII] emission imply that the heating in the quasar host is dominated by star formation, and not by the accreting black hole. The star-formation rate implied by the FIR continuum is 105-340 M_sun/yr, with a resulting star-formation rate surface density of ~100-350 M_sun/yr/kpc^2, well below the value for Eddington-accretion-limited star formation.
[48]  [pdf] - 1535738
Observations of the Lyman series forest towards the redshift 7.1 quasar ULAS J1120+0641
Comments: Accepted for publication in A&A. 11 pages, 10 figures
Submitted: 2017-02-13
We present a 30h integration Very Large Telescope X-shooter spectrum of the Lyman series forest towards the $z = 7.084$ quasar ULAS J1120+0641. The only detected transmission at ${\rm S/N}>5$ is confined to seven narrow spikes in the Ly$\alpha$ forest, over the redshift range $5.858<z<6.122$, just longward of the wavelength of the onset of the Ly$\beta$ forest. There is also a possible detection of one further unresolved spike in the Ly$\beta$ forest at $z=6.854$, with ${\rm S/N}=4.5$. We also present revised Hubble Space Telescope F814W photometry of the source. The summed flux from the transmission spikes is in agreement with the F814W photometry, so all the transmission in the Lyman series forest may have been detected. There is a Gunn-Peterson (GP) trough in the Ly$\alpha$ forest from $z=6.122$ all the way to the quasar near zone at $z=7.04$. The trough, of comoving length $240\,h^{-1}$Mpc, is over twice as long as the next longest known GP trough. We combine the spectroscopic and photometric results to constrain the evolution of the Ly$\alpha$ effective optical depth with redshift, $\tau_{\rm GP}^{\rm eff}$ extending a similar analysis by Simpson et al. We find $\tau_{\rm GP}^{\rm eff} \propto (1+z)^{\xi}$ where $\xi = 11.2^{+0.4}_{-0.6}$, for $z > 5.5$. The data nevertheless provide only a weak limit on the volume-weighted hydrogen intergalactic (IGM) neutral fraction at $z\sim 6.5$, $x_{{\rm HI}} > 10^{-4}$, similar to limits at redshift $z\sim6$ from less distant quasars. The new observations cannot extend measurements of the neutral fraction of the IGM to higher values because absorption in the Ly$\alpha$ forest is already saturated near $z\sim6$. For higher neutral fractions, other methods such as measuring the red damping wing of the IGM will be required.
[49]  [pdf] - 1542938
The HI Chronicles of LITTLE THINGS BCDS III: Gas Clouds in and around Mrk 178, VII Zw 403, AND NGC 3738
Comments: 39 pages; accepted for publication in AJ
Submitted: 2017-01-27
In most blue compact dwarf (BCD) galaxies, it remains unclear what triggers their bursts of star formation. We study the HI of three relatively isolated BCDs, Mrk 178, VII Zw 403, and NGC 3738, in detail to look for signatures of star formation triggers, such as gas cloud consumption, dwarf-dwarf mergers, and interactions with companions. High angular and velocity resolution atomic hydrogen (H I) data from the Very Large Array (VLA) dwarf galaxy HI survey, Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey (LITTLE THINGS), allows us to study the detailed kinematics and morphologies of the BCDs in HI. We also present high sensitivity HI maps from the NRAO Green Bank Telescope (GBT) of each BCD to search their surrounding regions for extended tenuous emission or companions. The GBT data do not show any distinct galaxies obviously interacting with the BCDs. The VLA data indicate several possible star formation triggers in these BCDs. Mrk 178 likely has a gas cloud impacting the southeast end of its disk or it is experiencing ram pressure stripping. VII Zw 403 has a large gas cloud in its foreground or background that shows evidence of accreting onto the disk. NGC 3738 has several possible explanations for its stellar morphology and H I morphology and kinematics: an advanced merger, strong stellar feedback, or ram pressure stripping. Although apparently isolated, the HI data of all three BCDs indicate that they may be interacting with their environments, which could be triggering their bursts of star formation.
[50]  [pdf] - 1552386
Evolution of cosmic star formation in the SCUBA-2 Cosmology Legacy Survey
Comments: MNRAS Accepted 2017 January 5. 29 pages, 18 figures
Submitted: 2016-07-14, last modified: 2017-01-05
We present a new exploration of the cosmic star-formation history and dust obscuration in massive galaxies at redshifts $0.5< z<6$. We utilize the deepest 450 and 850$\mu$m imaging from SCUBA-2 CLS, covering 230arcmin$^2$ in the AEGIS, COSMOS and UDS fields, together with 100-250$\mu$m imaging from Herschel. We demonstrate the capability of the T-PHOT deconfusion code to reach below the confusion limit, using multi-wavelength prior catalogues from CANDELS/3D-HST. By combining IR and UV data, we measure the relationship between total star-formation rate (SFR) and stellar mass up to $z\sim5$, indicating that UV-derived dust corrections underestimate the SFR in massive galaxies. We investigate the relationship between obscuration and the UV slope (the IRX-$\beta$ relation) in our sample, which is similar to that of low-redshift starburst galaxies, although it deviates at high stellar masses. Our data provide new measurements of the total SFR density (SFRD) in $M_\ast>10^{10}M_\odot$ galaxies at $0.5<z<6$. This is dominated by obscured star formation by a factor of $>10$. One third of this is accounted for by 450$\mu$m-detected sources, while one fifth is attributed to UV-luminous sources (brighter than $L^\ast_{UV}$), although even these are largely obscured. By extrapolating our results to include all stellar masses, we estimate a total SFRD that is in good agreement with previous results from IR and UV data at $z\lesssim3$, and from UV-only data at $z\sim5$. The cosmic star-formation history undergoes a transition at $z\sim3-4$, as predominantly unobscured growth in the early Universe is overtaken by obscured star formation, driven by the build-up of the most massive galaxies during the peak of cosmic assembly.
[51]  [pdf] - 1530633
Simulating cosmic ray physics on a moving mesh
Comments: Accepted by MNRAS, 31 pages, 10 figures
Submitted: 2016-04-25, last modified: 2016-11-24
We discuss new methods to integrate the cosmic ray (CR) evolution equations coupled to magneto-hydrodynamics (MHD) on an unstructured moving mesh, as realised in the massively parallel AREPO code for cosmological simulations. We account for diffusive shock acceleration of CRs at resolved shocks and at supernova remnants in the interstellar medium (ISM), and follow the advective CR transport within the magnetised plasma, as well as anisotropic diffusive transport of CRs along the local magnetic field. CR losses are included in terms of Coulomb and hadronic interactions with the thermal plasma. We demonstrate the accuracy of our formalism for CR acceleration at shocks through simulations of plane-parallel shock tubes that are compared to newly derived exact solutions of the Riemann shock tube problem with CR acceleration. We find that the increased compressibility of the post-shock plasma due to the produced CRs decreases the shock speed. However, CR acceleration at spherically expanding blast waves does not significantly break the self-similarity of the Sedov-Taylor solution; the resulting modifications can be approximated by a suitably adjusted, but constant adiabatic index. In first applications of the new CR formalism to simulations of isolated galaxies and cosmic structure formation, we find that CRs add an important pressure component to the ISM that increases the vertical scale height of disk galaxies, and thus reduces the star formation rate. Strong external structure formation shocks inject CRs into the gas, but the relative pressure of this component decreases towards halo centres as adiabatic compression favours the thermal over the CR pressure.
[52]  [pdf] - 1511464
Comments: 10 pages, 4 figures, proceedings of the SF2A conference, Lyon, 2016
Submitted: 2016-11-08
In these proceedings we highlight the primary scientific goals and design of the WEAVE-LOFAR survey, which will use the new WEAVE spectrograph on the 4.2m William Herschel Telescope to provide the primary source of spectroscopic information for the LOFAR Surveys Key Science Project. Beginning in 2018, WEAVE-LOFAR will generate more than 10$^6$ R=5000 365-960 nm spectra of low-frequency selected radio sources, across three tiers designed to efficiently sample the redshift-luminosity plane, and produce a data set of enormous legacy value. The radio frequency selection, combined with the high multiplex and throughput of the WEAVE spectrograph, make obtaining redshifts in this way very efficient, and we expect that the redshift success rate will approach 100 per cent at $z < 1$. This unprecedented spectroscopic sample - which will be complemented by an integral field component - will be transformational in key areas, including studying the star formation history of the Universe, the role of accretion and AGN-driven feedback, properties of the epoch of reionisation, cosmology, cluster haloes and relics, as well as the nature of radio galaxies and protoclusters. Each topic will be addressed in unprecedented detail, and with the most reliable source classifications and redshift information in existence.
[53]  [pdf] - 1488455
The impact of protocluster environments at z = 1.6
Comments: Accepted for publication in MNRAS.
Submitted: 2016-09-27
We investigate the effects of dense environments on galaxy evolution by examining how the properties of galaxies in the z = 1.6 protocluster Cl 0218.3-0510 depend on their location. We determine galaxy properties using spectral energy distribution fitting to 14-band photometry, including data at three wavelengths that tightly bracket the Balmer and 4000A breaks of the protocluster galaxies. We find that two-thirds of the protocluster galaxies, which lie between several compact groups, are indistinguishable from field galaxies. The other third, which reside within the groups, differ significantly from the intergroup galaxies in both colour and specific star formation rate. We find that the fraction of red galaxies within the massive protocluster groups is twice that of the intergroup region. These excess red galaxies are due to enhanced fractions of both passive galaxies (1.7 times that of the intergroup region) and dusty star-forming galaxies (3 times that of the intergroup region). We infer that some protocluster galaxies are processed in the groups before the cluster collapses. These processes act to suppress star formation and change the mode of star formation from unobscured to obscured.
[54]  [pdf] - 1513517
The SCUBA-2 Cosmology Legacy Survey: The clustering of submillimetre galaxies in the UKIDSS UDS field
Comments: Accepted for publication in MNRAS, 14 pages, 5 figures
Submitted: 2016-03-31, last modified: 2016-09-21
Submillimetre galaxies (SMGs) are among the most luminous dusty galaxies in the Universe, but their true nature remains unclear; are SMGs the progenitors of the massive elliptical galaxies we see in the local Universe, or are they just a short-lived phase among more typical star-forming galaxies? To explore this problem further, we investigate the clustering of SMGs identified in the SCUBA-2 Cosmology Legacy Survey. We use a catalogue of submillimetre ($850\mu$m) source identifications derived using a combination of radio counterparts and colour/IR selection to analyse a sample of 914 SMGs in the UKIDSS Ultra Deep Survey (UDS), making this the largest high redshift sample of these galaxies to date. Using angular cross-correlation techniques, we estimate the halo masses for this large sample of SMGs and compare them with passive and star-forming galaxies selected in the same field. We find that SMGs, on average, occupy high-mass dark matter halos (M$_{\text{halo}} >10^{13}$M$_{\odot}$) at redshifts $z > 2.5$, consistent with being the progenitors of massive quiescent galaxies in present-day galaxy clusters. We also find evidence of downsizing, in which SMG activity shifts to lower mass halos at lower redshifts. In terms of their clustering and halo masses, SMGs appear to be consistent with other star-forming galaxies at a given redshift.
[55]  [pdf] - 1510356
Faint Submillimeter Galaxies identified through their optical/near-infrared colours I: spatial clustering and halo masses
Comments: ApJ in press, 13 pages, 11 figures
Submitted: 2016-09-01
The properties of submillimeter galaxies (SMGs) that are fainter than the confusion limit of blank-field single-dish surveys ($S_{850} \lesssim$ 2 mJy) are poorly constrained. Using a newly developed color selection technique, Optical-Infrared Triple Color (OIRTC), that has been shown to successfully {select} such faint SMGs, we identify a sample of 2938 OIRTC-selected galaxies, dubbed Triple Color Galaxies (TCGs), in the UKIDSS-UDS field. We show that these galaxies have a median 850 $\mu$m flux of S$_{850} = 0.96\pm0.04$ mJy (equivalent to a star-formation rate SFR $\sim60-100$ M$_\odot$ yr$^{-1}$ based on SED fitting), representing the first large sample of faint SMGs that bridges the gap between bright SMGs and normal star-forming galaxies in S$_{850}$ and $L_{\rm IR}$. We assess the basic properties of TCGs and their relationship with other galaxy populations at $z\sim2$. We measure the two-point autocorrelation function for this population and derive a typical halo mass of log$_{10}$(M$_{\rm halo}$) $=12.9^{+0.2}_{-0.3}$, $12.7^{+0.1}_{-0.2}$, and $12.9^{+0.2}_{-0.3}$ $h^{-1}$M$_\odot$ at $z=1-2$, $2-3$, and $3-5$, respectively. Together with the bright SMGs (S$_{850} \gtrsim 2$ mJy) and a comparison sample of less far-infrared luminous star-forming galaxies, we find a lack of dependence between spatial clustering and S$_{850}$ (or SFR), suggesting that the difference between these populations may lie in their local galactic environment. Lastly, on the scale of $\sim8-17$ kpc at $1<z<5$ we find a tentative enhancement of the clustering of TCGs over the comparison star-forming galaxies, suggesting that some faint SMGs are physically associated pairs, perhaps reflecting a merging origin in their triggering.
[56]  [pdf] - 1460914
The role of cosmic ray pressure in accelerating galactic outflows
Comments: 8 pages, 4 figures, accepted for publication in ApJL; movie of simulated gas densities can be found here:
Submitted: 2016-06-07, last modified: 2016-08-16
We study the formation of galactic outflows from supernova explosions (SNe) with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SNe placement and energy feedback, including cosmic rays (CR), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.
[57]  [pdf] - 1457540
The evolution of post-starburst galaxies from z=2 to z= 0.5
Comments: 14 pages, 7 figures, accepted by MNRAS
Submitted: 2016-08-01
We present the evolution in the number density and stellar mass functions of photometrically selected post-starburst galaxies in the UKIDSS Deep Survey (UDS), with redshifts of 0.5<z<2 and stellar masses logM>10. We find that this transitionary species of galaxy is rare at all redshifts, contributing ~5% of the total population at z~2, to <1% by z~0.5. By comparing the mass functions of quiescent galaxies to post-starburst galaxies at three cosmic epochs, we show that rapid quenching of star formation can account for 100% of quiescent galaxy formation, if the post-starburst spectral features are visible for ~250Myr. The flattening of the low mass end of the quiescent galaxy stellar mass function seen at z~1 can be entirely explained by the addition of rapidly quenched galaxies. Only if a significant fraction of post-starburst galaxies have features that are visible for longer than 250Myr, or they acquire new gas and return to the star-forming sequence, can there be significant growth of the red sequence from a slower quenching route. The shape of the mass function of these transitory post-starburst galaxies resembles that of quiescent galaxies at z~2, with a preferred stellar mass of logM~10.6, but evolves steadily to resemble that of star-forming galaxies at z<1. This leads us to propose a dual origin for post-starburst galaxies: (1) at z>2 they are exclusively massive galaxies that have formed the bulk of their stars during a rapid assembly period, followed by complete quenching of further star formation, (2) at z<1 they are caused by the rapid quenching of gas-rich star-forming galaxies, independent of stellar mass, possibly due to environment and/or gas-rich major mergers.
[58]  [pdf] - 1470635
Semi-implicit anisotropic cosmic ray transport on an unstructured moving mesh
Comments: 16 pages, 16 figures, accepted by MNRAS
Submitted: 2016-04-28, last modified: 2016-07-16
In the interstellar medium of galaxies and the intracluster gas of galaxy clusters, the charged particles making up cosmic rays are moving almost exclusively along (but not across) magnetic field lines. The resulting anisotropic transport of cosmic rays in the form of diffusion or streaming not only affects the gas dynamics but also rearranges the magnetic fields themselves. The coupled dynamics of magnetic fields and cosmic rays can thus impact the formation and evolution of galaxies and the thermal evolution of galaxy clusters in critical ways. Numerically studying these effects requires solvers for anisotropic diffusion that are accurate, efficient, and robust, requirements that have proven difficult to satisfy in practice. Here, we present an anisotropic diffusion solver on an unstructured moving mesh that is conservative, does not violate the entropy condition, allows for semi-implicit time integration with individual timesteps, and only requires solving a single linear system of equations per timestep. We apply our new scheme to a large number of test problems and show that it works as well or better than previous implementations. Finally, we demonstrate for a numerically demanding simulation of the formation of an isolated disk galaxy that our local time-stepping scheme reproduces the results obtained with global time-stepping at a fraction of the computational cost.
[59]  [pdf] - 1439065
Galactic winds driven by isotropic and anisotropic cosmic ray diffusion in disk galaxies
Comments: 7 pages, 5 figures, accepted by ApJL
Submitted: 2016-05-02, last modified: 2016-07-16
The physics of cosmic rays (CR) is a promising candidate for explaining the driving of galactic winds and outflows. Recent galaxy formation simulations have demonstrated the need for active CR transport either in the form of diffusion or streaming to successfully launch winds in galaxies. However, due to computational limitations, most previous simulations have modeled CR transport isotropically. Here, we discuss high resolution simulations of isolated disk galaxies in a $10^{11}\rm{M_\odot}$ halo with the moving mesh code {\sc Arepo} that include injection of CRs from supernovae, advective transport, CR cooling, and CR transport through isotropic or anisotropic diffusion. We show that either mode of diffusion leads to the formation of strong bipolar outflows. However, they develop significantly later in the simulation with anisotropic diffusion compared to the simulation with isotropic diffusion. Moreover, we find that isotropic diffusion allows most of the CRs to quickly diffuse out of the disk, while in the simulation with anisotropic diffusion, most CRs remain in the disk once the magnetic field becomes dominated by its azimuthal component, which occurs after $\sim 300\,{\rm Myrs}$. This has important consequences for the gas dynamics in the disk. In particular, we show that isotropic diffusion strongly suppresses the amplification of the magnetic field in the disk compared to anisotropic or no diffusion models. We therefore conclude that reliable simulations which include CR transport inevitably need to account for anisotropic diffusion.
[60]  [pdf] - 1531028
The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts
Comments: Submitted to MNRAS, comments welcome. Catalogue and maps at
Submitted: 2016-07-13
We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.
[61]  [pdf] - 1444092
GMC Evolution in a Barred Spiral Galaxy with Star Formation and Thermal Feedback
Comments: 19 pages, 16 figures, accepted for publication in MNRAS
Submitted: 2016-06-15
We explore the impact of star formation and thermal stellar feedback on the giant molecular cloud (GMC) population forming in a M83-type barred spiral galaxy. We compare three high-resolution simulations (1.5 pc cell size) with different star formation/feedback models: one with no star formation, one with star formation but no feedback, and one with star formation and thermal energy injection. We analyze the resulting population of clouds, finding that we can identify the same population of massive, virialized clouds and transient, low-surface density clouds found in our previous work (that did not include star formation or feedback). Star formation and feedback can affect the mix of clouds we identify. In particular, star formation alone simply converts dense cloud gas into stars with only a small change to the cloud populations, principally resulting in a slight decrease in the transient population. Feedback, however, has a stronger impact: while it is not generally sufficient to entirely destroy the clouds, it does eject gas out of them, increasing the gas density in the inter-cloud region. This decreases the number of massive clouds, but substantially increases the transient cloud population. We also find that feedback tends to drive a net radial inflow of massive clouds, leading to an increase in the star formation rate in the bar region. We examine a number of possible reasons for this and conclude that it is possible that the drag force from the enhanced intercloud density could be responsible.
[62]  [pdf] - 1426311
The identification of post-starburst galaxies at z~1 using multiwavelength photometry: a spectroscopic verification
Comments: Accepted to MNRAS (29th March 2016); 5 pages, 3 figures
Submitted: 2016-03-29
Despite decades of study, we still do not fully understand why some massive galaxies abruptly switch off their star formation in the early Universe, and what causes their rapid transition to the red sequence. Post-starburst galaxies provide a rare opportunity to study this transition phase, but few have currently been spectroscopically identified at high redshift ($z>1$). In this paper we present the spectroscopic verification of a new photometric technique to identify post-starbursts in high-redshift surveys. The method classifies the broad-band optical-near--infrared spectral energy distributions (SEDs) of galaxies using three spectral shape parameters (super-colours), derived from a principal component analysis of model SEDs. When applied to the multiwavelength photometric data in the UKIDSS Ultra Deep Survey (UDS), this technique identified over 900 candidate post-starbursts at redshifts $0.5<z<2.0$. In this study we present deep optical spectroscopy for a subset of these galaxies, in order to confirm their post-starburst nature. Where a spectroscopic assessment was possible, we find the majority (19/24 galaxies; ~80 per cent) exhibit the strong Balmer absorption (H $\delta$ equivalent width $W_{\lambda}$ >5 Ang.) and Balmer break, characteristic of post-starburst galaxies. We conclude that photometric methods can be used to select large samples of recently-quenched galaxies in the distant Universe.
[63]  [pdf] - 1379076
The structure and evolution of a forming galaxy cluster at z = 1.62
Comments: Accepted for publication in MNRAS
Submitted: 2016-03-11
We present a comprehensive picture of the Cl0218.3-0510 protocluster at $z=1.623$ across 10 co-moving Mpc. Using filters that tightly bracket the Balmer and 4000 Angstrom breaks of the protocluster galaxies we obtain precise photometric redshifts resulting in a protocluster galaxy sample that is 89+/-5% complete and has a contamination of only 12+/-5%. Both star forming and quiescent protocluster galaxies are located allowing us to map the structure of the forming cluster for the first time. The protocluster contains 6 galaxy groups, the largest of which is the nascent cluster. Only a small minority of the protocluster galaxies are in the nascent cluster (11%) or in the other galaxy groups (22%), as most protocluster galaxies reside between the groups. Unobscured star forming galaxies predominantly reside between the protocluster's groups, whereas red galaxies make up a large fraction of the groups' galactic content, so observing the protocluster through only one of these types of galaxies results in a biased view of the protocluster's structure. The structure of the protocluster reveals how much mass is available for the future growth of the cluster and we use the Millennium Simulation, scaled to a Planck cosmology, to predict that Cl0218.3-0510 will evolve into a 2.7x 10^14 Msun cluster by the present day.
[64]  [pdf] - 1385445
The SCUBA-2 Cosmology Legacy Survey: Multi-wavelengths counterparts to 10$^3$ submillimeter galaxies in the UKIDSS-UDS field
Comments: ApJ in press, 23 pages, 17 figures
Submitted: 2016-01-11
We present multiwavelength identifications for the counterparts of 1088 submillimeter sources detected at 850$\mu$m in the SCUBA-2 Cosmology Legacy Survey study of the UKIDSS-UDS field. By utilising an ALMA pilot study on a subset of our bright SCUBA-2 sample as a training set, along with the deep optical-near-infrared data available in this field, we develop a novel technique, Optical-IR Triple Color (OIRTC), using $z-K$, $K-[3.6]$, $[3.6]-[4.5]$ colors to select the candidate submillimeter galaxy (SMG) counterparts. By combining radio identification and the OIRTC technique, we find counterpart candidates for 80% of the Class = 1 $\geq4\,\sigma$ SCUBA-2 sample, defined as those that are covered by both radio and OIR imaging and the base sample for our scientific analyses. Based on the ALMA training set, we expect the accuracy of these identifications to be $82\pm20$%, with a completeness of $69\pm16$%, essentially as accurate as the traditional $p$-value technique but with higher completeness. We find that the fraction of SCUBA-2 sources having candidate counterparts is lower for fainter 850$\mu$m sources, and we argue that for follow-up observations sensitive to SMGs with $S_{850}\gtrsim 1$ mJy across the whole ALMA beam, the fraction with multiple counterparts is likely to be $>40$% for SCUBA-2 sources at $S_{850} \gtrsim 4$ mJy. We find that the photometric redshift distribution for the SMGs is well fit by a lognormal distribution, with a median redshift of $z=2.3\pm0.1$. After accounting for the sources without any radio and/or OIRTC counterpart, we estimate the median redshift to be $z=2.6\pm0.1$ for SMGs with $S_{850} >1$ mJy. We also use this new large sample to study the clustering of SMGs and the the far-infrared properties of the unidentified submillimeter sources by stacking their Herschel SPIRE far-infrared emission.
[65]  [pdf] - 1296331
Evolution of Star Formation in the UKIDSS Ultra Deep Survey Field - II. Star Formation as a Function of Stellar Mass Between z=1.46 and z=0.63
Comments: Accepted for publication in MNRAS, 12 pages, 7 figures, 2 tables
Submitted: 2015-09-23
We present new results on the evolution of the cosmic star formation rate as a function of stellar mass in the SXDS-UDS field. We make use of narrow-band selected emission line galaxies in four redshift slices between z = 1.46 and z = 0.63, and compute stellar masses by fitting a series of templates to recreate each galaxy's star formation history. We determine mass-binned luminosity functions in each redshift slice, and derive the star formation rate density (rhoSFR) as a function of mass using the [OIII] or [OII] emission lines. We calculate dust extinction and metallicity as a function of stellar mass, and investigate the effect of these corrections on the shape of the overall rhoSFR(M). We find that both these corrections are crucial for determining the shape of the rhoSFR(M), and its evolution with redshift. The fully corrected rhoSFR(M) is a relatively flat distribution, with the normalisation moving towards lower values of rhoSFR with increasing cosmic time/decreasing redshift, and requiring star formation to be truncated across all masses studied here. The peak of rhoSFR(M) is found in the 10^10.5<Msun<10^11.0 mass bin at z = 1.46. In the lower redshift slices the location of the peak is less certain, however low mass galaxies in the range 10^7.0<Msun<10^8.0 play an important part in the overall rhoSFR(M) out to at least z ~ 1.2.
[66]  [pdf] - 1273235
The impact from survey depth and resolution on the morphological classification of galaxies
Comments: accepted for publication in MNRAS; 25 pages, 19 figures, 5 tables
Submitted: 2015-07-21
We consistently analyse for the first time the impact of survey depth and spatial resolution on the most used morphological parameters for classifying galaxies through non-parametric methods: Abraham and Conselice-Bershady concentration indices, Gini, M20 moment of light, asymmetry, and smoothness. Three different non-local datasets are used, ALHAMBRA and SXDS (examples of deep ground-based surveys), and COSMOS (deep space-based survey). We used a sample of 3000 local, visually classified galaxies, measuring their morphological parameters at their real redshifts (z ~ 0). Then we simulated them to match the redshift and magnitude distributions of galaxies in the non-local surveys. The comparisons of the two sets allow to put constraints on the use of each parameter for morphological classification and evaluate the effectiveness of the commonly used morphological diagnostic diagrams. All analysed parameters suffer from biases related to spatial resolution and depth, the impact of the former being much stronger. When including asymmetry and smoothness in classification diagrams, the noise effects must be taken into account carefully, especially for ground-based surveys. M20 is significantly affected, changing both the shape and range of its distribution at all brightness levels.We suggest that diagnostic diagrams based on 2 - 3 parameters should be avoided when classifying galaxies in ground-based surveys, independently of their brightness; for COSMOS they should be avoided for galaxies fainter than F814 = 23.0. These results can be applied directly to surveys similar to ALHAMBRA, SXDS and COSMOS, and also can serve as an upper/lower limit for shallower/deeper ones.
[67]  [pdf] - 1263860
Kinetic energy from supernova feedback in high-resolution galaxy simulations
Comments: 23 pages, 16 figures, accepted for publication in ApJ
Submitted: 2014-10-14, last modified: 2015-06-26
We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells' mass and velocity in a distributed fashion. The method is explored in the context of supernova feedback in high-resolution ($\sim 10$ pc) hydrodynamic simulations of galaxy formation. Resolution-dependence is a primary consideration in our application of the method and simulations of isolated explosions (performed at different resolutions) motivate a resolution-dependent scaling for the injected fraction of kinetic energy that we apply in cosmological simulations of a $10^9$ Msun dwarf halo. We find that in high density media ($\gtrsim$ 50 cm$^{-3}$) with coarse resolution ($\gtrsim 4$ pc per cell), results are sensitive to the initial kinetic energy fraction due to early and rapid cooling. In our galaxy simulations, the deposition of small amounts of supernova energy in kinetic form (as little as 1%) has a dramatic impact on the evolution of the system, resulting in an order of magnitude suppression of stellar mass. The overall behavior of the galaxy in the two highest resolution simulations we perform appears to converge. We discuss the resulting distribution of stellar metallicities, an observable sensitive to galactic wind properties, and find that while the new method demonstrates increased agreement with observed systems, significant discrepancies remain, likely due to simplistic assumptions that neglect contributions from Type Ia supernovae and stellar winds.
[68]  [pdf] - 1245792
The Subaru-XMM-Newton Deep Survey (SXDS) VIII.: Multi-wavelength Identification, Optical/NIR Spectroscopic Properties, and Photometric Redshifts of X-ray Sources
Comments: Accepted for publication in PASJ Subaru special issue. 42 pages, 22 figures. Entire contents of Tables 3, 8, 9, 10, and 11, and ASCII format tables are available from
Submitted: 2015-05-12
We report the multi-wavelength identification of the X-ray sources found in the Subaru-XMM-Newton Deep Survey (SXDS) using deep imaging data covering the wavelength range between the far-UV to the mid-IR. We select a primary counterpart of each X-ray source by applying the likelihood ratio method to R-band, 3.6micron, near-UV, and 24micron source catalogs as well as matching catalogs of AGN candidates selected in 1.4GHz radio and i'-band variability surveys. Once candidates of Galactic stars, ultra-luminous X-ray sources in a nearby galaxy, and clusters of galaxies are removed there are 896 AGN candidates in the sample. We conduct spectroscopic observations of the primary counterparts with multi-object spectrographs in the optical and NIR; 65\% of the X-ray AGN candidates are spectroscopically-identified. For the remaining X-ray AGN candidates, we evaluate their photometric redshift with photometric data in 15 bands. Utilising the multi-wavelength photometric data of the large sample of X-ray selected AGNs, we evaluate the stellar masses, M*, of the host galaxies of the narrow-line AGNs. The distribution of the stellar mass is remarkably constant from z=0.1 to 4.0. The relation between M* and 2--10 keV luminosity can be explained with strong cosmological evolution of the relationship between the black hole mass and M*. We also evaluate the scatter of the UV-MIR spectral energy distribution (SED) of the X-ray AGNs as a function of X-ray luminosity and absorption to the nucleus. The scatter is compared with galaxies which have redshift and stellar mass distribution matched with the X-ray AGN. The UV-NIR SEDs of obscured X-ray AGNs are similar to those of the galaxies in the matched sample. In the NIR-MIR range, the median SEDs of X-ray AGNs are redder, but the scatter of the SEDs of the X-ray AGN broadly overlaps that of the galaxies in the matched sample.
[69]  [pdf] - 1055775
High-resolution mass models of dwarf galaxies from LITTLE THINGS
Comments: 100 pages, 85 figures, Accepted for publication on AJ
Submitted: 2015-02-04
We present high-resolution rotation curves and mass models of 26 dwarf galaxies from LITTLE THINGS. LITTLE THINGS is a high-resolution Very Large Array HI survey for nearby dwarf galaxies in the local volume within 11 Mpc. The rotation curves of the sample galaxies derived in a homogeneous and consistent manner are combined with Spitzer archival 3.6 micron and ancillary optical U, B, and V images to construct mass models of the galaxies. We decompose the rotation curves in terms of the dynamical contributions by baryons and dark matter halos, and compare the latter with those of dwarf galaxies from THINGS as well as Lambda CDM SPH simulations in which the effect of baryonic feedback processes is included. Being generally consistent with THINGS and simulated dwarf galaxies, most of the LITTLE THINGS sample galaxies show a linear increase of the rotation curve in their inner regions, which gives shallower logarithmic inner slopes alpha of their dark matter density profiles. The mean value of the slopes of the 26 LITTLE THINGS dwarf galaxies is alpha =-0.32 +/- 0.24 which is in accordance with the previous results found for low surface brightness galaxies (alpha = -0.2 +/- 0.2) as well as the seven THINGS dwarf galaxies (alpha =-0.29 +/- 0.07). However, this significantly deviates from the cusp-like dark matter distribution predicted by dark-matter-only Lambda CDM simulations. Instead our results are more in line with the shallower slopes found in the Lambda CDM SPH simulations of dwarf galaxies in which the effect of baryonic feedback processes is included. In addition, we discuss the central dark matter distribution of DDO 210 whose stellar mass is relatively low in our sample to examine the scenario of inefficient supernova feedback in low mass dwarf galaxies predicted from recent Lambda SPH simulations of dwarf galaxies where central cusps still remain.
[70]  [pdf] - 910671
The Interplay between SF and AGN Activity, and its role in Galaxy Evolution
Comments: 16 pages, 3 figures to appear as part of "Continuum Science" in proceedings of "Advancing Astrophysics with the Square Kilometre Array" PoS(AASKA14)083
Submitted: 2014-12-18
It has become apparent that active galactic nuclei (AGN) may have a significant impact on the growth and evolution of their host galaxies and vice versa but a detailed understanding of the interplay between these processes remains elusive. Deep radio surveys provide a powerful, obscuration-independent tool for measuring both star formation and AGN activity in high-redshift galaxies. Multiwavelength studies of deep radio fields show a composite population of star-forming galaxies and AGN, with the former dominating at the lowest flux densities (S$_{1.4\mathrm{GHz}}<$100~$\mu$Jy). The sensitivity and resolution of the SKA will allow us to identify, and separately trace, the total star formation in the bulges of individual high-redshift galaxies, the related nuclear activity and any star formation occurring on larger scales within a disc. We will therefore gain a detailed picture of the apparently simultaneous development of stellar populations and black holes in the redshift range where both star-formation and AGN activity peak (1$\leq$z$\leq$4). In this chapter we discuss the role of the SKA in studying the connection between AGN activity and galaxy evolution, and the most critical technical requirements for such of studies
[71]  [pdf] - 910666
The star-formation history of the Universe with the SKA
Comments: 19 pages, to appear as part of 'Continuum Science' in Proceedings 'Advancing Astrophysics with the SKA (AASKA14)068'
Submitted: 2014-12-18
Radio wavelengths offer the unique possibility of tracing the total star-formation rate in galaxies, both obscured and unobscured. As such, they may provide the most robust measurement of the star-formation history of the Universe. In this chapter we highlight the constraints that the SKA can place on the evolution of the star-formation history of the Universe, the survey area required to overcome sample variance, the spatial resolution requirements, along with the multi-wavelength ancillary data that will play a major role in maximising the scientific promise of the SKA. The required combination of depth and resolution means that a survey to trace the star formation in the Universe should be carried out with a facility that has a resolution of at least ~0.5arcsec, with high sensitivity at < 1 GHz. We also suggest a strategy that will enable new parameter space to be explored as the SKA expands over the coming decade.
[72]  [pdf] - 936239
The spectral energy distribution of the redshift 7.1 quasar ULAS J1120+0641
Comments: Accepted for publication in A&A. 10 pages, 5 figures, 5 tables
Submitted: 2014-11-20
We present new observations of the highest-redshift quasar known, ULAS J1120+0641, redshift $z=7.084$, obtained in the optical, at near-, mid-, and far-infrared wavelengths, and in the sub-mm. We combine these results with published X-ray and radio observations to create the multiwavelength spectral energy distribution (SED), with the goals of measuring the bolometric luminosity $L_{\rm bol}$, and quantifying the respective contributions from the AGN and star formation. We find three components are needed to fit the data over the wavelength range $0.12-1000\,\mu$m: the unobscured quasar accretion disk and broad-line region, a dusty clumpy AGN torus, and a cool 47K modified black body to characterise star formation. Despite the low signal-to-noise ratio of the new long-wavelength data, the normalisation of any dusty torus model is constrained within $\pm40\%$. We measure a bolometric luminosity $L_{\rm bol}=2.6\pm0.6\times10^{47}\,$erg$\,$s$^{-1}=6.7 \pm 1.6\times10^{13}L_{\odot}$, to which the three components contribute $31\%,32\%,3\%$, respectively, with the remainder provided by the extreme UV $<0.12\,\mu$m. We tabulate the best-fit model SED. We use local scaling relations to estimate a star formation rate (SFR) in the range $60-270\,{\rm M}_\odot$/yr from the [C$\,{\scriptsize \rm II}$] line luminosity and the $158\,\mu$m continuum luminosity. An analysis of the equivalent widths of the [C$\,{\scriptsize \rm II}$] line in a sample of $z>5.7$ quasars suggests that these indicators are promising tools for estimating the SFR in high-redshift quasars in general. At the time observed the black hole was growing in mass more than 100 times faster than the stellar bulge, relative to the mass ratio measured in the local universe, i.e. compared to ${M_{\rm BH}}/{M_{\rm bulge}} \simeq 1.4\times10^{-3}$, for ULAS J1120+0641 we measure ${\dot{M}_{\rm BH}}/{\dot{M}_{\rm bulge}} \simeq 0.2$.
[73]  [pdf] - 1223235
The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the rest-frame far-infrared emission of sub-millimeter galaxies
Comments: 16 pages, 7 figures, accepted for publication in ApJ
Submitted: 2014-11-18
We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.
[74]  [pdf] - 1209644
Dwarf galaxies in CDM and SIDM with baryons: observational probes of the nature of dark matter
Comments: 16 pages, 11 figures. MNRAS accepted
Submitted: 2014-05-20, last modified: 2014-09-20
We present the first cosmological simulations of dwarf galaxies, which include dark matter self-interactions and baryons. We study two dwarf galaxies within cold dark matter, and four different elastic self-interacting scenarios with constant and velocity-dependent cross sections, motivated by a new force in the hidden dark matter sector. Our highest resolution simulation has a baryonic mass resolution of $1.8\times 10^2\,{\rm M}_\odot$ and a gravitational softening length of $34\,{\rm pc}$ at $z=0$. In this first study we focus on the regime of mostly isolated dwarf galaxies with halo masses $\sim10^{10}\,{\rm M}_\odot$ where dark matter dynamically dominates even at sub-kpc scales. We find that while the global properties of galaxies of this scale are minimally affected by allowed self-interactions, their internal structures change significantly if the cross section is large enough within the inner sub-kpc region. In these dark-matter-dominated systems, self-scattering ties the shape of the stellar distribution to that of the dark matter distribution. In particular, we find that the stellar core radius is closely related to the dark matter core radius generated by self-interactions. Dark matter collisions lead to dwarf galaxies with larger stellar cores and smaller stellar central densities compared to the cold dark matter case. The central metallicity within $1\,{\rm kpc}$ is also larger by up to $\sim 15\%$ in the former case. We conclude that the mass distribution, and characteristics of the central stars in dwarf galaxies can potentially be used to probe the self-interacting nature of dark matter.
[75]  [pdf] - 1217040
The HI Chronicles of LITTLE THINGS BCDs II: The Origin of IC 10's HI Structure
Comments: 36 pages, 17 figures, accepted for publication in The Astronomical Journal
Submitted: 2014-09-18
In this paper we analyze Very Large Array (VLA) telescope and Green Bank Telescope (GBT) atomic hydrogen (HI) data for the LITTLE THINGS(1) blue compact dwarf galaxy IC 10. The VLA data allow us to study the detailed HI kinematics and morphology of IC 10 at high resolution while the GBT data allow us to search the surrounding area at high sensitivity for tenuous HI. IC 10's HI appears highly disturbed in both the VLA and GBT HI maps with a kinematically distinct northern HI extension, a kinematically distinct southern plume, and several spurs in the VLA data that do not follow the general kinematics of the main disk. We discuss three possible origins of its HI structure and kinematics in detail: a current interaction with a nearby companion, an advanced merger, and accretion of intergalactic medium. We find that IC 10 is most likely an advanced merger or a galaxy undergoing accretion. 1:Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey;
[76]  [pdf] - 1216601
The cosmic evolution of radio-AGN feedback to z=1
Comments: Accepted for publication in MNRAS
Submitted: 2014-08-31
This paper presents the first measurement of the radio luminosity function of 'jet-mode' (radiatively-inefficient) radio-AGN out to z=1, in order to investigate the cosmic evolution of radio-AGN feedback. Eight radio source samples are combined to produce a catalogue of 211 radio-loud AGN with 0.5<z<1.0, which are spectroscopically classified into jet-mode and radiative-mode (radiatively-efficient) AGN classes. Comparing with large samples of local radio-AGN from the Sloan Digital Sky Survey, the cosmic evolution of the radio luminosity function of each radio-AGN class is independently derived. Radiative-mode radio-AGN show an order of magnitude increase in space density out to z~1 at all luminosities, consistent with these AGN being fuelled by cold gas. In contrast, the space density of jet-mode radio-AGN decreases with increasing redshift at low radio luminosities (L_1.4 < 1e24 W/Hz) but increases at higher radio luminosities. Simple models are developed to explain the observed evolution. In the best-fitting models, the characteristic space density of jet-mode AGN declines with redshift in accordance with the declining space density of massive quiescent galaxies, which fuel them via cooling of gas in their hot haloes. A time delay of 1.5-2 Gyr may be present between the quenching of star formation and the onset of jet-mode radio-AGN activity. The behaviour at higher radio luminosities can be explained either by an increasing characteristic luminosity of jet-mode radio-AGN activity with redshift (roughly as (1+z) cubed) or if the jet-mode radio-AGN population also includes some contribution of cold-gas-fuelled sources seen at a time when their accretion rate was low. Higher redshifts measurements would distinguish between these possibilities.
[77]  [pdf] - 903580
The SCUBA-2 Cosmology Legacy Survey: the submillimetre properties of Lyman break galaxies at z=3-5
Comments: 12 pages, 8 figures, submitted to MNRAS
Submitted: 2014-07-24
We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.
[78]  [pdf] - 841640
Galactic conformity and central / satellite quenching, from the satellite profiles of M$^{\ast}$ galaxies at $0.4<z<1.9$ in the UKIDSS UDS
Comments: Submitted to MNRAS. 25 pages, 13 figures, 1 table. Comments welcome. For readers with very little time, the central result of the paper is covered by Figures 9-11 (sections 5 and 6)
Submitted: 2014-06-23
We explore the redshift evolution of a curious correlation between the star-formation properties of central galaxies and their satellites (`galactic conformity') at intermediate to high redshift ($0.4<z<1.9$). Using an extremely deep near-infrared survey, we study the distribution and properties of satellite galaxies with stellar masses, ${\rm log} ({\rm M}_*/{\rm M}_{\odot})>9.7$, around central galaxies at the characteristic Schechter function mass, ${\rm M} \sim {\rm M}^{\ast}$. We fit the radial profiles of satellite number densities with simple power laws, finding slopes in the range -1.1 to -1.4 for mass-selected satellites, and -1.3 to -1.6 for passive satellites. We confirm the tendency for passive satellites to be preferentially located around passive central galaxies at $3\sigma$ significance and show that it exists to at least $z\sim2$. Meanwhile, the quenched fraction of satellites around star-forming galaxies is consistent with field galaxies of equal stellar masses. We find no convincing evidence for a redshift-dependent evolution of these trends. One simple interpretation of these results is that only passive central galaxies occupy an environment that is capable of independently shutting off star-formation in satellite galaxies. By examining the satellites of higher stellar mass star-forming galaxies (${\rm log} ({\rm M}_*/{\rm M}_{\odot}) > 11$), we conclude that the origin of galactic conformity is unlikely to be exclusively due to the host dark-matter halo mass. A halo-mass-independent correlation could be established by either formation bias or a more physical connection between central and satellite star-formation histories. For the latter, we argue that a star-formation (or AGN) related outburst event from the central galaxy could establish a hot halo environment which is then capable of quenching both central and satellite galaxies.
[79]  [pdf] - 1208479
Diffuse gas properties and stellar metallicities in cosmological simulations of disc galaxy formation
Comments: 18 pages, 10 figures, 1 table, accepted for publication in MNRAS. Added 1 table, updated figures 1, 2, 4 and 6 and a few text changes to match the accepted version
Submitted: 2014-03-19, last modified: 2014-06-06
We analyse the properties of the circum-galactic medium and the metal content of the stars comprising the central galaxy in eight hydrodynamical `zoom-in' simulations of disc galaxy formation. We use these properties as a benchmark for our model of galaxy formation physics implemented in the moving-mesh code AREPO, which succeeds in forming quite realistic late-type spirals in the set of `Aquarius' initial conditions of Milky Way-sized haloes. Galactic winds significantly influence the morphology of the circum-galactic medium and induce bipolar features in the distribution of heavy elements. They also affect the thermodynamic properties of the circum-galactic gas by supplying an energy input that sustains its radiative losses. Although a significant fraction of the heavy elements are transferred from the central galaxy to the halo, and even beyond the virial radius, enough metals are retained by stars to yield a peak in their metallicity distributions at about $Z_{\odot}$. All our default runs overestimate the stellar [O/Fe] ratio, an effect that we demonstrate can be rectified by an increase of the adopted SN type Ia rate. Nevertheless, the models have difficulty in producing stellar metallicity gradients of the same strength as observed in the Milky Way.
[80]  [pdf] - 1180698
Black hole mass estimates and emission-line properties of a sample of redshift z>6.5 quasars
Comments: 18 pages, 11 figures. Submitted to ApJ
Submitted: 2013-11-13, last modified: 2014-06-05
We present the analysis of optical and near-infrared spectra of the only four $z>6.5$ quasars known to date, discovered in the UKIDSS-LAS and VISTA-VIKING surveys. Our data-set consists of new VLT/X-Shooter and Magellan/FIRE observations. These are the best optical/NIR spectroscopic data that are likely to be obtained for the $z>6.5$ sample using current $6$ - $10$ m facilities. We estimate the black hole mass, the Eddington ratio, and the SiIV/CIV, CIII]/CIV, and FeII/MgII emission-line flux ratios. We perform spectral modeling using a procedure that allows us to derive a probability distribution for the continuum components and to obtain the quasar properties weighted upon the underlying distribution of continuum models. The $z>6.5$ quasars show the same emission properties as their counterparts at lower redshifts. The $z>6.5$ quasars host black holes with masses of $\sim 10^9$ M$_{\odot}$ that are accreting close to the Eddington luminosity ($\langle{\rm log} (L_{\rm Bol}/L_{\rm Edd})\rangle= -0.4\pm0.2$), in agreement with what has been observed for a sample of $4.0<z<6.5$ quasars. By comparing the SiIV/CIV and CIII]/CIV flux ratios with the results obtained from luminosity-matched samples at $z\sim6$ and $2\leq z\leq4.5$, we find no evidence of evolution of the line ratios with cosmic time. We compare the measured FeII/MgII flux ratios with those obtained for a sample of $4.0<z<6.4$ sources. The two samples are analyzed using a consistent procedure. There is no evidence that the FeII/MgII flux ratio evolves between $z=7$ and $z=4$. Under the assumption that the FeII/MgII traces the Fe/Mg abundance ratio, this implies the presence of major episodes of chemical enrichment in the quasar hosts in the first $\sim0.8$ Gyr after the Big Bang.
[81]  [pdf] - 1209924
No excess of bright galaxies around the redshift 7.1 quasar ULAS J1120+0641
Comments: Accepted for publication in MNRAS, 8 pages
Submitted: 2014-06-03
We present optical and near-infrared imaging of the field of the z=7.0842 quasar ULAS J112001.48+064124.3 taken with the Hubble Space Telescope. We use these data to search for galaxies that may be physically associated with the quasar, using the Lyman break technique, and find three such objects, although the detection of one in Spitzer Space Telescope imaging strongly suggests it lies at z~2. This is consistent with the field luminosity function and indicates that there is no excess of >L* galaxies within 1 Mpc of the quasar. A detection of the quasar shortward of the Ly-alpha line is consistent with the previously observed evolution of the intergalactic medium at z>5.5.
[82]  [pdf] - 1209407
A surprising consistency between the far-infrared galaxy luminosity functions of the field and Coma
Comments: 8 pages, 5 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society
Submitted: 2014-05-07
We present new deep images of the Coma Cluster from the ESA Herschel Space Observatory at wavelengths of 70, 100 and 160 microns, covering an area of 1.75 x 1.0 square degrees encompassing the core and southwest infall region. Our data display an excess of sources at flux densities above 100 mJy compared to blank-field surveys, as expected. We use extensive optical spectroscopy of this region to identify cluster members and hence produce cluster luminosity functions in all three photometric bands. We compare our results to the local field galaxy luminosity function, and the luminosity functions from the Herschel Virgo Cluster Survey (HeViCS). We find consistency between the shapes of the Coma and field galaxy luminosity functions at all three wavelengths, however we do not find the same level of agreement with that of the Virgo Cluster.
[83]  [pdf] - 1180567
X-rays from the redshift 7.1 quasar ULAS J1120+0641
Comments: Accepted for publication in MNRAS Letters
Submitted: 2013-11-07, last modified: 2014-02-03
We present X-ray imaging and spectroscopy of the redshift z=7.084 radio-quiet quasar ULAS J112001.48+064124.3 obtained with Chandra and XMM-Newton. The quasar is detected as a point source with both observatories. The Chandra observation provides a precise position, confirming the association of the X-ray source and the quasar, while a sufficient number of photons is detected in the XMM-Newton observation to yield a meaningful X-ray spectrum. In the XMM-Newton observation the quasar has a 2-10 keV luminosity of 4.7+-0.9 times 10^44 erg/s and a spectral slope alpha = 1.6+0.4/-0.3 (where f_nu is proportional to nu^-alpha). The quasar appears to have dimmed in the 15 months between the two observations, with a 2-10 keV luminosity of $1.8+1.0/-0.7 times 10^45 erg/s during the Chandra observation. We derive optical to X-ray spectral slopes alpha_OX of 1.76+-0.07 and 1.54+0.09/-0.08 at the times of the XMM-Newton and Chandra observations respectively, consistent with the range of alpha_OX found in other quasars of comparable ultraviolet luminosity. The very soft X-ray spectrum suggests that the quasar is accreting above the Eddington rate, L/L_Edd = 5+15/-4, compared to L/L_Edd = 1.2+0.6/-0.5 derived from the rest-frame ultraviolet. Super-Eddington accretion would help to reduce the discrepancy between the age of the quasar implied by the small size of the ionized near zone in which it sits (<10^7 years), and the characteristic e-folding time (2.5 times 10^7 years if L/L_Edd=2). Such super-Eddington accretion would also alleviate the challenging constraints on the seed black hole mass provided that the quasar has been rapidly accreting throughout its history. The remnant of an individual population III star is a plausible progenitor if an average L/L_Edd>1.46 has been maintained over the quasar's lifetime.
[84]  [pdf] - 771971
A JVLA 10~degree^2 deep survey
Comments: 18 pages, 6 figures. White Paper submitted to the call for Very Large Array Sky Surveys
Submitted: 2014-01-16
(Abridged)One of the fundamental challenges for astrophysics in the 21st century is finding a way to untangle the physical processes that govern galaxy formation and evolution. Given the importance and scope of this problem, the multi-wavelength astronomical community has used the past decade to build up a wealth of information over specific extragalactic deep fields to address key questions in galaxy formation and evolution. These fields generally cover at least 10square degrees to facilitate the investigation of the rarest, typically most massive, galaxies and AGN. Furthermore, such areal coverage allows the environments to be fully accounted for, thereby linking the single halo to the two-halo terms in the halo occupation distribution. Surveys at radio wavelengths have begun to lag behind those at other wavelengths, especially in this medium-deep survey tier. However, the survey speed offered by the JVLA means that we can now reach a point where we can begin to obtain commensurate data at radio wavelengths to those which already exists from the X-ray through to the far-infrared over ~10 square degrees. We therefore present the case for a 10 square degree survey to 1.5uJy at L-band in A or B Array, requiring ~4000 hours to provide census of star-formation and AGN-accretion activity in the Universe. For example, the observations will allow galaxies forming stars at 10Msolar/yr to be detected out to z~1 and luminous infrared galaxies (1000Msolar/yr to be found out to z~6. Furthermore, the survey area ensures that we will have enough cosmic volume to find these rare sources at all epochs. The bandwidth will allow us to determine the polarisation properties galaxies in the high-redshift Universe as a function of stellar mass, morphology and redshift.
[85]  [pdf] - 1202190
The SCUBA-2 Cosmology Legacy Survey: Ultraluminous star-forming galaxies in a z=1.6 cluster
Comments: ApJ in press, 16 pages, 8 figures, 1 table
Submitted: 2013-12-11
We analyse new SCUBA-2 submillimeter and archival SPIRE far-infrared imaging of a z=1.62 cluster, Cl0218.3-0510, which lies in the UKIDSS/UDS field of the SCUBA-2 Cosmology Legacy Survey. Combining these tracers of obscured star formation activity with the extensive photometric and spectroscopic information available for this field, we identify 31 far-infrared/submillimeter-detected probable cluster members with bolometric luminosities >1e12 Lo and show that by virtue of their dust content and activity, these represent some of the reddest and brightest galaxies in this structure. We exploit Cycle-1 ALMA submillimeter continuum imaging which covers one of these sources to confirm the identification of a SCUBA-2-detected ultraluminous star-forming galaxy in this structure. Integrating the total star-formation activity in the central region of the structure, we estimate that it is an order of magnitude higher (in a mass-normalised sense) than clusters at z~0.5-1. However, we also find that the most active cluster members do not reside in the densest regions of the structure, which instead host a population of passive and massive, red galaxies. We suggest that while the passive and active populations have comparable near-infrared luminosities at z=1.6, M(H)~-23, the subsequent stronger fading of the more active galaxies means that they will evolve into passive systems at the present-day which are less luminous than the descendants of those galaxies which were already passive at z~1.6 (M(H)~-20.5 and M(H)~-21.5 respectively at z~0). We conclude that the massive galaxy population in the dense cores of present-day clusters were already in place at z=1.6 and that in Cl0218.3-0510 we are seeing continuing infall of less extreme, but still ultraluminous, star-forming galaxies onto a pre-existing structure.
[86]  [pdf] - 1172214
First Systematic Search for Oxygen-Line Blobs at High Redshift: Uncovering AGN Feedback and Star-Formation Quenching
Comments: 13 pages, 11 figures, 4 tables, Accepted for Publication in ApJ
Submitted: 2013-06-21, last modified: 2013-12-06
We present the first systematic search for extended metal-line [OII]{\lambda}{\lambda}3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]{\lambda}3426 line at the 6{\sigma} level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeII{\lambda}2587 absorption and FeII*{\lambda}2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.
[87]  [pdf] - 1180296
Evidence for an Interaction in the Nearest Starbursting Dwarf Irregular Galaxy IC 10
Comments: 7 pages, 4 figures, 1 table. Accepted for publication in ApJ Letters
Submitted: 2013-10-28
Using deep 21-cm HI data from the Green Bank Telescope we have detected an ~18.3 kpc-long gaseous extension associated with the starbursting dwarf galaxy IC 10. The newly-found feature stretches 1.3 deg to the northwest and has a large radial velocity gradient reaching to ~65 km/s lower than the IC 10 systemic velocity. A region of higher column density at the end of the extension that possesses a coherent velocity gradient (~10 km/s across ~26 arcmin) transverse to the extension suggests rotation and may be a satellite galaxy of IC 10. The HI mass of IC 10 is 9.5x10^7 (d/805 kpc)^2 Msun and the mass of the new extension is 7.1x10^5 (d/805 kpc)^2 Msun. An IC 10-M31 orbit using known radial velocity and proper motion values for IC 10 show that the HI extension is inconsistent with the trailing portion of the orbit so that an M31-tidal or ram pressure origin seems unlikely. We argue that the most plausible explanation for the new feature is that it is the result of a recent interaction (and possible late merger) with another dwarf galaxy. This interaction could not only have triggered the origin of the recent starburst in IC 10, but could also explain the existence of previously-found counter-rotating HI gas in the periphery of the IC 10 which was interpreted as originating from primordial gas infall.
[88]  [pdf] - 1172667
Evidence for a correlation between the sizes of quiescent galaxies and local environment to z ~ 2
Comments: Accepted for publication in MNRAS. 16 pages, 11 figures, 6 tables
Submitted: 2013-07-11, last modified: 2013-07-15
We present evidence for a strong relationship between galaxy size and environment for the quiescent population in the redshift range 1 < z < 2. Environments were measured using projected galaxy overdensities on a scale of 400 kpc, as determined from ~ 96,000 K-band selected galaxies from the UKIDSS Ultra Deep Survey (UDS). Sizes were determined from ground-based K-band imaging, calibrated using space-based CANDELS HST observations in the centre of the UDS field, with photometric redshifts and stellar masses derived from 11-band photometric fitting. From the resulting size-mass relation, we confirm that quiescent galaxies at a given stellar mass were typically ~ 50 % smaller at z ~ 1.4 compared to the present day. At a given epoch, however, we find that passive galaxies in denser environments are on average significantly larger at a given stellar mass. The most massive quiescent galaxies (M_stellar > 2 x 10^11 M_sun) at z > 1 are typically 50 % larger in the highest density environments compared to those in the lowest density environments. Using Monte Carlo simulations, we reject the null hypothesis that the size-mass relation is independent of environment at a significance > 4.8 sigma for the redshift range 1 < z < 2. In contrast, the evidence for a relationship between size and environment is much weaker for star-forming galaxies.
[89]  [pdf] - 1172577
Enzo: An Adaptive Mesh Refinement Code for Astrophysics
Comments: 60 pages, 22 figures, submitted to ApJS. The code is available at
Submitted: 2013-07-08
This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.
[90]  [pdf] - 1172444
The HI Chronicles of LITTLE THINGS BCDs: Evidence for External Perturbations in the Morphology and Kinematics of Haro 29 and Haro 36
Comments: To be published in The Astronomical Journal
Submitted: 2013-07-02
We analyze high angular and velocity resolution HI-line data of two LITTLE THINGS (1) blue compact dwarfs (BCDs): Haro 29 and Haro 36. Both of these BCDs are disturbed morphologically and kinematically. Haro 29's HI data reveal a kinematic major axis that is offset from the optical major axis, and a disturbed outer HI component, indicating that Haro 29 may have had a past interaction. Position-velocity diagrams of Haro 36 indicate that it has two kinematically separate components at its center and a likely tidal tail in front of the galaxy. We find that Haro 36 most likely had an interaction in the past, is currently interacting with an unknown companion, or is a merger remnant. (1) "Local Irregulars That Trace Luminosity Extremes The HI Nearby Galaxy Survey"
[91]  [pdf] - 1172408
The prevalence of AGN feedback in massive galaxies at z~1
Comments: 10 pages, MNRAS, in press
Submitted: 2013-07-01
We use the optical--infrared imaging in the UKIDSS Ultra Deep Survey field, in combination with the new deep radio map of Arumugam et al., to calculate the distribution of radio luminosities among galaxies as a function of stellar mass in two redshift bins across the interval 0.4<z<1.2. This is done with the use of a new Bayesian method to classify stars and galaxies in surveys with multi-band photometry, and to derive photometric redshifts and stellar masses for those galaxies. We compare the distribution to that observed locally and find agreement if we consider only objects believed to be weak-lined radio-loud galaxies. Since the local distribution is believed to be the result of an energy balance between radiative cooling of the gaseous halo and mechanical AGN heating, we infer that this balance was also present as long ago as z~1. This supports the existence of a direct link between the presence of a low-luminosity ('hot-mode') radio-loud active galactic nucleus and the absence of ongoing star formation.
[92]  [pdf] - 1166443
Evolution of Star Formation in the UKIDSS Ultra Deep Survey Field - I. Luminosity Functions and Cosmic Star Formation Rate out to z=1.6
Comments: Accepted for publication in MNRAS. 16 pages, 11 figures (size of Figures 2c, 2d and 3b edited for Astro-ph)
Submitted: 2013-05-06
We present new results on the cosmic star formation history in the SXDS-UDS field out to z=1.6. We compile narrow-band data from the Subaru Telescope and the Visible and Infrared Survey Telescope for Astronomy (VISTA) in conjunction with broad-band data from the SXDS and UDS, to make a selection of 5725 emission-line galaxies in 12 redshift slices, spanning 10 Gyr of cosmic time. We determine photometric redshifts for the sample using 11-band photometry, and use a spectroscopically confirmed subset to fine tune the resultant redshift distribution. We use the maximum-likelihood technique to determine luminosity functions in each redshift slice and model the selection effects inherent in any narrow-band selection statistically, to obviate the retrospective corrections ordinarily required. The deep narrow-band data are sensitive to very low star formation rates (SFRs), and allow an accurate evaluation of the faint end slope of the Schechter function, alpha. We find that alpha is particularly sensitive to the assumed faintest broadband magnitude of a galaxy capable of hosting an emission line, and propose that this limit should be empirically motivated. For this analysis we base our threshold on the limiting observed equivalent widths of emission lines in the local Universe. We compute the characteristic SFR of galaxies in each redshift slice, and the integrated SFR density, rho_{SFR}. We find our results to be in good agreement with the literature and parametrize the evolution of the SFR density as rho_{SFR} proportional to (1 + z)^{4.58} confirming a steep decline in star formation activity since z = 1.6.
[93]  [pdf] - 662162
CANDELS Multiwavelength catalogs: Source Identification and Photometry in the CANDELS UKIDSS Ultra-Deep Survey Field
Comments: 22 pages, 20 figures, ApJ in press. The CANDELS UDS catalog (in ascii & fits format) and associated ReadMe file are made available on the CANDELS website ( as well as on the MAST CANDELS page (
Submitted: 2013-05-06
We present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources.
[94]  [pdf] - 1166255
High velocity outflows from young star-forming galaxies in the UKIDSS Ultra-Deep Survey
Comments: 16 pages, 15 figures, accepted by MNRAS
Submitted: 2013-04-26
We investigate galactic-scale outflows in the redshift range 0.71 < z < 1.63, using 413 K-band selected galaxies observed in the spectroscopic follow-up of the UKIDSS Ultra-Deep Survey (UDSz). The galaxies have an average stellar mass of ~10^9.5 solar masses and span a wide range in rest-frame colours, representing typical star-forming galaxies at this epoch. We stack the spectra by various galaxy properties, including stellar mass, [OII] equivalent width, star-formation rate, specific star-formation rate and rest-frame spectral indices. We find that outflows are present in virtually all spectral stacks, with velocities ranging from 100-1000 km s^-1, indicating that large-scale outflowing winds are a common property at these redshifts. The highest velocity outflows (>500 km s^-1) are found in galaxies with the highest stellar masses and the youngest stellar populations. Our findings suggest that high velocity galactic outflows are mostly driven by star-forming processes rather than AGN, with implied mass outflow rates comparable to the rates of star formation. Such behaviour is consistent with models required to reproduce the high-redshift mass-metallicity relation.
[95]  [pdf] - 640964
The effect of feedback and reionization on star formation in low-mass dwarf galaxy haloes
Comments: 26 pages, 14 figures, accepted for publication in MNRAS; additional explanation added to sections 2,3 and 4; Figure 4 updated; results unchanged
Submitted: 2012-11-05, last modified: 2013-03-18
We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc.
[96]  [pdf] - 1164987
Studying the emergence of the red sequence through galaxy clustering: host halo masses at z > 2
Comments: Accepted for publication in MNRAS. 16 pages, 9 figures, 1 table
Submitted: 2013-03-04, last modified: 2013-03-06
We use the UKIDSS Ultra-Deep Survey, the deepest degree-scale near-infrared survey to date, to investigate the clustering of star-forming and passive galaxies to z ~ 3.5. Our new measurements include the first determination of the clustering for passive galaxies at z > 2, which we achieve using a cross-correlation technique. We find that passive galaxies are the most strongly clustered, typically hosted by massive dark matter halos with M_halo > 5 x 10^12 M_sun irrespective of redshift or stellar mass. Our findings are consistent with models in which a critical halo mass determines the transition from star-forming to passive galaxies. Star-forming galaxies show no strong correlation between stellar mass and halo mass, but passive galaxies show evidence for an anti-correlation; low-mass passive galaxies appear, on average, to be located in the most massive halos. These results can be understood if the termination of star formation is most efficient for galaxies of low stellar mass in very dense environments.
[97]  [pdf] - 1157429
Black hole mass and Eddington ratio distribution functions of X-ray selected broad-line AGNs at z~1.4 in the Subaru XMM-Newton Deep Field
Comments: 23 pages with 18 figures. Accepted for publication in ApJ
Submitted: 2012-10-31
In order to investigate the growth of super-massive black holes (SMBHs), we construct the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of X-ray-selected broad-line AGNs at z~1.4 in the Subaru XMM-Newton Deep Survey field. In this redshift range, a significant part of the accretion growth of SMBHs is thought to be taking place. Black hole masses of X-ray-selected broad-line AGNs are estimated using the width of the broad MgII line and the 3000A monochromatic luminosity. We supplement the MgII FWHM values with the Ha FWHM obtained from our NIR spectroscopic survey. Using the black hole masses of broad-line AGNs at redshifts between 1.18 and 1.68, the binned broad-line AGN BHMF and ERDF are calculated using the Vmax method. To properly account for selection effects that impact the binned estimates, we derive the corrected broad-line AGN BHMF and ERDF by applying the Maximum Likelihood method, assuming that the ERDF is constant regardless of the black hole mass. We do not correct for the non-negligible uncertainties in virial BH mass estimates. If we compare the corrected broad-line AGN BHMF with that in the local Universe, the corrected BHMF at z~1.4 has a higher number density above 10^8 Msolar but a lower number density below that mass range. The evolution may be indicative of a down-sizing trend of accretion activity among the SMBH population. The evolution of broad-line AGN ERDF from z=1.4 to 0 indicates that the fraction of broad-line AGNs with accretion rate close to the Eddington-limit is higher at higher redshifts.
[98]  [pdf] - 1124227
The VISTA Deep Extragalactic Observations (VIDEO) Survey
Comments: 16 Pages, 12 figure and 5 tables, MNRAS in press. Figures have been degraded due to size constraints, Figure 10 also updated
Submitted: 2012-06-19, last modified: 2012-10-03
In this paper we describe the first data release of the the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey. VIDEO is a ~12degree^2 survey in the near-infrared Z,Y,J,H and K_s bands, specifically designed to enable the evolution of galaxies and large structures to be traced as a function of both epoch and environment from the present day out to z=4, and active galactic nuclei (AGN) and the most massive galaxies up to and into the epoch of reionization. With its depth and area, VIDEO will be able to fully explore the period in the Universe where AGN and starburst activity were at their peak and the first galaxy clusters were beginning to virialize. VIDEO therefore offers a unique data set with which to investigate the interplay between AGN, starbursts and environment, and the role of feedback at a time when it was potentially most crucial. We provide data over the VIDEO-XMM3 tile, which also covers the Canada-France-Hawaii-Telescope Legacy Survey Deep-1 field (CFHTLS-D1). The released VIDEO data reach a 5-sigma AB-magnitude depth of Z=25.7, Y=24.5, J=24.4, H=24.1 and K_s=23.8 in 2 arcsec diameter apertures (the full depth of Y=24.6 will be reached within the full integration time in future releases). The data are compared to previous surveys over this field and we find good astrometric agreement with the Two-Micron All Sky Survey, and source counts in agreement with the recently released UltraVISTA survey data. The addition of the VIDEO data to the CFHTLS-D1 optical data increases the accuracy of photometric redshifts and significantly reduces the fraction of catastrophic outliers over the redshift range 0<z<1 from 5.8 to 3.1 per cent in the absence of an i-band luminosity prior. (Truncated Abstract)
[99]  [pdf] - 968387
The Spitzer Extragalactic Representative Volume Survey (SERVS): survey definition and goals
Comments: 21 pages, 23 figures, published in PASP. This version differs slightly from PASP, mainly due to formatting issues. Figure 2 was also corrected from the previous version
Submitted: 2012-06-18, last modified: 2012-09-17
We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.
[100]  [pdf] - 1150952
Little Things
Comments: In press in AJ
Submitted: 2012-08-29
We present LITTLE THINGS (Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey) that is aimed at determining what drives star formation in dwarf galaxies. This is a multi-wavelength survey of 37 Dwarf Irregular and 4 Blue Compact Dwarf galaxies that is centered around HI-line data obtained with the National Radio Astronomy Observatory (NRAO) Very Large Array (VLA). The HI-line data are characterized by high sensitivity (less than 1.1 mJy/beam per channel), high spectral resolution (less than or equal to 2.6 km/s), and high angular resolution (~6 arcseconds. The LITTLE THINGS sample contains dwarf galaxies that are relatively nearby (less than or equal to 10.3 Mpc; 6 arcseconds is less than or equal to 300 pc), that were known to contain atomic hydrogen, the fuel for star formation, and that cover a large range in dwarf galactic properties. We describe our VLA data acquisition, calibration, and mapping procedures, as well as HI map characteristics, and show channel maps, moment maps, velocity-flux profiles, and surface gas density profiles. In addition to the HI data we have GALEX UV and ground-based UBV and Halpha images for most of the galaxies, and JHK images for some. Spitzer mid-IR images are available for many of the galaxies as well. These data sets are available on-line.
[101]  [pdf] - 1117575
Detection of atomic carbon [CII] 158 micron and dust emission from a z=7.1 quasar host galaxy
Comments: 7 pages, 3 figures. Published in ApJ Letters
Submitted: 2012-03-26, last modified: 2012-05-22
Using the IRAM Plateau de Bure Interferometer, we report the detection of the 158 micron [CII] emission line and underlying dust continuum in the host galaxy of the quasar ULAS J112001.48+064124.3 (hereafter J1120+0641) at z=7.0842+/-0.0004. This is the highest redshift detection of the [CII] line to date, and allows us to put first constraints on the physical properties of the host galaxy. The [CII] line luminosity is (1.2+/-0.2)x10^9 Lsun, which is a factor ~4 lower than observed in a luminous quasar at z=6.42 (SDSS J1148+5251). The underlying far-infrared (FIR) continuum has a flux density of 0.61+/-0.16 mJy, similar to the average flux density of z~6 quasars that were not individually detected in the rest-frame FIR. Assuming the FIR luminosity of L_FIR = 5.8x10^11-1.8x10^12 Lsun is mainly powered by star-formation, we derive a star-formation rate in the range 160-440 Msun/yr and a total dust mass in the host galaxy of 6.7x10^7-5.7x10^8 Msun (both numbers have significant uncertainties given the unknown nature of dust at these redshifts). The [CII] line width of sigma_V=100+/-15 km/s is among the smallest observed when compared to the molecular line widths detected in z~6 quasars. Both the [CII] and dust continuum emission are spatially unresolved at the current angular resolution of 2.0x1.7 arcsec^2 (corresponding to 10x9 kpc^2 at the redshift of J1120+0641).
[102]  [pdf] - 1092989
Radio imaging of the Subaru/XMM-Newton Deep Field - III. Evolution of the radio luminosity function beyond z=1
Comments: Accepted for publication in MNRAS: 36 pages, including 13 pages of figures to appear online only. In memory of Steve
Submitted: 2012-01-16
We present spectroscopic and eleven-band photometric redshifts for galaxies in the 100-uJy Subaru/XMM-Newton Deep Field radio source sample. We find good agreement between our redshift distribution and that predicted by the SKA Simulated Skies project. We find no correlation between K-band magnitude and radio flux, but show that sources with 1.4-GHz flux densities below ~1mJy are fainter in the near-infrared than brighter radio sources at the same redshift, and we discuss the implications of this result for spectroscopically-incomplete samples where the K-z relation has been used to estimate redshifts. We use the infrared--radio correlation to separate our sample into radio-loud and radio-quiet objects and show that only radio-loud hosts have spectral energy distributions consistent with predominantly old stellar populations, although the fraction of objects displaying such properties is a decreasing function of radio luminosity. We calculate the 1.4-GHz radio luminosity function (RLF) in redshift bins to z=4 and find that the space density of radio sources increases with lookback time to z~2, with a more rapid increase for more powerful sources. We demonstrate that radio-loud and radio-quiet sources of the same radio luminosity evolve very differently. Radio-quiet sources display strong evolution to z~2 while radio-loud AGNs below the break in the radio luminosity function evolve more modestly and show hints of a decline in their space density at z>1, with this decline occurring later for lower-luminosity objects. If the radio luminosities of these sources are a function of their black hole spins then slowly-rotating black holes must have a plentiful fuel supply for longer, perhaps because they have yet to encounter the major merger that will spin them up and use the remaining gas in a major burst of star formation.
[103]  [pdf] - 1083504
Ultra Steep Spectrum radio sources in the Lockman Hole: SERVS identifications and redshift distribution at the faintest radio fluxes
Comments: Accepted for publication in the Astrophysical Journal. 33 pages, 8 figures, 3 tables. v2: Corrected typos
Submitted: 2011-08-19, last modified: 2011-09-01
Ultra Steep Spectrum (USS) radio sources have been successfully used to select powerful radio sources at high redshifts (z>~2). Typically restricted to large-sky surveys and relatively bright radio flux densities, it has gradually become possible to extend the USS search to sub-mJy levels, thanks to the recent appearance of sensitive low-frequency radio facilities. Here a first detailed analysis of the nature of the faintest USS sources is presented. By using Giant Metrewave Radio Telescope and Very Large Array radio observations of the Lockman Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz integrated fluxes above 100 microJy, is assembled. Deep infrared data at 3.6 and 4.5 micron from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used to reliably identify counterparts for 48 (83%) of these sources, showing an average total magnitude of [3.6](AB)=19.8 mag. Spectroscopic redshifts for 14 USS sources, together with photometric redshift estimates, improved by the use of the deep SERVS data, for a further 19 objects, show redshifts ranging from z=0.1 to z=2.8, peaking at z~0.6 and tailing off at high redshifts. The remaining 25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, with 10 sources undetected at 3.6 and 4.5 micron (typically [3.6]>22-23 mag, from local measurements), which suggests the likely existence of higher redshifts among the sub-mJy USS population. The comparison with the Square Kilometre Array Design Studies Simulated Skies models indicate that Fanaroff-Riley type I radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of the faintest USS population, and raises the possibility that the high efficiency of the USS technique for the selection of high redshift sources remains even at the sub-mJy level.
[104]  [pdf] - 397808
VII Zw 403: HI Structure in a Blue Compact Dwarf Galaxy
Comments: 62 pages, 19 figures
Submitted: 2011-07-11
We present optical (UBVJ), ultraviolet (FUV, NUV), and high resolution atomic hydrogen (HI) observations of the nearby blue compact dwarf (BCD), VII Zw 403. We find that VII Zw 403 has a relatively high HI mass-to-light ratio for a BCD. The rotation velocity is nominally 10-15 km/s, but rises to ~20 km/s after correction for the ~8-10 km/s random motions present in the gas. The velocity field is complex; including a variation in the position angle of the major axis going from the NE to the SW parts of the galaxy. Our high resolution HI maps reveal structure in the central gas, including a large, low-density HI depression or hole between the southern and northern halves of the galaxy, coincident with an unresolved x-ray source. Although interactions have been proposed as the triggering mechanism for the vigorous star formation occurring in BCDs, VII Zw 403 does not seem to have been tidally triggered by an external interaction, as we have found no nearby possible perturbers. It also doesn't appear to fall in the set of galaxies that exhibit a strong central mass density concentration, as its optical scale length is large in comparison to similar systems. However, there are some features that are compatible with an accretion event: optical/HI axis misalignment, a change in position angle of the kinematic axis, and a complex velocity field.
[105]  [pdf] - 379812
A luminous quasar at a redshift of z = 7.085
Comments: 5 pages, 4 figures; source data available at
Submitted: 2011-06-29
The intergalactic medium was not completely reionized until approximately a billion years after the Big Bang, as revealed by observations of quasars with redshifts of less than 6.5. It has been difficult to probe to higher redshifts, however, because quasars have historically been identified in optical surveys, which are insensitive to sources at redshifts exceeding 6.5. Here we report observations of a quasar (ULAS J112001.48+064124.3) at a redshift of 7.085, which is 0.77 billion years after the Big Bang. ULAS J1120+0461 had a luminosity of 6.3x10^13 L_Sun and hosted a black hole with a mass of 2x10^9 M_Sun (where L_Sun and M_Sun are the luminosity and mass of the Sun). The measured radius of the ionized near zone around ULAS J1120+0641 is 1.9 megaparsecs, a factor of three smaller than typical for quasars at redshifts between 6.0 and 6.4. The near zone transmission profile is consistent with a Ly alpha damping wing, suggesting that the neutral fraction of the intergalactic medium in front of ULAS J1120+0641 exceeded 0.1.
[106]  [pdf] - 1077661
How neutral is the intergalactic medium surrounding the redshift z=7.085 quasar ULAS J1120+0641?
Comments: 5 pages, 4 figures, accepted to MNRAS letters
Submitted: 2011-06-29
The quasar ULAS J1120+0641 at redshift z=7.085 has a highly ionised near zone which is smaller than those around quasars of similar luminosity at z~6. The spectrum also exhibits evidence for a damping wing extending redward of the systemic Lya redshift. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionisation state of the inhomogeneous IGM surrounding this quasar. Our simulations show that the transmission profile is consistent with an IGM in the vicinity of the quasar with a volume averaged HI fraction of f_HI>0.1 and that ULAS J1120+0641 has been bright for 10^6--10^7 yr. The observed spectrum is also consistent with smaller IGM neutral fractions, f_HI ~ 10^-3--10-4, if a damped Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ~5 per cent of our simulated sight-lines for a bright phase of 10^6--10^7 yr. Unless ULAS J1120+0641 grows during a previous optically obscured phase, the low age inferred for the quasar adds to the theoretical challenge of forming a 2x10^9 M_sol black hole at this high redshift.
[107]  [pdf] - 1053390
AGN environments at z<1.5 in the UKIDSS Ultra-Deep Survey
Comments: 12 pages, 7 figures, Accepted for publication in MNRAS
Submitted: 2011-04-14, last modified: 2011-04-18
We investigate the environments of both X-ray and radio-loud AGN within the UKIDSS Ultra-deep Survey (UDS) using deep infrared selection to sample the galaxy density field in the redshift range 1.0 <= z <= 1.5. Using angular cross-correlation techniques we find that both X-ray and radio-loud AGN preferentially reside in overdense environments. We also find that both types of AGN cluster more strongly with those galaxies classified as `passive' rather than those that are actively star-forming. We infer clustering scale lengths comparable to those of passive red galaxies, suggesting that typical AGN at these epochs reside in dark-matter halos of mass M >~ 10^13 M_sun. A closer look at the small-scale environments of the AGN reveals that the neighbouring galaxies of radio-loud AGN have U-B colours more skewed towards the `green-valley' and the red sequence, whereas the neighbours of X-ray AGN show no difference to the general galaxy population. This suggests that although both AGN types live in overdense environments, the radio-loud AGN may be preferentially located in more evolved cluster cores, in a similar environment to low-powered radio AGN in the local Universe.
[108]  [pdf] - 958244
Probabilistic selection of high-redshift quasars
Comments: submitted to MNRAS; 20 pages, 13 figures
Submitted: 2011-01-25
High redshift quasars (HZQs) with redshifts of z >~ 6 are so rare that any photometrically-selected sample of sources with HZQ-like colours is likely to be dominated by Galactic stars and brown dwarfs scattered from the stellar locus. It is impractical to reobserve all such candidates, so an alternative approach was developed in which Bayesian model comparison techniques are used to calculate the probability that a candidate is a HZQ, P_q, by combining models of the quasar and star populations with the photometric measurements of the object. This method was motivated specifically by the large number of HZQ candidates identified by cross-matching the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS) to the Sloan Digital Sky Survey (SDSS): in the ~1900 deg^2 covered by the LAS in the UKIDSS Seventh Data Release (DR7) there are ~10^3 real astronomical point-sources with the measured colours of the target quasars, of which only ~10 are expected to be HZQs. Applying Bayesian model comparison to the sample reveals that most sources with HZQ-like colours have P_q <~ 0.1 and can be confidently rejected without the need for any further observations. In the case of the UKIDSS DR7 LAS, there were just 88 candidates with P_q >= 0.1; these object were prioritized for reobservation by ranking according to P_q (and their likely redshift, which was also inferred from the photometric data). Most candidates were rejected after one or two (moderate depth) photometric measurements by recalculating P_q using the new data. That left seven confirmed HZQs, three of which were previously identified in the SDSS and four of which were new UKIDSS discoveries. The high efficiency of this Bayesian selection method suggests that it could usefully be extended to other HZQ surveys (e.g. searches by Pan-STARRS or VISTA) as well as to other searches for rare objects.
[109]  [pdf] - 1051259
The Mass-Dependent Clustering History of K-selected Galaxies at z < 4 in the SXDS/UDS Field
Comments: 13 pages, 13 figures, 2 tables
Submitted: 2011-01-07
We investigate mass-dependent galaxy evolution based on a large sample of (more than 50,000) K-band selected galaxies in a multi-wavelength catalog of the Subaru/XMM-Newton Deep Survey (SXDS) and the UKIRT Infrared Deep Sky Survey (UKIDSS)/Ultra Deep Survey (UDS). We employ the optical to near-infrared photometry to determine photometric redshifts of these galaxies. Then, we estimate the stellar mass of our sample galaxies using a standard fitting procedure. From the sample galaxies, we obtain the stellar mass function of galaxies and the cosmic stellar mass density up to z<4. Our results are consistent with previous studies and we find a considerable number of low-mass galaxies (M<10^{10.5}) at the redshift range 3<z<4. The stellar-mass dependent correlation functions of our sample galaxies show clear evolution and they connect to that in the local universe consistently. Also, the massive galaxies show strong clustering throughout our studied redshift range. The correlation length of massive galaxies rapidly decreases from z=4 to 2. We also find some high mass density regions of massive galaxies at 1.4<z<2.5 in our sample, which may be candidate progenitors of the present-day clusters of galaxies. At this redshift range, massive star-forming galaxies are the dominant population making up the structures and the passively evolving galaxies show stronger clustering and they may have formed earlier than those star-forming galaxies.
[110]  [pdf] - 1042882
Galaxy Environments in the UKIDSS Ultra Deep Survey (UDS)
Comments: Accepted for publication in MNRAS. The paper contains 9 pages, 11 figures and 1 table
Submitted: 2011-01-04
We present a study of galaxy environments to z~2, based on a sample of over 33,000 K-band selected galaxies detected in the UKIDSS Ultra Deep Survey (UDS). The combination of infrared depth and area in the UDS allows us to extend previous studies of galaxy environment to z>1 without the strong biases associated with optical galaxy selection. We study the environments of galaxies divided by rest frame (U-B) colours, in addition to `passive' and `star-forming' subsets based on template fitting. We find that galaxy colour is strongly correlated with galaxy overdensity on small scales (<1Mpc diameter), with red/passive galaxies residing in significantly denser environments than blue/star-forming galaxies to z~1.5. On smaller scales (<0.5Mpc diameter) we also find a relationship between galaxy luminosity and environment, with the most luminous blue galaxies at z~1 inhabiting environments comparable to red, passive systems at the same redshift. Monte Carlo simulations demonstrate that these conclusions are robust to the uncertainties introduced by photometric redshift errors.
[111]  [pdf] - 254499
The Detection of a Population of Submillimeter-Bright, Strongly-Lensed Galaxies
Comments: This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science VOL 330, (Nov 5th, 2010), doi:10.1126/science.1193420
Submitted: 2010-11-04
Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.
[112]  [pdf] - 1033743
Statistics of 207 Lya Emitters at a Redshift Near 7: Constraints on Reionization and Galaxy Formation Models
Comments: 28 pages, 23 figures. Accepted for publication in ApJ
Submitted: 2010-07-17, last modified: 2010-09-03
We present Lya luminosity function (LF), clustering measurements, and Lya line profiles based on the largest sample, to date, of 207 Lya emitters (LAEs) at z=6.6 on the 1-deg^2 sky of Subaru/XMM-Newton Deep Survey (SXDS) field. Our z=6.6 Lya LF including cosmic variance estimates yields the best-fit Schechter parameters of phi*=8.5 +3.0/-2.2 x10^(-4) Mpc^(-3) and L*(Lya)=4.4 +/-0.6 x10^42 erg s^(-1) with a fixed alpha=-1.5, and indicates a decrease from z=5.7 at the >~90% confidence level. However, this decrease is not large, only =~30% in Lya luminosity, which is too small to be identified in the previous studies. A clustering signal of z=6.6 LAEs is detected for the first time. We obtain the correlation length of r_0=2-5 h^(-1) Mpc and bias of b=3-6, and find no significant boost of clustering amplitude by reionization at z=6.6. The average hosting dark halo mass inferred from clustering is 10^10-10^11 Mo, and duty cycle of LAE population is roughly ~1% albeit with large uncertainties. The average of our high-quality Keck/DEIMOS spectra shows an FWHM velocity width of 251 +/-16 km s^(-1). We find no large evolution of Lya line profile from z=5.7 to 6.6, and no anti-correlation between Lya luminosity and line width at z=6.6. The combination of various reionization models and our observational results about the LF, clustering, and line profile indicates that there would exist a small decrease of IGM's Lya transmission owing to reionization, but that the hydrogen IGM is not highly neutral at z=6.6. Our neutral-hydrogen fraction constraint implies that the major reionization process took place at z>~7.
[113]  [pdf] - 1032709
Measuring Transverse Motions for Nearby Galaxy Clusters
Comments: Accepted for publication in ApJ Letters
Submitted: 2010-05-24
Measuring the full three-dimensional motions of extra-galactic objects in the Universe presents a seemingly insurmountable challenge. In this paper we investigate the application of a technique to measure tangential motion that has previously only been applied nearby within the Local Group of galaxies, to clusters of galaxies far beyond its borders. We show that mapping the mean line-of-sight motion throughout a galaxy cluster could in principle be used to detect the "perspective rotation" induced by the projection of the cluster's tangential motion into the line-of-sight. The signal will be most prominent for clusters of the largest angular extent, most symmetric intrinsic velocity distribution and surveyed with the largest number of pointings possible. We investigate the feasibility of detecting this signal using three different approaches: measuring line-of-sight motions of individual cluster members; taking spectra of intracluster gas; and mapping distortions of the Cosmic Microwave Background radiation. We conclude that future spectroscopic surveys of 1000's of members of nearby galaxy clusters hold the most promise of measuring cluster tangential motions using this technique.
[114]  [pdf] - 317716
The evolution of galaxy clustering since z = 3 using the UKIDSS Ultra Deep Survey: the divergence of passive and star-forming galaxies
Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 11 pages containing 8 Figures and 1 table. The source contains additional results tables
Submitted: 2010-05-07
We use the UKIDSS Ultra-Deep Survey to trace the evolution of galaxy clustering to z = 3. Using photometric redshifts derived from data covering the wavelength range 0.3 - 4.5 um we examine this clustering as a function of absolute K-band luminosity, colour and star-formation rate. Comparing the deprojected clustering amplitudes, we find that red galaxies are more strongly clustered than blue galaxies out to at least z = 1.5, irrespective of rest-frame K-band luminosity. We then construct passive and star-forming samples based on stellar age, colour and star-formation histories calculated from the best fitting templates. The clustering strength of star-forming galaxies declines steadily from r_0 ~ 7 h^-1 Mpc at z ~ 2 to r_0 ~ 3 h^-1 Mpc at z ~ 0, while passive galaxies have clustering strengths up to a factor of two higher. Within the passive and star-forming subsamples, however, we find very little dependence of galaxy clustering on K-band luminosity. Galaxy `passivity' appears to be the strongest indicator of clustering strength. We compare these clustering measurements with those predicted for dark matter halos and conclude that passive galaxies typically reside in halos of mass M > 10^13 M_sun while luminous star-forming galaxies occupy halos an order of magnitude less massive over the range 0.5 < z < 1.5. The decline in the clustering strength of star-forming galaxies with decreasing redshift indicates a decline in the hosting halo mass for galaxies of a given luminosity. We find evidence for convergence of clustering in star-forming and passive galaxies around z ~ 2, which is consistent with this being the epoch at which the red sequence of galaxies becomes distinct.
[115]  [pdf] - 856214
Mid-infrared spectroscopy of candidate AGN-dominated submillimeter galaxies
Comments: Accepted for publication in ApJ. 17 pages, 10 figures
Submitted: 2010-03-02
Spitzer spectroscopy has revealed that ~80% of submm galaxies (SMGs) are starburst (SB) dominated in the mid-infrared. Here we focus on the remaining ~20% that show signs of harboring powerful active galactic nuclei (AGN). We have obtained Spitzer-IRS spectroscopy of a sample of eight SMGs which are candidates for harboring powerful AGN on the basis of IRAC color-selection (S8/S4.5>2; i.e. likely power-law mid-infrared SEDs). SMGs with an AGN dominating (>50%) their mid-infrared emission could represent `missing link' sources in an evolutionary sequence involving a major merger. First of all, we detect PAH features in all of the SMGs, indicating redshifts from 2.5-3.4, demonstrating the power of the mid-infrared to determine redshifts for these optically faint dusty galaxies. Secondly, we see signs of both star-formation (from the PAH features) and AGN activity (from continuum emission) in our sample: 62% of the sample are AGN-dominated in the mid-infrared with a median AGN content of 56%, compared with <30% on average for typical SMGs, revealing that our IRAC color selection has successfully singled out sources with proportionately more AGN emission than typical SB-dominated SMGs. However, we find that only about 10% of these AGN dominate the bolometric emission of the SMG when the results are extrapolated to longer infrared wavelengths, implying that AGN are not a significant power source to the SMG population overall, even when there is evidence in the mid-infrared for substantial AGN activity. When existing samples of mid-infrared AGN-dominated SMGs are considered, we find that S8/S4.5>1.65 works well at selecting mid-infrared energetically dominant AGN in SMGs, implying a duty cycle of ~15% if all SMGs go through a subsequent mid-infrared AGN-dominated phase in the proposed evolutionary sequence.
[116]  [pdf] - 1019147
When galaxies collide: understanding the broad absorption-line radio galaxy 4C +72.26
Comments: 11 pages, 10 figures, 1 table. Accepted for publication in MNRAS
Submitted: 2010-01-12
We present a range of new observations of the "broad absorption line radio galaxy" 4C +72.26 (z = 3.5), including sensitive rest-frame ultraviolet integral field spectroscopy using the Gemini/GMOS-N instrument and Subaru/CISCO K-band imaging and spectroscopy. We show that 4C +72.26 is a system of two vigorously star-forming galaxies superimposed along the line of sight separated by ~1300 +/- 200 km/s in velocity, with each demonstrating spectroscopically resolved absorption lines. The most active star-forming galaxy also hosts the accreting supermassive black hole which powers the extended radio source. We conclude that the star formation is unlikely to have been induced by a shock caused by the passage of the radio jet, and instead propose that a collision is a more probable trigger for the star formation. Despite the massive starburst, the UV-mid-infrared spectral energy distribution suggests that the pre-existing stellar population comprises ~10^12 Msolar of stellar mass, with the current burst only contributing a further ~2%, suggesting that 4C +72.26 has already assembled most of its final stellar mass.
[117]  [pdf] - 1018568
X-ray groups and clusters of galaxies in the Subaru-XMM Deep Field
Comments: 15 pages, MNRAS subm
Submitted: 2009-11-30
We present the results of a search for galaxy clusters in Subaru-XMM Deep Field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 2x10^{-15} ergs cm^{-2} s^{-1} over one of the widest XMM-Newton contiguous raster surveys, covering an area of 1.3 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. The red sequence technique allows us to identify 57 cluster candidates. We report on the progress with the cluster spectroscopic follow-up and derive their properties based on the X-ray luminosity and cluster scaling relations. In addition, 3 sources are identified as X-ray counterparts of radio lobes, and in 3 further sources, X-ray counterpart of radio lobes provides a significant fraction of the total flux of the source. In the area covered by NIR data, our identification success rate achieves 86%. We detect a number of radio galaxies within our groups and for a luminosity-limited sample of radio galaxies we compute halo occupation statistics using a marked cluster mass function. We compare the cluster detection statistics in the SXDF with the predictions of concordance cosmology and current knowledge of the X-ray cluster properties, concluding that a reduction of concordance sigma_8 value by 5% is required in order to match the prediction of the model and the data. This conclusion still needs verification through the completion of cluster follow-up.
[118]  [pdf] - 1001009
Discovery of a redshift 6.13 quasar in the UKIRT Infrared Deep Sky Survey
Comments: A&A in press; 8 pages, 7 figures; author list corrected
Submitted: 2008-10-22, last modified: 2009-09-17
Optical and near-infrared (NIR) spectra are presented for ULAS J131911.29+095951.4 (hereafter ULAS J1319+0950), a new redshift z = 6.127 +/- 0.004 quasar discovered in the Third Data Release (DR3) of the UKIRT Infrared Deep Sky Survey (UKIDSS). The source has Y_Vega = 19.10 +/- 0.03, corresponding to M_1450,AB = -27.12, which is comparable to the absolute magnitudes of the z ~= 6 quasars discovered in the SDSS. ULAS J1319+0950 was, in fact, registered by SDSS as a faint source with z_AB = 20.1 +/- 0.1, just below the signal--to--noise ratio limit of the SDSS high-redshift quasar survey. The faint z-band magnitude is a consequence of the weak Ly alpha/N V emission line, which has a rest-frame equivalent width of only ~20 A and provides only a small boost to the z-band flux. Nevertheless, there is no evidence from this UKIDSS-based search for a significant population of high-redshift quasars with weak emission lines. The Ly alpha optical depth to ULAS J1319+0950 is consistent with that measured towards similarly distant SDSS quasars, implying that results from optical- and NIR-selected quasars may be combined. Also presented is a new NIR-spectrum of the previously discovered UKIDSS quasar ULAS J020332.38+001229.2, which reveals the object to be a broad absorption line quasar. The new spectrum shows that the emission line previously identified as Ly alpha is actually N V, leading to a revised redshift of z = 5.72, rather than z = 5.86 as previously estimated.
[119]  [pdf] - 1002836
The discovery of a typical radio galaxy at z = 4.88
Comments: 5 pages, 3 figures, accepted for publication in MNRAS Letters
Submitted: 2009-07-09
In this letter we report the discovery of a z=4.88 radio galaxy discovered with a new technique which does not rely on pre-selection of a sample based on radio properties such as steep-spectral index or small angular size. This radio galaxy was discovered in the Elais-N2 field and has a spectral index of alpha = 0.75, i.e. not ultra-steep spectrum. It also has a luminosity consistent with being drawn from the break of the radio luminosity function and can therefore be considered as a typical radio galaxy. Using the Spitzer-SWIRE data over this field we find that the host galaxy is consistent with being similarly massive to the lower redshift powerful radio galaxies (~1-3L*). We note however, that at z=4.88 the H-alpha line is redshifted into the IRAC 3.6micron filter and some of the flux in this band may be due to this rather than stellar continuum emission. The discovery of such a distant radio source from our initial spectroscopic observations demonstrate the promise of our survey for finding the most distant radio sources.
[120]  [pdf] - 900726
Optical and near infra-red colours as a discriminant of the age and metallicity of stellar populations
Comments: 15 pages, 5 figures, accepted for publication in Monthly Notices of the Royal Astronomomical Society
Submitted: 2009-05-06
We present a comprehensive analysis of the ability of current stellar population models to reproduce the optical (ugriz) and near infra-red (JHK) colours of a small sample of well-studied nearby elliptical and S0 galaxies. We find broad agreement between the ages and metallicities derived using different population models, although different models show different systematic deviations from the measured broad-band fluxes. Although it is possible to constrain Simple Stellar Population models to a well defined area in age-metallicity space, there is a clear degeneracy between these parameters even with such a full range of precise colours. The precision to which age and metallicity can be determined independently, using only broad band photometry with realistic errors, is Delta{[Fe/H]} ~ 0.18 and Delta{log(Age)} ~ 0.25. To constrain the populations and therefore the star formation history further it will be necessary to combine broad-band optical-IR photometry with either spectral line indices, or else photometry at wavelengths outside of this range.
[121]  [pdf] - 15604
The star-formation history of K-selected galaxies
Comments: Published in MNRAS, 20 pages, colour
Submitted: 2008-08-22, last modified: 2009-03-16
We have studied the uJy radio properties of K-selected galaxies detected in the Ultra-Deep Survey portion of UKIDSS using 610- and 1,400-MHz images from the VLA and GMRT. These deep radio mosaics, combined with the largest and deepest K image currently available, allow high-S/N detections of many K-selected sub-populations, including sBzK and pBzK star-forming and passive galaxies. We find a strong correlation between the radio and K fluxes and a linear relationship between SFR and K luminosity. We find no evidence, from either radio spectral indices or a comparison with submm-derived SFRs, that the full sample is strongly contaminated by AGN. The sBzK and pBzK galaxies have similar levels of radio flux, SFR and specific SFR (SSFR) at z < 1.4, suggesting there is strong contamination of the pBzK sample by star-forming galaxies. At z > 1.4, pBzK galaxies become difficult to detect in the radio stack, though the implied SFRs are still much higher than expected for passively evolving galaxies. Their radio emission may come from low-luminosity AGN. EROs straddle the passive and star-forming regions of the BzK diagram and also straddle the two groups in terms of their radio properties. K-bright ERO samples are dominated by passive galaxies and faint ERO samples contain more star-forming galaxies. The star-formation history (SFH) from stacking all K sources in the UDS agrees well with that derived for other wavebands and other radio surveys, at least out to z ~ 2. The radio-derived SFH then appears to fall more steeply than that measured at other wavelengths. The SSFR for K-selected sources rises strongly with redshift at all stellar masses, and shows a weak dependence on stellar mass. High- and low-mass galaxies show a similar decline in SSFR since z ~ 2 (abridged).
[122]  [pdf] - 20419
Evidence of a link between the evolution of clusters and their AGN fraction
Comments: 19 pages, 19 figures, accepted for publication in MNRAS. Note that the figures in this version are of low resolution. A high-resolution version can be downloaded from
Submitted: 2009-01-16
We discuss the optical properties, X-ray detections, and Active Galactic Nucleus (AGN) populations of four clusters at z ~ 1 in the Subaru-XMM Deep Field (SXDF). The velocity distribution and plausible extended X-ray detections are examined, as well as the number of X-ray point sources and radio sources associated with the clusters. We find that the two clusters that appear virialised and have an extended X-ray detection contain few, if any, AGN, whereas the two pre-virialised clusters have a large AGN population. This constitutes evidence that the AGN fraction in clusters is linked to the clusters' evolutionary stage. The number of X-ray AGN in the pre-virialised clusters is consistent with an overdensity of factor ~ 200; the radio AGN appear to be clustered with a factor of three to six higher. The median K-band luminosities of L_K = 1.7 +/- 0.7 L* for the X-ray sources and L_K = 2.3 +/- 0.1 L* for the radio sources support the theory that these AGN are triggered by galaxy interaction and merging events in sub-groups with low internal velocity distributions, which make up the cluster environment in a pre-virialisation evolutionary stage.
[123]  [pdf] - 18428
An 80-kpc Lyman-alpha halo around a high redshift type-2 QSO
Comments: 8 pages, 7 figures. Accepted by MNRAS
Submitted: 2008-11-11
We announce the discovery of an extended emission line region associated with a high redshift type-2 QSO. The halo, which was discovered in our new wide-field narrow-band survey, resides at z = 2.85 in the Spitzer First Look Survey region and is extended over ~80 kpc. Deep VLBI observations imply that approximately 50 per cent of the radio emission is extended on scales > 200pc. The inferred AGN luminosity is sufficient to ionize the extended halo, and the optical emission is consistent with being triggered coevally with the radio source. The Lyman-alpha halo is as luminous as those found around high redshift radio galaxies, however the active nucleus is several orders of magnitude less luminous at radio wavelengths than those FRIIs more commonly associated with extended emission line regions. AMS05 appears to be a high-redshift analogue to the radio-quiet quasar E1821+643 which is core dominated but which also exhibits extended FRI-like structure and contains an optically powerful AGN. We also find evidence for more quiescent kinematics in the Lyman-alpha emission line in the outer regions of the halo, reminiscent of the haloes around the more powerful FRIIs. The optical to mid-infrared SED is well described by a combination of an obscured QSO (L_bol ~3.4 +/- 0.2 x 10^13 Lsolar) and a 1.4 Gyr old simple stellar population with mass ~3.9 +/- 0.3 x 10^11 Msolar.
[124]  [pdf] - 17682
The UKIRT Infrared Deep Sky Survey and the search for the most distant quasars
Comments: to appear in Proceedings of the XXth Moriond Meeting: Cosmology, J. Tran Thanh Van et al., eds.; 4 pages, 2 figures
Submitted: 2008-10-21
The UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS) has the necessary combination of filters (Y, J, H and K), depth (Y <~ 20.2) and area coverage (~4000 deg^2) to detect several redshift z >~ 6.4 quasars. The Third Data Release (DR3) included ~1000 deg^2 of LAS observations which have so far yielded two previously known z ~= 6 quasars and two new discoveries: ULAS J0203+0012, at z = 5.72; and ULAS J1319+0950, at z = 6.13.
[125]  [pdf] - 16371
The clustering and abundance of star-forming and passive galaxies at $z \sim 2$
Comments: 8 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2008-09-16
We use the UKIDSS Ultra-deep survey (UDS), currently the deepest panoramic near infra-red survey, together with deep Subaru optical imaging to measure the clustering, number counts and luminosity function of galaxies at $z\sim 2$ selected using the BzK selection technique. We find that both star-forming (sBzK) and passive (pBzK) galaxies, to a magnitude limit of $K_{AB} < 23$, are strongly clustered. The passive galaxies are the most strongly clustered population, with scale lengths of $r_0 = 15.0^{+1.9}_{-2.2}$h$^{-1}$Mpc compared with $r_0 = 6.75^{+0.34}_{-0.37}$h$^{-1}$Mpc for star-forming galaxies. The direct implication is that passive galaxies inhabit the most massive dark-matter halos, and are thus identified as the progenitors of the most massive galaxies at the present day. In addition, the pBzKs exhibit a sharp flattening and potential turn-over in their number counts, in agreement with other recent studies. This plateau cannot be explained by the effects of incompleteness. We conclude that only very massive galaxies are undergoing passive evolution at this early epoch, consistent with the downsizing scenario for galaxy evolution. Assuming a purely passive evolution for the pBzKs from their median redshift to the present day, their luminosity function suggests that only $\sim 2.5 %$ of present day massive ellipticals had a pBzK as a main progenitor.
[126]  [pdf] - 13642
The Subaru/XMM-Newton Deep Survey (SXDS): III. X-Ray Data
Comments: 17 pages, 32 figures, accepted for publication in ApJS
Submitted: 2008-06-17
We present the X-ray source catalog in the Subaru/XMM-Newton deep survey. A continuous area of 1.14 deg^2 centered at R.A. = 02h18m and Dec. = -05d is mapped by seven pointings with XMM-Newton covering the 0.2-10 keV band. From the combined images of the EPIC pn and MOS cameras, we detect 866, 1114, 645, and 136 sources with sensitivity limits of 6x10^{-16}, 8x10^{-16}, 3x10^{-15}, and 5x10^{-15} erg cm^{-2} s^{-1} in the 0.5-2, 0.5-4.5, 2-10, and 4.5-10 keV bands, respectively, with detection likelihood >= 7 (corresponding to a confidence level of 99.91%). The catalog consists of 1245 sources in total including 32 extended-source candidates. The averaged log N-log S relations are in good agreement with previous results, bridging the flux range between Chandra deep surveys and brighter surveys. The log N-log S relations show significant spatial variation among pointings on a scale of 0.2 deg^2. Analyzing the auto correlation function, we detect significant clustering signals from the 0.5-2 keV band sample, which can be fit with a power law form (\theta/\theta_c)^{-0.8} with a correlation length of \theta_c=5.9^{+1.0}_{-0.9} arcsec when the integral constraint term is included. In the 2-10 keV band, however, the clustering is not significant with a 90% upper limit of \theta_c < 1.5 arcsec.
[127]  [pdf] - 3304
The Subaru/XMM-Newton Deep Survey (SXDS). IV. Evolution of Lya Emitters from z=3.1 to 5.7 in the 1 deg^2 Field: Luminosity Functions and AGN
Comments: 75 pages, 27 figures; ApJS in press. High resolution version at
Submitted: 2007-07-20, last modified: 2008-06-02
We present luminosity functions (LFs) and various properties of Lya emitters (LAEs) at z=3.1, 3.7, and 5.7, in a 1 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field. We obtain a photometric sample of 858 LAE candidates based on deep Subaru/Suprime-Cam imaging data, and a spectroscopic sample of 84 confirmed LAEs from Subaru/FOCAS and VLT/VIMOS spectroscopy in a survey volume of ~10^6 Mpc^3 with a limiting Lya luminosity of ~3x10^42 erg/s. We derive the LFs of Lya and UV-continuum (~1500 \AA) for each redshift, taking into account the statistical error and the field-to-field variation. We find that the apparent Lya LF shows no significant evolution between z=3.1 and 5.7 within factors of 1.8 and 2.7 in L* and phi*, respectively. On the other hand, the UV LF of LAEs increases from z=3.1 to 5.7, indicating that galaxies with Lya emission are more common at earlier epochs. We identify six LAEs with AGN activities from our spectra combined with VLA, Spitzer, and XMM-Newton data. Among the photometrically selected LAEs at z=3.1 and 3.7, only ~1 % show AGN activities, while the brightest LAEs with logL(Lya) >~ 43.4-43.6 erg/s appear to always host AGNs. Our LAEs are bluer in UV-continuum color than dropout galaxies, suggesting lower extinction and/or younger stellar populations. Our stacking analyses provide upper limits to the radio luminosity and the f(HeII)/f(Lya) line fraction, and constrain the hidden star formation (+low-luminosity AGN) and the primordial population in LAEs.
[128]  [pdf] - 1000710
A new measurement of the evolving near-infrared galaxy luminosity function out to z~4: a continuing challenge to theoretical models of galaxy formation
Comments: 11 pages, 6 figures, submitted to MNRAS
Submitted: 2008-04-22
We present the most accurate measurement to date of cosmological evolution of the near-infrared galaxy luminosity function, from the local Universe out to z~4. The analysis is based on a large and highly complete sample of galaxies selected from the first data release of the UKIDSS Ultra Deep Survey. Exploiting a master catalogue of K- and z-band selected galaxies over an area of 0.7 square degrees, we analyse a sample of ~50,000 galaxies, all with reliable photometry in 16-bands from the far-ultraviolet to the mid-infrared. The unique combination of large area and depth provided by the Ultra Deep Survey allows us to trace the evolution of the K-band luminosity function with unprecedented accuracy. In particular, via a maximum likelihood analysis we obtain a simple parameterization for the luminosity function and its cosmological evolution, including both luminosity and density evolution, which provides an excellent description of the data from z =0 up to z~4. We find differential evolution for galaxies dependent on galaxy luminosity, revealing once again the ``down-sizing behaviour'' of galaxy formation. Finally, we compare our results with the predictions of the latest theoretical models of galaxy formation, based both on semi-analytical prescriptions, and on full hydrodynamical simulations.
[129]  [pdf] - 11073
The Subaru/XMM-Newton Deep Field - II. The 37 brightest radio sources
Comments: Accepted for publication by MNRAS; 39 pages, 7 figures, 4 tables
Submitted: 2008-03-19
We study the 37 brightest radio sources in the Subaru/XMM-Newton Deep Field (SXDF). We have spectroscopic redshifts for 24 of 37 objects and photometric redshifts for the remainder, yielding a median redshift z_med for the whole sample of z_med ~= 1.1 and a median radio luminosity close to the `FRI/FRII' luminosity divide. Using mid-IR (Spitzer MIPS 24 um) data we expect to trace nuclear accretion activity, even if it is obscured at optical wavelengths, unless the obscuring column is extreme. Our results suggest that above the FRI/FRII radio luminosity break most of the radio sources are associated with objects that have excess mid-IR emission, only some of which are broad-line objects, although there is one clear low-accretion-rate object with an FRI radio structure. For extended steep-spectrum radio sources, the fraction of objects with mid-IR excess drops dramatically below the FRI/FRII luminosity break, although there exists at least one high-accretion-rate `radio-quiet' QSO. We have therefore shown that the strong link between radio luminosity (or radio structure) and accretion properties, well known at z ~ 0.1, persists to z ~ 1. Investigation of mid-IR and blue excesses shows that they are correlated as predicted by a model in which, when significant accretion exists, a torus of dust absorbs ~30% of the light, and the dust above and below the torus scatters >~1% of the light.
[130]  [pdf] - 10645
The SCUBA Half-Degree Extragalactic Survey (SHADES) - VIII. The Nature of Faint Submm Galaxies in SHADES, SWIRE and SXDF Surveys
Comments: MNRAS in press
Submitted: 2008-03-04
We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.
[131]  [pdf] - 10233
Exploring the infrared/radio correlation at high redshift
Comments: The paper contains 7 figures and 1 table. In press at MNRAS
Submitted: 2008-02-19
We have analysed the 24um properties of a radio-selected sample in the Subaru-XMM/Newton Deep Field in order to explore the behaviour of the FIR/radio relation at high redshifts. Statistically, the correlation is described by q24, the ratio between the observed flux densities at 24um and 1.4GHz, respectively. Using 24um data results in considerably more scatter in the correlation than previous work using data at 60-70um. Nevertheless, we do observe a steady correlation as a function of redshift, up to z~3.5, suggesting its validity back to primeval times. We find q24 = 0.30 +/- 0.56 for the observed and q24 = 0.71 +/- 0.47 for the k-corrected radio sample, based on sources with 300uJy < S(1.4GHz) < 3.2mJy and 24um detections. A suitable k-correction given by a M82-like mid-IR template suggests no extreme silicate absorption in the bulk of our radio sample. Using thresholds in q24 to identify radio-excess sources, we have been able to characterise the transition from radio-loud AGN to star-forming galaxies and radio-quiet AGN at faint (<1mJy) radio flux densities. Our results are in broad agreement with previous studies which show a dominant radio-loud AGN population at >1mJy. The rest-frame U-B colours of the expected radio-excess population have redder distribution than those that follow the correlation. This is therefore a promising way to select obscured Type-2 AGN, with a radio loud nature, missed by deep X-ray observations. Spectroscopic follow-up of these sources is required to fully test this method.
[132]  [pdf] - 9748
The SCUBA HAlf Degree Extragalactic Survey (SHADES) - VII. Optical/IR photometry and stellar masses of sub-millimeter galaxies
Comments: 27 pages, 23 figures. Refereed resubmission to MNRAS. 5 figures containing postage stamp images have been placed in separate files due to upload limit
Submitted: 2008-02-04
We present estimates of the photometric redshifts, stellar masses and star formation histories of sources in the SCUBA HAlf Degree Extragalactic Survey (SHADES). This paper describes the 60 SCUBA sources detected in the Lockman Hole covering an area of ~320 square arcmin. Using photometry spanning the B band to 8um, we find that the average SCUBA source forms a significant fraction of its stars in an early period of star formation and that most of the remainder forms in a shorter more intense burst around the redshift it is observed. This trend does not vary significantly with source redshift but the exact ratio of early to late mass is quite sensitive to the way extinction is treated in the modelling. However, the sources show a clear increase in stellar mass with redshift, consistent with downsizing. In terms of SED types, only two out of the 51 sources we have obtained photometric redshifts for are best fit by a quasar-like spectral energy distribution, with approximately 80 per cent of the sources being best fit with late-type spectra (Sc, Im and starburst). By including photometry at 850um, we conclude that the average SCUBA source is forming stars at a rate somewhere between 6 and 30 times the rate implied from the rest-frame optical in a dust obscured burst and that this burst creates 15-65 per cent of the total stellar mass. Using a simplistic calculation, we estimate from the average star formation history that between one in five and one in 15 bright (L_* +2 mag < L_optical < L_* -1 mag) galaxies in the field over the interval 0 < z < 3 will at some point in their lifetime experience a similar energetic dusty burst of star formation. Finally, we compute the evolution of the star formation rate density and find it peaks around z=2.
[133]  [pdf] - 9406
The Subaru/XMM-Newton Deep Survey (SXDS) - II. Optical Imaging and Photometric Catalogs
Comments: 55 pages, 22 figures, 9 tables, accepted by ApJS, a higher-resolution version is available at
Submitted: 2008-01-25
We present multi-waveband optical imaging data obtained from observations of the Subaru/XMM-Newton Deep Survey (SXDS). The survey field, centered at R.A.=02:18:00, decl.=-05:00:00, has been the focus of a wide range of multi-wavelength observing programs spanning from X-ray to radio wavelengths. A large part of the optical imaging observations are carried out with Suprime-Cam on Subaru Telescope at Mauna Kea in the course of Subaru Telescope Observatory Projects. This paper describes our optical observations, data reduction and analysis procedures employed, and the characteristics of the data products. A total area of 1.22 sqdeg is covered in five contiguous sub-fields, each of which corresponds to a single Suprime-Cam field of view (34'x27'), in five broad-band filters B, V, Rc, i', z' to the depths of B=28.4, V=27.8, Rc=27.7, i'=27.7 and z'=26.6 (AB, 3-sigma, 2-arcsec aperture). The data are reduced and compiled into five multi-waveband photometric catalogs, separately for each Suprime-Cam pointing. The i'-band catalogs contain about 900,000 objects, making the SXDS catalogs one of the largest multi-waveband catalogs in corresponding depth and area coverage. The SXDS catalogs can be used for an extensive range of astronomical applications such as the number density of the Galactic halo stars to the large scale structures at the distant universe. The number counts of galaxies are derived and compared with those of existing deep extragalactic surveys. The optical data, the source catalogs, and configuration files used to create the catalogs are publicly available via the SXDS web page (
[134]  [pdf] - 260111
Accretion indicators for the 37 brightest radio sources in the Subaru/XMM-Newton Deep Field
Comments: 4 pages, 2 figures, to appear in conference proceedings of 'A Century of Cosmology; Past, Present and Future', Venice 2007
Submitted: 2007-10-26
We study the 37 brightest radio sources in the Subaru/XMM-Newton Deep Field (SXDF). Using mid-IR (Spitzer MIPS 24 micron) data we expect to trace nuclear accretion activity, even if it is obscured at optical wavelengths, unless the obscuring column is extreme. Our results suggest that above the `FRI/FRII' radio luminosity break most of the radio sources are associated with objects that have excess mid-IR emission, only some of which are broad-line objects, although there is one clear low-accretion-rate FRI. The fraction of objects with mid-IR excess drops dramatically below the FRI/FRII break, although there exists at least one high-accretion-rate QSO. Investigation of mid-IR and blue excesses shows that they are correlated as predicted by a model in which a torus of dust absorbs ~30% of the light, and the dust above and below the torus scatters >~1% of the light.
[135]  [pdf] - 6095
On the evolution of clustering of 24um-selected galaxies
Comments: 13 pages, 8 figures, to appear in MNRAS
Submitted: 2007-10-16
This paper investigates the clustering properties of a complete sample of 1041 24um-selected sources brighter than F[24um]=400 uJy in the overlapping region between the SWIRE and UKIDSS UDS surveys. We have concentrated on the two (photometric) interval ranges z=[0.6-1.2] (low-z sample) and z>1.6 (high-z sample) as it is in these regions were we expect the mid-IR population to be dominated by intense dust-enshrouded activity such as star formation and black hole accretion. Investigations of the angular correlation function produce a correlation length are r0~15.9 Mpc for the high-z sample and r0~8.5 Mpc for the low-z one. Comparisons with physical models reveal that the high-z sources are exclusively associated with very massive (M>~10^{13} M_sun)haloes, comparable to those which locally host groups-to-clusters of galaxies, and are very common within such (rare) structures. Conversely, lower-z galaxies are found to reside in smaller halos (M_min~10^{12} M_sun) and to be very rare in such systems. While recent studies have determined a strong evolution of the 24um luminosity function between z~2 and z~0, they cannot provide information on the physical nature of such an evolution. Our clustering results instead indicate that this is due to the presence of different populations of objects inhabiting different structures, as active systems at z<~1.5 are found to be exclusively associated with low-mass galaxies, while very massive sources appear to have concluded their active phase before this epoch. Finally, we note that the small-scale clustering data seem to require steep profiles for the distribution of galaxies within their halos. This is suggestive of close encounters and/or mergers which could strongly favour both AGN and star-formation activity.
[136]  [pdf] - 4336
Spectroscopic follow-up of a cluster candidate at z = 1.45
Comments: 15 pages, 7 figures, 2 tables, accepted for publication in MNRAS The figures in this version are low resolution; for a high-resolution version see (replaced to repair broken link for full-resolution version)
Submitted: 2007-08-28, last modified: 2007-09-03
We have obtained deep optical spectroscopic data of the highest-redshift cluster candidate (z ~ 1.4, CVB13) selected by Van Breukelen et al. (2006) in a photometric optical/infrared catalogue of the Subaru XMM-Newton Deep Field. The data, which comprise 104 targeted galaxies, were taken with the DEep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck 2 telescope and yielded 31 secure redshifts in the range 1.25 < z < 1.54 within a 7' x 4' field centred on CVB13. Instead of one massive cluster at z = 1.4, we find evidence for three projected structures at z = 1.40, z = 1.45, and z = 1.48. The most statistically robust of these structures, at z = 1.454, has six spectroscopically confirmed galaxies. Its total mass is estimated at > 10^14 M_sun and it may therefore be termed a cluster. There is an X-ray source at the cluster position which is marginally spatially resolved but whose X-ray spectrum is too hard to be thermal cluster emission. Its origin could be the summed X-ray emission from active galaxies in, and projected onto, the cluster. Serendipitously we have discovered a cluster at z = 1.28 with a mass of > 10^14 M_sun at the same position on the sky, comprising six spectroscopically confirmed cluster galaxies and at least one additional radio source. The selection of CVB13 for the cluster catalogue was evidently aided by the superposition of other, presumably lower-mass, structures, whereas the single cluster at z = 1.28 contained too few galaxies to be isolated by the same algorithm. Given the complicated nature of such structures, caution must be employed when measuring the mass function of putative high-redshift clusters with photometric techniques alone.
[137]  [pdf] - 81495
The UKIRT Infrared Deep Sky Survey (UKIDSS)
Comments: Final published version including significant revisions. Twenty four pages, fourteen figures. Original version April 2006; final version published in MNRAS August 2007
Submitted: 2006-04-20, last modified: 2007-08-15
We describe the goals, design, implementation, and initial progress of the UKIRT Infrared Deep Sky Survey (UKIDSS), a seven year sky survey which began in May 2005, using the UKIRT Wide Field Camera. It is a portfolio of five survey components covering various combinations of the filter set ZYJHK and H_2. The Large Area Survey, the Galactic Clusters Survey, and the Galactic Plane Survey cover approximately 7000 square degrees to a depth of K~18; the Deep Extragalactic Survey covers 35 square degrees to K~21, and the Ultra Deep Survey covers 0.77 square degrees to K~23. Summed together UKIDSS is 12 times larger in effective volume than the 2MASS survey. The prime aim of UKIDSS is to provide a long term astronomical legacy database; the design is however driven by a series of specific goals -- for example to find the nearest and faintest sub-stellar objects; to discover Population II brown dwarfs, if they exist; to determine the substellar mass function; to break the z=7 quasar barrier; to determine the epoch of re-ionisation; to measure the growth of structure from z=3 to the present day; to determine the epoch of spheroid formation; and to map the Milky Way through the dust, to several kpc. The survey data are being uniformly processed, and released in stages through the WFCAM Science Archive (WSA : Before the formal survey began, UKIRT and the UKIDSS consortium collaborated in obtaining and analysing a series of small science verification (SV) projects to complete the commissioning of the camera. We show some results from these SV projects in order to demonstrate the likely power of the eventual complete survey. Finally, using the data from the First Data Release we assess how well UKIDSS is meeting its design targets so far.
[138]  [pdf] - 84879
The evolution of the near-IR galaxy Luminosity Function and colour bimodality up to z ~ 2 from the UKIDSS Ultra Deep Survey Early Data Release
Comments: 12 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2006-09-11, last modified: 2007-08-09
We present new results on the cosmological evolution of the near-infrared galaxy luminosity function, derived from the analysis of a new sample of \~22,000 K(AB) < 22.5 galaxies selected over an area of 0.6 square degrees from the Early Data Release of the UKIDSS Ultra Deep Survey (UDS). Our study has exploited the multi-wavelength coverage of the UDS field provided by the new UKIDSS WFCAM K and J-band imaging, the Subaru/XMM-Newton Deep Survey and the Spitzer-SWIRE Survey. The unique combination of large area and depth provided by this new survey minimises the complicating effect of cosmic variance and has allowed us, for the first time, to trace the evolution of the brightest sources out to z~2 with good statistical accuracy. In agreement with previous studies we find that the characteristic luminosity of the near-infrared luminosity function brightens by 1 magnitude between z=0 and z~2, while the total density decreases by a factor ~2. Using the rest-frame (U-B) colour to split the sample into red and blue galaxies, we confirm the classic luminosity-dependent colour bimodality at z<1. However, the strength of the colour bimodality is found to be a decreasing function of redshift, and seems to disappear by z>1.5. Due to the large size of our sample we are able to investigate the differing cosmological evolution of the red and blue galaxy populations. It is found that the space density of the brightest red galaxies (M_K < -23) stays approximately constant with redshift, and that these sources dominate the bright-end of the luminosity function at redshifts z<1. In contrast, the brightening of the characteristic luminosity and mild decrease in space density displayed by the blue galaxy population leads them to dominate the bright-end of the luminosity function at redshifts z>1.
[139]  [pdf] - 3784
Low-power radio galaxy environments in the Subaru/XMM-Newton Deep Field at z~0.5
Comments: 12 pages, 7 figures, accepted for publication in MNRAS. High resolution version available upon request
Submitted: 2007-08-07
We present multi-object spectroscopy of galaxies in the immediate (Mpc-scale) environments of four low-power (L_1.4 GHz < 10^25 W/Hz) radio galaxies at z~0.5, selected from the Subaru/XMM-Newton Deep Field. We use the spectra to calculate velocity dispersions and central redshifts of the groups the radio galaxies inhabit, and combined with XMM-Newton (0.3-10 keV) X-ray observations investigate the L_X--sigma_v and T_X--sigma_v scaling relationships. All the radio galaxies reside in moderately rich groups -- intermediate environments between poor groups and rich clusters, with remarkably similar X-ray properties. We concentrate our discussion on our best statistical example that we interpret as a low-power (FRI) source triggered within a sub-group, which in turn is interacting with a nearby group of galaxies, containing the bulk of the X-ray emission for the system -- a basic scenario which can be compared to more powerful radio sources at both high (z>4) and low (z<0.1) redshifts. This suggests that galaxy-galaxy interactions triggered by group mergers may play an important role in the life-cycle of radio galaxies at all epochs and luminosities.
[140]  [pdf] - 142677
The UKIRT Infrared Deep Sky Survey First Data Release
Comments: Replaced with accepted version. 12 figures. 6 figures are low resolution. Version of paper with full-resolution figures available at
Submitted: 2006-10-06, last modified: 2007-07-30
The First Data Release (DR1) of the UKIRT Infrared Deep Sky Survey (UKIDSS) took place on 2006 July 21. UKIDSS is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. DR1 is the first large release of survey-quality data from UKIDSS and includes 320 sq degs of multicolour data to (Vega) K=18, complete (depending on the survey) in three to five bands from the set ZYJHK, together with 4 sq degs of deep JK data to an average depth K=21. In addition the release includes a similar quantity of data with incomplete filter coverage. In JHK, in regions of low extinction, the photometric uniformity of the calibration is better than 0.02 mag. in each band. The accuracy of the calibration in ZY remains to be quantified, and the same is true of JHK in regions of high extinction. The median image FWHM across the dataset is 0.82 arcsec. We describe changes since the Early Data Release in the implementation, pipeline and calibration, quality control, and archive procedures. We provide maps of the areas surveyed, and summarise the contents of each of the five surveys in terms of filters, areas, and depths. DR1 marks completion of 7 per cent of the UKIDSS 7-year goals.
[141]  [pdf] - 1807
The discovery of a massive supercluster at z=0.9 in the UKIDSS DXS
Comments: In press at MNRAS, 11 pages, 5 figures
Submitted: 2007-06-01
We analyse the first publicly released deep field of the UKIDSS Deep eXtragalactic Survey (DXS) to identify candidate galaxy over-densities at z~1 across ~1 sq. degree in the ELAIS-N1 field. Using I-K, J-K and K-3.6um colours we identify and spectroscopically follow-up five candidate structures with Gemini/GMOS and confirm they are all true over-densities with between five and nineteen members each. Surprisingly, all five structures lie in a narrow redshift range at z=0.89+/-0.01, although they are spread across 30Mpc on the sky. We also find a more distant over-density at z=1.09 in one of the spectroscopic survey regions. These five over-dense regions lying in a narrow redshift range indicate the presence of a supercluster in this field and by comparing with mock cluster catalogs from N-body simulations we discuss the likely properties of this structure. Overall, we show that the properties of this supercluster are similar to the well-studied Shapley and Hercules superclusters at lower redshift.
[142]  [pdf] - 89557
The SCUBA HAlf Degree Extragalactic Survey (SHADES) - IV: Radio-mm-FIR photometric redshifts
Comments: MNRAS, accepted (22 pages, 9 figures, 5 tables). Revised version, minor modifications and 2 extra figures
Submitted: 2007-02-19, last modified: 2007-05-31
We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.
[143]  [pdf] - 89598
The SCUBA HAlf Degree Extragalactic Survey (SHADES) - III. Identification of radio and mid-infrared counterparts to submillimetre galaxies
Comments: Accepted for MNRAS. 30 pages; 8 b&w figures; 1 colour figure; Appendix with 10 pages of b&w postage-stamp images. Full-resolution versions available at (PDF), (PS-A4) or (PS-Letter)
Submitted: 2007-02-20, last modified: 2007-05-30
Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.
[144]  [pdf] - 87125
Evidence for a large fraction of Compton-thick quasars at high redshift
Comments: 6 pages, 2 colour figures. Accepted by MNRAS. Full resolution version and supplementary figures can be found at:
Submitted: 2006-11-23, last modified: 2007-05-14
Using mid-infrared and radio selection criteria, we pre-select a sample of candidate high-redshift type-2 quasars in the Subaru XMM-Newton Deep Field (SXDF). To filter out starburst contaminants, we use a bayesian method to fit the spectral energy distributions (SEDs) between 24-microns and B-band, obtain photometric redshifts, and identify the best candidates for high-z type-2 quasars. This leaves us with 12 z_phot >= 1.7 type-2 quasar candidates in an area ~0.8 deg^2, of which only two have secure X-ray detections. The two detected sources have estimated column densities N_H~2 & 3x10^27 m^-2, i.e. heavily obscured but Compton-thin quasars. Given the large bolometric luminosities and redshifts of the undetected objects, the lack of X-ray detections suggests extreme absorbing columns N_H >= 10^28 m^-2 are typical. We have found evidence for a population of ``Compton-thick'' high-redshift type-2 quasars, at least comparable to, and probably larger than the type-1 quasar population, although spectroscopic confirmation of their AGN nature is important.
[145]  [pdf] - 471
The colour selection of distant galaxies in the UKIDSS Ultra-Deep Survey Early Data Release
Comments: Accepted for publication in MNRAS Letters. 5 pages, 5 figures, 1 table. Higher-resolution figures available from the authors on request
Submitted: 2007-04-17
We investigate colour selection techniques for high redshift galaxies in the UKIDSS Ultra Deep Survey Early Data Release (UDS EDR). Combined with very deep Subaru optical photometry, the depth (K_AB = 22.5) and area (0.62 deg^2) of the UDS EDR allows us to investigate optical/near-IR selection using a large sample of over 30,000 objects. By using the B-z, z-K colour-colour diagram (the BzK technique) we identify over 7500 candidate galaxies at z > 1.4, which can be further separated into passive and starforming systems (pBzK and sBzK respectively). Our unique sample allows us to identify a new feature not previously seen in BzK diagrams, consistent with the passively evolving track of early type galaxies at z < 1.4. We also compare the BzK technique with the R-K colour selection of Extremely Red Objects (EROs) and the J-K selection of Distant Red Galaxies (DRGs), and quantify the overlap between these populations. We find that the majority of DRGs, at these relatively bright magnitudes are also EROs. Since previous studies have found that DRGs at these magnitudes have redshifts of z ~ 1 we determine that these DRG/ERO galaxies have SEDs consistent with being dusty star-forming galaxies or AGN at z < 2. Finally we observe a flattening in the number counts of pBzK galaxies, similar to other studies, which may indicate that we are sampling the luminosity function of passive z > 1 galaxies over a narrow redshift range.
[146]  [pdf] - 82826
Number counts and clustering properties of bright Distant Red Galaxies in the UKIDSS Ultra Deep Survey Early Data Release
Comments: 5 pages, 4 figures, revised version accepted to MNRAS. Higher-resolution figures available from the authors on request
Submitted: 2006-06-15, last modified: 2006-12-19
We describe the number counts and spatial distribution of 239 Distant Red Galaxies (DRGs), selected from the Early Data Release of the UKIDSS Ultra Deep Survey. The DRGs are identified by their very red infrared colours with (J-K)AB>1.3, selected over 0.62 sq degree to a 90% completeness limit of KAB~20.7. This is the first time a large sample of bright DRGs has been studied within a contiguous area, and we provide the first measurements of their number counts and clustering. The population shows strong angular clustering, intermediate between those of K-selected field galaxies and optical/infrared-selected Extremely Red Galaxies. Adopting the redshift distributions determined from other recent studies, we infer a high correlation length of r0~11 h-1 Mpc. Such strong clustering could imply that our galaxies are hosted by very massive dark matter halos, consistent with the progenitors of present-day L>L* elliptical galaxies.
[147]  [pdf] - 85382
The 6C** sample of steep-spectrum radio sources: I - Radio data, near-infrared imaging and optical spectroscopy
Comments: 22 pages, 17 figures in .png format. Accepted for publication in MNRAS. Full size/resolution figures can be found at
Submitted: 2006-09-28
We present basic observational data on the 6C** sample. This is a new sample of radio sources drawn from the 151 MHz 6C survey, which was filtered with radio criteria chosen to optimize the chances of finding radio galaxies at z > 4. The filtering criteria are a steep-spectral index and a small angular size. The final sample consists of 68 sources from a region of sky covering 0.421 sr. We present VLA radio maps, and the results of K-band imaging and optical spectroscopy. Near-infrared counterparts are identified for 66 of the 68 sources, down to a 3-sigma limiting magnitude of K ~ 22 mag in a 3-arcsec aperture. Eight of these identifications are spatially compact, implying an unresolved nuclear source. The K-magnitude distribution peaks at a median K=18.7 mag, and is found to be statistically indistinguishable from that of the similarly selected 6C* sample, implying that the redshift distribution could extend to z > 4. Redshifts determined from spectroscopy are available for 22 (32 per cent) of the sources, over the range of 0.2 < z < 3.3 . We measure 15 of these, whereas the other 7 were previously known. Six sources are at z > 2.5. Four sources show broad emission lines in their spectra and are classified as quasars. Three of these show also an unresolved K-band identification. Eleven sources fail to show any distinctive emission and/or absorption features in their spectra. We suggest that these could be (i) in the so-called `redshift desert' region of 1.2 < z < 1.8, or (ii) at a greater redshift, but feature weak emission line spectra.
[148]  [pdf] - 85311
Probing accretion activity in radio sources using 24 micron Spitzer data
Comments: 2 pages, 2 figures, to be published by Crete University Press (conference proceedings of "Galaxy Evolution with Spitzer and Herschel", Crete, May 2006)
Submitted: 2006-09-26
We study the 36 brightest radio sources in the Subaru/XMM-Newton Deep Field SXDF). Using MIPS 24 micron data from Spitzer we expect to trace accretion activity, even if it is hidden at optical wavelengths, unless the obscuring column is extreme. Our results suggest that in the decade or so below the break in the radio luminosity function that at least half, and potentially nearly all, radio sources are associated with accreting quasar-like objects. This is not true at lower radio luminosities where the quasar-like fraction approaches zero once compact sources are excluded.
[149]  [pdf] - 85121
Radio imaging of the Subaru/XMM-Newton Deep Field - I. The 100-microJy catalogue, optical identifications, and the nature of the faint radio source population
Comments: 18 pages, MNRAS, in press, full catalogue and hi-res figures at
Submitted: 2006-09-19
We describe deep radio imaging at 1.4-GHz of the 1.3 square degree Subaru/XMM-Newton Deep Field (SXDF), made with the Very Large Array in B and C configurations. We present a radio map of the entire field, and a catalogue of 505 sources covering 0.8 square degrees to a peak flux density limit of 100 microJy. Robust optical identifications are provided for 90% of the sources, and suggested IDs are presented for all but 14 (of which 7 are optically blank, and 7 are close to bright contaminating objects). We show that the optical properties of the radio sources do not change with flux density, suggesting that AGNs continue to contribute significantly at faint flux densities. We test this assertion by cross-correlating our radio catalogue with the X-ray source catalogue and conclude that radio-quiet AGNs become a significant population at flux densities below 300 microJy, and may dominate the population responsible for the flattening of the radio source counts if a significant fraction of them are Compton-thick.
[150]  [pdf] - 142674
The SCUBA HAlf Degree Extragalactic Survey (SHADES) -- II. Submillimetre maps, catalogue and number counts
Comments: 36 pages including 21 figures and 8 tables. Accepted for publication in MNRAS. Plain-text source catalogue available on the SHADES website
Submitted: 2006-09-01
We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]
[151]  [pdf] - 84501
Galaxy clusters at 0.6 < z < 1.4 in the UKIDSS Ultra Deep Survey Early Data Release
Comments: 6 pages, 1 figure, accepted for publication in MNRAS. For full-resolution version see
Submitted: 2006-08-29
We present the first cluster catalogue extracted from the UKIRT Infrared Deep Sky Survey Early Data Release. The catalogue is created using UKIDSS Ultra Deep Survey infrared J and K data combined with 3.6 micro-m and 4.5 micro-m Spitzer bands and optical BVRi'z' imaging from the Subaru Telescope over 0.5 square degrees in the Subaru XMM-Newton Deep Field. We have created a new cluster-detection algorithm, based on the Friends-Of-Friends and Voronoi Tessellation methods, which utilises probability distribution functions derived from a photometric redshift analysis. We employ mock catalogues to understand the selection effects and contamination associated with the algorithm. The cluster catalogue contains 13 clusters at redshifts 0.61 <= z <= 1.39 with luminosities 10 L* <~ L_tot <~ 50 L*, corresponding to masses 5 x 10^13 M_sun <~ M_cluster <~ 3 x 10^14 M_sun for (M/M_sun) / (L/L_sun) = 75h. The measured sky surface density of ~ 10 deg^-2 for high-redshift (z=0.5-1.5), massive (>10^14 M_sun) clusters is precisely in line with theoretical predictions presented by Kneissl et al. (2001).
[152]  [pdf] - 84305
Extremely red objects in the UKIDSS Ultra Deep Survey Early Data Release
Comments: Accepted for publication in MNRAS Letters
Submitted: 2006-08-21
We construct a sample of extremely red objects (EROs) within the UKIDSS Ultra Deep Survey by combining the Early Data Release with optical data from the Subaru/XMM-Newton Deep Field. We find a total of 3715 objects over 2013 sq. arcmin with R-K>5.3 and K<=20.3, which is a higher surface density than found by previous studies. This is partly due to our ability to use a small aperture in which to measure colours, but is also the result of a genuine overdensity of objects compared to other fields. We separate our sample into passively-evolving and dusty star-forming galaxies using their RJK colours and investigate their radio properties using a deep radio map. The dusty population has a higher fraction of individually-detected radio sources and a higher mean radio flux density among the undetected objects, but the passive population has a higher fraction of bright radio sources, suggesting that AGNs are more prevalent among the passive ERO population.
[153]  [pdf] - 82556
The discovery of a significant sample of massive galaxies at redshifts 5 < z < 6 in the UKIDSS Ultra Deep Survey Early Data Release
Comments: 12 pages, 7 figures, revised version, accepted for publication in MNRAS
Submitted: 2006-06-06, last modified: 2006-08-04
We have exploited the large area coverage of the combined UKIDSS Ultra Deep Survey (UDS) and Subaru/XMM-Newton Deep Survey (SXDS) to search for bright Lyman-break galaxies (LBGs) at z >= 5. Using the available optical+near-infrared photometry to efficiently exclude low-redshift contaminants, we identify nine z >= 5 LBG candidates brighter than z'=25(AB) within the 0.6 square degree overlap region between the UDS early data release (EDR) and the optical coverage of the SXDS. Accounting for selection incompleteness, we estimate the corresponding surface density of z >= 5 LBGs with z'<=25(AB) to be 0.005+/-0.002 per square arcmin. Modelling of the optical+near-infrared photometry constrains the candidates' redshifts to lie in the range 5.1 < z < 5.9, and provides estimates for their stellar masses. Although the stellar mass estimates are individually uncertain, a stacking analysis suggests that the typical stellar mass of the LBG candidates is >~5x10^10 Msun which, if confirmed, places them amongst the most massive galaxies currently known at z >= 5. It is found that Lambda CDM structure formation can produce sufficient numbers of dark matter halos at z >= 5 to accommodate our estimated number density of massive LBGs for plausible values of sigma_8 and the ratio of stellar to dark matter. Moreover, it is found that recent galaxy formation models can also account for the existence of such massive galaxies at z >= 5. Finally, no evidence is found for the existence of LBGs with stellar masses in excess of ~3x10^11 Msun at this epoch, despite the large co-moving volume surveyed.
[154]  [pdf] - 81451
A population of high-redshift type-2 quasars-I. Selection Criteria and Optical Spectra
Comments: 22 Pages, 12 Figures. Accepted by MNRAS
Submitted: 2006-04-18, last modified: 2006-05-15
We discuss the relative merits of mid-infrared and X-ray selection of type-2 quasars. We describe the mid-infrared, near-infrared and radio selection criteria used to find a population of redshift z~2 type-2 quasars which we previously argued suggests that most supermassive black hole growth in the Universe is obscured (Martinez-Sansigre et al., 2005). We present the optical spectra obtained from the William Herschel Telescope, and we compare the narrow emission line luminosity, radio luminosity and maximum size of jets to those of objects from radio-selected samples. This analysis suggests that these are genuine radio-quiet type-2 quasars, albeit the radio-bright end of this population. We also discuss the possibility of two different types of quasar obscuration, which could explain how the ~2-3:1 ratio of type-2 to type-1 quasars preferred by modelling our population can be reconciled with the ~1:1 ratio predicted by unified schemes.
[155]  [pdf] - 133755
The first Swift X-ray Flash: The faint afterglow of XRF 050215B
Comments: 8 pages, 5 figures, accepted for publication in ApJ
Submitted: 2006-05-10
We present the discovery of XRF 050215B and its afterglow. The burst was detected by the Swift BAT during the check-out phase and observations with the X-ray telescope began approximately 30 minutes after the burst. These observations found a faint, slowly fading X-ray afterglow near the centre of the error box as reported by the BAT. Infrared data, obtained at UKIRT after 10 hours also revealed a very faint K-band afterglow. The afterglow appear unusual since it is very faint, especially in the infrared with K>20 only 9 hours post burst. The X-ray and infrared lightcurves exhibit a slow, monotonic decay with alpha=0.8 and no evidence for steepening associated with the jet break to 10 days post burst. We discuss possible explanations for the faintness and slow decay in the context of present models for the production of X-ray Flashes.
[156]  [pdf] - 74972
Definitive Identification of the Transition between Small- to Large-Scale Clustering for Lyman Break Galaxies
Comments: Accepted for publication in ApJL. 5 pages, 4 figures. The text and Figures 2-4 have been revised. There is no major change which affects to the main discussion shown in the original preprint. This paper with high resolution figures is available at (PDF)
Submitted: 2005-08-02, last modified: 2005-11-21
We report angular correlation function (ACF) of Lyman Break Galaxies (LBGs) with unprecedented statistical quality on the basis of 16,920 LBGs at z=4 detected in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field. The ACF significantly departs from a power law, and shows an excess on small scale. Particularly, the ACF of LBGs with i'<27.5 have a clear break between the small and large-scale regimes at the angular separation of ~7'' whose projected length corresponds to the virial radius of dark halos with a mass of 10^11-12 Mo, indicating multiple LBGs residing in a single dark halo. Both on small (2''<theta<3'') and large (40''<theta<400'') scales, clustering amplitudes monotonically increase with luminosity for the magnitude range of i'=24.5-27.5, and the small-scale clustering shows a stronger luminosity dependence than the large-scale clustering. The small-scale bias reaches b~10-50, and the outskirts of small-scale excess extend to a larger angular separation for brighter LBGs. The ACF and number density of LBGs can be explained by the cold dark matter model.
[157]  [pdf] - 75706
Most supermassive black hole growth is obscured by dust
Comments: 4 pages, 4 figures, "Granada Workshop on High Redshift Radiogalaxies", correction to elliptical galaxy in figure1
Submitted: 2005-09-06, last modified: 2005-09-18
We present an alternative method to X-ray surveys for hunting down the high-redshift type-2 quasar population, using Spitzer and VLA data on the Spitzer First Look Survey. By demanding objects to be bright at 24 microns but faint at 3.6 microns, and combining this with a radio criterion, we find 21 type-2 radio-quiet quasar candidates at the epoch at which the quasar activity peaked. Optical spectroscopy with the WHT confirmed 10 of these objects to be type-2s with 1.4 < z < 4.2 while the rest are blank. There is no evidence for contamination in our sample, and we postulate that our 11 blank-spectrum candidates are obscured by kpc-scale dust as opposed to dust from a torus around the accretion disk. By carefully modelling our selection criteria, we conclude that, at high redshift, 50-80 % of the supermassive black hole growth is obscured by dust.
[158]  [pdf] - 75483
The Number Density of Old Passively-Evolving Galaxies at z=1 in the Subaru/XMM-Newton Deep Survey Field
Comments: 28 pages, accepted for publication in Astrophysical Journal. The full version of the paper including Fig.3 and Fig.4 (large size) in full resolution is put at
Submitted: 2005-08-27
We obtained the number counts and the rest-frame B-band luminosity function of the color-selected old passively-evolving galaxies (OPEGs) at z=1 with very high statistical accuracy using a large and homogeneous sample of about 4000 such objects with z' <25 detected in the area of 1.03 deg^2 in the Subaru/XMM-Newton Deep Survey (SXDS) field. Our selection criteria are defined on the i'-z' and R-z' color-magnitude plane so that OPEGs at z=0.9-1.1 with formation redshift z_f=2-10 are properly sampled. The limiting magnitude corresponds to the luminosity of galaxies with M_*+3 at z=0. We made a pilot redshift observations for 99 OPEG candidates with 19 < z' < 22 and found that at least 78% (73/93) of the entire sample, or 95% (73/77) of these whose redshifts were obtained are indeed lie between z=0.87 and 1.12 and the most of their spectra show the continuum break and strong Ca H and K lines, indicating that these objects are indeed dominated by the old stellar populations. We then compare our results with the luminosity functions of the color- or the morphologically-selected early type galaxies at z=0 taking the evolutionary factor into account and found that the number density of old passive galaxies with sim M_* magnitude at z~1 averaged over the SXDS area is 40-60% of the equivalently red galaxies and 60-85% of the morphologically-selected E/S0 galaxies at z=0 depending on their luminosity evolution. It is revealed that more than half, but not all, of the present-day early-type galaxies had already been formed into quiescent passive galaxies at z=1.
[159]  [pdf] - 880853
The obscuration by dust of most of the growth of supermassive black holes
Comments: 10 pages, 2 figures, to appear in Nature on 4th August. Title and opening paragraph modified to Nature house style. Minor rearrangement and correction to table 1. Supplementary Information can be found at
Submitted: 2005-05-24, last modified: 2005-07-29
Supermassive black holes underwent periods of exponential growth during which we seem them as quasars in the distant Universe. The summed emission from these quasars generates the cosmic X-ray background, the spectrum of which has been used to argue that most black-hole growth is obscured. There are clear examples of obscured black-hole growth in the form of `type-2' quasars, but their numbers are fewer than expected from modelling of the X-ray background. Here we report on the direct detection of a population of distant type-2 quasars which is at least comparable in size to the well-known unobscured type-1 population. We selected objects that have mid-infrared and radio emissions characteristic of quasars, but which are faint at near-infrared and optical wavelengths. This population is responsible for most of the black hole growth in the young Universe and, throughout cosmic history, black-hole growth occurs in the dusty, gas-rich centres of active galaxies.
[160]  [pdf] - 73513
DDO 43: A Prototypical Dwarf Irregular Galaxy?
Comments: 50 pages, 22 figures, accepted by The Astronomical Journal, available at
Submitted: 2005-06-06
We present sensitive and high resolution 21 cm observations of the dwarf irregular (dIm) galaxy DDO 43, in conjunction with optical broad and narrow band images in U, B, V, and Halpha. The observations are used to examine the relationship between its HI morphology and kinematics to past and present star formation. Optically, it is a small (R_25 = 990 pc), faint (M_B of -14.0) Im with a slightly boxy shape. In HI, DDO 43 has an extended (R_HI/R_Ho = 2.8) gas envelope. There is a high density ridge associated with the optical body of the galaxy, containing several higher density knots and lower density holes. The largest hole is ~850x530 pc. No expansion is detected, so it must be relatively old. The largest and potentially oldest (7-70 Myr) of the six identified star clusters is located at the western edge of the hole. Four of the other clusters are located near high density peaks. There are several HII regions, most of which are associated with peaks in the HI surface density. The overall star formation rate is average for its type. In many ways, DDO 43 is a very typical dwarf irregular galaxy. Its HI morphology is consistent with a history of episodes of localized star ormation that create holes and shells in the ISM, some of which can overlap. These features are located within the area of solid-body rotation in the galaxy; the lack of shear in these small systems allows such structures to persist for long periods of time.
[161]  [pdf] - 72122
Vigorous star formation in a bulge-dominated ERO at z = 1.34
Comments: Accepted by MNRAS 1 April 2005. Ten pages, six figures
Submitted: 2005-04-01
We present near-IR spectroscopy of three Extremely Red Objects (EROs) using the OHS/CISCO spectrograph at Subaru telescope. One target exhibits a strong emission line which we identify as H-alpha at z = 1.34. Using new and existing ground-based optical and near-IR imaging, and archival HST imaging, we argue that this target is essentially an elliptical galaxy, with an old stellar population of around 4E11 solar masses, but having a dust-enshrouded star-forming component with a SFR of some 50--100 solar masses per yr. There is no evidence that the galaxy contains an AGN. Analysis of a further two targets, which do not exhibit any features in our near-IR spectra, suggests that one is a quiescent galaxy in the redshift range 1.2 < z < 1.6, but that the other cannot be conclusively categorised as either star-forming or quiescent. Even though our first target has many of the properties of an old elliptical, the ongoing star formation means that it cannot have formed_all_ of its stellar population at high redshift. While we cannot infer any robust values for the star formation rate in ellipticals at z > 1 from this one object, we argue that the presence of an object with such a high SFR in such a small sample suggests that a non-negligible fraction of the elliptical galaxy population may have formed a component of their stellar population at redshifts z ~ 1--2. We suggest that this is evidence for ongoing star formation in the history of elliptical galaxies.
[162]  [pdf] - 71882
The luminosity dependence of the Type 1 AGN fraction
Comments: 8 pages, accepted for publication by MNRAS
Submitted: 2005-03-23
Using a complete, magnitude-limited sample of active galaxies from the Sloan Digital Sky Survey (SDSS) we show that the fraction of broad-line (Type 1) active galactic nuclei increases with luminosity of the isotropically-emitted [O III] narrow emission line. Our results are quantitatively in agreement with, and far less uncertain than, similar trends found from studies of X-ray and radio-selected active galaxies. While the correlation between broad-line fraction and luminosity is qualitatively consistent with the receding torus model, its slope is shallower and we therefore propose a modification to this model where the height of the torus increases slowly with AGN luminosity. We demonstrate that the faint-end slope of the AGN luminosity function steepens significantly when a correction for `missing' Type 2 objects is made and that this can substantially affect the overall AGN luminosity density extrapolated from samples of more luminous objects.
[163]  [pdf] - 71675
The Connection between Ultraviolet and X-ray Outflows in AGN: the case of PDS 456
Comments: Accepted for publication in MNRAS. 6 pages. 6 figures
Submitted: 2005-03-14
High-velocity outflows from AGN are a well-known phenomena in terms of the Broad Absorption Lines seen in the UV/optical. More recently, similar, possibly related, outflows have been reported in the X-ray. The most extreme example is seen in the nearby, luminous QSO PDS 456, which displays a massive, high velocity (50000 km s-1), high-ionization X-ray outflow of 10 solar masses yr-1. Here we present the UV spectrum of PDS 456 as observed by the Hubble Space Telescope. We find the UV spectrum is also extreme, displaying very broad emission-lines, with CIV 1549 blueshifted by 5000 km s-1 and a broad Ly-alpha absorption trough blueshifted by 14000-24000 km s-1. No strong, broad high-ionization absorption features are seen. We interpret the combined UV and X-ray spectrum of PDS 456 as the signature of a decelerating, cooling outflow, which may be driven by radiation and/or magnetic field. This outflow may be the source of some of the broad emission and absorption-line gas.
[164]  [pdf] - 70052
The Discovery of Primeval Large-Scale Structures with Forming Clusters at Redshift 6
Comments: 8 pages, 4 figures, to be published in ApJL. The paper with high resolution figures is available at
Submitted: 2004-12-28
We report the discovery of primeval large-scale structures (LSSs) including two proto-clusters in a forming phase at z=5.7. We carried out extensive deep narrow-band imaging in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field, and obtained a cosmic map of 515 Ly alpha emitters (LAEs) in a volume with a transverse dimension of 180 Mpc x 180 Mpc and a depth of ~40 Mpc in comoving units. This cosmic map shows filamentary LSSs, including clusters and surrounding 10-40 Mpc scale voids, similar to the present-day LSSs. Our spectroscopic follow-up observations identify overdense regions in which two dense clumps of LAEs with a sphere of 1-Mpc diameter in physical units are included. These clumps show about 130 times higher star-formation rate density, mainly due to a large overdensity, ~80, of LAEs. These clumps would be clusters in a formation phase involving a burst of galaxy formation.
[165]  [pdf] - 68888
A Chandra observation of the z=2.285 galaxy FSC10214+4724: Evidence for a Compton-thick quasar?
Comments: 5 pages, 3 figures; MNRAS letters, in press
Submitted: 2004-11-11
We present a ~20 ks Chandra ACIS-S observation of the strongly lensed z=2.285 ultra-luminous infrared galaxy FSC10214+4724. Although this observation achieves the equivalent sensitivity of an up-to ~4 Ms Chandra exposure (when corrected for gravitational lensing), the rest-frame 1.6-26.3 keV emission from FSC10214+4724 is weak (L_X~2E42 erg/s for a lensing boost of ~100); a significant fraction of this X-ray emission appears to be due to vigorous star-formation activity. If FSC10214+4724 hosts a quasar, as previously suggested, then it must be obscured by Compton-thick material. We compare FSC10214+4724 to high-redshift SCUBA galaxies and discuss the X-ray identification of Compton-thick AGNs at high redshift.
[166]  [pdf] - 67820
DDO 88: A Galaxy-Sized Hole in the Interstellar Medium
Comments: 44 pages; 16 figures. To appear in AJ, January 2005. Available from, cd pub/dah/papers/d88 and
Submitted: 2004-09-29
We present an HI and optical study of the gas-rich dwarf irregular galaxy DDO 88. Although DDO 88's global optical and HI parameters are normal for its morphological type, it hosts a large (3 kpc diameter) and unusually complete ring of enhanced HI emission. The gas ring is located at approximately one-third of the total HI radius and one-half the optically-defined Holmberg radius, and contains 30% of the total HI of the galaxy. The ring surrounds a central depression in the HI distribution, so it may be a shell formed by a starburst episode. However, the UBV colors in the HI hole are not bluer than the rest of the galaxy as would be expected if an unusual star-forming event had taken place there recently, but there is an old (~1-3 Gyr), red cluster near the center of the hole that is massive enough to have produced the hole in the HI. An age estimate for the ring, however, is uncertain because it is not observed to be expanding. An expansion model produces a lower estimate of 0.5 Gyr, but the presence of faint star formation regions associated with the ring indicate a much younger age. We also estimate that the ring could have dispersed by now if it is older than 0.5 Gyr. This implies that the ring is younger than 0.5 Gyr. A younger age would indicate that the red cluster did not produce the hole and ring. If this ring and the depression in the gas which it surrounds were not formed by stellar winds and supernovae, this would indicate that some other, currently unidentified, mechanism is operating.
[167]  [pdf] - 64910
Deep near-infrared spectroscopy of submillimetre-selected galaxies
Comments: 10 pages, to appear in MNRAS
Submitted: 2004-05-20
We present the results of deep near-infrared spectroscopy of seven submillimetre-selected galaxies from the SCUBA 8-mJy and CUDSS surveys. These galaxies were selected because they are too faint to be accessible to optical spectrographs on large telescopes. We obtain a spectroscopic redshift for one object, and likely redshifts for two more, based on a combination of marginal emission line detections and the shape of the continuum. All three redshifts broadly agree with estimates from their radio/submm spectral energy distributions. From the emission line strengths of these objects, we infer star formation rates of 10-25 Msun/yr, while the lack of detections in the other objects imply even lower rates. By comparing our results with those of other authors, we conclude it is likely that the vast majority (more than 90 per cent) of the star formation in these objects is completely extinguished at rest-frame optical wavelengths, and the emission lines originate in a relatively unobscured region. Finally, we look at future prospects for making spectroscopic redshift determinations of submm galaxies.
[168]  [pdf] - 64925
Neutral Hydrogen in Arp 158
Comments: 11 pages, 16 figures. Accepted for publication in AJ
Submitted: 2004-05-20
We present 21 cm observations of Arp 158. We have performed a study of the neutral hydrogen (HI) to help us understand the overall formation and evolution of this system. This is a disturbed system with distinct optical knots connected by a linear structure embedded in luminous material. There is also a diffuse spray to the southeast. The HI seems to be made up of three distinct, kinematically separate systems. Arp 158 bears a certain optical resemblance to NGC 520 (Arp 157), which has been identified as a mid-stage merger. From our 21 cm observations of Arp 158, we also see a comparable HI content with NGC 520. These similarities suggest that Arp 158 is also an intermediate stage merger.
[169]  [pdf] - 63656
Dust and gas obscuration in ELAIS Deep X-ray Survey reddened quasars
Comments: 12 pages, 10 figures, ApJ, in press
Submitted: 2004-03-18
Hard X-ray surveys have uncovered a large population of heavily obscured AGN. They also reveal a population of quasars with moderate obscuration at both visible and X-ray wavelengths. We use Chandra selected samples of quasars from the ELAIS Deep X-ray Survey (EDXS) and CDF-N to investigate the obscuration towards the nuclei of moderately obscured AGN. We find an inverse correlation between the optical to X-ray flux ratio and the X-ray hardness ratio which can be interpreted as due to obscuration at visible and X-ray wavelengths. We present detailed optical and near-infrared data for a sample of optically-faint (R>23) quasars from the EDXS. These are used to constrain the amount of rest-frame UV/optical reddening towards these quasars. It is found that optically-faint quasars are mostly faint due to obscuration, not because they are intrinsically weak. After correcting for reddening, the optical magnitudes of most of these quasars are similar to the brighter quasars at these X-ray fluxes. Combining with gas column densities inferred from the X-ray observations we consider the gas-to-dust ratios of the obscuring matter. We find that the quasars generally have higher gas-to-dust absorption than that seen in the Milky Way - similar to what has been found for nearby Seyfert galaxies. We consider the possible existence of a large population of X-ray sources which have optical properties of Type 1 (unobscured) quasars, but X-ray properties of Type 2 (obscured) quasars. We show that the observed distribution of optical-to-X-ray flux ratios of quasars at z>1 is skewed to low values compared to the intrinsic distribution...
[170]  [pdf] - 51498
ERO R1 in CL0939+4713 field - Evidence for an S0-like galaxy at z \sim 1.5
Comments: 23 pages, 9 figures, Submitted to the Astrophysical Journal
Submitted: 2002-09-05, last modified: 2003-03-12
We present further observations of the extremely red object ERO J094258+4659.2, identified by \citet{iye00} as ERO R1 in their deep images of the cluster A851. We estimate its redshift independently by eight-band photometric redshift determination and cross-correlation of a new H-band spectrum with the optical spectra of local E/S0 galaxies, and conclude that it lies at $z \sim 1.5$. Although its colors are consistent both with an elliptical galaxy and an S0 galaxy at that redshift, its elongated shape and exponential luminosity profile suggest the presence of an evolved stellar disk component. We rule out the possibility that these properties are strongly influenced by gravitational lensing by the foreground cluster, and therefore conclude that this object is more likely to be an S0-like galaxy, rather than a lensed elliptical. The H-band spectrum does not show strong H$\alpha$ emission and the star formation rate therefore appears to be very modest. The presence of such a galaxy with an apparently relaxed disk of stars at this high redshift provides a new and strong constraint on theoretical models which aim to explain the formation and evolution of galaxies.
[171]  [pdf] - 54069
Quasars are more luminous than radio galaxies - so what?
Comments: 6 pages. To appear in "Radio galaxies: past, present and future", eds M. Jarvis et al., Leiden, Nov 2002
Submitted: 2003-01-06
Surveys to find high-redshift radio galaxies deliberately exclude optically-bright objects, which may be distant radio-loud quasars. In order to properly determine the space density of supermassive black holes, the fraction of such objects missed must be determined within a quantitative framework for AGN unification. I briefly describe the receding torus model, which predicts that quasars should have more luminous ionizing continua than radio galaxies of similar radio luminosity, and present evidence to support it. I also suggest two further tests of the model which should constrain some of its parameters.
[172]  [pdf] - 52286
Obscured AGN from the ELAIS Deep X-ray Survey
Comments: 14 pages, 11 figures, MNRAS in press
Submitted: 2002-10-10
The sources discovered in deep hard X-ray surveys with 2-8 keV fluxes of 10^-14 erg cm^-2 s^-1 make up the bulk of the X-ray background at these energies. We present here detailed multi-wavelength observations of three such sources from the ELAIS Deep X-ray Survey. The observations include sensitive near-infrared spectroscopy with the Subaru Telescope and X-ray spectral information from the Chandra X-ray Observatory. The sources observed all have optical-to-near-IR colours redder than an unobscured quasar and comprise a reddened quasar, a radio galaxy and an optically-obscured AGN. The reddened quasar is at a redshift z=2.61 and shows a very large X-ray absorbing column of N_H approx 3.10^23 cm^-2. This contrasts with the relatively small amount of dust reddening, implying a gas-to-dust ratio along the line-of-sight a hundred times greater than that of the Milky Way. The radio galaxy at z=1.57 shows only narrow emission lines, but has a surprisingly soft X-ray spectrum. The softness of this spectrum either indicates an unusually low gas-to-dust ratio for the absorbing medium or X-ray emission related to the young radio source. The host galaxy is extremely red (R-K=6.4) and its optical/near-IR spectrum is best fit by a strongly reddened (A_V~2) starburst. The third X-ray source discussed is also extremely red (R-K=6.1) and lies in a close grouping of three other R-K>6 galaxies. No emission or absorption lines were detected from this object, but its redshift (and that of one of the nearby galaxies) are constrained by SED-fitting to be just greater than z=1. The extremely red colours of these two galaxies can be accounted for by old stellar populations. These observations illustrate the diverse properties of hard X-ray selected AGN.
[173]  [pdf] - 48467
The radio galaxy 3C 356 and clues to the trigger mechanisms for powerful radio sources
Comments: 12 pages, accepted by MNRAS. Fig.1 has been reduced in resolution. Corrected citation to data in Fig.1
Submitted: 2002-03-25, last modified: 2002-04-04
We present deep near-infrared images, taken with Subaru Telescope, of the region around the z=1.08 radio source 3C 356 which show it to be associated with a poor cluster of galaxies. We discuss evidence that this cluster comprises two subclusters traced by the two galaxies previously proposed as identifications for 3C 356, which both seem to harbour AGN, and which have the disturbed morphologies expected if they underwent an interpenetrating collision at the time the radio jets were triggered. We explain the high luminosity and temperature of the diffuse X-ray emission from this system as the result of shock-heating of intracluster gas by the merger of two galaxy groups. Taken together with the results on other well-studied powerful radio sources, we suggest that the key ingredient for triggering a powerful radio source, at least at epochs corresponding to z~1, is a galaxy--galaxy interaction which can be orchestrated by the merger of their parent subclusters. This provides an explanation for the rapid decline in the number density of powerful radio sources since z~1. We argue that attempts to use distant radio-selected clusters to trace the formation and evolution of the general cluster population must address ways in which X-ray properties can be influenced by the radio source, both directly, by mechanisms such as inverse-Compton scattering, and indirectly, by the fact that the radio source may be preferentially triggered at a specific time during the formation of the cluster.
[174]  [pdf] - 45723
The Optical Counterpart of the Accreting Millisecond Pulsar SAX J1808.4-3658 in Outburst: Constraints on the Binary Inclination
Comments: 4 pages including 1 figure. Accepted for publication in ApJ Letters
Submitted: 2001-10-30
We present multiband optical/infrared photometry of V4580 Sgr, the optical counterpart of the accretion-powered millisecond pulsar SAX J1808.4-3658, taken during the 1998 X-ray outburst of the system. The optical flux is consistent with emission from an X-ray-heated accretion disk. Self-consistent modeling of the X-ray and optical emission during the outburst yields best-fit extinction Av=0.68 and inclination cos(i)=0.65, assuming a distance of 2.5 kpc. The allowed inclination range requires that the pulsar's stellar companion has extremely low mass, 0.05-0.10 Msun. Some of the infrared observations are not consistent with disk emission and are too bright to be from either the disk or the companion, even in the presence of X-ray heating.
[175]  [pdf] - 1232671
JHK Spectra of the z=2.39 Radio Galaxy 53W002
Comments: 10 pages, including 11 figures. Accepted for publication in PASJ(2001). Revised 5/15/2001
Submitted: 2001-04-29, last modified: 2001-05-16
We present low-resolution, near-IR JHK spectra of the weak z=2.39 radio galaxy 53W002, obtained with the OH-airglow Suppressor spectrograph (OHS) and Cooled Infrared Spectrograph and Camera for OHS (CISCO) on the Subaru Telescope. They cover rest-frame wavelengths of 3400-7200 A, and the emission lines of [O II]3727, Hb, [O III]4959, 5007, Ha, [N II]6548, 6583 and [S II]6716, 6731 were detected. Using the Ha/Hb line ratio, we find an extinction of E(B-V)=0.14. The emission-line ratios are reproduced by a cloud of electron density n_e=1x10^{3-4}(/cm3) with solar metallicity, ionized by an alpha=-0.7 power-law continuum with ionizing parameter U=1x10^-3. In addition to these emission lines, we make the first spectroscopic confirmation of the Balmer discontinuity in a high-z radio galaxy. Together with rest-frame UV photometry from the literature, we show that at least 1/3 of the present stellar mass was formed in the current starburst. The stellar mass was estimated to be (1-1.4)x10^11 M_sol by one-component model fitting, which is smaller than that of typical z~1 B2/6C radio galaxies. We suggest that 53W002 is currently assembling a large part of its stellar mass through merger events with the surrounding sub-galactic clumps, some of which can be identified with the Lya emitters detected in narrow-band imaging. After a few such events over the next few Gyr, 53W002 will evolve into a massive elliptical galaxy.
[176]  [pdf] - 39685
The non-stellar infrared continuum of Seyfert galaxies
Comments: Accepted for publication in AJ, March 2001 issue
Submitted: 2000-12-05
JHKL'M (1-5 micron) imaging of a sample of Seyfert (Sy) 2 galaxies is presented. We confirm that the 1-2.2 micron nuclear continuum of some Sy 2s is dominated by stellar emission, whereas the continuum emission at longer wavelengths (3-5 micron) is almost entirely non-stellar in origin. The IR non-stellar spectral energy distributions (SED) (up to 15 micron) of Sy galaxies show a variety of shapes, and they are well reproduced with the tapered disk models of Efstathiou & Rowan-Robinson (1995). We have used two models, one including an optically thin cone component and a coneless model. Although our modelling of the SEDs does not allow us to favor either model to account for all the observed SEDs, we find that the viewing angle towards the central source is well constrained by both models. We have also investigated non-stellar color-color diagrams. The colors of the Sy galaxies with viewing angles theta < 30 degree are better reproduced with the cone model. These diagrams provide a good means to separate Sy 2s with moderate obscuration (A_V < 20 mag, from hard X-ray observations) from those with high obscuration. The ground-based 4.8 micron and ISO 9.6 micron luminosities are well correlated with the hard X-ray luminosities of Sy 1s and 2s. These continuum emissions appear as a good indicator of the AGN luminosity, at least in the cases of hard X-ray Compton-thin Sy galaxies. We finally stress the finding that some Compton thick galaxies show bright non-stellar emission at 5 micron. This suggests that the near-IR emission in Sys is produced in an extended component illuminated by the central source, that is more visible from all viewing angles. We discuss possible implications of mid-IR surveys for the search of counterparts of highly obscured hard X-ray sources. (Abridged)
[177]  [pdf] - 37313
1-5 micron imaging of 3CRR galaxies: The K-z relation and the geometry of the torus
Comments: Accepted by MNRAS
Submitted: 2000-07-31
It has been claimed by Taylor et al. that the low-redshift end of the K-z relation for radio galaxies is too bright by about half a magnitude due to contributions from the obscured quasar nuclei. Such a result has major implications for the use of the K-band Hubble diagram in understanding the cosmological evolution of radio galaxies. In this paper we present 1-5 micron imaging data of a nearly-complete sample of low-redshift radio galaxies; this approach allows us to accurately determine the strengths of any unresolved nuclear components in the galaxies. We detect nuclear sources in five targets, whose broad-band colours are consistent with reddened quasar spectra. In all five cases the ratio of the inferred intrinsic near-infrared luminosity to the narrow-line luminosity is typical of quasars. We find a correlation between the inferred nuclear extinction and core-to-lobe ratio, which places constraints on the geometry of the torus. We find evidence for a shift of the K-z relation to fainter magnitudes, but by a much smaller amount (~0.1 mag) than Taylor et al. determined. Under the assumption that the nuclear sources in radio galaxies have the same intrinsic near-infrared spectra as quasars, our multi-wavelength images allow us to limit any possible shift to less than 0.3 magnitudes.
[178]  [pdf] - 36527
Hubble Space Telescope Imaging of the Circinus Galaxy
Comments: With 12 figures, to be published in the Astronomical Journal, September 2000 issue
Submitted: 2000-06-09
We present an HST imaging study of the nearby (4 Mpc distant) Circinus galaxy, which contains the nearest type 2 Seyfert nucleus and prominent circumnuclear star formation. Images have been obtained in the [OIII]$\lambda$5007, H$\alpha$ and H$_2$ v=1--0 S(1) emission lines, and in the green (5470 \AA), red (8140 \AA) and near infrared (2.04 and 2.15 $\mu$m) continua. An image in the [Fe II] $\lambda$1.644 $\mu$m line has been taken with a ground-based telescope. The [OIII] and H$\alpha$ images reveal the detailed structure of the complex of streamers and knots of high-excitation gas, which extends out of the galaxy disk. The morphology some 250 pc from the nucleus strongly suggests that the high-excitation gas is concentrated on the surface of a hollow cone with apex close to the nucleus. Such a structure may result through entrainment of dense gas from a circumnuclear torus in the galaxy disk by a low density, outflowing wind or jet. Within 40 pc of the nucleus, the high-excitation gas takes the form of a striking, filled V-shaped structure extending in the same direction as the larger scale high-excitation emission. This `V' can be described as an ionization cone, though a matter-bounded structure is also possible. The implied collimation of the ionizing photons or gaseous outflow must occur within 2 pc of the apex of the cone, presumed to be the location of the nucleus. The H$\alpha$ image shows a complex structure of HII regions, including the well known starburst ring of radius 150 - 270 pc. In addition, there is a more compact (40 pc radius), elliptical ring of HII regions around the `ionization cone'. We argue that this latter ring, which we call the nuclear ring, is intrinsically circular and located in the plane of the galaxy disk. Abstract truncated.
[179]  [pdf] - 36341
Infrared photometry of z~1 3C quasars
Comments: 6 pages, 3 figures, accepted by MNRAS Letters
Submitted: 2000-05-30
We present JHKL' photometry of a complete sample of steep-spectrum radio-loud quasars from the revised 3CR catalogue in the redshift range 0.65 < z < 1.20. After correcting for contributions from emission lines and the host galaxies, we investigate their spectral energy distributions (SEDs) around 1 micron. About 75% of the quasars are tightly grouped in the plane of optical spectral index, alpha_lo, versus near-infrared spectral index, alpha_hi, with the median value of alpha_lo close to the canonical value, and the median alpha_hi slightly flatter. We conclude that the fraction of moderately-obscured, red quasars decreases with increasing radio power, in accordance with the `receding torus' model which can also explain the relatively flat median near-infrared spectra of the 3CR quasars. Two of the red quasars have inverted infrared spectral indices, and we suggest that their unusual SEDs might result from a combination of dust-scattered and transmitted quasar light.
[180]  [pdf] - 35947
Studies of the high luminosity quasar, PDS 456
Comments: Poster contribution presented at the Joint MPE,AIP,ESO workshop on NLS1s, Bad Honnef, Dec. 1999, to appear in New Astronomy Reviews; also available at
Submitted: 2000-05-08
X-ray and multi-wavelength observations of the most luminous known local (z<0.3) AGN, the recently discovered radio-quiet quasar PDS 456, are presented. The spectral energy distribution shows that PDS 456 has a bolometric luminosity of 1e47 erg/s, peaking in the UV. The X-ray spectrum obtained by ASCA and RXTE shows considerable complexity. The most striking feature observed is a deep, highly-ionised, iron K edge (8.7 keV, rest-frame), originating via reprocessing from highly ionised material, possibly the inner accretion disk. PDS 456 was found to be remarkably variable for its luminosity; in one flare the X-ray flux doubled in just about 15 ksec. If confirmed this would be an unprecedented event in a high-luminosity source, with a light-crossing time corresponding to about 2RS. The implications are that either flaring occurs within the very central regions, or else that PDS 456 is a super-Eddington or relativistically beamed system.
[181]  [pdf] - 1232428
[OII]3727 Emission from the Companion to the Quasar BR 1202-0725 at z=4.7
Comments: 15 pages including 3 figures. Accepted for publication in PASJ
Submitted: 2000-03-08
Results of a narrow-band imaging for the redshifted [OII]3727 emission around a quasar at z=4.7 obtained with the Subaru telescope and CISCO (a Cassegrain near infrared camera) are presented. A significant emission line is detected in the narrow-band H_2 (v=1-0 S(1)) filter at a location 2.4" northwest from the quasar, where the presence of a companion has been reported in Lyman alpha emission and the rest-frame UV continuum. We identify this line as [OII]3727 emission and confirm that the source really is a companion at z=4.7. The [OII]3727 flux from the companion is estimated to be 2.5 x 10^{-17} erg s^{-1} cm^{-2}. If the companion is a star forming object, the inferred star formation rate is as high as 45-230 M_{solar} yr^{-1} even without assuming the extinction correction. This value is higher than those derived from the Lyman alpha emission or from the UV continuum. Thus, provided that the difference is caused by dust extinction, the extinction corrected star formation rate is calculated to be 45 to 2300 M_{solar} yr^{-1} depending on the assuming extinction curves.
[182]  [pdf] - 1232418
High-Resolution Near-Infrared Imaging of the Powerful Radio Galaxy 3C 324 at z = 1.21 with the Subaru Telescope
Comments: 21 pages (10 figures), accepted for publication in PASJ
Submitted: 2000-02-21
We have obtained high-resolution K'-band images of the powerful z=1.206 radio galaxy 3C 324 with the Subaru telescope under seeing conditions of 0.3--0.4 arcsec. We clearly resolved the galaxy and directly compared it to the optical images obtained with the Hubble Space Telescope. The host galaxy of 3C 324 is revealed to be a moderately luminous elliptical galaxy with a smooth light profile. The effective radius of the galaxy, as determined by profile fitting, is 1.3+-0.1 arcsec (1.2 kpc), which is significantly smaller than the value of 2.2 arcsec, published in Best et al. (1998, MNRAS, 292, 758). The peak of the K'-band light coincides with the position of the radio core, which implies that the powerful AGN lies at the nucleus of the host galaxy. The peak also coincides with the gap in the optical knotty structures which may be a dust lane hiding the UV-optical emission of the AGN from our line of sight; it is very likely that we are seeing the obscuring structure almost edge-on. We clearly detected the `aligned component' in the K'-band image by subtracting a model elliptical galaxy from the observed image. The red R_F702W-K color of the outer region of the galaxy avoiding the aligned component indicates that the near infrared light of the host galaxy is dominated by an old stellar population.
[183]  [pdf] - 1232415
Subaru Observations for the K-band Luminosity Distribution of Galaxies in Clusters near to 3C 324 at z$\sim$1.2
Comments: 12 pages, 9 figures, accepted for publication in PASJ
Submitted: 2000-02-17
We investigate the $K$-band luminosity distribution of galaxies in the region of clusters at $z\sim1.2$ near to the radio galaxy 3C 324. The imaging data were obtained during the commissioning period of the Subaru telescope. There is a significant excess of the surface number density of the galaxies with $K =$ 17--20 mag in the region within $\sim$ 40'' from 3C 324. At this bright end, the measured luminosity distribution shows a drop, which can be represented by the exponential cut off of the Schechter-function formula; the best-fitted value of the characteristic magnitude, $K^{*}$, is $\sim 18.4\pm0.8$. This measurement follows the evolutionary trend of the $K^*$ of the rich clusters observed at an intermediate redshift, which is consistent with passive evolution models with a formation redshift $z_f \gtsim 2$. At $K \gtsim 20$ mag, however, the excess of the galaxy surface density in the region of the clusters decreases abruptly, which may imply that the luminosity function of the cluster galaxies has a negative slope at the faint end. This may imply strong luminosity segregation between the inner and outer parts of the clusters, or some deficit of faint galaxies in the cluster central region of the cluster.
[184]  [pdf] - 1232416
Color-Magnitude Sequence in the Clusters at z$\sim$1.2 near the Radio Galaxy 3C 324
Comments: 16 pages, 10 figures, accepted for publication in PASJ
Submitted: 2000-02-17
We have investigated the optical and near-infrared colors of K'-selected galaxies in clusters at z ~ 1.2 near to the radio galaxy 3C 324 using images obtained with the Subaru telescope and archival HST data. The distribution of colors of the galaxies in the cluster region is found to be fairly broad, and it may imply significant scatter in their star-formation histories, although the effect of contamination of field galaxies is uncertain. The red sequence of galaxies whose R-K colors are consistent with passive evolution models for old galaxies is found to be truncated at K' ~ 20 mag, and there are few fainter galaxies with similar red colors in the cluster region. We find that the bulge-dominated galaxies selected by quantitative morphological classification form a broad sequence in the color-magnitude diagram, whose slope is much steeper than that expected from metallicity variations within a passively evolving coeval galaxy population. We argue that the observed color-magnitude sequence can be explained by metallicity and age variations, and the fainter galaxies with K' > 20 mag may be 1-2 Gyr younger than the brighter galaxies. Some spatial segregation of the color and K'-band luminosity is seen in the sky distribution; the redder and the brighter objects tend to be located near 3C 324.
[185]  [pdf] - 109603
PDS 456: an Extreme Accretion Rate Quasar?
Comments: 5 pages, 6 postscript figures. Accepted for publication in MNRAS Letters
Submitted: 1999-11-25
We present quasi-simultaneous ASCA and RXTE observations of the most luminous known AGN in the local (z<0.3) Universe, the recently discovered quasar PDS 456. Multiwavelength observations have been conducted which show that PDS 456 has a bolometric luminosity of 10^47 erg/s peaking in the UV part of the spectrum. In the X-ray band the 2-10 keV (rest-frame) luminosity is 10^45 erg/s. The broad-band X-ray spectrum obtained with ASCA and RXTE contains considerable complexity. The most striking feature observed is a very deep, ionised iron K edge, observed at 8.7 keV in the quasar rest-frame. We find that these features are consistent with reprocessing from highly ionised matter, probably the inner accretion disk. PDS 456 appeared to show a strong (factor of 2.1) outburst in just 17ksec, although non-intrinsic sources cannot be completely ruled out. If confirmed, this would be an unusual event for such a high-luminosity source, with the light-crossing-time corresponding to 2 Schwarzschild radii. The implication would be that flaring occurs within the very central regions, or else that PDS 456 is a `super-Eddington' or relativistically beamed system. Overall we conclude on the basis of the extreme blue/UV luminosity, the rapid X-ray variability and from the imprint of highly ionised material on the X-ray spectrum, that PDS 456 is a quasar with an unusually high accretion rate.
[186]  [pdf] - 107223
Young stars and non-stellar emission in the aligned radio galaxy 3C 256
Comments: 35 pages including 10 figures; to appear in Nov 10 ApJ
Submitted: 1999-06-30
We present ground-based images of the z=1.824 radio galaxy 3C 256 in the standard BVRIJHK filters and an interference filter centered at 8800A, a Hubble Space Telescope image in a filter dominated by Ly-alpha emission (F336W), and spectra covering rest-frame wavelengths from Ly-alpha to [O III] 5007. Together with published polarimetry observations, we use these to decompose the overall spectral energy distribution into nebular continuum emission, scattered quasar light, and stellar emission. The nebular continuum and scattered light together comprise half (one third) of the V-band (K-band) light within a 4-arcsec aperture, and are responsible for the strong alignment between the optical/near-infrared light and the radio emission. The stellar emission is dominated by a population estimated to be 100-200 Myr old (assuming a Salpeter IMF), and formed in a short burst with a peak star formation rate of 1-4x10^3 Msun/yr. The total stellar mass is estimated to be no more than 2x10^{11} Msun, which is far less than other luminous radio galaxies at similar redshifts, and suggests that 3C 256 will undergo further star formation or mergers.
[187]  [pdf] - 105472
The detection and photometric redshift determination of distant galaxies using SIRTF's Infrared Array Camera
Comments: 28 pages incl 12 figures; to appear in June 1999 PASP. Fig.12 replaced with corrected version
Submitted: 1999-03-03, last modified: 1999-04-23
We investigate the ability of the Space Infrared Telescope Facility's Infrared Array Camera to detect distant (z ~ 3)galaxies and measure their photometric redshifts. Our analysis shows that changing the original long wavelength filter specifications provides significant improvements in performance in this and other areas.
[188]  [pdf] - 105106
Thermal-infrared imaging of 3C radio galaxies at z~1
Comments: 17 pages incl 14 figures, accepted by MNRAS
Submitted: 1999-02-04
We present the results of a programme of thermal-IR imaging of nineteen z~1 radio galaxies from the 3CR and 3CRR samples. We detect emission at L' (3.8um) from four objects; in each case the emission is unresolved at 1" resolution. Fifteen radio galaxies remain undetected to sensitive limits of L'~15.5. Using these data in tandem with archived HST data and near-IR spectroscopy we show that 3 of the detected `radio galaxies' (3C22, 3C41, and 3C65) harbour quasars reddened by Av<5. Correcting for this reddening 3C22 and 3C41 are very similar to coeval 3C quasars, whilst 3C65 seems unusually underluminous. The fourth radio galaxy detection (3C265) is a more highly obscured (Av~15) but otherwise typical quasar which previously has been evident only in scattered light. We determine the fraction of dust-reddened quasars at z~1 to be 28(+25)(-13)% at 90% confidence. On the assumption that the undetected radio galaxies harbour quasars similar to those in 3C22, 3C41 and 3C265 (as seems reasonable given their similar narrow emission line luminosities) we deduce extinctions of Av>15 towards their nuclei. The contributions of reddened quasar nuclei to the total K-band light ranges from ~0 per cent for the non-detections, through ~10 per cent for 3C265 to ~80 per cent for 3C22 and 3C41. Correcting for these effects does not remove the previously reported differences between the K magnitudes of 3C and 6C radio galaxies, so contamination by reddened quasar nuclei is not a serious problem for drawing cosmological conclusions from the K-z relation for radio galaxies. We discuss these results in the context of the `receding torus' model which predicts a small fraction of lightly-reddened quasars in samples of high radio luminosity sources. We also examine the likely future importance of thermal-IR imaging in the study of distant powerful radio sources.
[189]  [pdf] - 104172
Optical and infrared observations of the luminous quasar PDS 456: a radio-quiet analogue of 3C 273?
Comments: Accepted by MNRAS (Letters), 5 pages incl. 5 figures
Submitted: 1998-12-02
We present infrared photometry and optical and infrared spectroscopy of the recently-discovered, extremely luminous nearby quasar PDS 456. A number of broad emission features are seen in the near-infrared which we are unable to identify. We measure a more accurate redshift from a narrow forbidden emission line and compare the optical-infrared spectrum to that of 3C 273. The close similarity suggests that PDS 456 is a radio-quiet analogue of 3C 273, although radio observations do not support this idea.
[190]  [pdf] - 102566
Beaming models and the correlations between the core dominance parameter and the core and extended powers of quasars
Comments: Accepted by A&A Letters: 4 pages, uses LaTeX2e and A&A macros
Submitted: 1998-08-20
We investigate the recent claim by Qin et al. that the observed correlations (or lack thereof) between the core dominance parameter and the core and extended powers of samples of lobe- and core-dominated quasars is in contradiction with beaming models. Contrary to their conclusion, we find that their results are in perfect agreement with such models, and support this assertion with Monte Carlo simulations.
[191]  [pdf] - 102307
The nuclear spectrum of the radio galaxy NGC 5128 (Centaurus A)
Comments: 12 pages including 2 figures, accepted by ApJ Letters
Submitted: 1998-07-29
We present near-infrared spectra of the nuclear disk in the nearby radio galaxy NGC 5128 (Centaurus A). On the basis of the observed strengths of the [S III] 0.9532um and [Fe II] 1.2567um lines, we classify NGC 5128 as a LINER. Modeling of the strengths of these and additional lines suggests that the nuclear region is powered by shocks rather than photoionization.
[192]  [pdf] - 102185
Near-infrared and X-ray obscuration to the nucleus of the Seyfert 2 galaxy NGC 3281
Comments: 23 pages including 9 figures
Submitted: 1998-07-17
We present the results of a near-infrared and X-ray study of the Seyfert 2 galaxy NGC 3281. Emission from the Seyfert nucleus is detected in both regions of the electromagnetic spectrum, allowing us to infer both the equivalent line of sight hydrogen column density, N_H = 71.0(+11.3,-12.3)e26/m^2 and the extinction due to dust, A_V = 22+/-11 magnitudes (90% confidence intervals). We infer a ratio of N_H/A_V which is an order of magnitude larger than that determined along lines of sight in the Milky Way and discuss possible interpretations. We consider the most plausible explanation to be a dense cloud in the foreground of both the X-ray and infrared emitting regions which obscures the entire X-ray source but only a fraction of the much larger infrared source.
[193]  [pdf] - 100558
A new look at the isotropy of narrow line emission in extragalactic radio sources
Comments: 4 pages, accepted by MNRAS
Submitted: 1998-03-03
We undertake a quantitative investigation, using Monte Carlo simulations, of the amount by which quasars are expected to exceed radio galaxies in optical luminosity in the context of the `receding torus' model. We compare these simulations with the known behaviour of the [O III] 5007 and [O II] 3727 emission lines and conclude that [O III] is the better indicator of the strength of the underlying non-stellar continuum.
[194]  [pdf] - 99935
HST and VLA Observations of Seyfert 2 Galaxies: The Relationship between Radio Ejecta and the Narrow Line Region
Comments: ApJ, in press, (AAS, emulateapj) LaTex, 22 pages, 8 PS figures, 3 GIF Figures, 2 JPEG color plates, preprint and ALL figures (at full resolution) as PS file available at
Submitted: 1998-01-12
We present HST/WFPC2 images, in narrow-band filters containing the [OIII] lambda5007 and Halpha+[NII] emission-lines and their adjacent continua, of a sample of seven Seyfert 2 galaxies selected on the basis of possessing either extended emission-line regions in ground-based observations or a hidden broad-line region in polarized light. Six of the galaxies have also been observed with the VLA to obtain radio maps with better quality and angular resolution than those in the literature. We find detailed correspondences between features in the radio and emission-line images that clearly indicate strong interactions between the radio jets and the interstellar medium. Such interactions play a major role in determining the morphology of the NLR, as the radio jets sweep up and compress ambient gas, producing ordered structures with enhanced surface brightness in line emission. In at least three galaxies, namely Mrk 573, ESO 428-G14, and Mrk 34, and perhaps also NGC 7212, off-nuclear radio lobes coincide with regions of low gaseous excitation (as measured by the [OIII]/(Halpha+[NII]) ratio). In Mrk 573 and NGC 4388, there is a clear trend for low brightness ionized gas to be of higher excitation. Those results may be understood if radio lobes and regions of high emission-line surface brightness are associated with high gas densities, reducing the ionization parameter. [OIII]/(Halpha+[NII]) excitation maps reveal bi-polar structures which can be interpreted as either the `ionization cones' expected in the unified scheme or widening, self-excited gaseous outflows. Only NGC 4388 and Mrk 573 show a clearly defined, straight-edged ionization cone.
[195]  [pdf] - 99736
Molecular Hydrogen and Paschen-alpha Emission in Cooling Flow Galaxies
Comments: ApJ Letters, in press, AAS LaTex, preprint also available at
Submitted: 1997-12-18
We present near-infrared spectra obtained to search for Pa-alpha and molecular hydrogen lines in edge-darkened (FR I-type) radio galaxies with bright Halpha emission in the redshift range 0.0535<z<0.15. We find that all three galaxies in our sample (PKS 0745-191, PKS 1346+26, & PKS2322-12) which are associated with strong cooling flows also have strong Pa-alpha and H_2 (1-0) S(1) through S(5) emission, while other radio galaxies do not. Together with earlier observations this confirms claims that cooling flow galaxies are copious emitters of molecular hydrogen with large H_2 (1-0) S(3)/Pa-alpha ratios in the range 0.5 to 2. The emission is centrally concentrated within the inner few kiloparsec and could come from warm (T ~ 1000-1500 K) molecular material which is being deposited by the cooling flow. We speculate that the H_2 emission could be related to the interaction between the jets and this molecular gas.
[196]  [pdf] - 95156
Helical Strands in the Jet-like Narrow Line Region of ESO 428-G14
Comments: ApJ Letters, accepted for publication, 10 pages, 2 PS Figures, AASTeX, also available at
Submitted: 1996-08-08
We present HST/WFPC2 images of the narrow line region (NLR) of the Seyfert 2 galaxy ESO 428-G14 (0714-2914, M4-1). The NLR consists of many individual, thin strands, which are very closely related to the radio jet and produce a highly complex yet ordered structure. We find that the jet is two-sided with a double-helix of emission-line gas apparently wrapped around the NW side. To the SE, the jet seems to be deflected at a ridge of highly excited gas. The strands to the SE may also wrap around the radio jet, but here complete helices are not seen. The overall structure is reminiscent of the jet seen in NGC 4258. Faint symmetric features aligned with the nucleus could indicate the presence of a highly collimated beam of photons or plasma from the center.