The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

We describe source code level parallelization for the {\tt kira} direct gravitational $N$-body integrator, the workhorse of the {\tt starlab} production environment for simulating dense stellar systems. The parallelization strategy, called ``j-parallelization'', involves the partition of the computational domain by distributing all particles in the system among the available processors. Partial forces on the particles to be advanced are calculated in parallel by their parent processors, and are then summed in a final global operation. Once total forces are obtained, the computing elements proceed to the computation of their particle trajectories. We report the results of timing measurements on four different parallel computers, and compare them with theoretical predictions. The computers employ either a high-speed interconnect, a NUMA architecture to minimize the communication overhead or are distributed in a grid. The code scales well in the domain tested, which ranges from 1024 - 65536 stars on 1 - 128 processors, providing satisfactory speedup. Running the production environment on a grid becomes inefficient for more than 60 processors distributed across three sites.

We investigate encounters between giant molecular clouds (GMCs) and star clusters. We propose a single expression for the energy gain of a cluster due to an encounter with a GMC, valid for all encounter distances and GMC properties. This relation is verified with N-body simulations of cluster-GMC encounters and excellent agreement is found. The fractional mass loss from the cluster is 0.25 times the fractional energy gain. This is because 75% of the injected energy goes to the velocities of escaping stars, that are higher than the escape velocity. We derive an expression for the cluster disruption time (t_dis) based on the mass loss from the simulations, taking into account the effect of gravitational focusing by the GMC. The disruption time depends on the cluster mass (M_c) and half-mass radius (r_h) as t_dis=2.0 S (M_c/10^4 M_sun)(3.75 pc/r_h)^3 Gyr, with S=1 for the solar neighbourhood and inversely proportional to the GMC density. The observed shallow relation between cluster radius and mass gives t_dis a power-law dependence on the mass with index 0.7, similar to that found from observations and from simulations of clusters dissolving in tidal fields (0.62). The constant of 2.0 Gyr is about a factor of 3.5 shorter than found from earlier simulations of clusters dissolving under the combined effect of galactic tidal field and stellar evolution. It is somewhat higher than the observationally determined value of 1.3 Gyr. It suggests, however, that the combined effect of tidal field and encounters with GMCs can explain the lack of old open clusters in the solar neighbourhood. GMC encounters can also explain the (very) short disruption time that was observed for star clusters in the central region of M51, since there rho_n is an order of magnitude higher than in the solar neighbourhood.

Hills (1988) predicted that runaway stars could be accelerated to velocities larger than 1000 km/s by dynamical encounters with the supermassive black hole (SMBH) in the Galactic center. The recently discovered hypervelocity star SDSS J090745.0+024507 (hereafter HVS) is escaping the Galaxy at high speed and could be the first object in this class. With the measured radial velocity and the estimated distance to the HVS, we trace back its trajectory in the Galactic potential. Assuming it was ejected from the center, we find that a $\sim$ 2 mas/yr proper motion is necessary for the star to have come within a few parsecs of the SMBH. We perform three-body scattering experiments to constrain the progenitor encounter which accelerated the HVS. As proposed by Yu & Tremaine (2003), we consider the tidal disruption of binary systems by the SMBH and the encounter between a star and a binary black hole, as well as an alternative scenario involving intermediate mass black holes. We find that the tidal disruption of a stellar binary ejects stars with a larger velocity compared to the encounter between a single star and a binary black hole, but has a somewhat smaller ejection rate due to the greater availability of single stars.

We investigate the possibility that mass transfer early in the evolution of a massive binary can effect a reversal of the end states of the two components, resulting in a neutron star which forms before a black hole. In this sense, such systems would comprise the high-mass analogues of white dwarf-neutron star systems such as PSR B2303+46. One consequence of this reversal is that a second episode of mass transfer from the black hole progenitor star can recycle the nascent neutron star, extending the life of the pulsar. An estimate of the formation rate through this channel is first performed via a simple analytical approach, and then refined using the results of the SeBa binary evolution package. The central role of kicks in determining the survival rate of these binaries is clearly demonstrated. The final result is expressed in terms of the number of field pulsars one can expect for every single neutron star-black hole (ns,bh) binary. We also calculate this figure for black hole-neutron star (bh,ns) systems formed through the usual channel. Assuming kicks drawn from the distribution of Cordes & Chernoff (1998), we find an expectation value of one (ns,bh) binary per 4 x 10^4 pulsars, and one (bh,ns) system for every 1500 pulsars. This helps to explain why neither system has been seen to date, though it suggests that detection of a (bh,ns) binary is imminent.

We present an X-ray binary population synthesis model, and use it to simulate the evolution of X-ray binaries formed in a burst of star formation of duration 20 Myr and star-formation rate 10 M_sun/yr. Our goal is to explain the hard (2-10 keV) X-ray properties of populations of extragalactic X-ray binaries recently observed by the Chandra X-Ray Observatory, especially those associated with recent or ongoing episodes of vigorous star formation. Our simulated X-ray binary population reaches a maximum 2-10 keV luminosity of ~4e40 erg/s after approximately 20 Myr. The X-ray luminous phase is sustained for a period of several hundreds of Myr by succeeding populations of systems with lighter secondary stars, i.e., it persists long after the star-formation episode has ended. These results are insensitive to the poorly-constrained values of the initial mass function and the average mass ratio between accreting and donor stars. The computed peak X-ray luminosity is consistent with observationally-derived correlations between the star-formation rate and the observed hard X-ray luminosity. Model cumulative luminosity functions at the earliest times have power-law indices in agreement with those derived from observations of actively star-forming galaxies. The model cumulative luminosity functions become increasingly steeper with time as the most luminous systems die off, which offers an explanation for the difference in the slopes of observed cumulative luminosity functions of young and old stellar populations.

Double neutron star binaries, analogous to the well known Hulse--Taylor pulsar PSR 1913+16, are guaranteed-to-exist sources of high frequency gravitational radiation detectable by LIGO. There is considerable uncertainty in the estimated rate of coalescence of such systems, with conservative estimates of ~1 per million years per galaxy, and optimistic theoretical estimates one or more magnitude larger. Formation rates of low-mass black hole-neutron star binaries may be higher than those of NS-NS binaries, and may dominate the detectable LIGO signal rate. We estimate the enhanced coalescence rate for BH-BH binaries due to weak asymmetric kicks during the formation of low mass black holes like Nova Sco, and find they may contribute significantly to the LIGO signal rate, possibly dominating the phase I detectable signals if the range of BH masses for which there is significant kick is broad enough. For a standard Salpeter IMF, assuming mild natal kicks, we project that the R6 merger rate of BH-BH systems is ~0.5, smaller than that of NS-NS systems. However, the higher chirp mass of these systems produces a signal nearly four times greater, on average, with a commensurate increase in search volume. The BH-BH coalescence channel considered here also predicts that a substantial fraction of BH-BH systems should have at least one component with near-maximal spin (a/M ~ 1).The waveforms produced by the coalescence of such a system should produce a clear spin signature, so this hypothesis could be directly tested by LIGO.

We develop a framework for constructing initial data sets for perturbations about spherically symmetric matter distributions. This framework facilitates setting initial data representing astrophysical sources of gravitational radiation involving relativistic stars. The procedure is based on Lichnerowicz-York's conformal approach to solve the constraints in Einstein's equations. The correspondence of these initial data sets in terms of the standard gauge perturbation variables in the Regge-Wheeler perturbation variables is established, and examples of initial data sets of merging neutron stars under the close-limit approximation are presented.

We report results of an HST "snapshot" survey aimed at finding resolved binary companions of the central stars of Galactic planetary nebulae (PNe). Using WF/PC and WFPC2, we searched the fields of 113 PNe for stars whose close proximity to the central star suggests a physical association. We find 10 binary nuclei that are very likely to be physically associated, and another six that are possible binary associations. By correcting for interstellar extinction and placing the central stars' companions on the main sequence, we derive distances to the objects, and thereby significantly increase the number of PNe with reliable distances. Comparison of our derived distances with those obtained from various statistical methods shows that all of the latter have systematically overestimated the distances, by factors ranging up to a factor of two or more. We show that this error is most likely due to the fact that the properties of our PNe with binary nuclei are systematically different from those of PNe used heretofore to calibrate statistical methods. Specifically, our PNe tend to have lower surface brightnesses at the same physical radius than the traditional calibration objects. This difference may arise from a selection effect: the PNe in our survey are typically nearby, old nebulae, whereas most of the objects that calibrate statistical techniques are low-latitude, high-surface-brightness, and more distant nebulae. As a result, the statistical methods that seem to work well with samples of distant PNe, e.g., those in the Galactic bulge or external galaxies, may not be applicable to the more diverse population of local PNe.