GWs from Eccentric Neutron Star Binaries

# Stochastic Background of Gravitational Waves Generated by Eccentric Neutron Star Binaries

## Abstract

Binary systems emit gravitational waves in a well-known pattern; for binaries in circular orbits, the emitted radiation has a frequency that is twice the orbital frequency. Systems in eccentric orbits, however, emit gravitational radiation in the higher harmonics too. In this paper, we are concerned with the stochastic background of gravitational waves generated by double neutron star systems of cosmological origin in eccentric orbits. We consider in particular the long-lived systems, that is, those binaries for which the time to coalescence is longer than the Hubble time (Gyr). Thus, we consider double neutron stars with orbital frequencies ranging from to Hz. Although in the literature some papers consider the spectra generated by eccentric binaries, there is still space for alternative approaches for the calculation of the backgrounds. In this paper, we use a method that consists in summing the spectra that would be generated by each harmonic separately in order to obtain the total background. This method allows us to clearly obtain the influence of each harmonic on the spectra. In addition, we consider different distribution functions for the eccentricities in order to investigate their effects on the background of gravitational waves generated. At last, we briefly discuss the detectability of this background by space-based gravitational wave antennas and pulsar timing arrays.

###### keywords:
gravitation – gravitational waves – binaries: close – stars: neutron.
23

## 1 Introduction

It is well known that binary systems emit gravitational waves (GWs) and, for circular orbits, the frequency of the emitted radiation is twice the orbital frequency. On the other hand, for eccentric orbits, the pattern of emission has a more complex form: the emitted signals are also formed by the higher harmonics of the orbital frequency. In fact, the higher the orbital eccentricity is, the greater the contribution of the higher harmonics to the emitted radiation (see e.g., Peters & Mathews, 1963; Douglass & Braginsky, 1979).

Since many binaries can have non-negligible orbital eccentricities (Kowaska et al., 2012), it is interesting to investigate the behaviour of the spectra when eccentricity is taken into account. Particularly, this paper is addressed to the stochastic background of GWs generated by long-lived double neutron star (DNS) systems of cosmological origin in eccentric orbits. We consider that a system is long-lived if its lifetime is longer than the Hubble time, namely Gyr (Kowaska et al., 2012). Quantitatively, this assumption implies that the binaries have orbital frequencies in the range Hz, as explained in Section 5. We should bear in mind that in this case we can neglect the time evolution of the orbital parameters of the systems.

It is worth noting that the inclusion of eccentricity in the calculation of stochastic backgrounds generated by compact binaries is not new. In fact, one can find in the literature interesting papers on the stochastic backgrounds generated by eccentric binaries, such as Ignatiev et al. (2001), Kowaska et al. (2012) and Enoki & Nagashima (2007). In Ignatiev et al. (2001), the authors study, for example, the GW background produced by coalescing neutron stars in the Galaxy in eccentric orbits; in Kowaska et al. (2012) the authors study the backgrounds generated by Population III binaries, giving special attention to the spectra generated by long-lived and short-lived binaries; in Enoki & Nagashima (2007), the authors studied the background generated by supermassive black hole binaries.

In general, the procedure to obtain the backgrounds generated by eccentric binaries includes, as a key component, the spectra of energy emitted by the source, which is written as , where is the energy and is, in our case, the orbital frequency. In the literature, the authors usually handle such a quantity following the steps: first, they sum the spectra corresponding to each harmonic, where is the number of the harmonic, in order to obtain the total spectrum ; subsequently, is included in the formula that determines the stochastic background. Here we use a different method: we calculate the backgrounds that would be generated by each harmonic of the emitted frequency and then we sum all these signals in order to obtain the total background. As we will see, this method allows us to clearly analyse the influence of each harmonic on the resulting backgrounds.

In addition, we consider in our calculations different distribution functions for the orbital eccentricity , and study the influence of such distributions on the resulting backgrounds. Particularly, we considered the case in which the distribution function of the eccentricities has the form of a Dirac delta function, that is, all the systems have the same eccentricity. We consider scenarios in which all systems have eccentricities of 0.25, 0.5, 0.75 and 0.9. Moreover, we consider the case of circular orbits (null eccentricity) and also another distribution that covers a wide range of eccentricities, namely, a uniform distribution ranging from to .

The paper is organized as follows: in Section 2, we derive the fundamental equations used to calculate the stochastic background of GWs generated by DNS systems of cosmological origin in eccentric orbits; in Section 3, we calculate the spectrum of energy for eccentric systems; in Section 4, we present and discuss the distribution functions of the eccentricities adopted; in Section 5, we discuss the maximum and minimum frequencies of the spectra and the population characteristics of the systems, such as the star formation rate density and the mass fraction of stars that is converted in DNSs; we present the results in Section 6 and the conclusions and perspectives in Section 7.

## 2 Fundamental Equation

For the calculation of the spectra generated by each harmonic, we follow the method described in Regimbau (2011) and Regimbau & Chauvineau (2007), who calculated the backgrounds in terms of the energy density parameter , namely

 Ωgw(ν0)=1ρcc3ν0F(ν0), (1)

where is the critical density of the universe, is the Hubble parameter, is the observed frequency and is the spectral density of the flux, which is given by (Ferrari et al., 1999a, b)

 F(ν0)=∫f(ν0,z)dRdzdz, (2)

where is the redshift and is the energy flux per unit frequency generated by a unique source, which has the form (Regimbau, 2011)

 f(ν0,z)=14πr(z)2dEdν∣∣∣ν=ν0(1+z). (3)

In equation (3), is the comoving distance, is the emitted frequency and is the spectrum of energy emitted by a binary system, calculated in the source frame. In addition, in equation (2), the event rate is calculated by

 dRdz=λ\scriptsize{dns}˙ρ∗(z)1+zdVdz, (4)

where is the star formation rate density (SFRD), is the mass fraction of stars that is converted into DNSs, and is given by

 dVdz=4π(cH0)r(z)2√ΩM(1+z)3+Ωk(1+z)2+ΩΛ. (5)

We use equation (5) in a scenario in which the density parameters obey the relations and , with the subindices M, DM, B, k and corresponding to matter, dark matter, baryonic matter, curvature and cosmological constant, respectively. Here, we consider a cosmological scenario in which , , and .

In addition, the eccentricity is introduced in the calculation of the background by means of a distribution function in equation (1). So, with the above considerations and performing all the substitutions, equation (1) assumes the following form

 Ω\scriptsize{gw}(ν0,ε)=8πG3H30c2λ\scriptsize{dns}ν0 (6) ×∫zf(ν0)zi(ν0)˙ρ∗(z)(1+z)E(z)dEdν∣∣∣ν=ν0(1+z)f(ε)dz,

where the meaning of the limits of integration, and , will be considered in Section 5.

## 3 Energy Emitted by Eccentric Binaries

In the previous section, we noticed that the spectrum of energy is a key ingredient for the calculation of the backgrounds, since it provides information concerning the patterns of emission of the systems as a function of the orbital frequency. We can now obtain an expression for the energy spectrum for eccentric systems starting from

 dEdνn=dEdtdtdΩdΩdνn, (7)

where we included the index in the notation, since equation (7) refers to the spectrum generated by each harmonic. In addition, is the angular orbital frequency and is given by

in which ’’ represents the orbital semimajor axis. In our case, and are given by (Peters & Mathews, 1963)

 dEdt=325G4c5m21m22(m1+m2)a5g(n,ε); (9)

and

where and are the masses of the components of the system. In addition, the functions and have the form

 g(n,ε)=n432{[Jn−2(nε)−2εJn−1(nε)+2nJn(nε)+ + (11) +2εJn+1(nε)−Jn+2(nε)]2+(1−ε2)[Jn−2(nε)++ −2εJn(nε)+Jn+2(nε)]2+43n2[Jn(nε)]2},

where are the Bessel functions of order , and

 F(ε)=1+73ε2/24+37ε4/96(1−ε2)7/2. (12)

Moreover, equations (11) and (12) are related to each other by

 F(ε)=∞∑n=2g(n,ε), (13)

and is obtained from Kepler’s third law.

Now, recall that systems in eccentric orbits emit radiation at harmonics of the orbital frequency, which is related to the angular orbital frequency by . Thus, we can write the emitted frequency of the harmonic as a function of and , namely

 νn=nΩ2πforn≥2, (14)

which enable us to calculate in equation (7).

Now, performing all the necessary substitutions, equation (7) is transformed into

 dEdνn=13(2πGn)2/3m1m2(m1+m2)1/3ν−1/3ng(n,ε)F(ε). (15)

Note that for , and , equation (15) comes down to the expression valid for circular orbits (see e.g., Ferrari et al., 2001).

In addition, has the following meaning: it provides the fraction of energy emitted at the harmonic by a system with orbital eccentricity .

## 4 The Distribution Functions of the Eccentricities

We are interested in studying the role of the distribution function of the orbital eccentricities on the stochastic background of GWs generated by DNS systems of cosmological origin. We consider distributions that, although could be considered unrealistic, can provide an insight into the influence of these functions on the spectra.

First, we consider a scenario where all the DNSs have the same orbital eccentricity, such that

 f(ε)=δ(ε−ε′), (16)

where is a Dirac delta function. Particularly, we consider the cases for (a population of circular DNSs), 0.25, 0.5, 0.75 and 0.9.

We also consider a distribution that covers a wide range of values of . For this purpose, we use a uniform and normalized function that is given by

 Missing dimension or its units for \hskip (17)

where and are the minimum and the maximum eccentricities, respectively. This distribution provides an interesting case to be compared with the distribution functions given in terms of the Dirac delta functions.

## 5 The Stochastic Background

Before proceeding, we pay attention to the choice of the minimum and maximum values of the orbital frequencies used in the calculations of the backgrounds. First, as already mentioned, we consider in this paper long-lived DNSs, which lead us to the determination of a maximum value for the initial frequency. We mean by a ’long-lived system’, a binary whose lifetime is at least of the order of the Hubble time, that is, Gyr.

In addition, recall that the greater the orbital eccentricity is, the shorter is the binary system lifetime. Therefore, we choose the maximum frequency such that the most eccentric systems considered in this paper have lifetimes equal to the Hubble time. The expression for , considering close to , is given by (Peters, 1964)

 T\scriptsize{life}≈768425a404β(1−ε20)7/2, (18)

where , is the initial orbital semimajor axis and is the initial eccentricity. Considering, for example, Gyr, and in equation (18) and using Kepler’s third law, we obtain a frequency of Hz.

On the other hand, according to Rosado (2011), there is a minimum frequency one must consider in the study of GWs generated by binary systems. He argues that the emitted radiation with frequencies lower than a given value should not be considered, since other mechanisms that take away energy from the systems are more effective (see Rosado 2011 for a detailed discussion.)

Still following the discussion of Rosado (2011), we should bear in mind two facts: it is difficult to define exactly the values of the minimum frequencies; on the other hand, such values are not important in practice, since the choice of a particular frequency will not have noticeable effects on the results. Considering these arguments, we adopted the rather arbitrary value of Hz in our calculations.

To proceed, the next step is to determine the SFRD and the value of .

Concerning the SFRD, many different proposals are considered in the literature, although they do not differ from each other very significantly. Here we adopt, as a fiducial one, that given by Springel & Hernquist (2003), namely

 ˙ρ∗(z)=ρmβeα(z−zm)β−α+αeβ(z−zm), (19)

where , , and fixing the normalization. It is worth mentioning that the authors considered a CDM cosmology, where the density parameters have the values given in Section 2. Moreover, they considered the Hubble parameter as with .

The mass fraction (see e.g., Regimbau & Pacheco, 2006, for a detailed discussion) is given by

 λ\scriptsize{dns}=β\scriptsize{ns}fpΦ\scriptsize{ns}, (20)

where is the fraction of binary systems that survive to the second supernova event, provides the fraction of massive binaries (that is, binary systems where both components could generate a supernova event) formed inside the whole population of stars, and is the mass fraction of neutron star progenitors. In our case, we assume that is calculated by

 Φ\scriptsize{ns}=∫258ϕ(m)dm, (21)

where is the Salpeter mass function (Salpeter, 1955), with , and . Moreover, note that we are considering that the progenitors of the neutron stars have masses in the range .

According to Regimbau & Pacheco (2006), the parameters and depend on different assumptions. First, in the evolutionary scenario of massive binaries the authors considered that none of the stars underwent recycle accretion. On the other hand, the velocity distribution of the natal kick affects the parameter , since an imparted kick could disrupt binaries or, less likely, prevent the disruption of the system. Particularly, the authors adopted a one-dimensional velocity dispersion of .

Still according to Regimbau & Pacheco (2006), the estimates of depend on itself and on the ratio between single and double NS systems in the Galaxy. As a result, there are uncertainties in the values of these parameters. Regimbau & Pacheco (2006) then estimate the relative uncertainties in and . Numerically, they obtained and , leading to a relative uncertain of in . The authors stated that these uncertainties are solely formal, resulting from the simulations they performed, which depend on the evolutionary scenario for the progenitors.

Concerning the range of masses of the neutron star progenitors, recall that the minimum and maximum values are subject of discussions. However, the limits we adopt in equation (21) are usual (see, e.g., Ferrari et al. 1999a, b). We refer the reader, for example, to Heger et al. (2003), Smartt (2009) and Carroll & Ostlie (2007), who discuss this subject in detail.

Numerically, we have , , and, consequently, .

For the limits of integration and in equation (6), we follow Regimbau (2011) and Regimbau & Chauvineau (2007), namely

 Missing dimension or its units for \hskip (22)

and

 Missing dimension or its units for \hskip (23)

where Hz, Hz, and are the minimum emitted frequency, the maximum emitted frequency, the minimum redshift and the maximum redshift, respectively.

Using the ingredients shown above, integrating over the eccentricity and substituting the numerical values, equation (6) takes the form

 Ωgw(n,ν0)=1.28×10−10(ν0n)2/3 (24) ×∫zf(ν0)zi(ν0)e0.6(z−5.4)0.333+0.6e0.933(z−5.4)dz(1+z)4/3E(z) ×∫ε\scriptsize{max}ε\scriptsize{min}g(n,ε)f(ε)F(ε)dε,

where and is given by either equation (16) or equation (17).

## 6 Results

We show in Figs 1-4 the spectra generated by some harmonics of the emitted signals, considering the distribution given by equation (16). It is worth mentioning that, although we show in these figures only some harmonics, we include in the calculations of the total spectra (Fig. 5) all relevant harmonics.

Note that for high eccentricities, the amplitudes corresponding to the higher harmonics become extremely relevant. Particularly, it is worth noticing the contrast between the cases corresponding to and .

We show in Fig.  5 the backgrounds generated when we sum the contributions of all relevant harmonics for the cases corresponding to Figs. 1-4. We also include in this very figure the spectra for: (a) the uniform distribution given by equation (17); (b) DNS systems in circular orbits.

Comparing the uniform curve in Fig.  5 with, for example, that by a population of DNSs in circular orbits, we can observe the effects of the orbital eccentricities: the frequency band of the background is broadened and shifted towards higher frequencies.

The frequency band of the spectrum for circular systems ranges from to Hz, whereas for the uniform distribution function the spectrum ranges from to Hz.

Also, the amplitudes of the spectra are slightly reduced in the region between and Hz whether the eccentricity of the population of DNSs is increased. This occurs due to the reduction of the amplitudes of the contributions corresponding to the lowest harmonics as compared to the highest ones.

This effect can be more clearly observed in the other spectra present in Fig. 5, namely for , the maximum amplitude occurs at a frequency of Hz; whereas for the maximum occurs at Hz. As expected, the higher the eccentricity is, the higher the frequency corresponding to the maximum amplitude of the spectrum. Moreover, the amplitudes have maximum values of and for and , respectively.

A relevant issue is to investigate the detectability of the backgrounds studied here. Given the frequency bands involved, we consider the detectability by pulsar timing array (PTA) techniques and by the planned (space-based) antenna eLISA.

We refer the reader to, for example, Moore et al. (2015), where the sensitivity curves of several detectors, including the ones we mentioned here, are provided.

Concerning the PTA technique, the frequency band of operation ranges from to Hz, with maximum sensitivity Hz. Particularly, in terms of the energy density parameter, EPTA (IPTA) has a maximum sensitivity of (), whereas SKA is planned to have maximum sensitivity of .

The planned eLISA (and LISA, its precursor) has maximum sensitivity, in terms of the energy density parameter, of () around Hz.

Given the sensitivity curves of PTA and eLISA (LISA), we conclude that the backgrounds shown in Fig. 5 would not be detected.

## 7 Conclusions

In this paper we show how to obtain the stochastic background of GWs generated by long-lived DNSs in eccentric orbits. We calculate the spectra by means of a different method: we determine the spectra that would be generated by each harmonic and then we sum these contributions in order to obtain the total backgrounds. In addition, we study the role of the eccentricity on the resulting spectra. This is achieved by using different distribution functions for the eccentricity. Particularly, we use two distribution functions: the Dirac delta function given by equation (16) and the uniform distribution given by equation (17).

From Fig. 5, we note that the distribution of eccentricities affect the amplitudes and the frequency bands of the backgrounds. For example, we note that using the distribution given by equation (17), the amplitudes decrease in the frequency band ranging from to Hz, relatively to the case of a population of systems in circular orbits (see Fig. 5). Concerning the frequency bands of the backgrounds, they are enlarged and shifted towards higher frequencies.

In general, we note that the higher the eccentricity of the population is the higher the frequency corresponding to the maximum amplitude of the spectrum. This behaviour is caused by the relative contribution given by the harmonics: for high eccentricities, higher harmonics affect more significantly the backgrounds.

Concerning the detectability of the spectra studied here, they would not be detected by the presently proposed space-based antennas, since the amplitudes involved are very small. Similar conclusion holds if we consider detections using the PTA technique.

In forthcoming papers we will address the study to systems formed by two black holes and those formed by a black hole and a neutron star. In addition, we will consider the high-frequency regimes, in which we cannot neglect the effects of the time evolution of the orbital parameters. Such cases are worth studying because, among other things, the spectra generated could form foregrounds for the planned space based antennas eLISA, DECIGO and BBO.

In addition, we expect that in the cases in which the time evolution of the parameters is considered, the form of the distribution will affect the results significantly, since in this case the frequency and the eccentricity will be related to each other.

## Acknowledgements

EFDE would like to thank the Brazilian agencies Capes, CNPq and FAPESP for financial support. JCNA would like to thank FAPESP and CNPq for the partial financial support. Last but not least, we thank the referee for the careful reading of the paper, the criticisms, and the very useful suggestions that greatly improved our paper.

### Footnotes

1. footnotemark:
2. pagerange: Stochastic Background of Gravitational Waves Generated by Eccentric Neutron Star BinariesReferences
3. pubyear: 2015

### References

1. Carroll B.W., Ostlie D.A., 2007, An Introduction to Modern Astrophysics, 2nd edn. Addison-Wesley, San Francisco, CA
2. Douglass D.H., Braginsky V.B., 1979, in Hawking S.W., Israel W., eds, General Relativity: An Einstein Centenary Survey. Cambridge Univ. Press, Cambridge, p. 99
3. Enoki M., Nagashima M., 2007, Prog. Theor. Phys., 117, 241
4. Ferrari V., Matarrese S. and Schneider R., 1999a, MNRAS, 303, 247
5. Ferrari V., Matarrese S. and Schneider R., 1999b, MNRAS, 303, 258
6. Ferrari V., Matarrese S., Schneider R. and Portegies Zwart S.F., 2001, MNRAS, 324, 797
7. Heger A., Fryer C. L., Woosley S. E., Langer N. and Hartmann D. H., 2003, ApJ, 591, 288
8. Ignatiev V. B., Kuranov A. G., Postnov K. A., Prokhorov M. E., 2001, MNRAS, 327, 531
9. Kowalska I., Bulik T., Belczynski K., 2012, A&A, 541, A120
10. Moore C. J., Cole R. H., Berry C. P. L., 2015, Class. Quantum Gravity, 32, 015014
11. Peters P.C., 1964, Phys. Rev., 136, B1224
12. Peters P.C., Mathews J., 1963, Phys. Rev., 131, 435
13. Regimbau T., 2011, Res. Astron. Astrophys., 11, 369
14. Regimbau T. and Chauvineau B., 2007, Class. Quantum Gravity, 24, S627
15. Regimbau T., de Freitas Pacheco J.A., 2006, ApJ, 642, 455
16. Rosado P.A., 2011, Phys. Rev. D, 84, 084004
17. Salpeter E.E., 1955, ApJ, 121, 161
18. Smartt S.J., 2009, ARA&A, 47, 63
19. Springel V. and Hernquist L., 2003, MNRAS, 339, 312
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters