Measuring the Physical Properties of Distant Galaxies and Black Holes in the Era of Surveys

In the era of deep and wide-field surveys (e.g. SDSS, LSST, LOFAR, SKA), we have access to an ever-increasing volume of multi-wavelength data for millions of galaxies both locally and at high redshifts. However, inferring the intrinsic properties of the whole population of galaxies requires robust statistical techniques and an understanding of observational bias. In this thesis, I present a study of the Far-Infrared Radio Correlation (FIRC) – a relation which is widely used to infer star-formation rates from otherwise featureless radio sources. Using LOFAR 150MHz, FIRST 1:4GHz, and Herschel infrared luminosities derived from the new LOFAR/H-ATLAS catalogue, we investigate possible variation in the monochromatic (250mm) FIRC at low and high radio frequencies. Although the average FIRC at high radio frequency is consistent with expectations based on a standard power-law radio spectrum, the average correlation at 150MHz is not. We see evidence for redshift evolution of the FIRC at 150MHz, and find that the FIRC varies with stellar mass, dust temperature and specific star formation rate, whether the latter is probed using MAGPHYS fitting, or using mid-infrared colour as a proxy. We can explain the variation, to within σ using a Bayesian partial correlation technique. This work was published as Read et al. (2018) in the Monthly Notices of the Royal Astronomical Society. Identifying an opportunity to increase in the efficiency of black-hole mass estimations, we perform photometric reverberation mapping using the Javelin photometric damped random walk model for the QSO SDSS J144645.44 +625304.0 at ɀ = 0:351 and estimate the Hβ lag of 72 +5-1 days and black hole mass of 108:28 +0.12-0.07 Mʘ. An analysis of the reliability of photometric reverberation mapping conducted using many thousands of simulated light curves shows that we can recover any input lag less than a third of the duration of our observing campaign to within 4per cent on average given our target’s observed signal-to-noise of > 20 and cadence of 14 days. We use our suite of simulated light curves to deconvolve artefacts from the QSO’s posterior lag distribution, increasing the signal-to-noise by a factor of ~3. We exceed the signal-to-noise of the Sloan Digital Sky Survey Reverberation Mapping Project (SDSS-RM) campaign with a quarter of the observing time resulting in a ~310 per cent per cent increase in SNR efficiency over SDSS-RM. Finally, I present a study of the radio luminosity star-formation rate relation directly with the LOFAR Two Metre Sky Survey (LoTSS) DR1, in an effort to understand the mass dependency of the L150MHz - SFR slope reported by Gürkan et al. (2018). Building on our previous study of the FIRC, we develop a fast, generalised algorithm to recover Complete And Noiseless Distributions from Incomplete Data (CANDID). We find that the mass dependency is real and in agreement with previous estimations in the literature when we include the effects of selection biases present in the LoTSS DR1 sample. We also propose that type-Ia supernovae may contribute to a L150MHz excess and construct a joint distribution of our LoTSS observations and the Horizon AGN simulation to test this

A Markov Chain Monte Carlo approach for measurement of jet precession in radio-loud active galactic nuclei

© 2020 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Jet precession can reveal the presence of binary systems of supermassive black holes. The ability to accurately measure the parameters of jet precession from radio-loud AGN is important for constraining the binary supermassive black hole population, which are expected as a result of hierarchical galaxy evolution. The age, morphology, and orientation along the line of sight of a given source often result in uncertainties regarding jet path. This paper presents a new approach for efficient determination of precession parameters using a 2D MCMC curve-fitting algorithm which provides us a full posterior probability distribution on the fitted parameters. Applying the method to Cygnus A, we find evidence for previous suggestions that the source is precessing. Interpreted in the context of binary black holes leads to a constraint of parsec scale and likely sub-parsec orbital separation for the putative supermassive binary.Peer reviewe

The performance of photometric reverberation mapping at high redshift and the reliability of damped random walk models

Accurate methods for reverberation mapping using photometry are highly sought after since they are inherently less resource intensive than spectroscopic techniques. However, the effectiveness of photometric reverberation mapping for estimating black hole masses is sparsely investigated at redshifts higher than z ≈ 0.04. Furthermore, photometric methods frequently assume a damped random walk (DRW) model, which may not be universally applicable

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

© 2023 The Author(s) . Published by Oxford University Press on behalf of Royal Astronomical Society. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY), https://creativecommons.org/licenses/by/4.0/WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.Peer reviewe

VizieR Online Data Catalog: H-ATLAS NGP LOFAR radio catalogue (Hardcastle+, 2016)

The NGP field was observed in four separate pointings, chosen to maximize sky covered, with the LOFAR HBA as part of the Surveys Key Science project. Observations used the HBADUALINNER mode, meaning that the station beams of core and remote stations roughly matched each other and giving the widest possible field of view. The first observation, which was made early on in LOFAR operations, was of slightly longer duration (~10h) than the others (~8h). International stations were included in some of the observations in 2014 but were not used in any of our analysis, which uses only the Dutch array. (1 data file)

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same regio

The wide-field, multiplexed, spectroscopic facility WEAVE : survey design, overview, and simulated implementation

Funding for the WEAVE facility has been provided by UKRI STFC, the University of Oxford, NOVA, NWO, Instituto de Astrofísica de Canarias (IAC), the Isaac Newton Group partners (STFC, NWO, and Spain, led by the IAC), INAF, CNRS-INSU, the Observatoire de Paris, Région Île-de-France, CONCYT through INAOE, Konkoly Observatory (CSFK), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Lund University, the Leibniz Institute for Astrophysics Potsdam (AIP), the Swedish Research Council, the European Commission, and the University of Pennsylvania.WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366-959 nm at R ∼ 5000, or two shorter ranges at R ∼ 20,000. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for ∼ 3 million stars and detailed abundances for ∼ 1.5 million brighter field and open-cluster stars; (ii) survey ∼ 0.4 million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey  ∼ 400 neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in z 1 million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at z > 2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.PostprintPeer reviewe

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.Comment: 41 pages, 27 figures, accepted for publication by MNRA

The development and validation of a scoring tool to predict the operative duration of elective laparoscopic cholecystectomy

Background: The ability to accurately predict operative duration has the potential to optimise theatre efficiency and utilisation, thus reducing costs and increasing staff and patient satisfaction. With laparoscopic cholecystectomy being one of the most commonly performed procedures worldwide, a tool to predict operative duration could be extremely beneficial to healthcare organisations. Methods: Data collected from the CholeS study on patients undergoing cholecystectomy in UK and Irish hospitals between 04/2014 and 05/2014 were used to study operative duration. A multivariable binary logistic regression model was produced in order to identify significant independent predictors of long (> 90 min) operations. The resulting model was converted to a risk score, which was subsequently validated on second cohort of patients using ROC curves. Results: After exclusions, data were available for 7227 patients in the derivation (CholeS) cohort. The median operative duration was 60 min (interquartile range 45–85), with 17.7% of operations lasting longer than 90 min. Ten factors were found to be significant independent predictors of operative durations > 90 min, including ASA, age, previous surgical admissions, BMI, gallbladder wall thickness and CBD diameter. A risk score was then produced from these factors, and applied to a cohort of 2405 patients from a tertiary centre for external validation. This returned an area under the ROC curve of 0.708 (SE = 0.013, p  90 min increasing more than eightfold from 5.1 to 41.8% in the extremes of the score. Conclusion: The scoring tool produced in this study was found to be significantly predictive of long operative durations on validation in an external cohort. As such, the tool may have the potential to enable organisations to better organise theatre lists and deliver greater efficiencies in care

Search for dark matter produced in association with bottom or top quarks in √s = 13 TeV pp collisions with the ATLAS detector

A search for weakly interacting massive particle dark matter produced in association with bottom or top quarks is presented. Final states containing third-generation quarks and miss- ing transverse momentum are considered. The analysis uses 36.1 fb−1 of proton–proton collision data recorded by the ATLAS experiment at √s = 13 TeV in 2015 and 2016. No significant excess of events above the estimated backgrounds is observed. The results are in- terpreted in the framework of simplified models of spin-0 dark-matter mediators. For colour- neutral spin-0 mediators produced in association with top quarks and decaying into a pair of dark-matter particles, mediator masses below 50 GeV are excluded assuming a dark-matter candidate mass of 1 GeV and unitary couplings. For scalar and pseudoscalar mediators produced in association with bottom quarks, the search sets limits on the production cross- section of 300 times the predicted rate for mediators with masses between 10 and 50 GeV and assuming a dark-matter mass of 1 GeV and unitary coupling. Constraints on colour- charged scalar simplified models are also presented. Assuming a dark-matter particle mass of 35 GeV, mediator particles with mass below 1.1 TeV are excluded for couplings yielding a dark-matter relic density consistent with measurements