Observational constraints on the neutron star mass distribution

Radio observations of neutron star binary pulsar systems have constrained strongly the masses of eight neutron stars. Assuming neutron star masses are uniformly distributed between lower and upper bounds $m_l$ and $m_u$, the observations determine with 95\% confidence that $1.01 < m_l/\text{M}_\odot < 1.34$ and $1.43 < m_u/\text{M}_\odot < 1.64$. These limits give observational support to neutron star formation scenarios that suggest that masses should fall predominantly in the range $1.3<m/\text{M}_\odot<1.6$, and will also be important in the interpretation of binary inspiral observations by the Laser Interferometer Gravitational-wave Observatory.Comment: Postscript, 4 pages, NU-GR-

Modeling the strangeness content of hadronic matter

The strangeness content of hadronic matter is studied in a string-flip model that reproduces various aspects of the QCD-inspired phenomenology, such as quark clustering at low density and color deconfinement at high density, while avoiding long range van der Waals forces. Hadronic matter is modeled in terms of its quark constituents by taking into account its internal flavor (u,d,s) and color (red, blue, green) degrees of freedom. Variational Monte-Carlo simulations in three spatial dimensions are performed for the ground-state energy of the system. The onset of the transition to strange matter is found to be influenced by weak, yet not negligible, clustering correlations. The phase diagram of the system displays an interesting structure containing both continuous and discontinuous phase transitions. Strange matter is found to be absolutely stable in the model.Comment: 14 pages, 1 table, 8 eps figures, revtex. Submitted to Phys. Rev. C, Presented at INPC2001 Berkeley, Ca. july 29-Aug

Measuring proper motions of isolated neutron stars with Chandra

The excellent spatial resolution of the Chandra observatory offers the unprecedented possibility to measure proper motions at X-ray wavelength with relatively high accuracy using as reference the background of extragalactic or remote galactic X-ray sources. We took advantage of this capability to constrain the proper motion of RX J0806.4-4123 and RX J0420.0-5022, two X-ray bright and radio quiet isolated neutron stars (INSs) discovered by ROSAT and lacking an optical counterpart. In this paper, we present results from a preliminary analysis from which we derive 2 sigma upper limits of 76 mas/yr and 138 mas/yr on the proper motions of RX J0806.4-4123 and RX J0420.0-5022 respectively. We use these values together with those of other ROSAT discovered INSs to constrain the origin, distance and evolutionary status of this particular group of objects. We find that the tangential velocities of radio quiet ROSAT neutron stars are probably consistent with those of 'normal' pulsars. Their distribution on the sky and, for those having accurate proper motion vectors, their possible birth places, all point to a local population, probably created in the part of the Gould Belt nearest to the earth.Comment: 8 pages, 3 figures, to appear in Astrophysics and Space Science, in the proceedings of "Isolated Neutron Stars: from the Interior to the Surface", edited by D. Page, R. Turolla and S. Zan

HST/STIS Observations of the Optical Counterpart to GRB 970228

We report on observations of the fading optical counterpart of the gamma-ray burst GRB 970228, made on 4~September~1997 using the STIS CCD on the Hubble Space Telescope. The unresolved counterpart is detected at V=28 +/- 0.25, consistent with a continued power-law decline with exponent -1.14 +/- 0.05. No proper motion is detected, in contradiction of some earlier claims. The counterpart is located within, but near the edge of, a faint extended source with diameter ~0."8 and integrated magnitude 25.7 +/- 0.25. Comparison with WFPC2 data taken one month after the initial burst and NTT data taken on March 13 shows no evidence for fading of the extended emission. After adjusting for the probable Galactic extinction in the direction of GRB 970228 of A_v=0.7, we find that the observed nebula is consistent with the sizes of galaxies of comparable magnitude found in the Hubble Deep Field and other deep HST images, and that only 2% of the sky is covered by galaxies of similar or greater surface brightness. Therefore, the extended source observed about GRB 970228 is most likely a galaxy at moderate redshift, and is almost certainly the host of the gamma-ray burst

Self-consistent description of nuclear compressional modes

Isoscalar monopole and dipole compressional modes are computed for a variety of closed-shell nuclei in a relativistic random-phase approximation to three different parametrizations of the Walecka model with scalar self-interactions. Particular emphasis is placed on the role of self-consistency which by itself, and with little else, guarantees the decoupling of the spurious isoscalar-dipole strength from the physical response and the conservation of the vector current. A powerful new relation is introduced to quantify the violation of the vector current in terms of various ground-state form-factors. For the isoscalar-dipole mode two distinct regions are clearly identified: (i) a high-energy component that is sensitive to the size of the nucleus and scales with the compressibility of the model and (ii) a low-energy component that is insensitivity to the nuclear compressibility. A fairly good description of both compressional modes is obtained by using a ``soft'' parametrization having a compression modulus of K=224 MeV.Comment: 28 pages and 10 figures; submitted to PR

Various features of quasiequilibrium sequences of binary neutron stars in general relativity

Quasiequilibrium sequences of binary neutron stars are numerically calculated in the framework of the Isenberg-Wilson-Mathews (IWM) approximation of general relativity. The results are presented for both rotation states of synchronized spins and irrotational motion, the latter being considered as the realistic one for binary neutron stars just prior to the merger. We assume a polytropic equation of state and compute several evolutionary sequences of binary systems composed of different-mass stars as well as identical-mass stars with adiabatic indices gamma=2.5, 2.25, 2, and 1.8. From our results, we propose as a conjecture that if the turning point of binding energy (and total angular momentum) locating the innermost stable circular orbit (ISCO) is found in Newtonian gravity for some value of the adiabatic index gamma_0, that of the ADM mass (and total angular momentum) should exist in the IWM approximation of general relativity for the same value of the adiabatic index.Comment: Text improved, some figures changed or deleted, new table, 38 pages, 31 figures, accepted for publication in Phys. Rev.

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.

Maximally incompressible neutron star matter

Relativistic kinetic theory, based on the Grad method of moments as developed by Israel and Stewart, is used to model viscous and thermal dissipation in neutron star matter and determine an upper limit on the maximum mass of neutron stars. In the context of kinetic theory, the equation of state must satisfy a set of constraints in order for the equilibrium states of the fluid to be thermodynamically stable and for perturbations from equilibrium to propagate causally via hyperbolic equations. Application of these constraints to neutron star matter restricts the stiffness of the most incompressible equation of state compatible with causality to be softer than the maximally incompressible equation of state that results from requiring the adiabatic sound speed to not exceed the speed of light. Using three equations of state based on experimental nucleon-nucleon scattering data and properties of light nuclei up to twice normal nuclear energy density, and the kinetic theory maximally incompressible equation of state at higher density, an upper limit on the maximum mass of neutron stars averaging 2.64 solar masses is derived.Comment: 8 pages, 2 figure

Topological Defects and CMB anisotropies : Are the predictions reliable ?

We consider a network of topological defects which can partly decay into neutrinos, photons, baryons, or Cold Dark Matter. We find that the degree-scale amplitude of the cosmic microwave background (CMB) anisotropies as well as the shape of the matter power spectrum can be considerably modified when such a decay is taken into account. We conclude that present predictions concerning structure formation by defects might be unreliable.Comment: 14 pages, accepted for publication in PR

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