arxiv: v2 [gr-qc] 1 Oct 2009

Transcription

1 A formalism for te construction of binary neutron stars wit arbitrary circulation Tomas W. Baumgarte 1, and Stuart L. Sapiro 2, 1 Department of Pysics and Astronomy, Bowdoin College, Brunswick, ME Department of Pysics, University of Illinois at Urbana-Campaign, Urbana, IL Most numerical models of binary stars in particular neutron stars in compact binaries assume te companions to be eiter corotational or irrotational. Eiter one of tese assumptions leads to a significant simplification in te ydrodynamic equations of stationary equilibrium. In tis paper we develop a new formalism for te construction of binary stars wit circulation intermediate between corotational and irrotational. Generalizing te equations for irrotational flow we cast te Euler equation, wic is an algebraic equation in te case of corotational or irrotational fluid flow, as an elliptic equation for a new auxiliary quantity. We also suggest a parameterized decomposition of te fluid flow tat allows for a variation of te stellar circulation. arxiv: v2 [gr-qc] 1 Oct 2009 PACS numbers: I. INTRODUCTION Significant effort as gone into te numerical modeling of compact binaries containing neutron stars see, e.g., [1] for a review). Quasiequilibrium models of suc binaries usually assume te existence of an approximate elical Killing vector tat generates te binary orbit. Te existence of suc an approximate Killing vector implies tat, in a coordinate system corotating wit te binary, te binary companions appear approximately stationary. However, tis notion leaves te rotation, or circulation, of te individual stars undetermined. For example, eac star could always sow te same side to its companion, making it corotational, or it could ave zero circulation as seen in an inertial frame, making it irrotational. In principle, of course, eac star could be rotating arbitrarily, about an arbitrary axis, and potentially even differentially. As long as neutron stars can be modeled as ideal fluids, constructing neutron star solutions amounts to finding solutions to te Euler equation and te continuity equation. For equilibrium solutions, te Euler equation simplifies in two special cases, namely for corotational and for irrotational binaries. In bot cases te Euler equation reduces to a total derivative and can be integrated once to yield an algebraic equation for te fluid entalpy. For corotational binaries te continuity equation is satisfied identically, wic simplifies te resulting equations significantly. Not surprisingly, tis assumption was adopted in te first quasiequilibrium studies of bot binary neutron stars see, e.g., [2, 3]) and black oleneutron star binaries e.g. [4]). Corotational flow does not seem very realistic, owever, since tidally locking a neutron star in a compact binary would require an un- Also at Department of Pysics, University of Illinois at Urbana- Campaign, Urbana, IL Also at Department of Astronomy and NCSA, University of Illinois at Urbana-Campaign, Urbana, IL pysically large viscosity see [5, 6]). Instead, it seems more realistic to assume te fluid to be irrotational, wic is a reasonable assumption as long as te orbital frequency is muc larger tan te spin frequency tat te neutron stars would ave in isolation. In tis case te velocity may be written as te gradient of a velocity potential. Te continuity equation ten becomes an elliptic equation for tis velocity potential. Tis formalism was developed independently by [7, 8, 9, 10] see also [11] for a comparison between te different approaces), and numerical solutions ave been constructed bot for binary neutron stars e.g. [12, 13, 14], see also [15] for Newtonian models) and black ole-neutron star binaries e.g. [16, 17, 18]). Constructing binaries containing rapidly rotating stars, e.g. millisecond pulsars, requires a formalism tat allows for arbitrary spin configurations. In tis case te Euler equation no longer yields a trivial first integral. So far, neutron stars wit arbitrary circulation ave been constructed only by [19], ereafter MS. MS make several assumptions in teir treatment of neutron stars wit arbitrary circulation. Instead of using te Euler equation, tey solve te Bernoulli equation, wic states tat a certain quantity is constant along fluid trajectories. MS ten find approximate solutions by assuming tat tis quantity is constant globally, wic turns out to be correct for corotational and irrotational fluid flow, but not in general. Wile tis approximation may be quite accurate, it is of interest, at te very least as a point of principle, to identify a formalism tat avoids tis approximation. MS furter assume tat a solenoidal part of te fluid flow is centered on te stellar centers, determined by te point of maximum baryonic density. Tis may also be a good approximation, and yields corotational and irrotational models in good agreement wit oter approaces, but again it would be desirable to develop a formalism tat avoids tis assumption. In tis paper we develop a formalism tat improves on tat of MS in several ways. We suggest a way in wic te Euler equation can be solved for circulation intermediate between corotational and irrotational, and specify te fluid flow in a geometrically motivated fasion

2 2 tat does not make any assumptions about te centers of rotation. After stating some preliminaries in Section II we will, for te sake of clarity, first develop tis formalism in a Newtonian framework in Section III. We will ten treat te problem in a relativistic framework in Section IV, and will briefly summarize in Section V. In Appendix A we will briefly discuss an alternative decomposition of te fluid flow. II. PRELIMINARIES We adopt a 3+1 decomposition of te spacetime metric and write te spacetime metric g ab in te form g ab dx a dx b = α 2 dt 2 +γ ij dx i +β i dt)dx j +β j dt). 1) Here α is te lapse function, β i is te sift vector and γ ij is te spatial metric induced on a spatial slice Σ of constant coordinate time t. Te normal on tis slice is n a = α 1 1, β i ). 2) Here and in te following we denote spacetime indices wit letters a,b,c,... and spatial indices wit letters i,j,k,... Te coordinate congruence t a, wic connects points wit te same spatial coordinate labels on different slices Σ, is ten given by t a = αn a +β a. We also assume te existence of a elical Killing vector ξel a, and assume it to be normalized in suc a way tat its time component is unity, ξel t = 1. We can write tis Killing vector as ξ a el = ta +k a, 3) were te vector k a is purely spatial and describes te orbital rotation ) a k a = Ω. 4) φ In cartesian coordinates, for a rotation about te z-axis, we may also write te spatial components of k a as k i = Ω y,x,0). 5) III. NEWTONIAN TREATMENT A. Newtonian equations of stationary equilibrium Te Newtonian equations of ydrodynamics are te Euler equation t v i +v j D j v i = ρ 1 0 P φ N 6) and te continuity equation t ρ 0 + ρ 0 v i ) = 0. 7) Here v i is te fluid velocity, ρ 0 te rest-mass density, P te pressure, and φ N te Newtonian potential, wic as to satisfy te Poisson equation D 2 φ N = 4πρ 0. 8) We denote a covariant derivative wit, and te Laplace operator wit D 2. It is also useful to define V i as te fluid velocity relative to k a, v i = k i +V i, 9) i.e. V i is te velocity as seen in a frame corotating wit te binary. For isentropic flow we also ave ρ 1 0 P =. 10) Here is te specific entalpy, wic, adopting a convention more common in a relativistic framework, is defined as = 1+ǫ+P/ρ 0, 11) were ǫ is te specific internal energy density. We also assumetat wearegivenanequationofstatetat relates te entalpy to te density ρ 0. For stationary solutions te Lie derivative of te ydrodynamic variables ρ 0 and v i as to vanis, i.e. L ξel v i = t v i +L k v i = t v i +k j D j v i v j k j = 0 12) and similar for ρ 0. Combining equation 12) wit te Euler equation 6) we ten obtain 1 2 v2 +v j V j +φ N )+V j D j v i v j ) = 0, wile te equation of continuity becomes 13) ρ 0 V i ) = 0. 14) Equations 13) and 14) form te Newtonian equations of stationary equilibrium. B. Corotational and irrotational flow It is easy to see ow te above equations simplify for corotational or irrotational fluid flow. For corotational flow, we must ave V i = 0, in wic case te continuity equation 14) is satisfied identically. Furtermore, te second term in te Euler equation 13) vanises identically, so tat te first term can be integrated to yield an algebraic equation for te entalpy, 1 2 v2 +φ N = C. 15) Here C is some constant of integration, and v 2 can be found from 9), wic now implies v i = k i. Here and

3 3 below we assume tat bot stars are identical for simplicity; te metod is easily generalized for non-identical stars. For irrotational flow, we require te fluid s vorticity to vanis, wic is automatically te case if we write te fluid velocity as a gradient of a velocity potential Φ, v i = Φ. 16) In tis case, te second term in te Euler equation 13) also vanises identically, so tat we again obtain an algebraic equation for te entalpy, 1 2 v2 +v j V j +φ N = C. 17) Te continuity equation 14) now becomes an elliptic equation for te velocity potential Φ, D 2 Φ = k i Φ) lnρ 0. 18) To ensure tat te rigt and side remains finite at te stellar surface, were ρ 0 becomes zero, we must ave k i Φ) ρ 0 surf = 0. 19) Given tat ρ 0 is normal to te stellar surface, tis relation furnises a Neumann boundary condition for Φ. We could ave derived te same condition by requiring tat, at te stellar surface, te fluid flow be tangent to te surface. C. Intermediate circulation For intermediate circulation, te second term in te Euler equation 13) no longer vanises, and we can no longer find a simple first integral. One approac to find a solution is to take a divergence of te equation, wic yields an elliptic equation D 2 H + V j D j v i v j ) ) = 0 20) for a new auxiliary quantity H 1 2 v2 +v j V j +φ N. 21) In te following we will interpret tis as an equation for te specific entalpy, = H v2 v j V j φ N. 22) Before proceeding we need to specify te circulation of te fluid. To do so, we coose to decompose te velocity v i according to wic, from 9), implies v i = Φ+ηk i, 23) V i = Φ+η 1)k i. 24) Here Φ is a generalized velocity potential, and η is an arbitrary constant. As we will discuss in Section III D we recover irrotational flow for η = 0 and corotational flow for η = 1. Wit our decomposition 23), η terefore parameterizes a sequence of stellar models connecting corotational and irrotational models. By generalizing te decomposition 23) it sould be possible to cange te axis of te internal rotation, or cange te degree of differential rotation. Inserting te decomposition 23) into equation 20) we first note tat te velocity potential Φ again drops out of te term D j v i v j = ηd j k i k j ) ηκ ij, 25) were we avedefined κ ij for convenience. Evidently, κ ij is antisymmetric by construction and can be computed directly from 5). Te second term in equation 20) ten becomes V j D j v i v j ) ) = η V j )κ ij +V j κ ij ). 26) Inserting 24) to evaluate te first term on te rigt and side we find V j )κ ij = D j Φ+η 1) k j )κ ij = η 1) k j )κ ij, 27) were te second equality olds because D j Φ is symmetric wile κ ij is antisymmetric. Tat means tat all second derivatives of Φ disappear in equation 20), leaving us wit D 2 H+ηη 1) k j )κ ij +η D j Φ+η 1)k j) κ ij = 0. 28) Tis is te form of te Euler equation tat generalizes equations 15) or 17) for intermediate circulation. Instead of an algebraic equation for te entalpy we now ave to solve an elliptic equation for te new auxiliary quantity H. Given solutions for te oter fields we can ten find te entalpy from 22). On te stellar surface te density vanises, meaning tat approaces unity. Tis requirement leads to te Diriclet boundary condition H surf = v2 +v j V j +φ N 29) on te stellar surface. To find an equation for te velocity potential Φ we insert 24) into te continuity equation 14), wic yields D 2 Φ = Φ+η 1)k i) lnρ 0. 30) To ensure tat te rigt and side remains finite on te stellar surface we must ave Φ+η 1)k i) ρ 0 surf = 0. 31) As in te case of irrotational flow, tis conditions forms a Neumann boundary condition.

4 4 We ave now assembled all equations and boundary conditions needed to construct binary stars wit circulation intermediate between corotational and irrotational flow in a Newtonian framework. In addition to te Poisson equation 8), we need to solve equation 28), subject to te boundary condition29), for H, and equation30), subject to boundary condition 31), for Φ. Given tese solutions we can find v i and V i from equations 23) and 24), and from 22). Finally, ρ 0 can be computed from given an equation of state, and Ω is determined by enforcing a circular orbit see, e.g., [15]). D. Limiting cases Before proceeding to a relativistic treatment of te problem it is useful to see ow te equations for corotational and irrotational flow appear as limiting cases of te equations for intermediate circulation. Clearly, our equations for intermediate circulation are consistent wit equations 13) and 14), and are terefore satisfied bot for corotational flow wit V i = 0 and irrotational flow wit v i = Φ. In te following we demonstrate tat tese solutions are also consistent wit te decomposition 23). We start wit irrotational flow, wic we obtain for η = 0. In tis limit, te decomposition 23) evidently reduces to16), and bot te equation30) and boundary condition 31) reduce to teir equivalents for irrotational flow, equations 18) and 19). For η = 0, equation 28) reduces to a Laplace equation for H, wic is uniquely solved by H = C H, were C H is a constant, wen H = C H on te boundary. Tis is consistent wit te oter equations, wic can be seen by evaluating equation 17) on te stellar surface. As a result, te equation for te entalpy reduces to equation 17), as expected. We obtain corotational flow by setting η = 1. To see tis, we first observe tat bot te equation 30) for te velocity potential Φ and its boundary condition 31) are solvedbyφ = C Φ, werec Φ isaconstant. Equation24) ten implies V i = 0, as expected for corotational flow. Equation28) again reduces to a Laplace equation wic, as in te irrotational case, is solved by H = C H. Te equation for te entalpy terefore reduces to equation 15). IV. RELATIVISTIC TREATMENT A. Relativistic equations of stationary equilibrium Te relativistic equations of ydrodynamics follow from te conservation of energy-momentum, b T ab = 0, 32) wose spatial components form te relativistic Euler equation, and te conservation of rest-mass a ρ 0 u a ) = 0. 33) Here T ab is te fluid s energy-momentum tensor, u a its four-velocity, and a denotes a covariant derivative associated wit te four-metric g ab. Following te approac of [10] we can write u a as u a = u t ξ a el +Va ), 34) were te elical Killing vector ξ a el was introduced in equation 3) and were V a is purely spatial, n a V a = 0. Just like its Newtonian counterpart in equation 9), V a now measures te fluid velocity in a frame comoving wit te binary. It is also convenient to define û i = γ i a u a. 35) In terms of û i te normalization of te four-velocity u a u a = 1 may be written as αu t = 1+ 2 û i û i) 1/2. 36) For stationary solutions te Lie derivative of te ydrodynamical variables along ξ a el again as to vanis. As sown in [10], te relativistic Euler equation ten becomes u t +û jv j )+V j D j û i û j ) = 0, 37) wile te continuity equation becomes αu t ρ 0 V i) = 0. 38) Here is te covariant derivative associated wit te spatial metric γ ij. We also assume tat, on te slice Σ, te spatial metric γ ij, te lapse α, and te sift β i are constructed by solving some subset of Einstein s equations, for example by using te conformal tin-sandwic decomposition. Tese equations ten play te role of te Poisson equation 8) in te Newtonian framework. As before we furter assume tat we are given an equation of state tat relates te specific entalpy to te rest density ρ 0. B. Corotational and irrotational flow As in te Newtonian case it is easy to see ow te equations simplify for bot corotational and irrotational flow. For corotational flow te four-velocity u a is aligned wit te elical Killing vector ξ a el, so tat V a = 0. Te only remaining term in te Euler equation 37) can ten be integrated to yield = C, 39) ut were C is some constant, wile te continuity equation 38) is satisfied identically.

5 5 To obtain irrotational flow, wit vanising relativistic fluid vorticity, we now assume û i = Φ, 40) were Φ is a velocity potential see [10] for a discussion). Wit tis coice te second term in te Euler equation 37) vanises again, and we can find te first integral u t +û iv i = C. 41) Combining equations 34), 36) and 41) we also find αu t = 1 α C +Bi Φ), 42) were we ave defined a rotational sift vector B a β a +k a. 43) Using 3) and t a = αn a +β a we may also write B a as B a = ξ a αn a. 44) We will derive a more general version of equation 42) for intermediate circulation below. Note tat V i can always be expressed as V i = γ i av a = γ i a ) 1 u tua ξel a For irrotational flow we can now insert 40), = 1 u t ûi B i. 45) V i = 1 u t Di Φ B i, 46) and insert tis result into te continuity equation 38) to obtain αρ 0 1 Φ) αu t ρ 0 B i ) = 0. 47) Wemayeliminateu t fromtisequationbyinserting42), wic yields C +B D 2 j ) D j Φ Φ B i α 2 C +B j ) D j Φ = α 2 B i Φ ln αρ 0. 48) As expected, tis equation reduces to 18) in te Newtonian limit. A Neumann boundary condition again results from requiring tat te rigt and side remain finite at te stellar surface, C +B j D j Φ α 2 ) B i Φ αρ 0 = 0. 49) surf To complete te argument, we could insert equations 40), 42) and 46) into 41) to obtain an equation for in terms of Φ, and independently of u t. C. Intermediate circulation As in te Newtonian case, we observetat for intermediate circulation te second term in te Euler equation 37) no longer vanises, so tat we can no longer find a first integral. Taking a divergence of te equation, we can again obtain an elliptic equation D 2 H + V j D j û i û j ) ) = 0 50) for a new auxiliary function H u t +û iv i. 51) We will again interpret tis as an algebraic equation for te specific entalpy. Generalizing te Newtonian expression23) we decompose û i according to û i = Φ+ηk i. 52) Repeating te steps in equation 46) we now find V i = 1 u t Di Φ+ηk i ) B i. 53) We point out tat tere are many different ways of generalizing te Newtonian decomposition 23), and we will discuss an attractive alternative to te decomposition 52) in Appendix A. Moreover, te Newtonian decomposition 23) represents only one particular sequence tat connects corotational and irrotational fluid flow; it sould also be possible to generate models tat rotate about an arbitrary axis of rotation. Similarly, te relativistic decomposition 52) represents only one possible sequence. Before proceeding we derive an equation for u t. Starting wit te definition 51) we insert equation 45) to obtain ) 1 u t +û i ûi u t B i = H, 54) or 1+ 2û u t i û i) = H +û i B i. 55) We now insert te normalization 36) as well as te decomposition 52) to find αu t = 1 H +B i ηk i + Φ) ). 56) α For te irrotational case, wen η = 0 and H = C, we recover equation 42), as expected. Returning to te Euler equation 50), we observe tat D j û i û j = ηd j k i k j ) ηκ ij, 57) in complete analogy to te Newtonian case. Expanding te terms in 50) we now find D 2 H + η D j Φ+ηk j ) 1 u t + η ) u t Di k j B j κ ij ) 1 + η u t Dj Φ+ηk j ) B j κ ij = 0. 58)

6 6 If desired, equation 56) can be inserted to eliminate te terms u t. Requiring tat te density vanis on te stellar surface, so tat = 1, again leads to a Diriclet boundary condition H surf = 1 u t +û iv i. 59) We now turn to te continuity equation 38). Inserting te decomposition 53), we obtain αρ0 ) Di Φ+ηk i ) αu t ρ 0 B i ) = 0, 60) or, eliminating u t wit te elp of equation 56), H +B D 2 j ) ηk j +D j Φ) Φ α 2 B i = η k i 61) H +B j ) ηk j +D j Φ) + α 2 B i Φ ηk i ln αρ 0 As before we require tat te rigt and side remain finite at te stellar surface, wic yields te boundary condition H +B j ηk j +D j Φ) α 2 ) B i Φ ηk i αρ 0 = 0. surf 62) We ave now assembled a complete set of equations and boundary conditions tat can be used to construct relativistic binary stars wit intermediate circulation. Given a coice for te parameter η we can solve equation 58) for H and equation 61) for Φ. We can ten find te combination u t from 56), and te velocity expressions V i and û i from equations 52) and 53). Finally, can be found from te definition 51), ρ 0 from te equation of state, and Ω can be determined by enforcing a circular orbit see [1] and references terein). Te above equations ave to be solved togeter wit gravitational field equations of te lapse α, te sift β i wic, togeter wit te elical Killing vector 3), determines te rotational sift vector B i ) and te spatial metric γ ij ; tis can be accomplised, for example, by solving te constraint equations of Einstein s field equations in te conformal tin-sandwic decomposition. All te above equations are coupled in complicated ways, and presumably te easiest way of constructing numerical solutions involves an iterative algoritm. D. Limiting cases Our equations are based on te Euler equation 37) and te continuity equation 38), wic are clearly satisfied for bot corotational flow wit V i = 0 and irrotational flow wit û i = Φ. However, it is less clear weter tese solutions are consistent wit our decomposition 52) and 53). As in te Newtonian case we can recover irrotational flow for η = 0. In tis limit, all equations, including a constant solution H = C H for H, are consistent wit te corresponding equations for irrotational flow in Section IVB. On te oter and, we do not seem to recover exactly corotational flow wit V i = 0 for η = 1 for te decomposition 53), since it is not clear ow, for η = 1, we would obtain a solution for Φ tat exactly cancels te oter terms in 53). As we ave seen in Section IIIB, we do obtain corotational flow for η = 1 in te Newtonian limit; for relativistic solutions, deviations from corotation can terefore be at most at te post-newtonian order. Our decomposition 52) terefore represents a sequence tat connects irrotational flow wit approximately corotational flow. In Appendix A we present an alternative decomposition of te velocity field, tat does reproduce exactly corotational flow for η = 1. However, tis decomposition also introduces several additional derivative terms tat may prove problematic, as we discuss in Appendix A. V. SUMMARY AND DISCUSSION In tis paper we propose a formalism for te construction of binary stars wit intermediate circulation. Adopting te approac of [10] we generalize te equations for irrotational fluid flow and sow tat stars wit intermediate circulation can be constructed by solving one additional elliptic equation. We also propose a decomposition of te fluid flow tat represents a sequence connecting irrotational stars wit corotational ones in te Newtonian case, and approximately corotational ones in te relativistic case. Exploring te pysical properties of tese sequences will require future numerical work. It sould also be possible to generalize te decomposition of te fluid flow to allow for rotation about oter axes of rotation, or to modify te degree of differential rotation. Our approac may also be used to construct neutron stars in quasiequilibrium orbit about black oles. Relaxing te assumption of eiter corotational or irrotational flow in binary stars is, at te very least, interesting as a point of principle. Suc a formalism is necessary, owever, for te construction of binaries containing rapidly rotating stars, for example millisecond pulsars, wen te binary s orbital frequency is no longer muc larger tan te spin frequency of te star. As suggested by MS, approximate inspiral sequences of suc binaries can ten be constructed by fixing te circulation along one particular fluid trajectory, for example along te stellar equator. Tis can be acieved by iterating, at a given binary separation, over te parameter η in te decomposition of te fluid flow until te desired value of te circulation as been realized.

7 7 Acknowledgments Tis work was supported in part by NSF grant PHY to Bowdoin College, NASA grants NNG04GK54G and NNX07AG96G to UIUC, and NSF grants PHY , PHY and PHY to UIUC. = 1 η) H ) +Bj D j Φ α 2 ηb k B i Φ ln αρ 0 B k. Taking a divergence of te Euler equation 37) results in D 2 H +η V j κ ij ) = 0, A5) APPENDIX A: AN ALTERNATIVE DECOMPOSITION OF THE FLUID VELOCITY were we ave now defined As we ave seen in Section IVD, te decomposition 52) does not necessarily lead to exactly corotational fluid flow in te limit η = 1. An attractive alternative to tis decomposition is wic leads to û i = Φ+ηu t B i, V i = 1 u t Di Φ+η 1)B i. A1) A2) Retracing te derivations of Section IVC we find tat equation 56) can now be written u t = H +Bi Φ α 2 ηb i B i, wile te continuity equation 38) becomes H +B D 2 j ) D j Φ Φ + η 1) α 2 ηb k B i B k A3) A4) κ ij D j u t B i ) u t B j ). A6) We first note tat, for η = 0, te above equations reduce to tose for irrotational flow as before. For η = 1, equation A4), togeter wit its boundary condition, is now solved by Φ = C Φ, were C Φ is a constant. Tat implies tat, unlike our original decomposition 52), te alternative decomposition A1) does reproduce corotational flow wit V i = 0 in te limit η = 1. Wile tis is a very attractive feature, tis decomposition also as disadvantages. Note tat te term u t contains first derivatives of Φ, meaning tat te term κ ij in equation A5) will introduce tird derivatives of Φ. Using te decomposition 52), on te oter and, we avoid tese tirdderivative terms. Te appearance of u t in equation 58) still introduces second derivatives of Φ. All of tese terms disappear in te Newtonian limit, owever, are terefore small, and opefully do not affect te well-beavedness of te system. [1] T. W. Baumgarte and S. L. Sapiro, Pys. Rept. 376, ). [2] T. W. Baumgarte, G. B. Cook, M. A. Sceel, S. L. Sapiro, and S. A. Teukolsky, Pys. Rev. Lett. 79, ). [3] T. W. Baumgarte, G. B. Cook, M. A. Sceel, S. L. Sapiro, and S. A. Teukolsky, Pys. Rev. D 57, ). [4] T. W. Baumgarte, M. L. Skoge, and S. L. Sapiro, Pys. Rev. D 70, ). [5] C. S. Kocanek, Astropys. J. 398, ). [6] L. Bildsten and C. Cutler, Astropys. J. 400, ). [7] S. Bonazzola, E. Gourgoulon, and J.-A. Marck, Pys. Rev. D 56, ). [8] H. Asada, Pys. Rev. D 57, ). [9] S. A. Teukolsky, Astropys. J. 504, ). [10] M. Sibata, Pys. Rev. D 58, ). [11] E. Gourgoulon 1998), gr-qc/ [12] S. Bonazzola, E. Gourgoulon, and J.-A. Marck, Pys. Rev. Lett. 82, ). [13] E. Gourgoulon, P. Grandclément, K. Taniguci, J.-A. Marck, and S. Bonazzola, Pys. Rev. D 63, ). [14] K. Taniguci and E. Gourgoulon, Pys. Rev. D 66, ). [15] K. Uryū andy. Eriguci, Mon. Not. R. Astron. Soc. 303, ). [16] K. Taniguci, T. W. Baumgarte, J. A. Faber, and S. L. Sapiro, Pys. Rev. D 72, ). [17] K. Taniguci, T. W. Baumgarte, J. A. Faber, and S. L. Sapiro, Pys. Rev. D 77, ). [18] P. Grandclément, Pys. Rev. D 74, ). [19] P. Marronetti ands.l. Sapiro, Pys.Rev.D68, ).

8 8 Erratum It was kindly brougt to our attention by E. Gourgoulon tat te formalism tat we propose to construct binary neutron stars wit arbitrary circulation solves only a subset of te equilibrium equations in general, but not all of tem. Specifically, our formalism considers only te divergence of te Euler equation for stationary equilibrium, Eq. 13) in te Newtonian case and Eq. 37) in te relativistic case. However, for te Euler equation to be satisfied, bot its divergence and its curl as to vanis. For te limiting cases of corotational or irrotational fluid flow te curl of tese equations does vanis identically, but tis is not necessarily te case for intermediate circulation. Wile our formalism still leads to valid binary initial data for fluid flow wit intermediate circulation, tese data are not necessarily in exact equilibrium for te Newtonian case or quasiequilibrium in te relativistic case. Given tat te curl vanises for bot limiting cases, it is possible tat te curl is small for intermediate configurations as well, meaning tat tese data may represent good approximations to equilibrium configurations. It may also be possible to construct equilibrium or quasiequilibrium) data by adding anoter contribution to te decomposition of te fluid velocity, Eq. 23) in te Newtonian case and Eq. 52) in te relativistic case, and solving for tis contribution by forcing te curl of te Euler equation to vanis.