Is pulsar timing a hopeful tool for detection of relic gravitational waves by using GW150914 data?

The inflation stage has a behaviour as power law expansion like S(η)∝η1+β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S(\eta )\propto \eta ^{1+\beta }$$\end{document} where β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} constrained on the 1+β<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1+\beta <0$$\end{document}. If the inflation were preceded by a radiation era, then there would be thermal spectrum of relic gravitational waves at the time of inflation. Based on this idea we find new upper bound on β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} by comparison the thermal spectrum with strain sensitivity of single pulsar timing. Also we show that sensitivity curve of single pulsar timing may be hopeful tool for detection of the spectrum in usual and thermal case by using the GW150914 data.


Introduction
The inflation stage has a behaviour as power law expansion like S(η) ∝ η 1+β . The S and η are scale factor and conformal time respectively and β constrained on the 1 + β < 0 1 [1,2]. If the inflation were preceded by a radiation era, then there would be thermal spectrum of gravitational waves at the time of inflation [3][4][5]. Based on this idea we have shown in [6] that there is some more chances for the detection of relic gravitational waves (RGWs) in usual and thermal spectrum by using the data of GW150914. On the other hand, there is an important procedure for direct detection of RGWs that called pulsar timing (PT). This procedure is based on the investigation of fluctuation of the pulses of the pulsars due to RGWs. The frequency range of PT is 10 −9 -10 −6 Hz [7]. The authors in [7] have found an upper bound on β and shown that there is so low chance for the detection of RGWs by using the mentioned method. But we think that the upper bound on β will modify if we use the thermal spectrum. Also the chance of detection of RGWs based on the strain sensitivity of PT will increase by using the obtained β for GW150914 [6] in usual and thermal case. Therefore the main purpose 1 Note: the upper bound on β is β ≤ −1.8 [1]. a e-mail: ba.ghayour@gmail.com of this work is investigation of this detection. In the present work, we use the unit c =h = k B = 1.

The spectrum of gravitational waves in usual and thermal case
The perturbed metric for a homogeneous isotropic flat Friedmann-Robertson-Walker universe can be written as where δ i j is the Kronecker delta symbol and h i j are metric perturbations with the transverse-traceless properties i.e; ∇ i h i j = 0, δ i j h i j = 0. The gravitational waves are described with the linearized field equation given by The tensor perturbations have two independent physical degrees of freedom (h + , h × ) that called polarization modes. We can write h + and h × in terms of the creation (a † ) and annihilation (a) operators then, where k is the comoving wave number with k = |k|, l pl = √ G is the Planck's length and p = +, × are polarization modes. The polarization tensors ε for each k and p. For a fixed wave number k and a fixed polarization state p the Eq. (2) gives coupled Klein-Gordon equation as follows: [1,2] and prime means derivative with respect to the conformal time. Since the polarization states are same, we consider f k (η) without the polarization index. The solution of the above equation for the different stages of the universe are given in [6]. There is another state that called 'thermal vacuum state' see [8] for more details.

The analysis of the spectrum
The expression of the noise spectral density of the relative frequency fluctuations for single pulsar timing is [7] where the ν is frequency. We consider the spectrum of RGWs in usual and thermal case in the range (10 −9 − 10 −6 ) Hz corresponding to the frequency range of PT. The spectrum of RGWs is related to characteristic strain spectrum like h mc (ν) ≡ h m (ν, η 0 )/ √ 2 [18,19]. We are interested to h mc (ν) √ F/ν for SNR = 1, where we consider the factor F = 1/3 in this work see [7] for more details. We plot h mc (ν) √ F/ν of RGWs compared to the strain sensitivity of PT based on Eqs. (6,9) in Fig. 1, respectively. We can see the obtained upper bound on β = −1.955 is modified compared to the β = −1.9 that obtained in [7] corresponds to the thermal spectrum for the lower frequency. Also we plot the spectrum h mc (ν) √ F/ν of RGWs for the obtained β based on the GW150914 data [6] compared to the strain sensitivity of PT in Fig. 2. The obtained results tell us that there is more Fig. 2 The spectrum h mc (ν) √ F/ν of RGWs in usual case (solid lines) and thermal case (dashed lines) for two obtained β based on GW150914 data [6] compared to the strain sensitivity of PT (green line) Therefore based on the our results, the PT may be more hopeful tool for the detection of RGWs.

Discussion and conclusion
The thermal spectrum of RGWs causes an enhanced amplitude that called 'modified amplitude'. The obtained results based on the thermal spectrum cause new upper bound on β which can give us more information about the inflation. Also using the obtained amounts of β based on the GW150914 data tell us that PT may be more hopeful tool for detection of RGWs in usual and thermal case.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendix A: extra-dimensional scenario and thermal gravitons
Cosmology with extra dimensions has been motivated by Kaluza and Klein (KK) [20][21][22]. The modern scenarios involving extra dimensions are being explored in particle physics, with most models possessing either a large volume or a large curvature. Although there exist different models of extra dimensions, there are some general features and signals common to all of them. In the presence of D extra spatial dimensions, the (3 + D + 1)-dimensional action for gravity can be written as and m pl is Planck mass, g is the four-dimensional metric, G N is Newton's constant, g D , G N and R D denote the higher dimensional counterparts of the metric, Newton's constant, and the Ricci scalar, respectively. The factor m D is the fundamental scale of the extra-dimensional theory.
Since the gravitational interactions are not strong enough to produce a thermal gravitons at temperatures below the Planck scale (m pl ∼ 1.22 × 10 19 GeV), the standard inflationary cosmology predicted the existence of the cosmic gravitational wave backgrounds which are non-thermal in nature. However if the universe contains extra dimensions that can generate the thermal gravitational waves, then its shape and amplitude of the stochastic cosmic background of gravitational waves (CGWB) may change significantly. This can happen when energies in the universe are higher than the fundamental scale m D . The gravitational coupling strength increases significantly as the gravitational field spreads out into the full spatial volume. Instead of freezing out at O(m pl ), as in 3+1 dimensions, gravitational interactions freeze-out at ∼ O(m D ). If the gravitational interactions become strong at an energy scale below the reheat temperature (m D < T R H ), gravitons get the opportunity to thermalize, creating a thermal CGWB. The creation of a thermal CGWB if (m D < T R H ), is unchanged by the type of extra dimensions chosen [9].
Thus, if extra dimensions do exist and the fundamental scale of those dimensions is below the reheat temperature, a relic thermal CGWB ought to exist today. Compared to the relic thermal photon background (CMB), a thermal CGWB would have the same shape, statistics and high degree of isotropy and homogeneity. The energy density (ρ g ) and fractional energy density (Ω g ) of a thermal CGWB are where ρ c is the critical energy density today, T C M B is the present temperature of the CMB and g * is the number of relativistic degrees of freedom at the scale of m D . g * is dependent on the particle content of the universe, i.e. whether (and at what scale) the universe is supersymmetric, has a KK tower, etc. Other quantities, such as the temperature (T ), peak frequency (ν), number density (n), and entropy density (s) of the thermal CGWB can be derived from the CMB if g * is known, as These quantities are not dependent on the number of extra dimensions, as the large discrepancy in size between the three large spatial dimensions and the D extra dimensions suppresses those corrections by at least a factor of ∼ 10 −29 . If m D is just barely above the scale of the standard model, then g * = 106.75. The thermal CGWB then has a temperature of T g = 1.52 × 10 25 Mpc −1 0.905 K and a peak frequency of 19 GHz [9].