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-

Gravitational collapse of a Hagedorn fluid in Vaidya geometry

The gravitational collapse of a high-density null charged matter fluid, satisfying the Hagedorn equation of state, is considered in the framework of the Vaidya geometry. The general solution of the gravitational field equations can be obtained in an exact parametric form. The conditions for the formation of a naked singularity, as a result of the collapse of the compact object, are also investigated. For an appropriate choice of the arbitrary integration functions the null radial outgoing geodesic, originating from the shell focussing central singularity, admits one or more positive roots. Hence a collapsing Hagedorn fluid could end either as a black hole, or as a naked singularity. A possible astrophysical application of the model, to describe the energy source of gamma-ray bursts, is also considered.Comment: 14 pages, 2 figures, 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

Early MIMD experience with a plasma physics simulation program on the CRAY X-MP

This paper describes some early experience with converting a plasma physics simulation program to the CRAY X-MP, a current multiple instruction, multiple data (MIMD) computer consisting of two processors with architecture similar to that of the CRAY-1. The computer program used in this study is an all Fortran version of SELF, a two species, one space, two velocity, electromagnetic, Newtonian, particle in cell, plasma simulation code. The approach to converting SELF to use both processors of the CRAY X-MP is described in some detail. The resulting multiprocessor version of SELF is nearly a factor of two faster in real time than the single processor version. The multiprocessor version obtains 58.2+-.1 seconds of central processor time in 30+-.5 seconds of real time. For comparison, the CRAY-1 execution time if 74.5 seconds. For SELF, which is mostly scalar coding, the CRAY X-MP is about 2.5 times faster overall than the CRAY-1

Understanding large software systems with the utility XREG

When a software system reaches a certain size, arteriosclerosis sets in. That is, the system becomes harder and harder to add new features and even more difficult to understand. Worse, the insertion of new features often introduces new programming errors as well as revealing ones already present but previously unmanifested. For Fortran programs, the 100,000 line size is critical. XREF fights hardening of the arteries by providing the user information on the organization, variable usage, and common block use within the software system. XREF is language independent. That is, it works, for example, with CAL, C, CIVIC, and CFT compiled subprograms. XREF performs its global symbol analysis from either the BUILD library file or the object (binary) file. Naturally more useful information is obtained when a symbol table is generated by the compiler. Specifically, eight global reports are produced. The most helpful report consists of a listing of all symbols contained in the BUILD library or binary file, the subprograms that use the symbol and the relocation basis (either local or common block name) of each symbol. Five secondary reports are provided for each subprogram. 3 refs