Post-Newtonian models of binary neutron stars

Using an energy variational method, we calculate quasiequilibrium configurations of binary neutron stars modeled as compressible triaxial ellipsoids obeying a polytropic equation of state. Our energy functional includes terms both for the internal hydrodynamics of the stars and for the external orbital motion. We add the leading post-Newtonian (PN) corrections to the internal and gravitational energies of the stars, and adopt hybrid orbital terms which are fully relativistic in the test-mass limit and always accurate to PN order. The total energy functional is varied to find quasiequilibrium sequences for both corotating and irrotational binaries in circular orbits. We examine how the orbital frequency at the innermost stable circular orbit depends on the polytropic index [Formula presented] and the compactness parameter [Formula presented]. We find that, for a given [Formula presented], the innermost stable circular orbit along an irrotational sequence is about 17% larger than the innermost secularly stable circular orbit along the corotating sequence when [Formula presented], and 20% larger when [Formula presented]. We also examine the dependence of the maximum neutron star mass on the orbital frequency and find that, if PN tidal effects can be neglected, the maximum equilibrium mass increases as the orbital separation decreases.