The recent detection of gravitational waves has generated interest in alternatives to the black hole interpretation of sources. One set of such alternatives involves a prediction of gravitational wave "echoes". We consider two aspects of possible echoes: First, general features of echoes coming from spacetime reflecting conditions. We find that the detailed nature of such echoes does not bear any clear relationship to quasi-normal frequencies. Second, we point out the pitfalls in the analysis of local reflecting "walls" near the horizon of rapidly rotating black holes.

Modelling of gravitational waves from binary black hole inspiral has played an important role in the recent observations of such signals. The late-stage ringdown phase of the gravitational waveform is often associated with the null particle orbit ("light ring") of the black hole spacetime. With simple models we show that this link between the light ring and spacetime ringing is based more on the history of specific models than on an actual constraining relationship. We also show, in particular, that a better understanding of the dissociation of the two may be relevant to the astrophysically interesting case of rotating (Kerr) black holes.

Modeling of gravitational waves (GWs) from binary black hole inspiral brings together early post-Newtonian waveforms and late quasinormal ringing waveforms. Attempts to bridge the two limits without recourse to numerical relativity involve predicting the time of the peak GW amplitude. This prediction will require solving the question of why the peak of the "source," i.e., the peak of the binary angular velocity, does not correspond to the peak of the GW amplitude. We show here that this offset can be understood as due to the existence two distinct components of the radiation: the "direct" radiation analogous to that in flat spacetime, and "scattered" radiation associated with curved spacetime. The time dependence of these two components, and of their relative phases determines the location of the peak amplitude. We use a highly simplified model to clarify the twocomponent nature of the source, then demonstrate that the explanation is valid also for an extreme mass ratio binary inspiral.

In this contribution Marcus Price gives a first hand account of some of the very first scientific observations undertaken with the "Giant Radio Telescope" at Parkes. It was clearly a very exciting time of discovery, enabled by superb engineering.

Some Galactic models predict a significant population of radio pulsars close to the our galactic center. Beams from these pulsars could get strongly deflected by the supermassive black hole (SMBH) believed to reside at the galactic center and reach the Earth. Earlier work assuming a Schwarzschild SMBH gave marginal chances of observing this exotic phenomenon with current telescopes and good chances with future telescopes. Here we calculate the odds of observability for a rotating SMBH. We find that the estimates of observation are not affected by the SMBH spin, but a pulsar timing analysis of deflected pulses might be able to provide an estimate of the spin of the central black hole.

In the field of gamma-ray astronomy, irregular and noisy datasets make difficult the characterization of light-curve features in terms of statistical significance while properly accounting for trial factors associated with the search for variability at different times and over different timescales. In order to address these difficulties, we propose a method based on the Haar wavelet decomposition of the data. It allows statistical characterization of possible variability, embedded in a white noise background, in terms of a confidence level. The method is applied to artificially generated data for characterization as well as to the the very high energy M87 light curve recorded with VERITAS in 2008 which serves here as a realistic application example.

According to some models, there may be a significant population of radio pulsars in the Galactic center. In principle, a beam from one of these pulsars could pass close to the supermassive black hole (SMBH) at the center, be deflected, and be detected by Earth telescopes. Such a configuration would be an unprecedented probe of the properties of spacetime in the moderate- to strong-field regime of the SMBH. We present here background on the problem, and approximations for the probability of detection of such beams. We conclude that detection is marginally probable with current telescopes, but that telescopes that will be operating in the near future, with an appropriate multiyear observational program, will have a good chance of detecting a beam deflected by the SMBH.

Experiments are in progress to prepare for intensity interferometry with arrays of air Cherenkov telescopes. At the Bonneville Seabase site, near Salt Lake City, a testbed observatory has been set up with two 3-m air Cherenkov telescopes on a 23-m baseline. Cameras are being constructed, with control electronics for either off- or online analysis of the data. At the Lund Observatory (Sweden), in Technion (Israel) and at the University of Utah (USA), laboratory intensity interferometers simulating stellar observations have been set up and experiments are in progress, using various analog and digital correlators, reaching 1.4 ns time resolution, to analyze signals from pairs of laboratory telescopes.

Massive gravitons are features of some alternatives to general relativity. This has motivated experiments and observations that, so far, have been consistent with the zero mass graviton of general relativity, but further tests will be valuable. A basis for new tests may be the high sensitivity gravitational wave experiments that are now being performed, and the higher sensitivity experiments that are being planned. In these experiments it should be feasible to detect low levels of dispersion due to nonzero graviton mass. One of the most promising techniques for such a detection may be the pulsar timing program that is sensitive to nano-Hertz gravitational waves. Here we present some details of such a detection scheme. The pulsar timing response to a gravitational wave background with the massive graviton is calculated, and the algorithm to detect the massive graviton is presented. We conclude that, with 90% probability, massles gravitons can be distinguished from gravitons heavier than $3\times 10^{-22}$ eV (Compton wave length $\lambda_{\rm g}=4.1 \times 10^{12}$ km), if biweekly observation of 60 pulsars are performed for 5 years with pulsar RMS timing accuracy of 100 ns. If 60 pulsars are observed for 10 years with the same accuracy, the detectable graviton mass is reduced to $5\times 10^{-23}$ eV ($\lambda_{\rm g}=2.5 \times 10^{13}$ km); for 5-year observations of 100 or 300 pulsars, the sensitivity is respectively $2.5\times 10^{-22}$ ($\lambda_{\rm g}=5.0\times 10^{12}$ km) and $10^{-22}$ eV ($\lambda_{\rm g}=1.2\times 10^{13}$ km). Finally, a 10-year observation of 300 pulsars with 100 ns timing accuracy would probe graviton masses down to $3\times 10^{-23}$ eV ($\lambda_{\rm g}=4.1\times 10^{13}$ km).

A pulsar beam passing close to a black hole can provide a probe of very strong gravitational fields even if the pulsar itself is not in a strong field region. In the case that the spin of the hole can be ignored, we have previously shown that all strong field effects on the beam can be understood in terms of two "universal" functions, $F(\phi_{\rm in})$ and $T(\phi_{\rm in})$ of the angle of beam emission $\phi_{\rm in}$; these functions are universal in that they depend only on a single parameter, the pulsar/black hole distance from which the beam is emitted. Here we apply this formalism to general pulsar-hole-observer geometries, with arbitrary alignment of the pulsar spin axis and arbitrary pulsar beam direction and angular width. We show that the analysis of the observational problem has two distinct elements: (i) the computation of the location and trajectory of an observer-dependent "keyhole" direction of emission in which a signal can be received by the observer; (ii) the determination of an annulus that represents the set of directions containing beam energy. Examples of each are given along with an example of a specific observational scenario.

For an equation of state in which pressure is a function only of density, the analysis of Newtonian stellar structure is simple in principle if the system is axisymmetric, or consists of a corotating binary. It is then required only to solve two equations: one stating that the "injection energy", $\kappa$, a potential, is constant throughout the stellar fluid, and the other being the integral over the stellar fluid to give the gravitational potential. An iterative solution of these equations generally diverges if $\kappa$ is held fixed, but converges with other choices. We investigate the mathematical reason for this convergence/divergence by starting the iteration from an approximation that is perturbatively different from the actual solution. A cycle of iteration is then treated as a linear "updating" operator, and the properties of the linear operator, especially its spectrum, determine the convergence properties. For simplicity, we confine ourselves to spherically symmetric models in which we analyze updating operators both in the finite dimensional space corresponding to a finite difference representation of the problem, and in the continuum, and we find that the fixed-$\kappa$ operator is self-adjoint and generally has an eigenvalue greater than unity; in the particularly important case of a polytropic equation of state with index greater than unity, we prove that there must be such an eigenvalue. For fixed central density, on the other hand, we find that the updating operator has only a single eigenvector, with zero eigenvalue, and is nilpotent in finite dimension, thereby giving a convergent solution.

We analyze the effect of weak field gravitational waves on the timing of pulsars, with particular attention to gauge invariance, that is, to the effects that are independent of the choice of coordinates. We find: (i) the Doppler shift cannot be separated into gauge invariant gravitational wave and kinetic contributions; (ii) a gauge invariant separation can be made for the time derivative of the Doppler shift in which the gravitational wave contribution is directly related to the Riemann tensor, and the kinetic contribution is that for special relativity; (iii) the gauge dependent effects in the Doppler shift play no role in the program of gravitational wave detection via pulsar timing. The direct connection shown between pulsar timing and the Riemann tensor of the gravitational waves will be of importance in discussions of gravitational waves from alternative (non-Einsteinian) theories of gravitation.

If a pulsar orbits a supermassive black hole, the timing of pulses that pass close to the hole will show a variety of strong field effects. To compute the intensity and timing of pulses that have passed close to a nonrotating black hole we introduce here a simple formalism based on two "universal functions," one for the bending of photon trajectories and the other for the photon travel time on these trajectories. We apply this simple formalism to the case of a pulsar in circular orbit that beams its pulses into the orbital plane. In addition to the "primary" pulses that reach the receiver by a more-or-less direct path, we find that there are secondary and higher order pulses. These are usually much dimmer than the primary pulses, but they can be of comparable or even greater intensity if they are emitted when pulsar is on the side of the hole furthest from the receiver. We show that there is a phase relationship of the primary and secondary pulses that is a probe of the strongly curved spacetime geometry. Analogs of these phenomena are expected in more general configurations, in which a pulsar in orbit around a hole emits pulses that are not confined to the orbital plane.

We present the largest catalogue to date of optical counterparts for HI radio-selected galaxies, Hopcat. Of the 4315 HI radio-detected sources from the HI Parkes All Sky Survey (Hipass) catalogue, we find optical counterparts for 3618 (84%) galaxies. Of these, 1798 (42%) have confirmed optical velocities and 848 (20%) are single matches without confirmed velocities. Some galaxy matches are members of galaxy groups. From these multiple galaxy matches, 714 (16%) have confirmed optical velocities and a further 258 (6%) galaxies are without confirmed velocities. For 481 (11%), multiple galaxies are present but no single optical counterpart can be chosen and 216 (5%) have no obvious optical galaxy present. Most of these 'blank fields' are in crowded fields along the Galactic plane or have high extinctions. Isolated 'Dark galaxy' candidates are investigated using an extinction cut of ABj < 1 mag and the blank fields category. Of the 3692 galaxies with an ABj extinction < 1 mag, only 13 are also blank fields. Of these, 12 are eliminated either with follow-up Parkes observations or are in crowded fields. The remaining one has a low surface brightness optical counterpart. Hence, no isolated optically dark galaxies have been found within the limits of the Hipass survey.

Transient microstructure in the diffuse interstellar medium has been observed towards galactic and extragalactic sources for decades, usually in lines of atoms and ions, and, more recently, in molecular lines. Evidently, there is a molecular component to the transient microstructure. In this paper we explore the chemistry that may arise in such microstructure. We use a PDR code to model the conditions of relatively high density, low temperature, very low visual extinction and very short elapsed time that are appropriate for these objects. We find that there is a well defined region of parameter space where detectable abundances of molecular species might be found. The best matching models are those where the interstellar microstructure is young (<100 yr), small (~100 AU), and dense (>10^4 cm^-3).

We have placed limits on the cosmological significance of gas-rich low surface-brightness (LSB) galaxies as a proportion of the total population of gas-rich galaxies by carrying out a very deep survey (HIDEEP) for neutral hydrogen (HI) with the Parkes multibeam system. Such a survey avoids the surface-brightness selection effects that limit the usefulness of optical surveys for finding LSB galaxies. To complement the HIDEEP survey we have digitally stacked eight 1-hour R-band Tech Pan films from the UK Schmidt Telescope covering 36 square degrees of the survey area to reach a very deep isophotal limit of 26.5 R mag/sq. arcsec. At this level, we find that all of the 129 HI sources within this area have optical counterparts and that 107 of them can be identified with individual galaxies. We have used the properties of the galaxies identified as the optical counterparts of the HI sources to estimate the significance of LSB galaxies (defined to be those at least 1.5 magnitudes dimmer in effective surface-brightness than the peak in the observed distribution seen in optical surveys). We calculate the contribution of LSB galaxies to the total number, neutral hydrogen density, luminosity density, baryonic mass density, dynamical mass density and cross-sectional area of gas-rich galaxies. We do not find any `Crouching Giant' LSB galaxies such as Malin 1, nor do we find a population of extremely low surface-brightness galaxies not previously found by optical surveys. Such objects must either be rare, gas-poor or outside the survey detection limits.

We have carried out an extremely long integration-time (9000 s/beam) 21-cm blind survey of 60 square degrees in Centaurus using the Parkes multibeam system. We find that the noise continues to fall as sqrt(t) throughout, enabling us to reach an HI column-density limit of 4.2 x 10^18 cm^-2 for galaxies with a velocity width of 200 km/s in the central 32 square degree region, making this the deepest survey to date in terms of column density sensitivity. The HI data are complemented by very deep optical observations from digital stacking of multi-exposure UK Schmidt Telescope R-band films, which reach an isophotal level of 26.5 R mag arcsec^-2 (~27.5 B mag arcsec^-2). 173 HI sources have been found, 96 of which have been uniquely identified with optical counterparts in the overlap area. There is not a single source without an optical counterpart. Although we have not measured the column-densities directly, we have inferred them from the optical sizes of their counterparts. All appear to have a column-density of N_HI = 10^(20.65 +/- 0.38). This is at least an order of magnitude above our sensitivity limit, with a scatter only marginally larger than the errors on N_HI. This needs explaining. If confirmed it means that HI surveys will only find low surface brightness (LSB) galaxies with high M_HI/L_B. Gas-rich LSB galaxies with lower HI mass to light ratios do not exist. The paucity of low column-density galaxies also implies that no significant population will be missed by the all-sky HI surveys being carried out at Parkes and Jodrell Bank.

We present a new accurate measurement of the HI mass function of galaxies from the HIPASS Bright Galaxy Catalog, a sample of 1000 galaxies with the highest HI peak flux densities in the southern hemisphere (Koribalski et al. 2003). This sample spans nearly four orders of magnitude in HI mass (from log M_HI/M_sun=6.8 to 10.6, H0=75) and is the largest sample of HI selected galaxies to date. We develop a bivariate maximum likelihood technique to measure the space density of galaxies, and show that this is a robust method, insensitive to the effects of large scale structure. The resulting HI mass function can be fitted satisfactorily with a Schechter function with faint-end slope alpha=-1.30. This slope is found to be dependent on morphological type, with later type galaxies giving steeper slopes. We extensively test various effects that potentially bias the determination of the HI mass function, including peculiar motions of galaxies, large scale structure, selection bias, and inclination effects, and quantify these biases. The large sample of galaxies enables an accurate measurement of the cosmological mass density of neutral gas: Omega_HI=(3.8 +/- 0.6) x 10^{-4}. Low surface brightness galaxies contribute only 15% to this value, consistent with previous findings.

Using HIPASS data, we have searched for HI in a ~25x25 sq.deg. region centred on the Fornax cluster. Within a velocity search range of 300 - 3700 km/s and a lower flux limit of ~40 mJy, 110 galaxies with HI emission were detected, one of which is previously uncatalogued. None of the detections has early-type morphology. Previously unknown velocities for 14 galaxies have been determined, with a further 4 velocity measurements being significantly dissimilar to published values. Identification of an optical counterpart is relatively unambiguous for more than ~90% of our HI galaxies. The galaxies appear to be embedded in a sheet at the cluster velocity which extends for more than 30 deg across the search area. At the nominal cluster distance of ~20 Mpc, this corresponds to an elongated structure more than 10 Mpc in extent. A velocity gradient across the structure is detected, with radial velocities increasing by \~500 km/s from SE to NW. The clustering of galaxies evident in optical surveys is only weakly suggested in the spatial distribution of our HI detections. Our results suggest a considerable deficit of HI-rich galaxies in the centre of the cluster. However, relative to the field, there is a 3(+/-1)-fold excess of HI-rich galaxies in the outer parts of the cluster where galaxies may be infalling towards the cluster for the first time.

We report the discovery, from the HI Parkes All-Sky Survey (HIPASS), of a gas cloud associated with the asymmetric spiral galaxy NGC 2442. This object, designated HIPASS J0731-69, contains ~10^9 M_sun of HI, or nearly one-third as much atomic gas as NGC 2442 itself. No optical counterpart to any part of HIPASS J0731-69 has yet been identified, consistent with the gas being diffuse, and with its stream-like kinematics. If the gas in HIPASS J0731-69 was once part of NGC 2442, then it was most likely a fairly recent tidal encounter with a moderately massive companion which tore it loose, although the possibility of ram-pressure stripping cannot be ruled out. This discovery highlights the potential of the HIPASS data for yielding new clues to the nature of some of the best-known galaxies in the local universe.

We report the discovery, from the HI Parkes All-Sky Survey (HIPASS), of an isolated cloud of neutral hydrogen which we believe to be extragalactic. The HI mass of the cloud (HIPASS J1712-64) is very low, 1.7 x 10^7 Msun, using an estimated distance of ~3.2 Mpc. Most significantly, we have found no optical companion to this object to very faint limits (mu(B)~ 27 mag arcsec^-2). HIPASS J1712-64 appears to be a binary system similar to, but much less massive than, HI 1225+01 (the Virgo HI Cloud) and has a size of at least 15 kpc. The mean velocity dispersion, measured with the Australia Telescope Compact Array (ATCA), is only 4 km/s for the main component and because of the weak or non-existent star-formation, possibly reflects the thermal linewidth (T<2000 K) rather than bulk motion or turbulence. The peak column density for HIPASS J1712-64, from the combined Parkes and ATCA data, is only 3.5 x 10^19 cm^-2, which is estimated to be a factor of two below the critical threshold for star formation. Apart from its significantly higher velocity, the properties of HIPASS J1712-64 are similar to the recently recognised class of Compact High Velocity Clouds. We therefore consider the evidence for a Local Group or Galactic origin, although a more plausible alternative is that HIPASS J1712-64 was ejected from the interacting Magellanic Cloud/Galaxy system at perigalacticon ~ 2 x 10^8 yr ago.

A blind survey for HI bright galaxies in the southern Zone of Avoidance, (212 deg < l < 36 deg; |b| < 5 deg), has been made with the 21 cm multibeam receiver on the Parkes 64 m radiotelescope. The survey, sensitive to normal spiral galaxies to a distance of about 40 Mpc and more nearby dwarfs, detected 110 galaxies. Of these, 67 have no counterparts in the NASA/IPAC Extragalactic Database. In general, the uncataloged galaxies lie behind thicker obscuration than do the cataloged objects. All of the newly-discovered galaxies have HI flux integrals more than an order of magnitude lower than the Circinus galaxy. The survey recovers the Puppis cluster and foreground group (Kraan-Korteweg & Huchtmeier 1992), and the Local Void remains empty. The HI mass function derived for the sample is satisfactorily fit by a Schechter function with parameters alpha = 1.51 +- 0.12, Phi* = 0.006 +- 0.003, and log M* = 9.7 +- 0.10.

Using approximate techniques we study the final moments of the collision of two (individually non-spinnning) black holes which inspiral into each other. The approximation is based on treating the whole space-time as a single distorted black hole. We obtain estimates for the radiated energy, angular momentum and waveforms for the gravitational waves produced in such a collision. The results can be of interest for analyzing the data that will be forthcoming from gravitational wave interferometric detectors, like the LIGO, GEO, LISA, VIRGO and TAMA projects.

The gravitational radiation generated by a particle in a close unbounded orbit around a neutron star is computed as a means to study the importance of the $w$ modes of the neutron star. For simplicity, attention is restricted to odd parity (``axial'') modes which do not couple to the neutron star's fluid modes. We find that for realistic neutron star models, particles in unbounded orbits only weakly excite the $w$ modes; we conjecture that this is also the case for astrophysically interesting sources of neutron star perturbations. We also find that for cases in which there is significant excitation of quadrupole $w$ modes, there is comparable excitation of higher multipole modes.

We have commenced a 21-cm survey of the entire southern sky (\delta < 0 degrees, -1200 km/s < v < 12700 km/s) which is ''blind'', i.e. unbiased by previous optical information. In the present paper we report on the results of a pilot project which is based on data from this all-sky survey. The project was carried out on an area of 600 square degrees centred on the nearby Centaurus A (Cen A) group of galaxies at a mean velocity of v ~ 500 km/s. This was recently the subject of a separate and thorough optical survey. We found 10 new group members to add to the 21 galaxies already known in the Cen A group: five of these are previously uncatalogued galaxies, while five were previously catalogued but not known to be associated with the group. We found optical counterparts for all the HI detections, most of them intrinsically very faint low surface brightness dwarfs. The new group members add approximately 6% to the HI mass of the group and 4% to its light. The HI mass function, derived from all the known group galaxies in the interval 10^7 \Msun of HI to 10^9 \Msun of HI, has a faint-end slope of 1.30 +/- 0.15, allowing us to rule out a slope of 1.7 at 95% confidence. Even if the number in the lowest mass bin is increased by 50%, the slope only increases to 1.45 +/- 0.15.

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.

The perturbation theory of black holes has been useful recently for providing estimates of gravitational radiation from black hole collisions. Second order perturbation theory, relatively undeveloped until recently, has proved to be important both for providing refined estimates and for indicating the range of validity of perturbation theory. Here we present the second order formalism for perturbations of Schwarzschild spacetimes. The emphasis is on practical methods for carrying out second order computations of outgoing radiation. General issues are illustrated throughout with examples from ``close-limit'' results, perturbation calculations in which black holes start from small separation.

We present the first results of a blind HI survey for galaxies in the southern Zone of Avoidance with a multibeam receiver on the Parkes telescope. This survey is eventually expected to catalog several thousand galaxies within Galactic latitude |b|<5 degrees, mostly unrecognised before due to Galactic extinction and confusion. We present here results of the first three detections to have been imaged with the Australia Telescope Compact Array (ATCA). The galaxies all lie near Galactic longitude 325 degrees and were selected because of their large angular sizes, up to 1.3 degrees. Linear sizes range from 53 to 108 kpc. The first galaxy is a massive 5.7x10^11 solar mass disk galaxy with a faint optical counterpart, SGC 1511.1--5249. The second is probably an interacting group of galaxies straddling the Galactic equator. No optical identification is possible. The third object appears to be an interacting pair of low column density galaxies, possibly belonging to an extended Circinus or Centaurus A galaxy group. No optical counterpart has been seen despite the predicted extinction (A(B) = 2.7 - 4.4 mag) not being excessive. We discuss the implications of the results, in particular the low HI column densities (~10^19 atoms/sq.cm) found for two of the three galaxies.

Interactions between galaxies are common and are an important factor in determining their physical properties such as position along the Hubble sequence and star-formation rate. There are many possible galaxy interaction mechanisms, including merging, ram-pressure stripping, gas compression, gravitational interaction and cluster tides. The relative importance of these mechanisms is often not clear, as their strength depends on poorly known parameters such as the density, extent and nature of the massive dark halos that surround galaxies. A nearby example of a galaxy interaction where the mechanism is controversial is that between our own Galaxy and two of its neighbours -- the Large and Small Magellanic Clouds. Here we present the first results of a new HI survey which provides a spectacular view of this interaction. In addition to the previously known Magellanic Stream, which trails 100 degrees behind the Clouds, the new data reveal a counter-stream which lies in the opposite direction and leads the motion of the Clouds. This result supports the gravitational model in which leading and trailing streams are tidally torn from the body of the Magellanic Clouds.

We study the head-on collision of two equal-mass momentarily stationary black holes, using black hole perturbation theory up to second order. Compared to first-order results, this significantly improves agreement with numerically computed waveforms and energy. Much more important, second-order results correctly indicate the range of validity of perturbation theory. This use of second-order, to provide ``error bars,'' makes perturbation theory a viable tool for providing benchmarks for numerical relativity in more generic collisions and, in some range of collision parameters, for supplying waveform templates for gravitational wave detection.

We compute gravitational radiation waveforms, spectra and energies for a point particle of mass $m_0$ falling from rest at radius $r_0$ into a Schwarzschild hole of mass $M$. This radiation is found to lowest order in $(m_0/M)$ with the use of a Laplace transform. In contrast with numerical relativity results for head-on collisions of equal-mass holes, the radiated energy is found not to be a monotonically increasing function of initial separation; there is a local radiated-energy maximum at $r_0\approx4.5M$. The present results, along with results for infall from infinity, provide a complete catalog of waveforms and spectra for particle infall. We give a representative sample from that catalog and an interesting observation: Unlike the simple spectra for other head-on collisions (either of particle and hole, or of equal mass holes) the spectra for $\infty>r_0>\sim5M$ show a series of evenly spaced bumps. A simple explanation is given for this. Lastly, our energy vs. $r_0$ results are compared with approximation methods used elsewhere, for small and for large initial separation.

We study the even-parity $\ell=2$ perturbations of a Schwarzschild black hole to second order. The Einstein equations can be reduced to a single linear wave equation with a potential and a source term. The source term is quadratic in terms of the first order perturbations. This provides a formalism to address the validity of many first order calculations of interest in astrophysics.

The problem of the mutual attraction and joining of two black holes is of importance as both a source of gravitational waves and as a testbed of numerical relativity. If the holes start out close enough that they are initially surrounded by a common horizon, the problem can be viewed as a perturbation of a single black hole. We take initial data due to Misner for close black holes, apply perturbation theory and evolve the data with the Zerilli equation. The computed gravitational radiation agrees with and extends the results of full numerical computations.

Problem with the figures should be corrected. Apparently a broken uuencoder was the cause.

We consider a compact source of gravitational waves of frequency $\omega$, in or near a massive spherically symmetric distribution of matter or a black hole. Recent calculations have led to apparently contradictory results for the influence of the massive body on the propagation of the waves. We show here that the results are in fact consistent and in agreement with the ``standard'' viewpoint in which the high frequency compact source produces the radiation as if in a flat background, and the background curvature affects the propagation of these waves.