Source author record

R. J. Smith

R. J. Smith appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

69works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

69 published item(s)

preprint2026arXiv

PHANGS-JWST: the largest extragalactic molecular cloud catalog traced by polycyclic aromatic hydrocarbon emission

High-resolution JWST images of nearby spiral galaxies reveal polycyclic aromatic hydrocarbon (PAH) emission structures that trace molecular gas, including CO-dark regions. We identify ISM cloud structures in PHANGS-JWST 7.7 $μ$m PAH maps for 66 galaxies, smoothed to 30 pc and at native resolution, extracting 108,466 and 146,040 clouds, respectively. Molecular properties were inferred using a linear conversion from PAH to CO. Given the tendency for clouds in galaxy centers to overlap in velocity space, we opted to flag these and omit them from the analysis in this work. The remaining clouds correspond to giant molecular clouds, such as those detected in CO(2-1) emission by ALMA, or lower surface density clouds that either fall below the ALMA detection limits of existing maps or genuinely have no molecular counterpart. Cross-matching with ALMA CO maps at 90 pc in 27 galaxies shows that 41 % of PAH clouds have CO associations. The converted molecular properties vary little across environments, but the most massive clouds are preferentially found in spiral arms. Fitting lognormal mass distributions down to $2\times10^{3} M_{\odot}$ shows that spiral arms host the highest-mass clouds, consistent with enhanced formation in arm gravitational potentials. Cloud molecular surface densities decline by a factor of $\sim 1.5-2$ toward $2 - 3 R_{e}$. However, the trend largely varies in individual galaxies, with flat, decreasing, and even no trend as a function of galactocentric radius. Factors like large-scale processes and morphologies might influence the observed trends. We publish two catalogs online, one at the common resolution of 30 pc and another at the native resolution. We expect them to have broad utility for future PAH clouds, molecular clouds, and star formation studies.

preprint2025arXiv

Anisotropy in the carbon monoxide (CO) line emission across the Milky Way's disk

We present a study of the $^{12}$CO(1-0) line emission anisotropy across the Milky Way's disk to examine the effect of stellar feedback and Galactic dynamics on the distribution of the dense interstellar medium. The Hessian matrix method is used to characterize the CO line emission distribution and identify the preferential orientation across line-of-sight velocity channels in the Dame et al. 2001 composite Galactic plane survey, which covers the Galactic latitude range $|b|<5^{\circ}$. The structures sampled with this tracer are predominantly parallel to the Galactic plane toward the inner Galaxy, in clear contrast with the predominantly perpendicular orientation of the structures traced by neutral atomic hydrogen (HI) emission toward the same regions. The analysis of the Galactic plane portions sampled at higher angular resolution with other surveys reveals that the alignment with the Galactic plane is also prevalent at smaller scales. We find no preferential orientation in the CO emission toward the outer Galaxy, in contrast with the preferential alignment with the Galactic plane displayed by HI in that portion of the Milky Way. We interpret these results as the combined effect of the decrease in mid-plane pressure with increasing Galactocentric radius and SN feedback lifting diffuse gas more efficiently than dense gas off the Galactic plane.

preprint2022arXiv

A backing detector for order-keV neutrons

We have designed and tested a large-area (0.15~m$^2$) neutron detector based on neutron capture on \ce{^{6}Li}. The neutron detector design has been optimized for the purpose of tagging the scattering angle of keV-scale neutrons. These neutron detectors would be employed to calibrate the low-energy ($<$100 eV) nuclear recoil in detectors for dark matter and coherent elastic neutrino nucleus scattering (CE$ν$NS). We describe the design, construction, and characterization of a prototype. The prototype is designed to have a tagging efficiency of $\sim$25\% at the relevant $\mathcal{O}$(keV) neutron energies, and with a mean capture time of $\sim$17$~μ$s. The prototype was characterized using a \ce{^{252}Cf} neutron source and agreement with the simulation was observed within a few percent level.

preprint2022arXiv

Scintillation yield from electronic and nuclear recoils in superfluid $^4$He

Superfluid $^4$He is a promising target material for direct detection of light ($<$ 1 GeV) dark matter. Possible signal channels available for readout in this medium include prompt photons, triplet excimers, and roton and phonon quasiparticles. The relative yield of these signals has implications for the sensitivity and discrimination power of a superfluid $^4$He dark matter detector. Using a 16~cm$^3$ volume of 1.75~K superfluid $^4$He read out by six immersed photomultiplier tubes, we measured the scintillation from electronic recoils ranging between 36.3 and 185 keV$_\mathrm{ee}$, yielding a mean signal size of $1.25^{+0.03}_{-0.03}$~phe/keV$_\mathrm{ee}$, and nuclear recoils from 53.2 to 1090 keV$_\mathrm{nr}$. We compare the results of our relative scintillation yield measurements to an existing semiempirical model based on helium-helium and electron-helium interaction cross sections. We also study the behavior of delayed scintillation components as a function of recoil type and energy, a further avenue for signal discrimination in superfluid $^4$He.

preprint2022arXiv

The LUX-ZEPLIN (LZ) radioactivity and cleanliness control programs

LUX-ZEPLIN (LZ) is a second-generation direct dark matter experiment with spin-independent WIMP-nucleon scattering sensitivity above $1.4 \times 10^{-48}$ cm$^{2}$ for a WIMP mass of 40 GeV/c$^{2}$ and a 1000 d exposure. LZ achieves this sensitivity through a combination of a large 5.6 t fiducial volume, active inner and outer veto systems, and radio-pure construction using materials with inherently low radioactivity content. The LZ collaboration performed an extensive radioassay campaign over a period of six years to inform material selection for construction and provide an input to the experimental background model against which any possible signal excess may be evaluated. The campaign and its results are described in this paper. We present assays of dust and radon daughters depositing on the surface of components as well as cleanliness controls necessary to maintain background expectations through detector construction and assembly. Finally, examples from the campaign to highlight fixed contaminant radioassays for the LZ photomultiplier tubes, quality control and quality assurance procedures through fabrication, radon emanation measurements of major sub-systems, and bespoke detector systems to assay scintillator are presented.

preprint2022arXiv

The width of Herschel filaments varies with distance

Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the distance from the observer, by using previously published results from the Herschel Gould Belt Survey. Methods: We obtained updated estimates on the distances to nearby molecular clouds observed with Herschel by using recent results based on 3D dust extinction mapping and Gaia. We examined the widths of filaments from individual clouds separately, as opposed to treating them as a single population. We used these per-cloud filament widths to search for signs of variation amongst the clouds of the previously published study. Results: We find a significant dependence of the mean per-cloud filament width with distance. The distribution of mean filament widths for nearby clouds is incompatible with that of farther away clouds. The mean per-cloud widths scale with distance approximately as 4--5 times the beam size. We examine the effects of resolution by performing a convergence study of a filament profile in the Herschel image of the Taurus Molecular Cloud. We find that resolution can severely affect the shapes of radial profiles over the observed range of distances. Conclusions: We conclude that the data are inconsistent with 0.1 pc being the universal characteristic width of filaments.

preprint2020arXiv

Cloud formation in the atomic and molecular phase: HI self absorption (HISA) towards a Giant Molecular Filament

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

preprint2020arXiv

First discovery of an ultra-cool white dwarf benchmark in common proper motion with an M dwarf

Ultra-cool white dwarfs are among the oldest stellar remnants in the Universe. Their efficient gravitational settling and low effective temperatures are responsible for the smooth spectra they exhibit. For that reason, it is not possible to derive their radial velocities or to find the chemistry of the progenitors. The best that can be done is to infer such properties from associated sources, which are coeval. The simplest form of such a system is a common proper motion pair where one star is an evolved degenerate and the other a main-sequence star. In this work, we present the discovery of the first of such a system, the M dwarf LHS 6328 and the ultra-cool white dwarf PSO J1801+625, from the Pan-STARRS 1 3$π$ survey and the $Gaia$ Data Release 2. Follow-up spectra were collected covering a usable wavelength range of 3750-24500 Å. Their spectra show that the white dwarf has an effective temperature of 3550K and surface gravity of $\log g = 7.45 \pm 0.13$ or $\log g = 7.49 \pm 0.13$ for a CO or He core, respectively, when compared against synthetic spectra of ultra-cool white dwarf atmosphere models. The system has slightly subsolar metallicity with -0.25 < [Fe/H] < 0.0, and a spatial velocity of (U, V, W) = ($-114.26 \pm 0.24$, $222.94 \pm 0.60$, $10.25 \pm 0.34$) km s $^{-1}$ the first radial velocity and metallicity measurements of an ultra-cool white dwarf. This makes it the first and only benchmark of its kind to date.

preprint2020arXiv

Lowly polarized light from a highly magnetized jet of GRB 190114C

We report multi-color optical imaging and polarimetry observations of the afterglow of the first TeV- detected gamma-ray burst, GRB 190114C, using RINGO3 and MASTER II polarimeters. Observations begin 31 s after the onset of the GRB and continue until $\sim 7000\,$s post-burst. The light curves reveal a chromatic break at $\sim 400- 500\,$s, with initial temporal decay $α= 1.669 \pm 0.013$ flattening to $α\sim 1$ post-break, which we model as a combination of reverse and forward-shock components, with magnetization parameter $R_{\rm B} \sim 70$. The observed polarization degree decreases from $7.7 \pm 1.1\%$ to $2-4\%$ during $52-109\,$s post-burst and remains steady at this level for the subsequent $\sim 2000$-s, at constant position angle. Broadband spectral energy distribution modeling of the afterglow confirms GRB 190114C is highly obscured (A$_{\rm v, HG} = 1.49 \pm 0.12 \,$mag; N$_{\rm H, HG}= (9.0 \pm 0.3) \times 10^{22}\,$cm$^{-2}$). We interpret the measured afterglow polarization as intrinsically low and dominated by dust, in contrast to ${\rm P} >10\%$ measured previously for other GRB reverse shocks, with a small contribution from polarized prompt photons in the first minute. We test whether 1st and higher-order inverse Compton scattering in a magnetized reverse shock can explain the low optical polarization and the sub-TeV emission but conclude neither is explained in the reverse shock Inverse Compton model. Instead, the unexpectedly low intrinsic polarization degree in GRB 190114C can be explained if large-scale jet magnetic fields are distorted on timescales prior to reverse shock emission.

preprint2020arXiv

Measurements of Angle-Resolved Reflectivity of PTFE in Liquid Xenon with IBEX

Liquid xenon particle detectors rely on excellent light collection efficiency for their performance. This depends on the high reflectivity of polytetrafluoroethylene (PTFE) at the xenon scintillation wavelength of 178 nm, but the angular dependence of this reflectivity is not well-understood. IBEX is designed to directly measure the angular distribution of xenon scintillation light reflected off PTFE in liquid xenon. These measurements are fully described by a microphysical reflectivity model with few free parameters. Dependence on PTFE type, surface finish, xenon pressure, and wavelength of incident light is explored. Total internal reflection is observed, which results in the dominance of specular over diffuse reflection and a reflectivity near 100% for high angles of incidence.

preprint2019arXiv

Feedback in W49A diagnosed with radio recombination lines and models

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

preprint2019arXiv

Strong excess Faraday rotation on the Inside of the Sagittarius spiral arm

We present first results for Faraday rotation of compact polarized sources (1 to 2 GHz continuum) in The HI/OH/Recombination line (THOR) survey of the inner Galaxy. In the Galactic longitude range 39 degr < l < 52 degr, we find rotation measures in the range -310 rad/m2 < RM < +4219 rad/m2, with the highest values concentrated within a degree of l = 48 degrees at the Sagittarius arm tangent. Most of the high RMs arise in diffuse plasma, along lines of sight that do not intersect HII regions. For l > 49 degr, RM drops off rapidly, while at l < 47 degr, the mean RM is higher with a larger standard deviation than at l > 49 degr. We attribute the RM structure to the compressed diffuse Warm Ionized Medium in the spiral arm, upstream of the major star formation regions. The Sagittarius arm acts as a significant Faraday screen inside the Galaxy. This has implications for models of the Galactic magnetic field and the expected amount of Faraday rotation of Fast Radio Bursts from their host galaxies. We emphasize the importance of sensitivity to high Faraday depth in future polarization surveys.

preprint2019arXiv

The HI/OH/Recombination line survey of the inner Milky Way (THOR): data release 2 and HI overview

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$nα$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $τ$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

preprint2016arXiv

LOTUS: A low cost, ultraviolet spectrograph

We describe the design, construction and commissioning of LOTUS; a simple, low-cost long-slit spectrograph for the Liverpool Telescope. The design is optimized for near-UV and visible wavelengths and uses all transmitting optics. It exploits the instrument focal plane field curvature to partially correct axial chromatic aberration. A stepped slit provides narrow (2.5x95 arcsec) and wide (5x25 arcsec) options that are optimized for spectral resolution and flux calibration respectively. On sky testing shows a wavelength range of 3200-6300 Angstroms with a peak system throughput (including detector quantum efficiency) of 15 per cent and wavelength dependant spectral resolution of R=225-430. By repeated observations of the symbiotic emission line star AG Peg we demonstrate the wavelength stability of the system is less than 2 Angstroms rms and is limited by the positioning of the object in the slit. The spectrograph is now in routine operation monitoring the activity of comet 67P/Churyumov-Gerasimenko during its current post-perihelion apparition.

preprint2015arXiv

A remarkable recurrent nova in M31: Discovery and optical/UV observations of the predicted 2014 eruption

The Andromeda Galaxy recurrent nova M31N 2008-12a had been caught in eruption eight times. The inter-eruption period of M31N 2008-12a is ~1 year, making it the most rapidly recurring system known, and a strong single-degenerate Type Ia Supernova progenitor candidate. Following the 2013 eruption, a campaign was initiated to detect the predicted 2014 eruption and to then perform high cadence optical photometric and spectroscopic monitoring using ground-based telescopes, along with rapid UV and X-ray follow-up with the Swift satellite. Here we report the results of a high cadence multicolour optical monitoring campaign, the spectroscopic evolution, and the UV photometry. We also discuss tantalising evidence of a potentially related, vastly-extended, nebulosity. The 2014 eruption was discovered, before optical maximum, on October 2, 2014. We find that the optical properties of M31N 2008-12a evolve faster than all Galactic recurrent novae known, and all its eruptions show remarkable similarity both photometrically and spectroscopically. Optical spectra were obtained as early as 0.26 days post maximum, and again confirm the nova nature of the eruption. A significant deceleration of the inferred ejecta expansion velocity is observed which may be caused by interaction of the ejecta with surrounding material, possibly a red giant wind. We find a low ejected mass and low ejection velocity, which are consistent with high mass-accretion rate, high mass white dwarf, and short recurrence time models of novae. We encourage additional observations, especially around the predicted time of the next eruption, towards the end of 2015.

preprint2015arXiv

Hierarchical fragmentation and collapse signatures in a high-mass starless region

Aims: Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. Methods: We observed the massive (~800M_sun) starless gas clump IRDC18310-4 with the Plateau de Bure Interferometer (PdBI) at sub-arcsecond resolution in the 1.07mm continuum andN2H+(3-2) line emission. Results: Zooming from a single-dish low-resolution map to previous 3mm PdBI data, and now the new 1.07mm continuum observations, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70mum, we identify weak 70mum emission toward one core with a comparably low luminosity of ~16L_sun, re-enforcing the general youth of the region. The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the individual components (in the 0.3 to 1.0km/s regime). Based on single-dish C18O(2-1) data we estimate a low virial-to-gas-mass ratio <=0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear to have very low levels of internal turbulence.

preprint2015arXiv

Kinematic and Thermal Structure at the onset of high-mass star formation

We want to understand the kinematic and thermal properties of young massive gas clumps prior to and at the earliest evolutionary stages of high-mass star formation. Do we find signatures of gravitational collapse? Do we find temperature gradients in the vicinity or absence of infrared emission sources? Do we find coherent velocity structures toward the center of the dense and cold gas clumps? To determine kinematics and gas temperatures, we used ammonia, because it is known to be a good tracer and thermometer of dense gas. We observed the NH$_3$(1,1) and (2,2) lines within seven very young high-mass star-forming regions with the VLA and the Effelsberg 100m telescope. This allows us to study velocity structures, linewidths, and gas temperatures at high spatial resolution of 3-5$"$, corresponding to $\sim$0.05 pc. We find on average cold gas clumps with temperatures in the range between 10 K and 30 K. The observations do not reveal a clear correlation between infrared emission peaks and ammonia temperature peaks. We report an upper limit for the linewidth of $\sim$1.3 km s$^{-1}$, at the spectral resolution limit of our VLA observation. This indicates a relatively low level of turbulence on the scale of the observations. Velocity gradients are present in almost all regions with typical velocity differences of 1 to 2 km s$^{-1}$ and gradients of 5 to 10 km s$^{-1}$ pc$^{-1}$. These velocity gradients are smooth in most cases, but there is one exceptional source (ISOSS23053), for which we find several velocity components with a steep velocity gradient toward the clump centers that is larger than 30 km s$^{-1}$ pc$^{-1}$. This steep velocity gradient is consistent with recent models of cloud collapse. Furthermore, we report a spatial correlation of ammonia and cold dust, but we also find decreasing ammonia emission close to infrared emission sources.

preprint2015arXiv

Limits on optical polarization during the prompt phase of GRB 140430A

Gamma-ray burst GRB 140430A was detected by the Swift satellite and observed promptly with the imaging polarimeter RINGO3 mounted on the Liverpool Telescope, with observations beginning while the prompt $γ$-ray emission was still ongoing. In this paper, we present densely sampled (10-second temporal resolution) early optical light curves in 3 optical bands and limits to the degree of optical polarization. We compare optical, X-ray and gamma-ray properties and present an analysis of the optical emission during a period of high-energy flaring. The complex optical light curve cannot be explained merely with a combination of forward and reverse shock emission from a standard external shock, implying additional contribution of emission from internal shock dissipation. We estimate an upper limit for time averaged optical polarization during the prompt phase to be as low as P < 12% (1$σ$). This suggests that the optical flares and early afterglow emission in this GRB are not highly polarized. Alternatively, time averaging could mask the presence of otherwise polarized components of distinct origin at different polarization position angles.

preprint2015arXiv

Liverpool Telescope 2: a new robotic facility for rapid transient follow-up

The Liverpool Telescope is one of the world's premier facilities for time domain astronomy. The time domain landscape is set to radically change in the coming decade, with surveys such as LSST providing huge numbers of transient detections on a nightly basis; transient detections across the electromagnetic spectrum from other facilities such as SVOM, SKA and CTA; and the era of `multi-messenger astronomy', wherein events are detected via non-electromagnetic means, such as gravitational wave emission. We describe here our plans for Liverpool Telescope 2: a new robotic telescope designed to capitalise on this new era of time domain astronomy. LT2 will be a 4-metre class facility co-located with the LT at the Observatorio del Roque de Los Muchachos on the Canary island of La Palma. The telescope will be designed for extremely rapid response: the aim is that the telescope will take data within 30 seconds of the receipt of a trigger from another facility. The motivation for this is twofold: firstly it will make it a world-leading facility for the study of fast fading transients and explosive phenomena discovered at early times. Secondly, it will enable large-scale programmes of low-to-intermediate resolution spectral classification of transients to be performed with great efficiency. In the target-rich environment of the LSST era, minimising acquisition overheads will be key to maximising the science gains from any follow-up programme. The telescope will have a diverse instrument suite which is simultaneously mounted for automatic changes, but it is envisaged that the primary instrument will be an intermediate resolution, optical/infrared spectrograph for scientific exploitation of transients discovered with the next generation of synoptic survey facilities. In this paper we outline the core science drivers for the telescope, and the requirements for the optical and mechanical design.

preprint2015arXiv

THOR - The HI, OH, Recombination Line Survey of the Milky Way - The pilot study: HI observations of the giant molecular cloud W43

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

preprint2014arXiv

A broadband silicon quarter-wave retarder for far-infrared spectroscopic circular dichroism

The high brightness, broad spectral coverage and pulsed characteristics of infrared synchrotron radiation enable time-resolved spectroscopy under throughput-limited optical systems, as can occur with the high-field magnet cryostat systems used to study electron dynamics and cyclotron resonance by far-infrared techniques. A natural extension for magnetospectroscopy is to sense circular dichroism, i.e. the difference in a material's optical response for left and right circularly polarized light. A key component for spectroscopic circular dichroism is an achromatic 1/4 wave retarder functioning over the spectral range of interest. We report here the development of an in-line retarder using total internal reflection in high-resistivity silicon. We demonstrate its performance by distinguishing electronic excitations of different handednesses for GaAs in a magnetic field. This 1/4 wave retarder is expected to be useful for far-infrared spectroscopy of circular dichroism in many materials.

preprint2014arXiv

A surprising consistency between the far-infrared galaxy luminosity functions of the field and Coma

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.

preprint2014arXiv

Liverpool Telescope 2: a new robotic facility for time domain astronomy in 2020+

The robotic 2m Liverpool Telescope, based on the Canary island of La Palma, has a diverse instrument suite and a strong track record in time domain science, with highlights including early time photometry and spectra of supernovae, measurements of the polarization of gamma-ray burst afterglows, and high cadence light curves of transiting extrasolar planets. In the next decade the time domain will become an increasingly prominent part of the astronomical agenda with new facilities such as LSST, SKA, CTA and Gaia, and promised detections of astrophysical gravitational wave and neutrino sources opening new windows on the transient universe. To capitalise on this exciting new era we intend to build Liverpool Telescope 2: a new robotic facility on La Palma dedicated to time domain science. The next generation of survey facilities will discover large numbers of new transient sources, but there will be a pressing need for follow-up observations for scientific exploitation, in particular spectroscopic follow-up. Liverpool Telescope 2 will have a 4-metre aperture, enabling optical/infrared spectroscopy of faint objects. Robotic telescopes are capable of rapid reaction to unpredictable phenomena, and for fast-fading transients like gamma-ray burst afterglows. This rapid reaction enables observations which would be impossible on less agile telescopes of much larger aperture. We intend Liverpool Telescope 2 to have a world-leading response time, with the aim that we will be taking data with a few tens of seconds of receipt of a trigger from a ground- or space-based transient detection facility. We outline here our scientific goals and present the results of our preliminary optical design studies.

preprint2014arXiv

Measurement of the intrinsic electron neutrino component in the T2K neutrino beam with the ND280 detector

The T2K experiment has reported the first observation of the appearance of electron neutrinos in a muon neutrino beam. The main and irreducible background to the appearance signal comes from the presence in the neutrino beam of a small intrinsic component of electron neutrinos originating from muon and kaon decays. In T2K, this component is expected to represent 1.2% of the total neutrino flux. A measurement of this component using the near detector (ND280), located 280 m from the target, is presented. The charged current interactions of electron neutrinos are selected by combining the particle identification capabilities of both the time projection chambers and electromagnetic calorimeters of ND280. The measured ratio between the observed electron neutrino beam component and the prediction is 1.01+-0.10 providing a direct confirmation of the neutrino fluxes and neutrino cross section modeling used for T2K neutrino oscillation analyses. Electron neutrinos coming from muons and kaons decay are also separately measured, resulting in a ratio with respect to the prediction of 0.68+-0.30 and 1.10+-0.14, respectively.

preprint2014arXiv

Measurement of the neutrino-oxygen neutral-current interaction cross section by observing nuclear deexcitation $γ$ rays

We report the first measurement of the neutrino-oxygen neutral-current quasielastic (NCQE) cross section. It is obtained by observing nuclear deexcitation $γ$-rays which follow neutrino-oxygen interactions at the Super-Kamiokande water Cherenkov detector. We use T2K data corresponding to $3.01 \times 10^{20}$ protons on target. By selecting only events during the T2K beam window and with well-reconstructed vertices in the fiducial volume, the large background rate from natural radioactivity is dramatically reduced. We observe 43 events in the $4-30$ MeV reconstructed energy window, compared with an expectation of 51.0, which includes an estimated 16.2 background events. The background is primarily nonquasielastic neutral-current interactions and has only 1.2 events from natural radioactivity. The flux-averaged NCQE cross section we measure is $1.55 \times 10^{-38}$ cm$^2$ with a 68\% confidence interval of $(1.22, 2.20) \times 10^{-38}$ cm$^2$ at a median neutrino energy of 630 MeV, compared with the theoretical prediction of $2.01 \times 10^{-38}$ cm$^2$.

preprint2014arXiv

Precise Measurement of the Neutrino Mixing Parameter θ_{23} from Muon Neutrino Disappearance in an Off-axis Beam

New data from the T2K neutrino oscillation experiment produce the most precise measurement of the neutrino mixing parameter theta_{23}. Using an off-axis neutrino beam with a peak energy of 0.6 GeV and a data set corresponding to 6.57 x 10^{20} protons on target, T2K has fit the energy-dependent nu_mu oscillation probability to determine oscillation parameters. Marginalizing over the values of other oscillation parameters yields sin^2 (theta_{23}) = 0.514 +0.055/-0.056 (0.511 +- 0.055), assuming normal (inverted) mass hierarchy. The best-fit mass-squared splitting for normal hierarchy is Delta m^2_{32} = (2.51 +- 0.10) x 10^{-3} eV^2/c^4 (inverted hierarchy: Delta m^2_{13} = (2.48 +- 0.10) x 10^{-3} eV^2/c^4). Adding a model of multinucleon interactions that affect neutrino energy reconstruction is found to produce only small biases in neutrino oscillation parameter extraction at current levels of statistical uncertainty.

preprint2014arXiv

The nature of the late achromatic bump in GRB 120326A

The long ${\it Swift}$ gamma-ray burst GRB 120326A at redshift $z=1.798$ exhibited a multi-band light curve with a striking feature: a late-time, long-lasting achromatic rebrightening, rarely seen in such events. Peaking in optical and X-ray bands $\sim 35$ ks ($\sim 12.5$ ks in the GRB rest frame) after the 70-s GRB prompt burst, the feature brightens nearly two orders of magnitude above the underlying optical power-law decay. Modelling the multiwavelength light curves, we investigate possible causes of the rebrightening in the context of the standard fireball model. We exclude a range of scenarios for the origin of this feature: reverse-shock flash, late-time forward shock peak due to the passage of the maximal synchrotron frequency through the optical band, late central engine optical/X-ray flares, interaction between the expanding blast wave and a density enhancement in the circumburst medium and gravitational microlensing. Instead we conclude that the achromatic rebrightening may be caused by a refreshed forward shock or a geometrical effect. In addition, we identify an additional component after the end of the prompt emission, that shapes the observed X-ray and optical light curves differently, ruling out a single overall emission component to explain the observed early time emission.

preprint2013arXiv

GRB 090727 and gamma-ray bursts with early time optical emission

We present a multi-wavelength analysis of gamma-ray burst GRB 090727, for which optical emission was detected during the prompt gamma-ray emission by the 2-m autonomous robotic Liverpool Telescope and subsequently monitored for a further two days with the Liverpool and Faulkes telescopes. Within the context of the standard fireball model, we rule out a reverse shock origin for the early time optical emission in GRB 090727 and instead conclude that the early time optical flash likely corresponds to emission from an internal dissipation processes. Putting GRB 090727 into a broader observational and theoretical context, we build a sample of 36 gamma-ray bursts (GRBs) with contemporaneous early time optical and gamma-ray detections. From these GRBs, we extract a sub-sample of 18 GRBs, which show optical peaks during prompt gamma-ray emission, and perform detailed temporal and spectral analysis in gamma-ray, X-ray, and optical bands. We find that in most cases early time optical emission shows sharp and steep behaviour, and notice a rich diversity of spectral properties. Using a simple internal shock dissipation model, we show that the emission during prompt GRB phase can occur at very different frequencies via synchrotron radiation. Based on the results obtained from observations and simulation, we conclude that the standard external shock interpretation for early time optical emission is disfavoured in most cases due to sharp peaks ($Δt/t < 1$) and steep rise/decay indices, and that internal dissipation can explain the properties of GRBs with optical peaks during gamma-ray emission.

preprint2013arXiv

Measurement of Neutrino Oscillation Parameters from Muon Neutrino Disappearance with an Off-axis Beam

The T2K collaboration reports a precision measurement of muon neutrino disappearance with an off-axis neutrino beam with a peak energy of 0.6 GeV. Near detector measurements are used to constrain the neutrino flux and cross section parameters. The Super-Kamiokande far detector, which is 295 km downstream of the neutrino production target, collected data corresponding to $3.01 \times 10^{20}$ protons on target. In the absence of neutrino oscillations, $205 \pm 17$ (syst.) events are expected to be detected and only 58 muon neutrino event candidates are observed. A fit to the neutrino rate and energy spectrum assuming three neutrino flavors, normal mass hierarchy and $θ_{23}\leq π/4$ yields a best-fit mixing angle $\sin^2(2θ_{23})=1.000$ and mass splitting $|Δm^2_{32}| =2.44 \times 10^{-3}$ eV$^2$/c$^4$. If $θ_{23}\geq π/4$ is assumed, the best-fit mixing angle changes to $\sin^2(2θ_{23})=0.999$ and the mass splitting remains unchanged.

preprint2013arXiv

Measurement of the Inclusive NuMu Charged Current Cross Section on Carbon in the Near Detector of the T2K Experiment

T2K has performed the first measurement of νμ inclusive charged current interactions on carbon at neutrino energies of ~1 GeV where the measurement is reported as a flux-averaged double differential cross section in muon momentum and angle. The flux is predicted by the beam Monte Carlo and external data, including the results from the NA61/SHINE experiment. The data used for this measurement were taken in 2010 and 2011, with a total of 10.8 x 10^{19} protons-on-target. The analysis is performed on 4485 inclusive charged current interaction candidates selected in the most upstream fine-grained scintillator detector of the near detector. The flux-averaged total cross section is <σ_CC>_ϕ=(6.91 +/- 0.13 (stat) +/- 0.84 (syst)) x10^{-39} cm^2/nucleon for a mean neutrino energy of 0.85 GeV.

preprint2013arXiv

New constraints on GRB jet geometry and relativistic shock physics

We use high--quality, multi-band observations of Swift GRB120404A, from gamma-ray to radio frequencies, together with the new hydrodynamics code of van Eerten et al. (2012) to test the standard synchrotron shock model. The evolution of the radio and optical afterglow, with its prominent optical rebrightening at t_rest 260-2600 s, is remarkably well modelled by a decelerating jet viewed close to the jet edge, combined with some early re-energization of the shock. We thus constrain the geometry of the jet with half-opening and viewing angles of 23 and 21 deg respectively and suggest that wide jets viewed off-axis are more common in GRBs than previously thought. We also derive the fireball microphysics parameters epsilon_B=2.4e-4 and epsilon_e=9.3e-2 and a circumburst density of n=240 cm^-3. The ability to self-consistently model the microphysics parameters and jet geometry in this way offers an alternative to trying to identify elusive canonical jet breaks at late times. The mismatch between the observed and model-predicted X-ray fluxes is explained by the local rather than the global cooling approximation in the synchrotron radiation model, constraining the microphysics of particle acceleration taking place in a relativistic shock and, in turn, emphasising the need for a more realistic treatment of cooling in future developments of theoretical models. Finally, our interpretation of the optical peak as due to the passage of the forward shock synchrotron frequency highlights the importance of high quality multi-band data to prevent some optical peaks from being erroneously attributed to the onset of fireball deceleration.

preprint2013arXiv

Polyhedrality in Pieces

The aim of this paper is to present two tools, Theorems 4 and 7, that make the task of finding equivalent polyhedral norms on certain Banach spaces easier and more transparent. The hypotheses of both tools are based on countable decompositions, either of the unit sphere S_X or of certain subsets of the dual ball of a given Banach space X. The sufficient conditions of Theorem 4 are shown to be necessary in the separable case. Using Theorem 7, we can unify two known results regarding the polyhedral renorming of certain C(K) spaces, and spaces having an (uncountable) unconditional basis. New examples of spaces having equivalent polyhedral norms are given in the final section.

preprint2013arXiv

STILT: System Design & Performance

The Small Telescopes Installed at the Liverpool Telescope (STILT) have been in operation since March 2009, collecting wide field data from their position, mounted to the Liverpool Telescope. The two instruments; SkycamT and SkycamZ have been used to create a variability search of the skies visible at La Palma with the limits of 12th and 18th R band magnitude with fields of view of 21x21 and 1x1 degrees. We provide here a description of the hardware and software setup and the performance of the system to date.

preprint2013arXiv

The T2K Neutrino Flux Prediction

The Tokai-to-Kamioka (T2K) experiment studies neutrino oscillations using an off-axis muon neutrino beam with a peak energy of about 0.6 GeV that originates at the J-PARC accelerator facility. Interactions of the neutrinos are observed at near detectors placed at 280 m from the production target and at the far detector -- Super-Kamiokande (SK) -- located 295 km away. The flux prediction is an essential part of the successful prediction of neutrino interaction rates at the T2K detectors and is an important input to T2K neutrino oscillation and cross section measurements. A FLUKA and GEANT3 based simulation models the physical processes involved in the neutrino production, from the interaction of primary beam protons in the T2K target, to the decay of hadrons and muons that produce neutrinos. The simulation uses proton beam monitor measurements as inputs. The modeling of hadronic interactions is re-weighted using thin target hadron production data, including recent charged pion and kaon measurements from the NA61/SHINE experiment. For the first T2K analyses the uncertainties on the flux prediction are evaluated to be below 15% near the flux peak. The uncertainty on the ratio of the flux predictions at the far and near detectors is less than 2% near the flux peak.

preprint2012arXiv

ACCESS - V. Dissecting ram-pressure stripping through integral-field spectroscopy and multi-band imaging

We study the case of a bright (L>L*) barred spiral galaxy from the rich cluster A3558 in the Shapley supercluster core (z=0.05) undergoing ram-pressure stripping. Integral-field spectroscopy, complemented by multi-band imaging, allows us to reveal the impact of ram pressure on the interstellar medium. We study in detail the kinematics and the physical conditions of the ionized gas and the properties of the stellar populations. We observe one-sided extraplanar ionized gas along the full extent of the galaxy disc. Narrow-band Halpha imaging resolves this outflow into a complex of knots and filaments. The gas velocity field is complex with the extraplanar gas showing signature of rotation. In all parts of the galaxy, we find a significant contribution from shock excitation, as well as emission powered by star formation. Shock-ionized gas is associated with the turbulent gas outflow and highly attenuated by dust. All these findings cover the whole phenomenology of early-stage ram-pressure stripping. Intense, highly obscured star formation is taking place in the nucleus, probably related to the bar, and in a region 12 kpc South-West from the centre. In the SW region we identify a starburst characterized by a 5x increase in the star-formation rate over the last ~100 Myr, possibly related to the compression of the interstellar gas by the ram pressure. The scenario suggested by the observations is supported and refined by ad hoc N-body/hydrodynamical simulations which identify a rather narrow temporal range for the onset of ram-pressure stripping around t~60 Myr ago, and an angle between the galaxy rotation axis and the intra-cluster medium wind of ~45 deg. Taking into account that the galaxy is found ~1 Mpc from the cluster centre in a relatively low-density region, this study shows that ram-pressure stripping still acts efficiently on massive galaxies well outside the cluster cores.

preprint2012arXiv

First Muon-Neutrino Disappearance Study with an Off-Axis Beam

We report a measurement of muon-neutrino disappearance in the T2K experiment. The 295-km muon-neutrino beam from Tokai to Kamioka is the first implementation of the off-axis technique in a long-baseline neutrino oscillation experiment. With data corresponding to 1.43 10**20 protons on target, we observe 31 fully-contained single muon-like ring events in Super-Kamiokande, compared with an expectation of 104 +- 14 (syst) events without neutrino oscillations. The best-fit point for two-flavor nu_mu -> nu_tau oscillations is sin**2(2 theta_23) = 0.98 and |Δm**2_32| = 2.65 10**-3 eV**2. The boundary of the 90 % confidence region includes the points (sin**2(2 theta_23),|Δm**2_32|) = (1.0, 3.1 10**-3 eV**2), (0.84, 2.65 10**-3 eV**2) and (1.0, 2.2 10**-3 eV**2).

preprint2011arXiv

A faint optical flash in dust-obscured GRB 080603A - implications for GRB prompt emission mechanisms

We report the detection of a faint optical flash by the 2-m Faulkes Telescope North simultaneously with the second of two prompt gamma-ray pulses in INTEGRAL gamma-ray burst (GRB) 080603A, beginning at t_rest = 37 s after the onset of the GRB. This optical flash appears to be distinct from the subsequent emerging afterglow emission, for which we present comprehensive broadband radio to X-ray light curves to 13 days post-burst and rigorously test the standard fireball model. The intrinsic extinction toward GRB 080603A is high (A_V,z = 0.8 mag), and the well-sampled X-ray-to-near-infrared spectral energy distribution is interesting in requiring an LMC2 extinction profile, in contrast to the majority of GRBs. Comparison of the gamma-ray and extinction-corrected optical flux densities of the flash rules out an inverse-Compton origin for the prompt gamma-rays; instead, we suggest that the optical flash could originate from the inhomogeneity of the relativistic flow. In this scenario, a large velocity irregularity in the flow produces the prompt gamma-rays, followed by a milder internal shock at a larger radius that would cause the optical flash. Flat gamma-ray spectra, roughly F propto nu^-0.1, are observed in many GRBs. If the flat spectrum extends down to the optical band in GRB 080603A, the optical flare could be explained as the low-energy tail of the gamma-ray emission. If this is indeed the case, it provides an important clue to understanding the nature of the emission process in the prompt phase of GRBs and highlights the importance of deep (R> 20 mag), rapid follow-up observations capable of detecting faint, prompt optical emission.

preprint2011arXiv

A fully automated data reduction pipeline for the FRODOSpec integral field spectrograph

A fully autonomous data reduction pipeline has been developed for FRODOSpec, an optical fibre-fed integral field spectrograph currently in use at the Liverpool Telescope. This paper details the process required for the reduction of data taken using an integral field spectrograph and presents an overview of the computational methods implemented to create the pipeline. Analysis of errors and possible future enhancements are also discussed.

preprint2011arXiv

ACCESS IV: The quenching of star formation in a cluster population of dusty S0s

We present an analysis of the mid-infrared (MIR) colours of 165 70um-detected galaxies in the Shapley supercluster core (SSC) at z=0.048 using panoramic Spitzer/MIPS 24 and 70um imaging. While the bulk of galaxies show f70/f24 colours typical of local star-forming galaxies, we identify a significant sub-population of 23 70micron-excess galaxies, whose MIR colours (f70/f24>25) are much redder and cannot be reproduced by any of the standard model infrared SEDs. These galaxies are found to be strongly concentrated towards the cores of the five clusters that make up the SSC, and also appear rare among local field galaxies, confirming them as a cluster-specific phenomenon. Their optical spectra and lack of significant UV emission imply little or no ongoing star formation, while fits to their panchromatic SEDs require the far-IR emission to come mostly from a diffuse dust component heated by the general interstellar radiation field rather than ongoing star formation. Most of these 70micron-excess galaxies are identified as ~L* S0s with smooth profiles. We find that almost every cluster galaxy in the process of star-formation quenching is already either an S0 or Sa, while we find no passive galaxies of class Sb or later. Hence the formation of passive early-type galaxies in cluster cores must involve the prior morphological transformation of late-type spirals into Sa/S0s, perhaps via pre-processing or the impact of cluster tidal fields, before a subsequent quenching of star formation once the lenticular encounters the dense environment of the cluster core. In the cases of many cluster S0s, this phase of star-formation quenching is characterised by an excess of 70um emission, indicating that the cold dust content is declining at a slower rate than star formation.

preprint2011arXiv

Constraining GRB Emission Physics with Extensive Early-Time, Multiband Follow-up

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

preprint2011arXiv

Indication of Electron Neutrino Appearance from an Accelerator-produced Off-axis Muon Neutrino Beam

The T2K experiment observes indications of $ν_μ\rightarrow ν_e$ appearance in data accumulated with $1.43\times10^{20}$ protons on target. Six events pass all selection criteria at the far detector. In a three-flavor neutrino oscillation scenario with $|Δm_{23}^2|=2.4\times10^{-3}$ eV$^2$, $\sin^2 2θ_{23}=1$ and $\sin^2 2θ_{13}=0$, the expected number of such events is 1.5$\pm$0.3(syst.). Under this hypothesis, the probability to observe six or more candidate events is 7$\times10^{-3}$, equivalent to 2.5$σ$ significance. At 90% C.L., the data are consistent with 0.03(0.04)$<\sin^2 2θ_{13}<$ 0.28(0.34) for $δ_{\rm CP}=0$ and a normal (inverted) hierarchy.

preprint2011arXiv

Measurements of the T2K neutrino beam properties using the INGRID on-axis near detector

Precise measurement of neutrino beam direction and intensity was achieved based on a new concept with modularized neutrino detectors. INGRID (Interactive Neutrino GRID) is an on-axis near detector for the T2K long baseline neutrino oscillation experiment. INGRID consists of 16 identical modules arranged in horizontal and vertical arrays around the beam center. The module has a sandwich structure of iron target plates and scintillator trackers. INGRID directly monitors the muon neutrino beam profile center and intensity using the number of observed neutrino events in each module. The neutrino beam direction is measured with accuracy better than 0.4 mrad from the measured profile center. The normalized event rate is measured with 4% precision.

preprint2011arXiv

XRF 100316D / SN 2010bh and the Nature of Gamma Ray Burst Supernovae

We present ground-based and HST optical and infrared observations of XRF 100316D / SN 2010bh. It is seen that the optical light curves of SN 2010bh evolve at a faster rate than the archetype GRB-SN 1998bw, but at a similar rate to SN 2006aj, a supernova that was spectroscopically linked with XRF 060218, and at a similar rate to non-GRB associated type Ic SN 1994I. We estimate the rest-frame extinction of this event from our optical data to be E(B-V)=0.18 +/- 0.08 mag. We find the V-band absolute magnitude of SN 2010bh to be M_{V}=-18.62 +/- 0.08, which is the faintest peak V-band magnitude observed to-date for a spectroscopically-confirmed GRB-SNe. When we investigate the origin of the flux at t-t_{o}=0.598 days, it is shown that the light is not synchrotron in origin, but is likely coming from the supernova shock break-out. We then use our optical and infrared data to create a quasi-bolometric light curve of SN 2010bh which we model with a simple analytical formula. The results of our modeling imply that SN 2010bh synthesized a nickel mass of M_{Ni} \approx 0.10 M_{sun}, ejected M_{ej} \approx 2.2 M_{sun} and has an explosion energy of E_{k} \approx 1.4 x 10^{52} erg. Finally, for a sample 22 GRB-SNe we check for a correlation between the stretch factors and luminosity factors in the R band and conclude that no statistically-significant correlation exists.

preprint2010arXiv

A tale of two GRB-SNe at a common redshift of z = 0.54

We present ground-based and HST optical observations of the optical transients (OTs) of long-duration Gamma Ray Bursts (GRBs) 060729 and 090618, both at a redshift of z = 0.54. For GRB 060729, bumps are seen in the optical light curves (LCs), and the late-time broadband spectral energy distributions (SEDs) of the OT resemble those of local type Ic supernovae (SNe). For GRB 090618, the dense sampling of our optical observations has allowed us to detect well-defined bumps in the optical LCs, as well as a change in colour, that are indicative of light coming from a core-collapse SN. The accompanying SNe for both events are individually compared with SN1998bw, a known GRB-supernova, and SN1994I, a typical type Ic supernova without a known GRB counterpart, and in both cases the brightness and temporal evolution more closely resemble SN1998bw. We also exploit our extensive optical and radio data for GRB 090618, as well as the publicly-available SWIFT -XRT data, and discuss the properties of the afterglow at early times. In the context of a simple jet-like model, the afterglow of GRB 090618 is best explained by the presence of a jet-break at t-to > 0.5 days. We then compare the rest-frame, peak V -band absolute magnitudes of all of the GRB and X-Ray Flash (XRF)-associated SNe with a large sample of local type Ibc SNe, concluding that, when host extinction is considered, the peak magnitudes of the GRB/XRF-SNe cannot be distinguished from the peak magnitudes of non-GRB/XRF SNe.

preprint2010arXiv

GRB 090313 and the Origin of Optical Peaks in GRB Light Curves: Implications for Lorentz Factors and Radio Flares

We use a sample of 19 Gamma Ray Bursts (GRBs) that exhibit single-peaked optical light curves to test the standard fireball model by investigating the relationship between the time of the onset of the afterglow and the temporal rising index. Our sample includes GRBs and X-ray flashes for which we derive a wide range of initial Lorentz factors ($40 < Γ< 450$). Using plausible model parameters the typical frequency of the forward shock is expected to lie close to the optical band; within this low typical frequency framework, we use the optical data to constrain $ε_e$ and show that values derived from the early time light curve properties are consistent with published typical values derived from other afterglow studies. We produce expected radio light curves by predicting the temporal evolution of the expected radio emission from forward and reverse shock components, including synchrotron self-absorption effects at early time. Although a number of the GRBs in this sample do not have published radio measurements, we demonstrate the effectiveness of this method in the case of {\it Swift} GRB 090313, for which millimetric and centrimetric observations were available, and conclude that future detections of reverse-shock radio flares with new radio facilities such as the EVLA and ALMA will test the low frequency model and provide constraints on magnetic models.

preprint2010arXiv

GRB 090902B: afterglow observations and implications

The optical-infrared afterglow of the LAT-detected long duration burst, GRB 090902B, has been observed by several instruments. The earliest detection by ROTSE-IIIa occurred 80 minutes after detection by the GBM instrument onboard the Fermi Gamma-Ray Space Telescope, revealing a bright afterglow and a decay slope suggestive of a reverse shock origin. Subsequent optical-IR observations followed the light curve for 6.5 days. The temporal and spectral behavior at optical-infrared frequencies is consistent with synchrotron fireball model predictions; the cooling break lies between optical and XRT frequencies ~ 1.9 days after the burst. The inferred electron energy index is $p = 1.8 \pm 0.2$, which would however imply an X-ray decay slope flatter than observed. The XRT and LAT data have similar spectral indices and the observed steeper value of the LAT temporal index is marginally consistent with the predicted temporal decay in the radiative regime of the forward shock model. Absence of a jet break during the first 6 days implies a collimation-corrected $γ$-ray energy $E_γ > 2.2\times10^{52}\rm$ ergs, one of the highest ever seen in a long-duration GRBs. More events combining GeV photon emission with multi-wavelength observations will be required to constrain the nature of the central engine powering these energetic explosions and to explore the correlations between energetic quanta and afterglow emission.

preprint2010arXiv

Ten per cent polarized optical emission from GRB 090102

The nature of the jets and the role of magnetic fields in gamma-ray bursts (GRB) remains unclear. In a baryon-dominated jet only weak, tangled fields generated in situ through shocks would be present. In an alternative model, jets are threaded with large scale magnetic fields that originate at the central engine and which accelerate and collimate the jets. The way to distinguish between the models is to measure the degree of polarization in early-time emission, however previous claims of gamma-ray polarization have been controversial. Here we report that the early optical emission from GRB 090102 was polarized at the level of P=10+/-1%, indicating the presence of large-scale fields originating in the expanding fireball. If the degree of polarization and its position angle were variable on timescales shorter than our 60-s exposure, then the peak polarization may have been larger than 10 per cent.

preprint2010arXiv

UV-IR luminosity functions and stellar mass functions of galaxies in the Shapley supercluster core

We present a panchromatic study of luminosity functions (LFs) and stellar mass functions (SMFs) of galaxies in the core of the Shapley supercluster at z=0.048, in order to investigate how the dense environment affects the galaxy properties, such as star formation (SF) or stellar masses. We find that while faint-end slopes of optical and NIR LFs steepen with decreasing density, no environment effect is found in the slope of the SMFs. This suggests that mechanisms transforming galaxies in different environments are mainly related to the quench of SF rather than to mass-loss. The Near-UV (NUV) and Far-UV (FUV) LFs obtained have steeper faint-end slopes than the local field population, while the 24$μ$m and 70$μ$m galaxy LFs for the Shapley supercluster have shapes fully consistent with those obtained for the local field galaxy population. This apparent lack of environmental dependence for the infrared (IR) LFs suggests that the bulk of the star-forming galaxies that make up the observed cluster IR LF have been recently accreted from the field and have yet to have their SF activity significantly affected by the cluster environment.

preprint2009arXiv

ACCESS: NIR Luminosity Function and Stellar Mass Function of Galaxies in the Shapley Supercluster Environment

We present the NIR luminosity (LF) and stellar mass functions (SMF) of galaxies in the core of the Shapley supercluster at z=0.048, based on new K-band observations in conjunction with B- and R-band photometry and a subsample of ~650 galaxies spectroscopically confirmed supercluster members, allowing to investigate the galaxies down to M_K^*+6 and M=10^8.75 M_sun. For the 3 deg^2 field the K-band LF is described by a Schechter function with M_K^*=-24.96+-0.10 and α=-1.42+-0.03. We investigate the effect of environment by deriving the LF in three regions selected according to the local galaxy density, and observe a significant increase in the faint-end slope going from the high- (α=-1.33) to the low-density (α=-1.49) environments. The SMF is fitted well by a Schechter function with log_10(M^*)=11.16+-0.04 and α=-1.20+-0.02. The SMF of supercluster galaxies is also characterised by an excess of massive galaxies that are associated to the cluster BCGs. While the value of M* depends on environment the slope of the galaxy SMF does not vary with the environment. By comparing our findings with cosmological simulations, we conclude that the environmental dependences of the LF are not primary due to variations in the merging histories, but to processes which are not treated in the semi-analytical models, such as tidal stripping or harassment. In field regions the SMF shows a sharp upturn below M=10^9 M_sun, close to our mass limit, suggesting that the upturns seen in our K-band LFs, but not in the SMF, are due to this dwarf population. The environmental variations seen in the faint-end of the K-band LF suggests that dwarf galaxies, which are easier to strip than their more massive counterparts, are affected by tidal/gas stripping upon entering the supercluster.

preprint2009arXiv

LoCuSS: The mid-infrared Butcher-Oemler effect

We study the mid-infrared (MIR) properties of galaxies in 30 massive galaxy clusters at 0.02<z<0.40, using panoramic Spitzer/MIPS 24micron and NIR data. This is the largest sample of clusters to date with MIR data covering not only the cluster cores, but extending into the infall regions. We revisit the Butcher-Oemler effect, measuring the fraction of massive infrared-luminous galaxies (K<K*+1.5, L_IR>5x10^10L_sun) within r_200, finding a steady increase in the fraction with redshift from ~3% at z=0.02 to ~10% by z=0.30, and an rms cluster-to-cluster scatter about this trend of 0.03. The best-fit redshift evolution model is of the form f_SF ~ (1+z)^5.7, which is stronger redshift evolution than that of L*_IR in both clusters and the field. We find that, statistically, this excess is associated with galaxies found at large cluster-centric radii, implying that the MIR Butcher-Oemler effect can be explained by a combination of both the global decline in star-formation in the universe since z~1 and enhanced star formation in the infall regions of clusters at intermediate redshifts. This picture is supported by a simple infall model based on the Millennium Simulation semi-analytic galaxy catalogs, whereby star-formation in infalling galaxies is instantaneously quenched upon their first passage through the cluster, in that the observed radial trends of f_SF trace those inferred from the simulations. We also find that f_SF does not depend on simple indicators of the dynamical state of clusters, including the offset between the brightest cluster galaxy and the peak of the X-ray emission. This is consistent with the picture described above in that most new star-formation in clusters occurs in the infall regions, and is thus not sensitive to the details of cluster-cluster mergers in the core regions.

preprint2009arXiv

On the origin of the scatter around the Fundamental Plane: correlations with stellar population parameters

We present a fundamental plane (FP) analysis of 141 early-type galaxies in the Shapley supercluster at z=0.049 based on spectroscopy from the AAOmega spectrograph at the AAT and photometry from the WFI on the ESO/MPI 2.2m telescope. The key feature of the survey is its coverage of low-mass galaxies down to sigma_0~50km/s. We obtain a best-fitting FP relation log r_e=1.06 log sigma_0 - 0.82 log < I >_e + cst in the R band. The shallow exponent of sigma_0 is a result of the extension of our sample to low velocity dispersions. We investigate the origin of the intrinsic FP scatter, using estimates of age, metallicity and alpha/Fe. We find that the FP residuals anti-correlate (>3sigma) with the mean stellar age in agreement with previous work. However, a stronger (>4sigma) correlation with alpha/Fe is also found. These correlations indicate that galaxies with effective radii smaller than those predicted by the FP have stellar populations systematically older and with alpha over-abundances larger than average, for their sigma_0. Including alpha/Fe as a fourth parameter in the FP, the total scatter decreases from 0.088 dex to 0.075 dex and the estimated intrinsic scatter decreases from 0.068 dex to 0.049 dex. Thus, variations in alpha/Fe account for >> 30% of the total variance around the FP, and >> 50% of the estimated intrinsic variance. This result indicates that the distribution of galaxies around the FP are tightly related to the enrichment, and hence to the timescale of star-formation. Our results appear to be consistent with the merger hypothesis for the formation of ellipticals which predicts that a significant fraction of the scatter is due to variations in the importance of dissipation in forming merger remnants of a given mass.

preprint2009arXiv

The HST/ACS Coma Cluster Survey: V - Compact Stellar Systems in the Coma Cluster

The HST ACS Coma Cluster Treasury Survey is a deep two passband imaging survey of the nearest very rich cluster of galaxies, covering a range of galaxy density environments. The imaging is complemented by a recent wide field redshift survey of the cluster conducted with Hectospec on the 6.5m MMT. Among the many scientific applications for this data are the search for compact galaxies. In this paper, we present the discovery of seven compact (but quite luminous) stellar systems, ranging from M32-like galaxies down to ultra-compact dwarfs (UCDs)/dwarf to globular transition objects (DGTOs). We find that all seven compact galaxies require a two-component fit to their light profile and have measured velocity dispersions that exceed those expected for typical early-type galaxies at their luminosity. From our structural parameter analysis we conclude that three of the sample should be classified as compact ellipticals or M32-like galaxies, the remaining four being less extreme systems. The three compact ellipticals are all found to have old luminosity weighted ages (> 12 Gyr), intermediate metallicities (-0.6 < [Fe/H] < -0.1) and high [Mg/Fe] (> 0.25). Our findings support a tidal stripping scenario as the formation mode of compact galaxies covering the luminosity range studied here. We speculate that at least two early-type morphologies may serve as the progenitor of compact galaxies in clusters.

preprint2006arXiv

A Review of Early-Time Optical Follow-ups with 2-m Robotic Telescopes

We summarise recent deep, rapid GRB follow-up observations using the RoboNet-1.0 network which comprises three fully-robotic 2-m telescopes, the Liverpool Telescope and the Faulkes Telescopes North and South. Observations begin automatically within minutes of receipt of a GRB alert and may continue for hours or days to provide well-sampled multi-colour light curves or deep upper limits. Our light curves show a variety of early afterglow behaviour, from smooth, simple or broken power laws to 'bumpy', for a wide range of optical brightness (from the unprecedented faint detections of GRB 060108 and GRB 060510B to classical bright ones). We discuss GRB 051111 as an example of how the combination of optical and X-ray light curves can provide insight into the circumburst environment, in particular the role played by intrinsic extinction soon after the burst.

preprint2006arXiv

High Quality Early Time Light Curves of GRB 060206: Implications for Gamma Ray Burst Environments and Energetics

The 2-m robotic Liverpool Telescope (LT) reacted promptly to the high-redshift (z=4.048) gamma-ray burst GRB 060206. The afterglow was identified automatically and multicolor r'i'z' imaging was triggered without human intervention. Combining our data with those obtained from later follow-ups provides a well-sampled optical light-curve from 5 minutes to >2d after the gamma event. The light-curve is highly structured with at least three bumps evident in the first 75 minutes, including a major rebrightening (Delta_r'=-1.6 at t=3000s), interpreted as late energy injection. At early time (t=440s), we find evidence for fast (Delta t_(rest)<4s<<t) variability, indicating on-going internal-engine activity. We emphasise that a low redshift GRB (z<1) with similar intrinsic properties would have been interpreted completely differently due to undersampling of the light curve in the rest frame at early times; the light-curve behaviour of GRB 060206 should therefore not be considered peculiar. Finally, although the observed late-time steepening of the optical light curve resembles a jet break if taken in isolation, the lack of a corresponding change in the X-ray slope rules out a jet-break interpretation. Traditionally, GRB jet breaks have been inferred from optical data in the absence of simultaneous X-ray data. We suggest therefore that current estimates of the jet opening angle distribution might be biased by events like GRB060206. Consequently, the GRB explosion energy distribution and event rates may have to be revised.

preprint2005arXiv

Early GRB Optical and Infrared Afterglow Observations with the 2-m Robotic Liverpool Telescope

We present the first optical observations of a Gamma Ray Burst (GRB) afterglow using the 2-m robotic Liverpool Telescope (LT), which is owned and operated by Liverpool John Moores University and situated on La Palma. We briefly discuss the capabilities of LT and its suitability for rapid follow-up observations of early optical and infrared GRB light curves. In particular, the combination of aperture, site, instrumentation and rapid response (robotic over-ride mode aided by telescope's rapid slew and fully-opening enclosure) makes the LT ideal for investigating the nature of short bursts, optically-dark bursts, and GRB blast-wave physics in general. We briefly describe the LT's key position in the RoboNet-1.0 network of robotic telescopes. We present the LT observations of GRB041006 and use its gamma-ray properties to predict the time of the break in optical light curve, a prediction consistent with the observations.

preprint2005arXiv

On Statistical Lensing and the Anti-Correlation Between 2dF QSOs and Foreground Galaxies

We cross-correlate APM and SDSS galaxies with background QSOs from the 2dF QSO Redshift Survey, and detect a significant (2.8sigma) anti-correlation. The lack of a signal between 2dF stars and our galaxy samples suggests the anti-correlation is not due to a systematic error. The possibility that dust in the foreground galaxies could produce the anti-correlation is marginally rejected, at the 2sigma level through consideration of QSO colours. It is possible that dust that obscures QSOs without reddening them, or preferentially discards reddened QSOs from the 2QZ sample, could produce such an anti-correlation, however, such models are at odds with the positive QSO-galaxy correlations found at bright magnitudes by other authors. Our detection of a galaxy-QSO anti-correlation is consistent with statistical lensing theory. When combined with earlier results that have reported a positive galaxy-QSO correlation, a consistent, compelling picture emerges that spans faint and bright QSO samples showing positive or negative correlations according to the QSO N(m) slope. We find that galaxies are highly anti-biased on small scales. We consider two models that use different descriptions of the lensing matter and find they yield consistent predictions for the strength of galaxy bias on 0.1Mpc/h scales of b~0.1 (for LCDM). Whilst the slope of our power-law fit to the QSO-galaxy cross-correlation does not rule out linear bias, when we compare our measurement of b on 100 kpc/h scales to independent methods that determine b~1 on Mpc/h scales, we conclude that bias, on these small scales, is scale-dependent. These results indicate more mass, at least on the 100 kpc/h scales probed, than predicted by simple LCDM biasing prescriptions, and can thus constrain halo occupation models of the galaxy distribution.

preprint2005arXiv

The Liverpool Telescope Automatic Pipeline for Real-time GRB Afterglow Detection

The 2-m robotic Liverpool Telescope (LT) is ideally suited to the rapid follow-up of unpredictable and transient events such as GRBs. Our GRB follow-up strategy is designed to identify optical/IR counterparts in real time; it involves the automatic triggering of initial observations, on receipt of an alert from Gamma Ray Observatories HETE-2, INTEGRAL and Swift, followed by automated data reduction, analysis, OT identification and subsequent observing mode choice. The lack of human intervention in this process requires robustness at all stages of the procedure. Here we describe the telescope, its instrumentation and GRB pipeline.

preprint2004arXiv

The 2dF QSO Redshift Survey - XII. The spectroscopic catalogue and luminosity function

We present the final catalogue of the 2dF QSO Redshift Survey (2QZ), based on Anglo-Australian Telescope 2dF spectroscopic observations of 44576 colour-selected (u b_J r) objects with 18.25<b_J<20.85 selected from APM scans of UK Schmidt Telescope (UKST) photographic plates. The 2QZ comprises 23338 QSOs, 12292 galactic stars (including 2071 white dwarfs) and 4558 compact narrow-emission-line galaxies. We obtained a reliable spectroscopic identification for 86 per cent of objects observed with 2dF. We also report on the 6dF QSO Redshift Survey (6QZ), based on UKST 6dF observations of 1564 brighter 16<b_J<18.25 sources selected from the same photographic input catalogue. In total, we identified 322 QSOs spectroscopically in the 6QZ. The completed 2QZ is, by more than a factor 50, the largest homogeneous QSO catalogue ever constructed at these faint limits (b_J<20.85) and high QSO surface densities (35 QSOs deg^-2). As such it represents an important resource in the study of the Universe at moderate-to-high redshifts. As an example of the results possible with the 2QZ, we also present our most recent analysis of the optical QSO luminosity function and its cosmological evolution with redshift. For a flat, Omega_m=0.3 and Omega_lam=0.7, Universe, we find that a double power law with luminosity evolution that is exponential in look-back time, t, of the form L*(z) exp(6.15t), equivalent to an e-folding time of 2Gyr, provides an acceptable fit to the redshift dependence of the QSO luminosity function over the range 0.4 < z < 2.1 and M_bJ<-22.5. Evolution described by a quadratic in redshift is also an acceptable fit, with L*(z)~10^(1.39z-0.29z^2).

preprint2004arXiv

The 2dF QSO Redshift Survey - XIV. Structure and evolution from the two-point correlation function

We present a clustering analysis of QSOs using over 20000 objects from the final catalogue of the 2dF QSO Redshift Survey (2QZ), measuring the z-space correlation function, xi(s). When averaged over the range 0.3<z<2.2 we find that xi(s) is flat on small scales, steepening on scales above ~25h-1Mpc. In a WMAP/2dF cosmology we find a best fit power law with s_0=5.48+0.42-0.48h-1Mpc and gamma=1.20+-0.10 on scales s=1-25h-1Mpc. A CDM model assuming WMAP/2dF cosmological parameters is a good description of the QSO xi(s) after accounting for non-linear clustering and z-space distortions, and a linear bias of b_qso(z=1.35)=2.02+-0.07. We subdivide the 2QZ into 10 redshift intervals from z=0.53 to 2.48 and find a significant increase in clustering amplitude at high redshift in the WMAP/2dF cosmology. We derive the bias of the QSOs which is a strong function of redshift with b_qso(z=0.53)=1.13+-0.18 and b_qso(z=2.48)=4.24+-0.53. We use these bias values to derive the mean dark matter halo (DMH) mass occupied by the QSOs. At all redshifts 2QZ QSOs inhabit approximately the same mass DMHs with M_DH=(3.0+-1.6)x10^12h-1M_sun, which is close to the characteristic mass in the Press-Schechter mass function, M*, at z=0. If the relation between black hole (BH) mass and M_DH or host velocity dispersion does not evolve, then we find that the accretion efficiency (L/L_edd) for L* QSOs is approximately constant with redshift. Thus the fading of the QSO population from z~2 to 0 appears to be due to less massive BHs being active at low redshift. We apply different methods to estimate, t_qso, the active lifetime of QSOs and constrain this to be in the range 4x10^6-6x10^8 years at z~2. (Abridged).

preprint2003arXiv

Emission line widths and QSO black hole mass estimates from the 2dF QSO Redshift Survey

We have used composite spectra generated from more than 22000 QSOs observed in the course of the 2dF and 6dF QSO Redshift Surveys to investigate the relationship between the velocity width of emission lines and QSO luminosity. We find that the velocity width of the broad emission lines Hbeta, Hgamma, MgII, CIII] and CIV are correlated with the continuum luminosity, with a significance of more than 99 per cent. Of the major narrow emission lines ([OIII] 5007, [OII] 3727, NeIII 3870 and NeV 3426) only [OIII] exhibits a significant correlation between line width and luminosity. Assuming that the gas is moving in Keplerian orbits and that the radius of the broad line region is related to the QSO continuum luminosity, we use the velocity widths of the broad lines to derive average black hole masses for the QSOs contributing to the composite spectra. The resultant QSO mass-luminosity relationship is consistent with M ~ L^0.97+-0.16. We find that the correlation between line width and redshift, if present, must be weak, and only CIV shows significant evidence of evolution. This enables us to constrain the redshift evolution of the black hole mass-luminosity ratio to be ~(1+z)^beta with beta ~< 1, much less than the ~(1+z)^3 evolution seen in QSO luminosity evolution. Assuming that the motion of the broad line region gas is Keplerian and that its radius depends on the QSO luminosity, our models indicate that the observed weak redshift dependence is too small for the observed QSO luminosity function to be due to the evolution of a single long-lived population of sources.

preprint2003arXiv

Ghosts of the Milky Way: a search for topology in new quasar catalogues

We revisit the possibility that we inhabit a compact multi-connected flat, or nearly-flat, Universe. Analysis of COBE data has shown that, for such a case, the size of the fundamental domain must be a substantial fraction of the horizon size. Nevertheless, there could be several copies of the Universe within the horizon. If the Milky Way was once a quasar we might detect its `ghost' images. Using new large quasar catalogues we repeat the search by Fagundes & Wichoski for antipodal quasar pairs. By applying linear theory to account for the peculiar velocity of the local group, we are able to narrow the search radius to 134 arcsec. We find seven candidate antipodal quasar pairs within this search radius. However, a similar number would be expected by chance. We argue that, even with larger quasar catalogues, and more accurate values of the cosmological parameters, it is unlikely to be possible to identify putative ghost pairs unambiguously, because of the uncertainty of the correction for peculiar motion of the Milky Way.

preprint2003arXiv

Possible arcminute-separation gravitational lensed QSOs in the 2dF QSO survey

We report the possible discovery of multiple gravitationally-lensed images of QSOs with angular separations on arcminute scales. The QSOs were selected from the completed 2dF QSO survey as having redshifts and optical colours which are statistically consistent. In this paper we present higher-quality optical spectra of the candidates and discuss the likelihood of their genuinely being systems lensed by massive clusters of galaxies. From a comparison of the spectra it appears that up to six pairs of QSOs may be lensed multiple images, although the true number may be less than that and further observations should be undertaken to amass more evidence and to detect the lensing clusters. Two of the candidates may be associated with low redshift clusters of galaxies.

preprint2003arXiv

The 2dF QSO Redshift Survey - X. Lensing of Background QSOs by Galaxy Groups

We cross-correlate QSOs from the 2dF Survey with galaxy groups. The galaxy samples are limited to B < 20.5. We use an objective algorithm to detect galaxy groups. A 3sigma anti-correlation is observed between QSOs and galaxy groups. This paucity of faint QSOs around groups is neither a selection effect nor due to restrictions on the placement of 2dF fibres. By observing the colours of QSOs on the scales of the anti-correlation, we limit dust in galaxy groups, finding a maximum reddening of E(b_j-r) < 0.012 at the 95% level. The small amount of dust thus inferred is insufficient to cause the anti-correlation, supporting the suggestion by Croom & Shanks that the signal is due to gravitational lensing. The possibility remains that tailored dust models, such as grey dust, heavy patches of dust or a combination of dust and lensing, could explain the anti-correlation. Assuming the signal is caused by lensing rather than dust, we measure the average velocity dispersion of a Singular Isothermal Sphere that would cause the anti-correlation as around 1150 km/s. Simulations reject 600 km/s at the 5% significance level. We also model foreground lenses as NFW haloes and measure the typical mass within 1.5 Mpc/h of the halo centre as M_{1.5} = (1.2 +/- 0.9) x 10^{15} solarmasses/h. Regardless of whether we utilise a SIS or NFW dark matter profile, our model favours more mass in groups than accounted for in a universe with density parameter Omega_m = 0.3. Detailed simulations and galaxy group redshifts will significantly reduce the current systematic uncertainties in these $Ω_m$ estimates. Reducing the remaining uncertainty will require larger QSO and galaxy group surveys (abridged).

preprint2003arXiv

The 2dF QSO Redshift Survey - XI. The QSO Power Spectrum

We present a power spectrum analysis of the final 2dF QSO Redshift Survey catalogue containing 22652 QSOs. Utilising the huge volume probed by the QSOs, we can accurately measure power out to scales of ~500Mpc and derive new constraints, at z~1.4, on the matter and baryonic contents of the Universe. Importantly, these new cosmological constraints are derived at an intermediate epoch between the CMB observations at z~1000, and local (z~0) studies of large-scale structure; the average QSO redshift corresponds to a look-back time of approximately two-thirds of the age of the Universe. We find that the amplitude of clustering of the QSOs at z~1.4 is similar to that of present day galaxies. The power spectra of the QSOs at high and low redshift are compared and we find little evidence for any evolution in the amplitude. Assuming a lambda cosmology to derive the comoving distances, r(z), to the QSOs, the power spectrum derived can be well described by a model with shape parameter Gamma=0.13+-0.02. If an Einstein-de Sitter model r(z) is instead assumed, a slightly higher value of Gamma=0.16+-0.03 is obtained. A comparison with the Hubble Volume LCDM simulation shows very good agreement over the whole range of scales considered. A standard (Omega_m=1) CDM model, however, predicts a much higher value of Gamma than is observed, and it is difficult to reconcile such a model with these data. We fit CDM model power spectra (assuming scale-invariant initial fluctuations), convolved with the survey window function, and corrected for redshift space distortions, and find that models with baryon oscillations are slightly preferred, with the baryon fraction Omega_b/Omega_m=0.18+-0.10. The overall shape of the power spectrum provides a strong constraint on Omega_m*h (where h is the Hubble parameter), with Omega_m*h=0.19+-0.05.

preprint1997arXiv

Galaxy clusters in the Perseus-Pisces region - I. Spectroscopic and photometric data for early-type galaxies

We present new spectroscopic and photometric data for 137 early-type galaxies in nine clusters, and for a set of nearby standard galaxies. The clusters studied are Perseus (A0426), Pisces, A0262, A0347, J8, HMS0122+3305, 7S21, A2199 and A2634. Our spectroscopic data comprise radial velocities (cz), central velocity dispersions (sigma) and magnesium line strength indices (Mg). Internal errors (derived from repeat observations) are 7.6% on each measurement of sigma, and 0.010 mag on each Mg measurement. We correct the major published datasets to a standard aperture size, and define a new `standard system' by the aperture-corrected Lick data of Davies et al. Through extensive intercomparisons with data from the literature, we present the corrections required to bring the major published datasets onto the standard system. We demonstrate that our new velocity dispersion data can be brought into consistency with the standard system, to an uncertainty of <0.01 dex. From R-band CCD photometry, we derive effective diameter (A_e), mean surface brightness within effective diameter (mu_e) and an R-band diameter equivalent to the D_n parameter of Dressler et al. Internal comparisons indicate an average error of 0.005 in each measurement of log(D_n). The combination log(A_e) - 0.3 mu_e, approximately the quantity used in the Fundamental Plane distance indicator, has an uncertainty of 0.006 per measurement. The photometric data can be brought onto a system consistent with external data at the level of 0.5% in distance. These data will be used in a companion paper to derive distance and peculiar velocity estimates for the nine clusters studied. (Abridged)

preprint1997arXiv

Galaxy clusters in the Perseus-Pisces region - II. The peculiar velocity field

We have measured the mean peculiar motions of 103 early-type galaxies in 7 clusters in the Perseus-Pisces (PP) ridge or PP background, and a further 249 such galaxies in 9 calibrating clusters from the literature, using the inverse Fundamental Plane relation. This relation is found to have a distance error of 20% per galaxy. None of the 6 clusters in the PP ridge has a significant motion with respect to the CMB frame, but the PP background cluster J8 shows marginal evidence of `backside infall' into the PP supercluster. The full 16 cluster sample has a mean CMB-frame bulk motion of 420 +/- 280 km/s, towards l=262, b=-25. This result is consistent both with no bulk motion in the CMB frame and with the ~350 km/s bulk motion found by Courteau et al. It is inconsistent at the 98% confidence level with the ~700 km/s bulk flow found by Lauer & Postman (LP). The PP ridge clusters are found to have a small and statistically insignificant mean radial motion with respect to the CMB frame: -60 +/- 220 km/s. A comparison between our cluster peculiar velocities and the predicted peculiar velocities from the IRAS 1.2 Jy density field, smoothed on a 500 km/s scale, yields beta_I = Omega^0.6/b_I = 0.95 +/- 0.48, consistent with previous results. We find agreement between our peculiar motions and published Tully-Fisher results for the same clusters. The disagreement between the 11 clusters common to our sample and that of LP, based on BCGs, is statistically significant at the >~ 99.7% confidence level indicating that the errors of one or both of these data sets are underestimated. When the BCG distances corrected for the X-ray luminosity of the host cluster are used, the disagreement is reduced to the ~94% confidence level. (Abridged)

preprint1995arXiv

The environments of z<0.3 QSOs

We have carried out an investigation of the galaxy environments of low redshift ($z<0.3$) QSOs by cross-correlating the positions on the sky of X-ray-selected QSOs/AGN identified in the {\it Einstein} Medium Sensitivity Survey (EMSS) with those of $B_J<20.5$ galaxies in the APM galaxy catalogues. At $<5\,$arcmin, we find a significant ($5σ$) galaxy excess around $z<0.3$ QSOs. The amplitude of the low redshift ($z<0.3$) QSO-galaxy angular cross-correlation function is identical to that of the APM galaxy angular correlation function, implying that these (predominantly radio-quiet) QSOs inhabit environments similar to those of normal galaxies. No significant galaxy excess was found around a `control' sample of $z>0.3$ QSOs. Coupled with previous observations, these results imply that the environment of radio-quiet QSOs undergoes little evolution over a wide range in redshift ($0 less than z less than 1.5$). This is in marked contrast to the rapid increase in the richness of the environments associated with radio-loud QSOs over the same redshift range. The similarity between QSO-galaxy clustering and galaxy-galaxy clustering also suggests that QSOs are unbiased with respect to galaxies and make useful tracers of large-scale structure in the Universe.