We consider the observation of stellar-mass black holes binaries with the Laser Interferometer Space Antenna (LISA). Preliminary results based on Fisher information matrix analyses have suggested that gravitational waves from those sources could be very sensitive to possible deviations from the theory of general relativity and from the strong equivalence principle during the low-frequency binary inspiral. We perform a full Markov Chain Monte Carlo Bayesian analysis to quantify the sensitivity of these signals to two phenomenological modifications of general relativity, namely a putative gravitational dipole emission and a non-zero mass for the graviton, properly accounting for the detector's response. Moreover, we consider a scenario where those sources could be observed also with Earth-based detectors, which should measure the coalescence time with precision better than $1 \ {\rm ms}$. This constraint on the coalescence time further improves the bounds that we can set on those phenomenological deviations from general relativity. We show that tests of dipole radiation and the graviton's mass should improve respectively by seven and half an order(s) of magnitude over current bounds. Finally, we discuss under which conditions one may claim the detection of a modification to general relativity.

The space-based gravitational wave detector LISA will observe mergers of massive black hole binary systems (MBHBs) to cosmological distances, as well as inspiralling stellar-origin (or stellar-mass) binaries (SBHBs) years before they enter the LIGO/Virgo band. Much remains to be explored for the parameter recovery of both classes of systems. Previous MBHB analyses relied on inspiral-only signals and/or a simplified Fisher matrix analysis, while SBHBs have not yet been extensively analyzed with Bayesian methods. We accelerate likelihood computations by (i) using a Fourier-domain response of the LISA instrument, (ii) using a reduced-order model for non-spinning waveforms that include a merger-ringdown and higher harmonics, (iii) setting the noise realization to zero and computing overlaps in the amplitude/phase representation. We present the first simulations of Bayesian inference for the parameters of massive black hole systems including consistently the merger and ringdown of the signal, as well as higher harmonics. We clarify the roles of LISA response time and frequency dependencies in breaking degeneracies and illustrate how degeneracy breaking unfolds over time. We also find that restricting the merger-dominated signal to its dominant harmonic can make the extrinsic likelihood very degenerate. Including higher harmonics proves to be crucial to break degeneracies and considerably improves the localization of the source, with a surviving bimodality in the sky position. We also present simulations of Bayesian inference for the extrinsic parameters of SBHBs, and show that although unimodal, their posterior distributions can have non-Gaussian features.

The signal processing firmware that has been developed for the Low Frequency Aperture Array component of the Square Kilometre Array is described. The firmware is implemented on a dual FPGA board, that is capable of processing the streams from 16 dual polarization antennas. Data processing includes channelization of the sampled data for each antenna, correction for instrumental response and for geometric delays and formation of one or more beams by combining the aligned streams. The channelizer uses an oversampling polyphase filterbank architecture, allowing a frequency continuous processing of the input signal without discontinuities between spectral channels. Each board processes the streams from 16 antennas, as part of larger beamforming system, linked by standard Ethernet interconnections. There are envisaged to be 8192 of these signal processing platforms in the first phase of the Square Kilometre array so particular attention has been devoted to ensure the design is low cost and low power.

Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field.

The zodiacal dust complex, a population of dust and small particles that pervades the Solar System, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. Here we present a new set of data obtained using a novel technique: direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible by the extreme precision of the instruments flown on the LISA Pathfinder spacecraft, a technology demonstrator for a future space-based gravitational wave observatory that operated near the first Sun-Earth Lagrange point from early 2016 through Summer of 2017. Using a simple model of the impacts and knowledge of the control system, we show that it is possible to detect impacts and measure properties such as the transferred momentum (related to the particle's mass and velocity), direction of travel, and location of impact on the spacecraft. In this paper, we present the results of a systematic search for impacts during 4348 hours of Pathfinder data. We report a total of 54 candidates with momenta ranging from 0.2$\,\mu\textrm{Ns}$ to 230$\,\mu\textrm{Ns}$. We furthermore make a comparison of these candidates with models of micrometeoroid populations in the inner solar system including those resulting from Jupiter-family comets, Oort-cloud comets, Hailey-type comets, and Asteroids. We find that our measured population is consistent with a population dominated by Jupiter-family comets with some evidence for a smaller contribution from Hailey-type comets. This is in agreement with consensus models of the zodiacal dust complex in the momentum range sampled by LISA Pathfinder.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps.

The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge of gravitational-wave astronomy will strongly compel a subsequent mission to further explore the frequency bands of the GW spectrum that can only be accessed from space. The 2020s is the time to start developing technology and studying mission concepts for a large-scale mission to be launched in the 2040s. The mission concept would then be proposed to Astro2030. Only space based missions can access the GW spectrum between 10 nHz and 1 Hz because of the Earths seismic noise. This white paper surveys the science in this band and mission concepts that could accomplish that science. The proposed small scale activity is a technology development program that would support a range of concepts and a mission concept study to choose a specific mission concept for Astro2030. In this white paper, we will refer to a generic GW mission beyond LISA as bLISA.

Experience suggests that structural issues in how institutional Astrophysics approaches data-driven science and the development of discovery technology may be hampering the community's ability to respond effectively to a rapidly changing environment in which increasingly complex, heterogeneous datasets are challenging our existing information infrastructure and traditional approaches to analysis. We stand at the confluence of a new epoch of multimessenger science, remote co-location of data and processing power and new observing strategies based on miniaturized spacecraft. Significant effort will be required by the community to adapt to this rapidly evolving range of possible discovery moduses. In the suggested creation of a new Astrophysics element, Advanced Astrophysics Discovery Technology, we offer an affirmative solution that places the visibility of discovery technologies at a level that we suggest is fully commensurate with their importance to the future of the field.

By listening to gravity in the low frequency band, between 0.1 mHz and 1 Hz, the future space-based gravitational-wave observatory LISA will be able to detect tens of thousands of astrophysical sources from cosmic dawn to the present. The detection and characterization of all resolvable sources is a challenge in itself, but LISA data analysis will be further complicated by interruptions occurring in the interferometric measurements. These interruptions will be due to various causes occurring at various rates, such as laser frequency switches, high-gain antenna re-pointing, orbit corrections, or even unplanned random events. Extracting long-lasting gravitational-wave signals from gapped data raises problems such as noise leakage and increased computational complexity. We address these issues by using Bayesian data augmentation, a method that reintroduces the missing data as auxiliary variables in the sampling of the posterior distribution of astrophysical parameters. This provides a statistically consistent way to handle gaps while improving the sampling efficiency and mitigating leakage effects. We apply the method to the estimation of galactic binaries parameters with different gap patterns, and we compare the results to the case of complete data.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

Binary black hole (BBH) mergers provide a prime source for current and future interferometric GW observatories. Massive BBH mergers may often take place in plasma-rich environments, leading to the exciting possibility of a concurrent electromagnetic (EM) signal observable by traditional astronomical facilities. However, many critical questions about the generation of such counterparts remain unanswered. We explore mechanisms that may drive EM counterparts with magnetohydrodynamic simulations treating a range of scenarios involving equal-mass black-hole binaries immersed in an initially homogeneous fluid with uniform, orbitally aligned magnetic fields. We find that the time development of Poynting luminosity, which may drive jet-like emissions, is relatively insensitive to aspects of the initial configuration. In particular, over a significant range of initial values, the central magnetic field strength is effectively regulated by the gas flow to yield a Poynting luminosity of $10^{45}-10^{46} \rho_{-13} M_8^2 \, {\rm erg}\,{\rm s}^{-1}$, with BBH mass scaled to $M_8 \equiv M/(10^8 M_{\odot})$ and ambient density $\rho_{-13} \equiv \rho/(10^{-13} \, {\rm g} \, {\rm cm}^{-3})$. We also calculate the direct plasma synchrotron emissions processed through geodesic ray-tracing. Despite lensing effects and dynamics, we find the observed synchrotron flux varies little leading up to merger.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

The study of galaxy mergers and supermassive binary black holes (SMBBHs) is central to our understanding of the galaxy and black hole assembly and (co-)evolution at the epoch of structure formation and throughout cosmic history. Galaxy mergers are the sites of major accretion episodes, they power quasars, grow supermassive black holes (SMBHs), and drive SMBH-host scaling relations. The coalescing SMBBHs at their centers are the loudest sources of gravitational waves (GWs) in the universe, and the subsequent GW recoil has a variety of potential astrophysical implications which are still under exploration. Future GW astronomy will open a completely new window on structure formation and galaxy mergers, including the direct detection of coalescing SMBBHs, high-precision measurements of their masses and spins, and constraints on BH formation and evolution in the high-redshift universe.

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift/BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien 2014 is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of $\gtrsim97\%$ ($\lesssim 3\%$ error), which is a significant improvement on a cut in GRB flux which has an accuracy of $89.6\%$ ($10.4\%$ error). These models are then used to measure the detection efficiency of Swift as a function of redshift $z$, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of $n_0 \sim 0.48^{+0.41}_{-0.23} \ {\rm Gpc}^{-3} {\rm yr}^{-1}$ with power-law indices of $n_1 \sim 1.7^{+0.6}_{-0.5}$ and $n_2 \sim -5.9^{+5.7}_{-0.1}$ for GRBs above and below a break point of $z_1 \sim 6.8^{+2.8}_{-3.2}$. This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting. The code used in this is analysis is publicly available online (https://github.com/PBGraff/SwiftGRB_PEanalysis).

Robust gauge conditions are critically important to the stability and accuracy of numerical relativity (NR) simulations involving compact objects. Most of the NR community use the highly robust---though decade-old---moving-puncture (MP) gauge conditions for such simulations. It has been argued that in binary black hole (BBH) evolutions adopting this gauge, noise generated near adaptive-mesh-refinement (AMR) boundaries does not converge away cleanly with increasing resolution, severely limiting gravitational waveform accuracy at computationally feasible resolutions. We link this noise to a sharp (short-wavelength), initial outgoing gauge wave crossing into progressively lower resolution AMR grids, and present improvements to the standard MP gauge conditions that focus on stretching, smoothing, and more rapidly settling this outgoing wave. Our best gauge choice greatly reduces gravitational waveform noise during inspiral, yielding less fluctuation in convergence order and $\sim 40%$ lower waveform phase and amplitude errors at typical resolutions. Noise in other physical quantities of interest is also reduced, and constraint violations drop by more than an order of magnitude. We expect these improvements will carry over to simulations of all types of compact binary systems, as well as other $N$+1 formulations of gravity for which MP-like gauge conditions can be chosen.

Advanced LIGO and Advanced Virgo will be all-sky monitors for merging compact objects within a few hundred Mpc. Finding the electromagnetic counterparts to these events will require an understanding of the transient sky at low red-shift (z<0.1). We performed a systematic search for extragalactic, low red-shift, transient events in the XMM-Newton Slew Survey. In a flux limited sample, we found that highly-variable objects comprised 10% of the sample, and that of these, 10% were spatially coincident with cataloged optical galaxies. This led to 4x10^-4 transients per square degree above a flux threshold of 3x10^-12 erg cm-2 s-1 [0.2-2 keV] which might be confused with LIGO/Virgo counterparts. This represents the first extragalactic measurement of the soft X-ray transient rate within the Advanced LIGO/Virgo horizon. Our search revealed six objects that were spatially coincident with previously cataloged galaxies, lacked evidence for optical AGNs, displayed high luminosities around 10^43 erg s-1, and varied in flux by more than a factor of ten when compared with the ROSAT All-Sky Survey. At least four of these displayed properties consistent with previously observed tidal disruption events.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this Letter, we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe a total amplification of the magnetic field of ~2 orders of magnitude which is driven by the accretion onto the binary and that leads to much stronger electromagnetic signals, more than a factor of 10^4 larger than comparable calculations done in the force-free regime where such amplifications are not possible.

We have used the Taurus Tunable Filter to search for Lyman alpha emitters in the fields of three high redshift quasar fields: two at z~2.2 (MRC B1256-243 and MRC B2158-206) and one at z~4.5 (BR B0019-1522). Our observations had a field of view of around 35 square arcminutes, and reached AB magnitudes of magnitudes of ~21 (MRC B1256-243), ~22 (MRC B2158-206), and ~22.6 (BR B0019-1522), dependent on wavelength. We have identified candidate emission line galaxies in all three of the fields, with the higher redshift field being by far the richest. By combining our observations with simulations of the instrumental response, we estimate the total density of emission line galaxies in each field. Seventeen candidate emission line galaxies were found in within 1.5 Mpc of BR0019-1522, a number density of 4.9 +/- 1.2 x 10^-3 Mpc^-3, suggesting a significant galaxy overdensity at z~4.5.

We conduct a descriptive analysis of the multipolar structure of gravitational-radiation waveforms from equal-mass aligned-spin mergers, following an approach first presented in the complementary context of nonspinning black holes of varying mass ratio [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)]. We find that, as with the nonspinning mergers, the dominant waveform mode phases evolve together in lock-step through inspiral and merger, supporting the previous waveform description in terms of an adiabatically rigid rotator driving gravitational-wave emission-an implicit rotating source. We further apply the late-time merger-ringdown model for the rotational frequency introduced in [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)], along with an improved amplitude model appropriate for the dominant (2, \pm2) modes. This provides a quantitative description of the merger-ringdown waveforms, and suggests that the major features of these waveforms can be described with reference only to the intrinsic parameters associated with the state of the final black hole formed in the merger. We provide an explicit model for the merger-ringdown radiation, and demonstrate that this model agrees to fitting factors better than 95% with the original numerical waveforms for system masses above \sim150M\odot. This model may be directly applicable to gravitational-wave detection of intermediate-mass black-hole mergers.

We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning black hole binaries, combining for the first time the entire inspiral-merger-ringdown signal, the effect of the LISA orbits, and the complete three-channel LISA response. We consider an ensemble of systems near the peak of LISA's sensitivity band, with total rest mass of 2\times10^6 M\odot, a redshift of z = 1, and randomly chosen orientations and sky positions. We find median sky localization errors of approximately \sim3 arcminutes. This is comparable to the field of view of powerful electromagnetic telescopes, such as the James Webb Space Telescope, that could be used to search for electromagnetic signals associated with merging massive black holes. We investigate the way in which parameter errors decrease with measurement time, focusing specifically on the additional information provided during the merger-ringdown segment of the signal. We find that this information improves all parameter estimates directly, rather than through diminishing correlations with any subset of well- determined parameters. Although we have employed the baseline LISA design for this study, many of our conclusions regarding the information provided by mergers will be applicable to alternative mission designs as well.

Recent breakthroughs in the field of numerical relativity have led to dramatic progress in understanding the predictions of General Relativity for the dynamical interactions of two black holes in the regime of very strong gravitational fields. Such black-hole binaries are important astrophysical systems and are a key target of current and developing gravitational-wave detectors. The waveform signature of strong gravitational radiation emitted as the black holes fall together and merge provides a clear observable record of the process. After decades of slow progress, these mergers and the gravitational-wave signals they generate can now be routinely calculated using the methods of numerical relativity. We review recent advances in understanding the predicted physics of events and the consequent radiation, and discuss some of the impacts this new knowledge is having in various areas of astrophysics.

Understanding the predictions of general relativity for the dynamical interactions of two black holes has been a long-standing unsolved problem in theoretical physics. Black-hole mergers are monumental astrophysical events, releasing tremendous amounts of energy in the form of gravitational radiation, and are key sources for both ground- and space-based gravitational-wave detectors. The black-hole merger dynamics and the resulting gravitational waveforms can only be calculated through numerical simulations of Einstein's equations of general relativity. For many years, numerical relativists attempting to model these mergers encountered a host of problems, causing their codes to crash after just a fraction of a binary orbit could be simulated. Recently, however, a series of dramatic advances in numerical relativity has allowed stable, robust black-hole merger simulations. This remarkable progress in the rapidly maturing field of numerical relativity, and the new understanding of black-hole binary dynamics that is emerging is chronicled. Important applications of these fundamental physics results to astrophysics, to gravitational-wave astronomy, and in other areas are also discussed.

We investigate the precision with which the parameters describing the characteristics and location of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna (LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the signals, we find that LISA will have far greater precision than previous estimates for nonspinning mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with moderate mass ratios, q >= 1/10, and total masses 10^5 < M/M_{Sun} < 10^7. We compare the parameter uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry and higher-order content to the predicted parameter uncertainties resulting from inclusion of the merger. In real-time observations, the later parts of the signal lead to significant improvements in sky-position precision in the last hours and even the final minutes of observation. For comparable mass systems with total mass M/M_{Sun} = ~10^6, we find that the increased precision resulting from including the merger is comparable to the increase in signal-to-noise ratio. For the most precise systems under investigation, half can be localized to within O(10 arcmin), and 10% can be localized to within O(1 arcmin).

Although the gravitational wave kick velocity in the orbital plane of coalescing black holes has been understood for some time, apparently conflicting formulae have been proposed for the dominant out-of-plane kick, each a good fit to different data sets. This is important to resolve because it is only the out-of-plane kicks that can reach more than 500 km/s and can thus eject merged remnants from galaxies. Using a different ansatz for the out-of-plane kick, we show that we can fit almost all existing data to better than 5 %. This is good enough for any astrophysical calculation, and shows that the previous apparent conflict was only because the two data sets explored different aspects of the kick parameter space.

Infrared spectroscopy of the H-alpha emission lines of a sub-sample of 19 high-redshift (0.8 < z < 2.3) Molonglo quasars, selected at 408 MHz, is presented. These emission lines are fitted with composite models of broad and narrow emission, which include combinations of classical broad-line regions of fast-moving gas clouds lying outside the quasar nucleus, and/or a theoretical model of emission from an optically-thick, flattened, rotating accretion disk. All bar one of the nineteen sources are found to have emission consistent with the presence of an optically-emitting accretion disk, with the exception appearing to display complex emission including at least three broad components. Ten of the quasars have strong Bayesian evidence for broad-line emission arising from an accretion disk together with a standard broad-line region, selected in preference to a model with two simple broad lines. Thus the best explanation for the complexity required to fit the broad H-alpha lines in this sample is optical emission from an accretion disk in addition to a region of fast-moving clouds. We derive estimates of the angle between the rotation axis of the accretion disk and the line of sight. A weak correlation is found between the accretion disk angle and the logarithm of the low-frequency radio luminosity. This is direct, albeit tenuous, evidence for the receding torus model. Velocity shifts of the broad H-alpha components are analysed and the results found to be consistent with a two-component model comprising one single-peaked broad line emitted at the same redshift as the narrow lines, and emission from an accretion disk which appears to be preferentially redshifted with respect to the narrow lines for high-redshift sources and blueshifted relative to the narrow lines for low-redshift sources.

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

Recent advances in numerical relativity provide a detailed description of the waveforms of coalescing massive black hole binaries (MBHBs), expected to be the strongest detectable LISA sources. We present a preliminary study of LISA's sensitivity to MBHB parameters using a hybrid numerical/analytic waveform for equal-mass, non-spinning holes. The Synthetic LISA software package is used to simulate the instrument response and the Fisher information matrix method is used to estimate errors in the parameters. Initial results indicate that inclusion of the merger signal can significantly improve the precision of some parameter estimates. For example, the median parameter errors for an ensemble of systems with total redshifted mass of one million Solar masses at a redshift of one were found to decrease by a factor of slightly more than two for signals with merger as compared to signals truncated at the Schwarzchild ISCO.

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning black holes has a maximum ``recoil kick'' of up to $\sim 4000 \kms$. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick perpendicular to the orbital plane does not scale as $\sim\eta^2$ (where $\eta$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\sim\eta^3$, at least for systems with low orbital precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters.

We present a multipolar analysis of the gravitational recoil computed in recent numerical simulations of binary black hole (BH) coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including l=4 are sufficient to accurately reproduce the final recoil velocity (within ~2%) and that only a few dominant modes contribute significantly to it (within ~5%). We describe how the relative amplitudes, and more importantly, the relative phases, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ringdown phases. We also find that the numerical results can be reproduced by an ``effective Newtonian'' formula for the multipole moments obtained by replacing the radial separation in the Newtonian formulae with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes (QNMs). Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of ``anti-kick'' for each of the simulations. Lastly, we apply this multipolar analysis to help explain the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Several groups have recently computed the gravitational radiation recoil produced by the merger of two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black holes with mass ratio $q\equiv m_1/m_2=2/3$ and dimensionless spin parameters of $a_1/m_1$ and $a_2/m_2$ equal to 0 or 0.2, with directions aligned or anti-aligned with the orbital angular momentum. We also run an equal-mass $a_1/m_1=-a_2/m_2=0.2$ case, and find good agreement with previous results. We find that, for currently reported kicks from aligned or anti-aligned spins, a simple kick formula inspired by post-Newtonian analyses can reproduce the numerical results to better than $\sim$10%.

Coalescing black-hole binaries are expected to be the strongest sources of gravitational waves for ground-based interferometers as well as the space-based interferometer LISA. Recent progress in numerical relativity now makes it possible to calculate the waveforms from the strong-field dynamical merger and is revolutionizing our understanding of these systems. We review these dramatic developments, emphasizing applications to issues in gravitational wave observations. These new capabilities also make possible accurate calculations of the recoil or kick imparted to the final remnant black hole when the merging components have unequal masses, or unequal or unaligned spins. We highlight recent work in this area, focusing on results of interest to astrophysics.

Models of late-time neutrino mass generation contain new interactions of the cosmic background neutrinos with supernova relic neutrinos (SRNs) through exchange of the on-shell light boson, leading to significant modification of the differential SRN flux observed at earth. We consider Abelian U(1) model for generating neutrino masses at low scales and we show that there is a large parameter space in this model for which the changes induced in the flux by the exchange of the light bosons might allow one to distinguish between neutrinos being Majorana or Dirac particles, the type of neutrino mass hierarchy (normal or inverted or quasi-degenerate), and could also possibly determine the absolute values of the neutrino masses. Measurements of the presence of these effects would be possible at the next-generation water Cerenkov detectors enriched with Gadolinium, or a large 100 kton liquid argon detector.

This a particularly exciting time for gravitational wave physics. Ground-based gravitational wave detectors are now operating at a sensitivity such that gravitational radiation may soon be directly detected, and recently several groups have independently made significant breakthroughs that have finally enabled numerical relativists to solve the Einstein field equations for coalescing black-hole binaries, a key source of gravitational radiation. The numerical relativity community is now in the position to begin providing simulated merger waveforms for use by the data analysis community, and it is therefore very important that we provide ways to validate the results produced by various numerical approaches. Here, we present a simple comparison of the waveforms produced by two very different, but equally successful approaches--the generalized harmonic gauge and the moving puncture methods. We compare waveforms of equal-mass black hole mergers with minimal or vanishing spins. The results show exceptional agreement for the final burst of radiation, with some differences attributable to small spins on the black holes in one case.

Recent developments in numerical relativity have made it possible to follow reliably the coalescence of two black holes from near the innermost stable circular orbit to final ringdown. This opens up a wide variety of exciting astrophysical applications of these simulations. Chief among these is the net kick received when two unequal mass or spinning black holes merge. The magnitude of this kick has bearing on the production and growth of supermassive black holes during the epoch of structure formation, and on the retention of black holes in stellar clusters. Here we report the first accurate numerical calculation of this kick, for two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on analytic considerations to give a significant fraction of the maximum possible recoil. We have performed multiple runs with different initial separations, orbital angular momenta, resolutions, extraction radii, and gauges. The full range of our kick speeds is 86--116 km s$^{-1}$, and the most reliable runs give kicks between 86 and 97 km s$^{-1}$. This is intermediate between the estimates from two recent post-Newtonian analyses and suggests that at redshifts $z\gtrsim 10$, halos with masses $\lesssim 10^9 M_\odot$ will have difficulty retaining coalesced black holes after major mergers.

We study dynamics and radiation generation in the last few orbits and merger of a binary black hole system, applying recently developed techniques for simulations of moving black holes. Our analysis of the gravitational radiation waveforms and dynamical black hole trajectories produces a consistent picture for a set of simulations with black holes beginning on circular-orbit trajectories at a variety of initial separations. We find profound agreement at the level of one percent among the simulations for the last orbit, merger and ringdown. We are confident that this part of our waveform result accurately represents the predictions from Einstein's General Relativity for the final burst of gravitational radiation resulting from the merger of an astrophysical system of equal-mass non-spinning black holes. The simulations result in a final black hole with spin parameter a/m=0.69. We also find good agreement at a level of roughly 10 percent for the radiation generated in the preceding few orbits.

We present new techniqes for evolving binary black hole systems which allow the accurate determination of gravitational waveforms directly from the wave zone region of the numerical simulations. Rather than excising the black hole interiors, our approach follows the "puncture" treatment of black holes, but utilzing a new gauge condition which allows the black holes to move successfully through the computational domain. We apply these techniques to an inspiraling binary, modeling the radiation generated during the final plunge and ringdown. We demonstrate convergence of the waveforms and good conservation of mass-energy, with just over 3% of the system's mass converted to gravitional radiation.

We present a detailed description of techniques developed to combine 3D numerical simulations and, subsequently, a single black hole close-limit approximation. This method has made it possible to compute the first complete waveforms covering the post-orbital dynamics of a binary black hole system with the numerical simulation covering the essential non-linear interaction before the close limit becomes applicable for the late time dynamics. To determine when close-limit perturbation theory is applicable we apply a combination of invariant a priori estimates and a posteriori consistency checks of the robustness of our results against exchange of linear and non-linear treatments near the interface. Once the numerically modeled binary system reaches a regime that can be treated as perturbations of the Kerr spacetime, we must approximately relate the numerical coordinates to the perturbative background coordinates. We also perform a rotation of a numerically defined tetrad to asymptotically reproduce the tetrad required in the perturbative treatment. We can then produce numerical Cauchy data for the close-limit evolution in the form of the Weyl scalar $\psi_4$ and its time derivative $\partial_t\psi_4$ with both objects being first order coordinate and tetrad invariant. The Teukolsky equation in Boyer-Lindquist coordinates is adopted to further continue the evolution. To illustrate the application of these techniques we evolve a single Kerr hole and compute the spurious radiation as a measure of the error of the whole procedure. We also briefly discuss the extension of the project to make use of improved full numerical evolutions and outline the approach to a full understanding of astrophysical black hole binary systems which we can now pursue.

LISA will be able to detect gravitational waves from inspiralling massive black hole (MBH) binaries out to redshifts z > 10. If the binary masses and luminosity distances can be extracted from the LISA data stream, this information can be used to reveal the merger history of MBH binaries and their host galaxies in the evolving universe. Since this parameter extraction generally requires that LISA observe the inspiral for a significant fraction of its yearly orbit, carrying out this program requires adequate sensitivity at low frequencies, f < 10^{-4} Hz. Using several candidate low frequency sensitivities, we examine LISA's potential for characterizing MBH binary coalescences at redshifts z > 1.

We have scanned the fields of six radio-loud quasars using the Taurus Tunable Filter to detect redshifted [OII] 3727 line-emitting galaxies at redshifts 0.8 < z < 1.3. Forty-seven new emission-line galaxy (ELG) candidates are found. This number corresponds to an average space density about 100 times that found locally and, at L([OII]) < 10^{42} erg s^{-1} cm^{-2}, is 2 - 5 times greater than the field ELG density at similar redshifts, implying that radio-loud quasars inhabit sites of above-average star formation activity. The implied star-formation rates are consistent with surveys of field galaxies at z ~ 1. However, the variation in candidate density between fields is large and indicative of a range of environments, from the field to rich clusters. The ELG candidates also cluster -- both spatially and in terms of velocity -- about the radio sources. In fields known to contain rich galaxy clusters, the ELGs lie at the edges and outside the concentrated cores of red, evolved galaxies, consistent with the morphology-density relation seen in low-redshift clusters. This work, combined with other studies, suggests that the ELG environments of powerful AGN look very much the same from moderate to high redshifts, i.e. 0.8 < z < 4.

We compute the gravitational radiation generated in the evolution of a family of close binary black hole configurations, using a combination of numerical and perturbative approximation methods. We evolve the binaries with spins, $s$, aligned or counter-aligned with the orbital angular momentum from near the innermost stable circular orbit (ISCO) down to the final single rotating black hole. For the moderately spinning holes studied here the remnant Kerr black hole formed at the end of an inspiral process have a rotation parameter $a/M\approx0.72+0.32(s/m_H)$, suggesting it is difficult (though not excluded) to end up with near maximally rotating holes from such scenarios.

We have carried out multicolour imaging of a complete sample of radio-loud quasars at 0.6 < z < 1.1 and find groups or clusters of galaxies in the fields of at least 8 and possibly 13 of the 21 sources. There is no evidence for an evolution in the richness of the environments of radio-loud quasars from other low-redshift studies to z >~ 0.9. The quasars associated with groups and clusters in our sample do not necessarily reside in the centre of the galaxy distribution which rarely displays a spherical geometry. Clustering is preferentially associated with small or asymmetric steep-spectrum radio sources. The quasars with the largest projected angular size are, in nearly all cases, found in non-clustered environments. Radio-based selection (including source size) of high-redshift groups and clusters can be a very efficient method of detecting rich environments at these redshifts. We find that in optical searches for galaxy overdensities above z ~ 0.6 multiple filters must be used. If the single-filter counting statistics used by groups at lower redshift are applied to our data, uncertainties are too large to make accurate quantifications of cluster richness. This means that genuine clustering of galaxies about quasars will be missed and, in ~10% of cases, putative clusters turn out to be false detections. The statistics are further diluted by the fact that galaxy overdensities are generally not centred on the quasar.

With the goal of bringing theory, particularly numerical relativity, to bear on an astrophysical problem of critical interest to gravitational wave observers we introduce a model for coalescence radiation from binary black hole systems. We build our model using the "Lazarus approach", a technique that bridges far and close limit approaches with full numerical relativity to solve Einstein equations applied in the truly nonlinear dynamical regime. We specifically study the post-orbital radiation from a system of equal-mass non-spinning black holes, deriving waveforms which indicate strongly circularly polarized radiation of roughly 3% of the system's total energy and 12% of its total angular momentum in just a few cycles. Supporting this result we first establish the reliability of the late-time part of our model, including the numerical relativity and close-limit components, with a thorough study of waveforms from a sequence of black hole configurations varying from previously treated head-on collisions to representative target for ``ISCO'' data corresponding to the end of the inspiral period. We then complete our model with a simple treatment for the early part of the spacetime based on a standard family of initial data for binary black holes in circular orbit. A detailed analysis shows strong robustness in the results as the initial separation of the black holes is increased from 5.0 to 7.8M supporting our waveforms as a suitable basic description of the astrophysical radiation from this system. Finally, a simple fitting of the plunge waveforms is introduced as a first attempt to facilitate the task of analyzing data from gravitational wave detectors.

We assess the global properties of associated C IV 1548,1550 absorption lines measured in the spectra of radio-loud quasars drawn from a near-complete, low-frequency selected sample. The observations span restframe C IV in two redshift ranges -- 0.7<z<1.0 and 1.5<z<3.0 -- which were targetted in the UV with HST/STIS and in the optical with ground-based telescopes, respectively. First, we corroborate trends for C IV associated absorption to be found preferentially in steep-spectrum and lobe-dominated quasars, implying the absorbing material tends to lie away from the radio-jet axis. Furthermore, we find a clear anticorrelation between C IV absorption strength and the projected linear size of steep-spectrum quasars, indicative of an evolutionary sequence. We also find that heavily-absorbed quasars are systematically redder, implying dust is mixed in with the C IV-absorbing gas. No redshift dependence was found in any of the trends considered. These new results show that radio sources are triggered in galaxies which are exceptionally rich in gas and dust, which then dissipates on a timescale comparable with, but less than, that of the radio source. This observational sequence, together with the lack of redshift-dependence, points to a direct causal link between the event which triggered the radio source and the build-up of absorbing gas and dust, whose make-up is tantalisingly similar to the products of a post-merger starburst. Thus, these new results provide direct evidence for the clearing of absorbing material around quasars with time, as well as the probable association of starburst activity with the onset of the radio activity in AGN.

We study the coalescence of non-spinning binary black holes from near the innermost stable circular orbit down to the final single rotating black hole. We use a technique that combines the full numerical approach to solve Einstein equations, applied in the truly non-linear regime, and linearized perturbation theory around the final distorted single black hole at later times. We compute the plunge waveforms which present a non negligible signal lasting for $t\sim 100M$ showing early non-linear ringing, and we obtain estimates for the total gravitational energy and angular momentum radiated.

We have carried out multi-colour imaging of the fields of a statistically complete sample of low-frequency selected radio loud quasars at 0.6<z<1.1, in order to determine the characteristics of their environments. The largest radio sources are located in the field, and smaller steep-spectrum sources are more likely to be found in richer environments, from compact groups through to clusters. This radio-based selection (including source size) of high redshift groups and clusters is a highly efficient method of detecting rich environments at these redshifts. Although our single filter clustering measures agree with those of other workers, we show that these statistics cannot be used reliably on fields individually, colour information is required for this.

Results are presented from a survey for CIV associated absorption in a complete low-frequency-selected sample of quasars. In agreement with previous work, associated absorbers are most common in steep-spectrum and lobe-dominated quasars, indicative of an anisotropic cloud distribution. Furthermore, we find the strongest CIV absorption occurs in sources of small radio size, suggesting that the absorbing clouds are destroyed or displaced as radio sources expand. Evidence for dust in the clouds is also found, such that quasars with strong absorption are systematically redder. Finally, we find no evidence for evolution in the frequency or properties of the absorbers from z~0.7 to z~3.

We have obtained multi-colour imaging of a representative, statistically complete sample of low-frequency selected (S_408MHz > 0.95Jy) radio loud quasars at intermediate (0.6 < z < 1.1) redshifts. These sources are found in a variety of environments, from the field through to rich clusters. We show that statistical measures of environmental richness, based upon single-band observations are inadequate at these redshifts for a variety of reasons. Environmental richness seems correlated with the size and morphology of the radio source, as expected if the energy density in the radio lobes is approximately the equipartition value and the lobes are in pressure equilbrium with a surrounding intragroup/cluster medium. Selecting on radio size therefore efficiently selects dense galactic sytems at these redshifts.

Using a combination of multicolour broad- and narrow-band imaging techniques and follow-up spectroscopy, we have detected an overdensity of galaxies in the field of quasar MRC B0450-221, whose properties are consistent with a cluster at the quasar redshift z=0.9. An excess of red galaxies (V-I>2.2, I-K'>3.8) is evident within 1' of the quasar, with the colours expected for galaxies at z=0.9 that have evolved passively for 3 Gyr or more. A number of line-emitting galaxies (nine candidates with equivalent widths EW>70A) are also detected in the field using the TAURUS Tunable Filter (TTF). Three have been confirmed spectroscopically to indeed lie at z=0.9. The TTF candidates with the strongest [O II] line emission cluster in a group which lies 200-700 kpc away from the quasar and the red galaxy excess, and therefore most likely on the outskirts of the cluster. These observations are the first in a series probing quasar environments at z~1 with TTF.

We present the first results in a new program intended to make the best use of all available technologies to provide an effective understanding of waves from inspiralling black hole binaries in time for imminent observations. In particular, we address the problem of combining the close-limit approximation describing ringing black holes and full numerical relativity, required for essentially nonlinear interactions. We demonstrate the effectiveness of our approach using general methods for a model problem, the head-on collision of black holes. Our method allows a more direct physical understanding of these collisions indicating clearly when non-linear methods are important. The success of this method supports our expectation that this unified approach will be able to provide astrophysically relevant results for black hole binaries in time to assist gravitational wave observations.

We compare the fully nonlinear and perturbative evolution of nonrotating black holes with odd-parity distortions utilizing the perturbative results to interpret the nonlinear results. This introduction of the second polarization (odd-parity) mode of the system, and the systematic use of combined techniques brings us closer to the goal of studying more complicated systems like distorted, rotating black holes, such as those formed in the final inspiral stage of two black holes. The nonlinear evolutions are performed with the 3D parallel code for Numerical Relativity, {Cactus}, and an independent axisymmetric code, {Magor}. The linearized calculation is performed in two ways: (a) We treat the system as a metric perturbation on Schwarzschild, using the Regge-Wheeler equation to obtain the waveforms produced. (b) We treat the system as a curvature perturbation of a Kerr black hole (but here restricted to the case of vanishing rotation parameter a) and evolve it with the Teukolsky equation The comparisons of the waveforms obtained show an excellent agreement in all cases.

We have carried out a study to simulate distant clusters of galaxies in deep ground-based optical images. We find that when model galaxies are added to deep images obtained with the William Herschel Telescope, there is considerable scatter of the recovered galaxy colours away from the model values; this scatter is larger than that expected from photometric errors and is significantly affected by confusion, due to ground-based seeing, between objects in the field. In typical conditions of $\approx$ 1-arcsec seeing, the combination of confusion and incompleteness causes a considerable underestimation of the true surface density of $z \approx 3$ galaxies. We argue that the actual surface density of $z \approx 3$ galaxies may be several times greater than that estimated by previous ground-based studies, consistent with the surface density of such objects found in the HDF.

We describe a new statistic for measuring the small-scale velocity dispersion of galaxies directly from redshift surveys. This statistic is based on the object-weighted statistic proposed by Davis, Miller, & White (1997). Compared with the traditional pair-weighted velocity dispersion, our statistic is less sensitive to the presence or absence of rare, rich clusters of galaxies. This measure of the thermal energy of the galaxy distribution is ideally suited for use with a filtered version of the cosmic energy equation. We discuss the application of the statistic to the Las Campanas Redshift Survey. The low observed dispersion strongly favors cosmological models with low matter density, Omega_m ~ 0.2.

The pair-weighted relative velocity dispersion of galaxies provides a measure of the thermal energy of fluctuations of the observed galaxy distribution, but the measure is difficult to interpret and is very sensitive to the existence of rare, rich clusters of galaxies. Several alternative statistical procedures have recently been suggested to relieve these problems. We apply a variant of the object-weighted statistical method of Davis, Miller, & White (1997) to the Las Campanas Redshift Survey (LCRS), which is the largest and deepest existing redshift survey that is nearly fully sampled. The derived one-dimensional dispersion on scales ~ 1 Mpc/h is quite low: sigma_1 = 126 +/- 10 km/s, with a modest decrease at larger scales. The statistic is very stable; the six independent slices of the LCRS all yield consistent results. We apply the same statistical procedure to halos in numerical simulations of an open cosmological model and flat models with and without a cosmological constant. In contrast to the LCRS, all the models show a dispersion which increases for scales > 1 Mpc/h; it is uncertain whether this is a numerical artifact or a real physical effect. The standard cluster-normalized Cold Dark Matter model with Omega_m = 1 as well as a tilted variant with n = 0.8 yield dispersions substantially hotter than the LCRS value, while models with low matter density (Omega_m = 0.3) are broadly consistent with the LCRS data. Using a filtered cosmic energy equation, we measure Omega_m ~ 0.2, with small-scale bias factors b = 1.0--1.5 for high-density models and b = 0.7--1.1 for low-density models.

Using several approximations, we calculate an estimate of the gravitational radiation emitted when two equal mass black holes coalesce at the end of their binary inspiral. We find that about 1% of the mass energy of the pair will emerge as gravitational waves during the final ringdown and a negligible fraction of the angular momentum will be radiated.

We carry out N-body simulations of several non-Gaussian structure formation models, including Peebles' isocurvature cold dark matter model, cosmic string models, and a model with primordial voids. We compare the evolution of the cluster mass function in these simulations with that predicted by a modified version of the Press-Schechter formalism. We find that the Press-Schechter formula can accurately fit the cluster evolution over a wide range of redshifts for all of the models considered, with typical errors in the mass function of less than 25%, considerably smaller than the amount by which predictions for different models may differ. This work demonstrates that the Press-Schechter formalism can be used to place strong model independent constraints on non-Gaussianity in the universe.

We describe observations at frequencies near 15 GHz of the second 2x2 degree field imaged with the Cambridge Cosmic Anisotropy Telescope (CAT). After the removal of discrete radio sources, structure is detected in the images on characteristic scales of about half a degree, corresponding to spherical harmonic multipoles in the approximate range l= 330--680. A Bayesian analysis confirms that the signal arises predominantly from the cosmic microwave background (CMB) radiation for multipoles in the lower half of this range; the average broad-band power in a bin with centroid l=422 (theta = 51') is estimated to be Delta_T/T=2.1^{+0.4}_{-0.5} x 10^{-5}. For multipoles centred on l=615 (theta =35'), we find contamination from Galactic emission is significant, and constrain the CMB contribution to the measured power in this bin to be Delta_T/T <2.0 x 10^{-5} (1-sigma upper limit). These new results are consistent with the first detection made by CAT in a completely different area of sky. Together with data from other experiments, this new CAT detection adds weight to earlier evidence from CAT for a downturn in the CMB power spectrum on scales smaller than 1 degree. Improved limits on the values of H_0 and Omega are determined using the new CAT data.

The Wavelength-Oriented Microwave Background Analysis Team (WOMBAT) is constructing microwave maps which will be more realistic than previous simulations. Our foreground models represent a considerable improvement: where spatial templates are available for a given foreground, we predict the flux and spectral index of that component at each place on the sky and estimate uncertainties. We will produce maps containing simulated CMB anisotropy combined with expected foregrounds. The simulated maps will be provided to the community as the WOMBAT Challenge, so such maps can be analyzed to extract cosmological parameters by scientists who are unaware of their input values. This will test the efficacy of foreground subtraction, power spectrum analysis, and parameter estimation techniques and help identify the areas most in need of progress. These maps are also part of the FORECAST project, which allows web-based access to the known foreground maps for the planning of CMB missions.

Optical spectroscopic data are presented here for quasars from the Molonglo Quasar Sample (MQS), which forms part of a complete survey of 1-Jy radio sources from the Molonglo Reference Catalogue. The combination of low-frequency selection and complete identifications means that the MQS is relatively free from the orientation biases which affect most other quasar samples. To date, the sample includes 105 quasars and 6 BL Lac objects, 106 of which have now been confirmed spectroscopically. This paper presents a homogenous set of low-resolution optical spectra for 79 MQS quasars, the majority of which have been obtained at the Anglo-Australian Telescope. Full observational details are given and redshifts, continuum and emission-line data tabulated for all confirmed quasars.

The Wavelength-Oriented Microwave Background Analysis Team (WOMBAT) is constructing microwave skymaps which will be more realistic than previous simulations. Our foreground models represent a considerable improvement: where spatial templates are available for a given foreground, we predict the flux and spectral index of that component at each place on the sky and estimate the uncertainties in these quantities. We will produce maps containing simulated Cosmic Microwave Background anisotropies combined with all major expected foreground components. The simulated maps will be provided to the cosmology community as the WOMBAT Challenge, a "hounds and hares" exercise where such maps can be analyzed to extract cosmological parameters by scientists who are unaware of their input values. This exercise will test the efficacy of current foreground subtraction, power spectrum analysis, and parameter estimation techniques and will help identify the areas most in need of progress.

We present deep optical images of the PC1643+4631 field obtained at the WHT. This field contains two quasars at redshifts z=3.79 & 3.83 and a cosmic microwave background (CMB) decrement detected with the Ryle Telescope. The images are in U,G,V,R and I filters, and are complete to 25th magnitude in R and G and to 25.5 in U. The isophotal galaxy counts are consistent with the results of Metcalde et al. (1996), Hogg et al. (1997), and others. We find an excess of robust high-redshift Ly-break galaxy candidates with 25.0<R<25.5 compared with the mean number found in the fields studied by Steidel et al. -we expect 7 but find 16 - but we do not find that the galaxies are concentrated in the direction of the CMB decrement. However, we are still not sure of the distance to the system causing the CMB decrement. We have also used our images to compare the commonly used object-finding algorithms of FOCAS and SExtractor: we find FOCAS the more efficient at detecting faint objects and the better at dealing with composite objects, whereas SExtractor's morphological classification is more reliable, especially for faint objects near the resolution limit. More generally, we have also compared the flux lost using isophotal apertures on a real image with that on a noise-only image: recovery of artificial galaxies from the noise-only image significantly overestimates the flux lost from the galaxies, and we find that the corrections made using this technique suffer a systematic error of some 0.4 magnitudes.

In an investigation of the cause of the cosmic microwave background decrement in the field of the z = 3.8 quasar pair PC1643+4631, we have carried out a study to photometrically estimate the redshifts of galaxies in deep multi-colour optical images of the field taken with the WHT. To examine the possibility that a massive cluster of galaxies lies in the field, we have attempted to recover simulated galaxies with intrinsic colours matching those of the model galaxies used in the photometric redshift estimation. We find that when such model galaxies are added to our images, there is considerable scatter of the recovered galaxy redshifts away from the model value; this scatter is larger than that expected from photometric errors and is the result of confusion, simply due to ground-based seeing, between objects in the field. We have also compared the likely efficiency of the photometric redshift technique against the colour criteria used to select z>3 galaxies via the strong colour signature of the Lyman-limit break. We find that these techniques may significantly underestimate the true surface density of z>3, due to confusion between the high-redshift galaxies and other objects near the line of sight. We argue that the actual surface density of z=3 galaxies may be as much as 6 times greater than that estimated by previous ground-based studies, and note that this conclusion is consistent with the surface density of high-redshift objects found in the HDF. Finally, we conclude that all ground-based deep field surveys are inevitably affected by confusion, and note that reducing the effective seeing in ground-based images will be of paramount importance in observing the distant universe.

We have used the Optical Redshift Survey (ORS; Santiago et al. 1995) to construct the gravity field due to fluctuations in the galaxy density field out to distances of 8000 km/s. At large scales where linear theory applies, the comparison of this gravity field with the observed peculiar velocity field offers a powerful cosmological probe, because the predicted flow field is proportional to the parameter Omega^{0.6}/b, where Omega is the matter density and b is the bias of the galaxy distribution. The more densely sampled ORS gravity field, to excellent approximation, matches that of the earlier IRAS 1.2-Jy redshift survey (Fisher et al. 1995), provided beta is reduced by a factor b(opt)/b(IRAS)=1.4. Apart from this scaling, the most significant difference between the ORS and IRAS fields is induced by differing estimates of the over-density of the Virgo cluster. Neither of these gravity fields is consistent with the peculiar velocity field constructed from the full Mark III (Willick et al. 1997) sample. We find that a simple but plausible non-linear bias algorithm for the galaxy distribution relative to the mass has a negligible effect on the derived fields. We conclude that the substitution of optical for IRAS catalogues cannot alone resolve the discrepancies between the IRAS gravity field and the Mark III peculiar velocity field.

We solve the problem of expressing the Weyl scalars $\psi $ that describe gravitational perturbations of a Kerr black hole in terms of Cauchy data. To do so we use geometrical identities (like the Gauss-Codazzi relations) as well as Einstein equations. We are able to explicitly express $\psi $ and $\partial _t\psi $ as functions only of the extrinsic curvature and the three-metric (and geometrical objects built out of it) of a generic spacelike slice of the spacetime. These results provide the link between initial data and $\psi $ to be evolved by the Teukolsky equation, and can be used to compute the gravitational radiation generated by two orbiting black holes in the close limit approximation. They can also be used to extract waveforms from spacetimes completely generated by numerical methods.

Preliminary results are described for emission-line imaging of the field of a quasar at z=0.9 with a unique new instrument --- the TAURUS Tunable Filter (TTF). At least fourteen [O II]-emitting galaxy candidates are found within +/- 750 km/s of the quasar redshift with a significance of >3 sigma. Another eight candidates (>3 sigma) are also found with relative velocities +/- 750-1500 km/s from the quasar, indicative of a large velocity dispersion in the field. Together with the existence of a group of red galaxies in the field identified through broad-band imaging, we suggest that the TTF data are probing cluster merging which is ongoing at z=0.9.

Spectroscopy of high-redshift compact steep-spectrum (CSS) quasars shows a high incidence of Ly-alpha absorption close to the quasar redshift. Associated absorption systems like these are rare in other types of radio sources, with the notable exception of high-redshift radio galaxies (HZRGs). CSS quasars and HZRGs share many properties; the hypothesis that they are intrinsically similar objects is presented and discussed.

Using complete samples of steep-spectrum quasars, we present evidence for a correlation between radio and optical luminosity which is not caused by selection effects, nor caused by an orientation dependence (such as relativistic beaming), nor a byproduct of cosmic evolution. We argue that this rules out models of jet formation in which there are no parameters in common with the production of the optical continuum. This is arguably the most direct evidence to date for a close link between accretion onto a black hole and the fuelling of relativistic jets. The correlation also provides a natural explanation for the presence of aligned optical/radio structures in only the most radio luminous high-redshift galaxies.

We compare the peculiar velocities of nearby SNe Ia with those predicted by the gravity fields of full sky galaxy catalogs. The method provides a powerful test of the gravitational instability paradigm and strong constraints on the density parameter beta = Omega^0.6/b. For 24 SNe Ia within 10,000 km/s we find the observed SNe Ia peculiar velocities are well modeled by the predictions derived from the 1.2 Jy IRAS survey and the Optical Redshift Survey (ORS). Our best $\beta$ is 0.4 from IRAS, and 0.3 from the ORS, with beta>0.7 and beta<0.15 ruled out at 95% confidence levels from the IRAS comparison. Bootstrap resampling tests show these results to be robust in the mean and in its error. The precision of this technique will improve as additional nearby SNe Ia are discovered and monitored.

This paper reports the measurement of the Sunyaev-Zeldovich effect (SZE) in the nearby clusters A478, A2142 and A2256 at 32 GHz. These clusters are part of a complete X-ray flux limited sample, which includes the previously measured Coma cluster. We find surface baryonic mass densities for the hot IGM within the 7.35' FWHM Gaussian beam of the OVRO 5.5-m telescope projected on the cluster centers. This is used to place lower limits on the ratio of the SZE determined baryonic mass to the gravitational binding mass. When combined with available X-ray models, the data from these three clusters and Coma yields a sample average value for the Hubble constant of 54 +/- 14 km/s/Mpc.

In a companion paper (Paper I: Jones et al. 1996) we report the discovery of a cosmic microwave background decrement, indicative of a distant cluster with mass ~10^{15} solar masses, towards the quasar pair PC1643+4631A&B (z=3.79,3.83, separation 198''). To search for the cluster responsible, we have obtained R-, J- and K-band images of the field and have also carried out optical spectroscopy of selected objects in it. No such cluster is evident in these images. Assuming the cluster causing the decrement is similar to massive clusters already known, our magnitude limits imply that it must lie about or beyond z=1. This provides independent support for the X-ray-based distance argument of Paper I. The cluster must gravitationally lens objects behind it; for a cluster z around 1-2, the Einstein-ring radius for sources at z ~= 3.8 is ~ 100''. Simple modelling, producing simultaneously the S-Z effect and the lensing, shows that the source positions of quasars A and B lie within ~10'' of each other and may indeed be coincident. The two quasar spectra are found to be remarkably similar apart from their one-percent redshift difference. Assuming A and B are images of a single quasar, we present a possible explanation of this difference.

In a 15-GHz Ryle-Telescope observation of PC1643+4631 A & B, a pair of quasars at redshifts z = 3.79 and 3.83 separated by 198 arcsec on the sky, we find a decrement in the cosmic microwave background (CMB) of -380 +/- 64 uJy in a 110 arcsec x 175 arcsec beam. Assuming this to be a Sunyaev-Zel'dovich effect due to an intervening cluster, the minimum magnitude of the central temperature decrement is 560 uK. A serendipitous ROSAT observation shows that there is no X-ray-luminous cluster in the direction of the decrement at z < 1. The implied gas mass is >~ 2 x 10^{14} solar masses (assuming a temperature of ~ 5 keV), indicating a total mass of > 10^{15} solar masses. This result demonstrates the existence of a massive system too distant to be detected by its emission, but which can be found via its imprint on the CMB.

The orientation-dependence of the optical continuum emission in radio-loud quasars is investigated using a new, complete sample of low-frequency-selected quasars, the Molonglo Quasar Sample (MQS). The optical continuum is found to be highly anisotropic, brightening continuously from lobe- to core-dominated quasars by 3-5 mag. It is argued that aspect-dependent extinction, rather than relativistic boosting as has been previously proposed, provides the simplest explanation consistent with the data. The reddening hypothesis is supported by both the steeper optical slopes and the larger Balmer decrements found in lobe-dominated quasars, as well as the stronger anisotropy seen at blue wavelengths. The dust responsible is shown to be physically associated with the quasar, lying mostly at radii between the broad and narrow-line regions in a clumpy distribution. Such a geometry is reminiscent of a torus. However, substantial numbers of dust clouds must lie within the torus opening angle, contributing to an increasing average optical depth with increasing viewing angle away from the jet axis.