Description of the Scenario Machine

We present here an updated description of the "Scenario Machine" code. This tool is used to carry out a population synthesis of binary stars. Previous version of the description can be found at http://xray.sai.msu.ru/~mystery//articles/review/contents.htmlComment: 32 pages, 3 figures. Corrected typo

Magnetic fileds of coalescing neutron stars and the luminosity function of short gamma-ray burst

Coalescing neutron star binaries are believed to be the most reliable sources for ground-based detectors of gravitational waves and likely progenitors of short gamma-ray bursts. In the process of coalescence, magnetic fields of neutron stars can induce interesting observational manifestations and affect the form of gravitational wave signal. In this papaer we use the population synthesis method to model the expected distribution of neutron star magnetic fields during the coalescence under different assumptions on the initial parameters of neutron stars and their magnetic field evolution. We discuss possible elecotrmagnetic phenomena preceding the coalescence of magnetized neutron star binaries and the effect of magnetic field on the gravitational wave signal. We find that a log-normal (Gaussian in logarithms) distribution of the initial magnetic fields of neutron stars, which agrees with observed properties of radio pulsars, produces the distribution of the magnetic field energy during the coalescence that adequately describes the observed luminosity function of short gamma-ray bursts under different assumptions on the field evolution and initial parameters of neutron stars. This agreement lends further support to the model of coalescing neutron star binaries as progenitors of gamma-ray bursts.Comment: v.2, LATEX, 25 pages, inc. 7 ps figures, Astron. Lett., in press. Typos corrected, reference adde

Radio precursors to neutron star binary mergings

We discuss a possible generation of radio bursts preceding final stages of binary neutron star mergings which can be accompanied by short gamma-ray bursts. Detection of such bursts appear to be advantageous in the low-frequency radio band due to a time delay of ten to several hundred seconds required for radio signal to propagate in the ionized intergalactic medium. This delay makes it possible to use short gamma-ray burst alerts to promptly monitor specific regions on the sky by low-frequency radio facilities, especially by LOFAR. To estimate the strength of the radio signal, we assume a power-law dependence of the radio luminosity on the total energy release in a magnetically dominated outflow, as found in millisecond pulsars. Based on the planned LOFAR sensitivity at 120 MHz, we estimate that the LOFAR detection rate of such radio transients could be about several events per month from redshifts up to $z\sim1.3$ in the most optimistic scenario. The LOFAR ability to detect such events would crucially depend on exact efficiency of low-frequency radio emission mechanism.Comment: 6 pages, 2 figures, Accepted for publication in Astrophysics & Space Science. Largely extended version of ArXiv:0912.521

A Comparison of search templates for gravitational waves from binary inspiral

We compare the performances of the templates defined by three different types of approaches: traditional post-Newtonian templates (Taylor-approximants), ``resummed'' post-Newtonian templates assuming the adiabatic approximation and stopping before the plunge (P-approximants), and further ``resummed'' post-Newtonian templates going beyond the adiabatic approximation and incorporating the plunge with its transition from the inspiral (Effective-one-body approximants). The signal to noise ratio is significantly enhanced (mainly because of the inclusion of the plunge signal) by using these new effective-one-body templates relative to the usual post-Newtonian ones for binary masses greater than $30 M_\odot$, the most likely sources for initial laser interferometers. Independently of the question of the plunge signal, the comparison of the various templates confirms the usefulness of using resummation methods. The paper also summarizes the key elements of the construction of various templates and thus can serve as a resource for those involved in writing inspiral search software.Comment: eta-dependent tail terms corrected after related errata by Blanchet (2005

Magnetars, Gamma-ray Bursts, and Very Close Binaries

We consider the possible existence of a common channel of evolution of binary systems, which results in a gamma-ray burst during the formation of a black hole or the birth of a magnetar during the formation of a neutron star. We assume that the rapid rotation of the core of a collapsing star can be explained by tidal synchronization in a very close binary. The calculated rate of formation of rapidly rotating neutron stars is qualitatively consistent with estimates of the formation rate of magnetars. However, our analysis of the binarity of newly-born compact objects with short rotational periods indicates that the fraction of binaries among them substantially exceeds the observational estimates. To bring this fraction into agreement with the statistics for magnetars, the additional velocity acquired by a magnetar during its formation must be primarily perpendicular to the orbital plane before the supernova explosion, and be large.Comment: 16 pages, 4 figure

Circular orbits of corotating binary black holes: comparison between analytical and numerical results

We compare recent numerical results, obtained within a ``helical Killing vector'' (HKV) approach, on circular orbits of corotating binary black holes to the analytical predictions made by the effective one body (EOB) method (which has been recently extended to the case of spinning bodies). On the scale of the differences between the results obtained by different numerical methods, we find good agreement between numerical data and analytical predictions for several invariant functions describing the dynamical properties of circular orbits. This agreement is robust against the post-Newtonian accuracy used for the analytical estimates, as well as under choices of resummation method for the EOB ``effective potential'', and gets better as one uses a higher post-Newtonian accuracy. These findings open the way to a significant ``merging'' of analytical and numerical methods, i.e. to matching an EOB-based analytical description of the (early and late) inspiral, up to the beginning of the plunge, to a numerical description of the plunge and merger. We illustrate also the ``flexibility'' of the EOB approach, i.e. the possibility of determining some ``best fit'' values for the analytical parameters by comparison with numerical data.Comment: Minor revisions, accepted for publication in Phys. Rev. D, 19 pages, 6 figure

Post-Newtonian SPH calculations of binary neutron star coalescence. II. Binary mass ratio, equation of state, and spin dependence

Using our new Post-Newtonian SPH (smoothed particle hydrodynamics) code, we study the final coalescence and merging of neutron star (NS) binaries. We vary the stiffness of the equation of state (EOS) as well as the initial binary mass ratio and stellar spins. Results are compared to those of Newtonian calculations, with and without the inclusion of the gravitational radiation reaction. We find a much steeper decrease in the gravity wave peak strain and luminosity with decreasing mass ratio than would be predicted by simple point-mass formulae. For NS with softer EOS (which we model as simple $\Gamma=2$ polytropes) we find a stronger gravity wave emission, with a different morphology than for stiffer EOS (modeled as $\Gamma=3$ polytropes as in our previous work). We also calculate the coalescence of NS binaries with an irrotational initial condition, and find that the gravity wave signal is relatively suppressed compared to the synchronized case, but shows a very significant second peak of emission. Mass shedding is also greatly reduced, and occurs via a different mechanism than in the synchronized case. We discuss the implications of our results for gravity wave astronomy with laser interferometers such as LIGO, and for theoretical models of gamma-ray bursts (GRBs) based on NS mergers.Comment: RevTeX, 38 pages, 24 figures, Minor Corrections, to appear in Phys. Rev.

Post-Newtonian SPH calculations of binary neutron star coalescence. I. Method and first results

We present the first results from our Post-Newtonian (PN) Smoothed Particle Hydrodynamics (SPH) code, which has been used to study the coalescence of binary neutron star (NS) systems. The Lagrangian particle-based code incorporates consistently all lowest-order (1PN) relativistic effects, as well as gravitational radiation reaction, the lowest-order dissipative term in general relativity. We test our code on sequences of single NS models of varying compactness, and we discuss ways to make PN simulations more relevant to realistic NS models. We also present a PN SPH relaxation procedure for constructing equilibrium models of synchronized binaries, and we use these equilibrium models as initial conditions for our dynamical calculations of binary coalescence. Though unphysical, since tidal synchronization is not expected in NS binaries, these initial conditions allow us to compare our PN work with previous Newtonian results. We compare calculations with and without 1PN effects, for NS with stiff equations of state, modeled as polytropes with $\Gamma=3$. We find that 1PN effects can play a major role in the coalescence, accelerating the final inspiral and causing a significant misalignment in the binary just prior to final merging. In addition, the character of the gravitational wave signal is altered dramatically, showing strong modulation of the exponentially decaying waveform near the end of the merger. We also discuss briefly the implications of our results for models of gamma-ray bursts at cosmological distances.Comment: RevTeX, 37 pages, 17 figures, to appear in Phys. Rev. D, minor corrections onl

Multimessenger astronomy with the Einstein Telescope

Gravitational waves (GWs) are expected to play a crucial role in the development of multimessenger astrophysics. The combination of GW observations with other astrophysical triggers, such as from gamma-ray and X-ray satellites, optical/radio telescopes, and neutrino detectors allows us to decipher science that would otherwise be inaccessible. In this paper, we provide a broad review from the multimessenger perspective of the science reach offered by the third generation interferometric GW detectors and by the Einstein Telescope (ET) in particular. We focus on cosmic transients, and base our estimates on the results obtained by ET's predecessors GEO, LIGO, and Virgo.Comment: 26 pages. 3 figures. Special issue of GRG on the Einstein Telescope. Minor corrections include

Accretion, Outflows, and Winds of Magnetized Stars

Many types of stars have strong magnetic fields that can dynamically influence the flow of circumstellar matter. In stars with accretion disks, the stellar magnetic field can truncate the inner disk and determine the paths that matter can take to flow onto the star. These paths are different in stars with different magnetospheres and periods of rotation. External field lines of the magnetosphere may inflate and produce favorable conditions for outflows from the disk-magnetosphere boundary. Outflows can be particularly strong in the propeller regime, wherein a star rotates more rapidly than the inner disk. Outflows may also form at the disk-magnetosphere boundary of slowly rotating stars, if the magnetosphere is compressed by the accreting matter. In isolated, strongly magnetized stars, the magnetic field can influence formation and/or propagation of stellar wind outflows. Winds from low-mass, solar-type stars may be either thermally or magnetically driven, while winds from massive, luminous O and B type stars are radiatively driven. In all of these cases, the magnetic field influences matter flow from the stars and determines many observational properties. In this chapter we review recent studies of accretion, outflows, and winds of magnetized stars with a focus on three main topics: (1) accretion onto magnetized stars; (2) outflows from the disk-magnetosphere boundary; and (3) winds from isolated massive magnetized stars. We show results obtained from global magnetohydrodynamic simulations and, in a number of cases compare global simulations with observations.Comment: 60 pages, 44 figure