BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

John L. Friedman

John L. Friedman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-phastro-ph.SRhep-phmath-phmath.MPnucl-th

Catalog footprint

What is connected

19works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Astrophysical Implications of Neutron Star Inspiral and Coalescence

The first inspiral of two neutron stars observed in gravitational waves was remarkably close, allowing the kind of simultaneous gravitational wave and electromagnetic observation that had not been expected for several years. Their merger, followed by a gamma-ray burst and a kilonova, was observed across the spectral bands of electromagnetic telescopes. These GW and electromagnetic observations have led to dramatic advances in understanding short gamma-ray bursts; determining the origin of the heaviest elements; and determining the maximum mass of neutron stars. From the imprint of tides on the gravitational waveforms and from observations of X-ray binaries, one can extract the radius and deformability of inspiraling neutron stars. Together, the radius, maximum mass, and causality constrain the neutron-star equation of state, and future constraints can come from observations of post-merger oscillations. We selectively review these results, filling in some of the physics with derivations and estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Combined gravitational and electromagnetic self-force on charged particles in electrovac spacetimes

We consider the self-force on a charged particle moving in a curved spacetime with a background electromagnetic field, extending previous studies to situations in which gravitational and electromagnetic perturbations are comparable. The formal expression $f^{ret}_α$ for the self-force on a particle, written in terms of the retarded perturbed fields, is divergent, and a renormalization is needed to find the particle's acceleration at linear order in its mass $m$ and charge $e$. We assume that, as in previous work in a Lorenz gauge, the renormalization for accelerated motion comprises an angle average and mass renormalization. Using the short distance expansion of the perturbed electromagnetic and gravitational fields, we show that the renormalization is equivalent to that obtained from a mode sum regularization in which one subtracts from the expression for the self-force in terms of the retarded fields a singular part field comprising only the leading and subleading terms in the mode sum. The most striking part of our result, arising from a remarkable cancellation, is that the renormalization involves no mixing of electromagnetic and gravitational fields. In particular, the renormalized mass is obtained by subtracting (1) the purely electromagnetic contribution from a point charge moving along an accelerated trajectory and (2) the purely gravitational contribution from a point mass moving along the same trajectory. In a mode-sum regularization, the same cancellation implies that the required regularization parameters are sums of their purely electromagnetic and gravitational values.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational Self-force in a Radiation Gauge

In this, the first of two companion papers, we present a method for finding the gravitational self-force in a modified radiation gauge for a particle moving on a geodesic in a Schwarzschild or Kerr spacetime. An extension of an earlier result by Wald is used to show the spin-weight $\pm 2$ perturbed Weyl scalar ($ψ_0$ or $ψ_4$) determines the metric perturbation outside the particle up to a gauge transformation and an infinitesimal change in mass and angular momentum. A Hertz potential is used to construct the part of the retarded metric perturbation that involves no change in mass or angular momentum from $ψ_0$ in a radiation gauge. The metric perturbation is completed by adding changes in the mass and angular momentum of the background spacetime outside the radial coordinate $r_0$ of the particle in any convenient gauge. The resulting metric perturbation is singular on the trajectory of the particle and discontinuous across the sphere $r=r_0$. A mode-sum method can be used to renormalize the self-force, but the justification given in the published version of this paper \cite{sf2} referred to work by Sam Gralla \cite{gralla10} to justify the use of the renormalized self-force, and the radiation gauge we use does not satisfy the regularity conditions required by Gralla. Instead we show that the renormalized self-force, computed either from the retarded field for $r>r_0$ or for $r<r_0$ gives the correct equations of motion in a gauge smoothly related to a Lorenz gauge; and Pound et al. \cite{pmb13} argue that the average of the self-force obtained in the way described in our paper for $r>r_0$ and for $r<r_0$ gives the correct equation of motion for our gauge (what Pound et al. call the no-string gauge).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self force on an accelerated particle

We calculate the singular field of an accelerated point particle (scalar charge, electric charge or small gravitating mass) moving on an accelerated (non-geodesic) trajectory in a generic background spacetime. Using a mode-sum regularization scheme, we obtain explicit expressions for the self-force regularization parameters. In the electromagnetic and gravitational case, we use a Lorenz gauge. This work extends the work of Barack and Ori [1] who demonstrated that the regularization parameters for a point particle in geodesic motion in a Schwarzschild spacetime can be described solely by the leading and subleading terms in the mode-sum (commonly known as the $A$ and $B$ terms) and that all terms of higher order in $\ell$ vanish upon summation (later they showed the same behavior for geodesic motion in Kerr [2], [3]). We demonstrate that these properties are universal to point particles moving through any smooth spacetime along arbitrary (accelerated) trajectories. Our renormalization scheme is based on, but not identical to, the Quinn-Wald axioms. As we develop our approach, we review and extend work showing that that different definitions of the singular field used in the literature are equivalent to our approach. Because our approach does not assume geodesic motion of the perturbing particle, we are able use our mode-sum formalism to explicitly recover a well-known result: The self-force on static scalar charges near a Schwarzschild black hole vanishes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extracting equation of state parameters from black hole-neutron star mergers: aligned-spin black holes and a preliminary waveform model

Information about the neutron-star equation of state is encoded in the waveform of a black hole-neutron star system through tidal interactions and the possible tidal disruption of the neutron star. During the inspiral this information depends on the tidal deformability Lambda of the neutron star, and we find that Lambda is the best measured parameter during the merger and ringdown as well. We performed 134 simulations where we systematically varied the equation of state as well as the mass ratio, neutron star mass, and aligned spin of the black hole. Using these simulations we have developed an analytic representation of the full inspiral-merger-ringdown waveform calibrated to these numerical waveforms, and we use this analytic waveform to estimate the accuracy to which Lambda can be measured with gravitational-wave detectors. We find that although the inspiral tidal signal is small, coherently combining this signal with the merger-ringdown matter effect improves the measurability of Lambda by a factor of ~3 over using just the merger-ringdown matter effect alone. However, incorporating correlations between all the waveform parameters then decreases the measurability of Lambda by a factor of ~3. The uncertainty in Lambda increases with the mass ratio, but decreases as the black hole spin increases. Overall, a single Advanced LIGO detector can measure Lambda for mass ratios Q = 2--5, black hole spins J_BH/M_BH^2 = -0.5--0.75, neutron star masses M_NS = 1.2M_sun--1.45M_sun, and an optimally oriented distance of 100Mpc to an uncertainty of ~10%--100%. For the proposed Einstein Telescope, the uncertainty in Lambda is an order of magnitude smaller.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Matter effects on binary neutron star waveforms

Using an extended set of equations of state and a multiple-group multiple-code collaborative effort to generate waveforms, we improve numerical-relativity-based data-analysis estimates of the measurability of matter effects in neutron-star binaries. We vary two parameters of a parameterized piecewise-polytropic equation of state (EOS) to analyze the measurability of EOS properties, via a parameter Λ that characterizes the quadrupole deformability of an isolated neutron star. We find that, to within the accuracy of the simulations, the departure of the waveform from point-particle (or spinless double black-hole binary) inspiral increases monotonically with Λ, and changes in the EOS that did not change Λ are not measurable. We estimate with two methods the minimal and expected measurability of Λ in second- and third- generation gravitational-wave detectors. The first estimate, using numerical waveforms alone, shows two EOS which vary in radius by 1.3km are distinguishable in mergers at 100Mpc. The second estimate relies on the construction of hybrid waveforms by matching to post-Newtonian inspiral, and estimates that the same EOS are distinguishable in mergers at 300Mpc. We calculate systematic errors arising from numerical uncertainties and hybrid construction, and we estimate the frequency at which such effects would interfere with template-based searches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

EMRI corrections to the angular velocity and redshift factor of a mass in circular orbit about a Kerr black hole

This is the first of two papers on computing the self-force in a radiation gauge for a particle moving in circular, equatorial orbit about a Kerr black hole. In the EMRI (extreme-mass-ratio inspiral) framework, with mode-sum renormalization, we compute the renormalized value of the quantity $h_{αβ}u^αu^β$, gauge-invariant under gauge transformations generated by a helically symmetric gauge vector; and we find the related order $\frak{m}$ correction to the particle's angular velocity at fixed renormalized redshift (and to its redshift at fixed angular velocity). The radiative part of the perturbed metric is constructed from the Hertz potential which is extracted from the Weyl scalar by an algebraic inversion\cite{sf2}. We then write the spin-weighted spheroidal harmonics as a sum over spin-weighted spherical harmonics and use mode-sum renormalization to find the renormalization coefficients by matching a series in $L=\ell+1/2$ to the large-$L$ behavior of the expression for $H := \frac12 h_{αβ}u^αu^β$. The non-radiative parts of the perturbed metric associated with changes in mass and angular momentum are calculated in the Kerr gauge.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Extracting equation of state parameters from black hole-neutron star mergers. I. Nonspinning black holes

The late inspiral, merger, and ringdown of a black hole-neutron star (BHNS) system can provide information about the neutron-star equation of state (EOS). Candidate EOSs can be approximated by a parametrized piecewise-polytropic EOS above nuclear density, matched to a fixed low-density EOS; and we report results from a large set of BHNS inspiral simulations that systematically vary two parameters. To within the accuracy of the simulations, we find that, apart from the neutron-star mass, a single physical parameter Lambda, describing its deformability, can be extracted from the late inspiral, merger, and ringdown waveform. This parameter is related to the radius, mass, and l=2 Love number, k_2, of the neutron star by Lambda = 2k_2 R^5/3M_{NS}^5, and it is the same parameter that determines the departure from point-particle dynamics during the early inspiral. Observations of gravitational waves from BHNS inspiral thus restrict the EOS to a surface of constant Lambda in the parameter space, thickened by the measurement error. Using various configurations of a single Advanced LIGO detector, we find that neutron stars are distinguishable from black holes of the same mass and that Lambda^{1/5} or equivalently R can be extracted to 10-40% accuracy from single events for mass ratios of Q=2 and 3 at a distance of 100 Mpc, while with the proposed Einstein Telescope, EOS parameters can be extracted to accuracy an order of magnitude better.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Neutron star equation of state via gravitational wave observations

Gravitational wave observations can potentially measure properties of neutron star equations of state by measuring departures from the point-particle limit of the gravitational waveform produced in the late inspiral of a neutron star binary. Numerical simulations of inspiraling neutron star binaries computed for equations of state with varying stiffness are compared. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly if the neutron stars are more compact. This suggests that gravitational wave observations at frequencies around 1 kHz will be able to measure a compactness parameter and place stringent bounds on possible neutron star equations of state. Advanced laser interferometric gravitational wave observatories will be able to tune their frequency band to optimize sensitivity in the required frequency range to make sensitive measures of the late-inspiral phase of the coalescence.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Iteration Stability for Simple Newtonian Stellar Systems

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", $κ$, 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 $κ$ 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-$κ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Measuring the neutron star equation of state with gravitational wave observations

We report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matter described by a parametrized equation of state. Previous work suggested the use of an effective cutoff frequency to place constraints on the equation of state. We find, however, that greater accuracy is obtained by measuring departures from the point-particle limit of the gravitational waveform produced during the late inspiral. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly for smaller values of the neutron star compactness (the ratio of the mass of the neutron star to its radius). We estimate that realistic equations of state will lead to gravitational waveforms that are distinguishable from point particle inspirals at an effective distance (the distance to an optimally oriented and located system that would produce an equivalent waveform amplitude) of 100 Mpc or less. As Lattimer and Prakash observed, neutron-star radius is closely tied to the pressure at density not far above nuclear. Our results suggest that broadband gravitational wave observations at frequencies between 500 and 1000 Hz will constrain this pressure, and we estimate the accuracy with which it can be measured. Related first estimates of radius measurability show that the radius can be determined to an accuracy of ~1 km at 100 Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-conformally flat initial data for binary compact objects

A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all metric components, including the spatially non-conformally flat potentials, and for irrotational flow. A numerical method applicable to both formulations is explained with an emphasis on the imposition of a spatial gauge condition. Results are shown for solution sequences of irrotational BNS with matter approximated by parametrized equations of state that use a few segments of polytropic equations of state. The binding energy and total angular momentum of solution sequences computed within the conformally flat -- Isenberg-Wilson-Mathews (IWM) -- formulation are closer to those of the third post-Newtonian (3PN) two point particles up to the closest orbits, for the more compact stars, whereas sequences resulting from the WL/NHS formulations deviate from the 3PN curve even more for the sequences with larger compactness. We think it likely that this correction reflects an overestimation in the IWM formulation as well as in the 3PN formula, by $\sim 1$ cycle in the gravitational wave phase during the last several orbits. The work suggests that imposing spatial conformal flatness results in an underestimate of the quadrupole deformation of the components of binary neutron-star systems in the last few orbits prior to merger.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Deriving formulations for numerical computation of binary neutron stars in quasicircular orbits

Two relations, the virial relation $M_{\rm ADM}=M_{\rm K}$ and the first law in the form $δM_{\rm ADM}=ΩδJ$, should be satisfied by a solution and a sequence of solutions describing binary compact objects in quasiequilibrium circular orbits. Here, $M_{\rm ADM}$, $M_{\rm K}$, $J$, and $Ω$ are the ADM mass, Komar mass, angular momentum, and orbital angular velocity, respectively. $δ$ denotes an Eulerian variation. These two conditions restrict the allowed formulations that we may adopt. First, we derive relations between $M_{\rm ADM}$ and $M_{\rm K}$ and between $δM_{\rm ADM}$ and $ΩδJ$ for general asymptotically flat spacetimes. Then, to obtain solutions that satisfy the virial relation and sequences of solutions that satisfy the first law at least approximately, we propose a formulation for computation of quasiequilibrium binary neutron stars in general relativity. In contrast to previous approaches in which a part of the Einstein equation is solved, in the new formulation, the full Einstein equation is solved with maximal slicing and in a transverse gauge for the conformal three-metric. Helical symmetry is imposed in the near zone, while in the distant zone, a waveless condition is assumed. We expect the solutions obtained in this formulation to be excellent quasiequilibria as well as initial data for numerical simulations of binary neutron star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Nonlinear r-Modes in Neutron Stars: Instability of an unstable mode

We study the dynamical evolution of a large amplitude r-mode by numerical simulations. R-modes in neutron stars are unstable growing modes, driven by gravitational radiation reaction. In these simulations, r-modes of amplitude unity or above are destroyed by a catastrophic decay: A large amplitude r-mode gradually leaks energy into other fluid modes, which in turn act nonlinearly with the r-mode, leading to the onset of the rapid decay. As a result the r-mode suddenly breaks down into a differentially rotating configuration. The catastrophic decay does not appear to be related to shock waves at the star's surface. The limit it imposes on the r-mode amplitude is significantly smaller than that suggested by previous fully nonlinear numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Thermodynamics of binary black holes and neutron stars

We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid spacetimes with a helical Killing vector k^α, heuristically, the generator of time-translations in a corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black-hole binaries, a rigidity theorem implies that the Killing vector lies along the horizon's generators, and from this one can deduce the zeroth law (constant surface gravity of the horizon). Remarkably, although the mass and angular momentum of such a system are not defined, there is an exact first law, relating the change in the asymptotic Noether charge to the changes in the vorticity, baryon mass, and entropy of the fluid, and in the area of black holes. Binary systems with MΩsmall have an approximate asymptopia in which one can write the first law in terms of the asymptotic mass and angular momentum. Asymptotic flatness is precise in two classes of solutions used to model binary systems: spacetimes satisfying the post-Newtonian equations, and solutions to a modified set of field equations that have a spatially conformally flat metric. (The spatial conformal flatness formalism with helical symmetry, however, is consistent with maximal slicing only if replaces the extrinsic curvature in the field equations by an artificially tracefree expression in terms of the shift vector.) For these spacetimes, nearby equilibria whose stars have the same vorticity obey the relation δM = ΩδJ, from which one can obtain a turning point criterion that governs the stability of orbits.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Hamiltonian spacetime dynamics with a spherical null-dust shell

We consider the Hamiltonian dynamics of spherically symmetric Einstein gravity with a thin null-dust shell, under boundary conditions that fix the evolution of the spatial hypersurfaces at the two asymptotically flat infinities of a Kruskal-like manifold. The constraints are eliminated via a Kuchar-type canonical transformation and Hamiltonian reduction. The reduced phase space $\tildeΓ$ consists of two disconnected copies of $R^4$, each associated with one direction of the shell motion. The right-moving and left-moving test shell limits can be attached to the respective components of $\tildeΓ$ as smooth boundaries with topology $R^3$. Choosing the right-hand-side and left-hand-side masses as configuration variables provides a global canonical chart on each component of $\tildeΓ$, and renders the Hamiltonian simple, but encodes the shell dynamics in the momenta in a convoluted way. Choosing the shell curvature radius and the "interior" mass as configuration variables renders the shell dynamics transparent in an arbitrarily specifiable stationary gauge "exterior" to the shell, but the resulting local canonical charts do not cover the three-dimensional subset of $\tildeΓ$ that corresponds to a horizon-straddling shell. When the evolution at the infinities is freed by introducing parametrization clocks, we find on the unreduced phase space a global canonical chart that completely decouples the physical degrees of freedom from the pure gauge degrees of freedom. Replacing one infinity by a flat interior leads to analogous results, but with the reduced phase space $R^2 \cup R^2$. The utility of the results for quantization is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Reduced phase space formalism for spherically symmetric geometry with a massive dust shell

We perform a Hamiltonian reduction of spherically symmetric Einstein gravity with a thin dust shell of positive rest mass. Three spatial topologies are considered: Euclidean (R^3), Kruskal (S^2 x R), and the spatial topology of a diametrically identified Kruskal (RP^3 - {a point at infinity}). For the Kruskal and RP^3 topologies the reduced phase space is four-dimensional, with one canonical pair associated with the shell and the other with the geometry; the latter pair disappears if one prescribes the value of the Schwarzschild mass at an asymptopia or at a throat. For the Euclidean topology the reduced phase space is necessarily two-dimensional, with only the canonical pair associated with the shell surviving. A time-reparametrization on a two-dimensional phase space is introduced and used to bring the shell Hamiltonians to a simpler (and known) form associated with the proper time of the shell. An alternative reparametrization yields a square-root Hamiltonian that generalizes the Hamiltonian of a test shell in Minkowski space with respect to Minkowski time. Quantization is briefly discussed. The discrete mass spectrum that characterizes natural minisuperspace quantizations of vacuum wormholes and RP^3-geons appears to persist as the geometrical part of the mass spectrum when the additional matter degree of freedom is added.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Quantum field theory in Lorentzian universes-from-nothing

We examine quantum field theory in spacetimes that are time nonorientable but have no other causal pathology. These are Lorentzian universes-from-nothing, spacetimes with a single spacelike boundary that nevertheless have a smooth Lorentzian metric. Classically, such spacetimes are locally indistinguishable from their globally hyperbolic covering spaces. However, the construction of a quantum field theory (QFT) is more problematic. One can define a family of local algebras on an atlas of globally hyperbolic subspacetimes. But one cannot extend a generic positive linear function from a single algebra to the collection of all local algebras without violating positivity, while satisfying the physically appropriate overlap conditions. This difficulty can be overcome by restricting the size of neighborhoods so that the union of any pair is time- orientable. The structure of local algebras and states is then locally indistinguishable from that of QFT on a globally hyperbolic spacetime. But this size restriction on neighborhoods makes the structure unsatisfactory as a global field theory. The theory allows less information than QFT in a globally hyperbolic spacetime, because correlations between field operators at a pair of points are defined only if a curve joining the points lies in a single neigh- borhood. Moreover, to extend a local state to a collection of states, we use an antipodally symmetric state on the covering space, a state that would not yield a sensible state on the spacetime if all correlations could be measured.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Comparing Models of Rapidly Rotating Relativistic Stars Constructed by Two Numerical Methods

We present the first direct comparison of codes based on two different numerical methods for constructing rapidly rotating relativistic stars. A code based on the Komatsu-Eriguchi-Hachisu (KEH) method (Komatsu et al. 1989), written by Stergioulas, is compared to the Butterworth-Ipser code (BI), as modified by Friedman, Ipser and Parker. We compare models obtained by each method and evaluate the accuracy and efficiency of the two codes. The agreement is surprisingly good. A relatively large discrepancy recently reported (Eriguchi et al. 1994) is found to arise from the use of two different versions of the equation of state. We find, for a given equation of state, that equilibrium models with maximum values of mass, baryon mass, and angular momentum are (generically) all distinct and either all unstable to collapse or are all stable. Our implementation of the KEH method will be available as a public domain program for interested users.

0 citations0 reviewsNo code linkedOpen access