Single field inflation in the light of Pulsar Timing Array Data: Quintessential interpretation of blue tilted tensor spectrum through Non-Bunch Davies initial condition

In this work, we present a quintessential interpretation of having a blue-tilted tensor power spectrum for canonical single-field slow-roll inflation to explain the recently observed Pulsar Timing Array (NANOGrav 15-year and EPTA) signal of Gravitational Waves (GW). We formulate the complete semi-classical description of cosmological perturbation theory in terms of scalar and tensor modes using the Non-Bunch Davies initial condition. We found that the existence of the blue tilt ($n_t$) within the favoured range $1.2<n_t<2.5$ can be explained in terms of a newly derived consistency relation. Further, we compute a new field excursion formula using the Non-Bunch Davies initial condition, that validates the requirement of Effective Field Theory in the sub-Planckian regime, $|\Delta\phi|\ll M_{\rm pl}$ for the allowed value of the tensor-to-scalar ratio, $r<0.06$ from CMB observations. In our study, we refer to this result as Anti Lyth bound as it violates the well-known Lyth bound originally derived for Bunch Davies initial condition. Further, we study the behaviour of the spectral density of GW and the associated abundance with the frequency, which shows that within the frequency domain $10^{-9}{\rm Hz}<f<10^{-7}{\rm Hz}$ the outcome obtained from our analysis is completely consistent with the Pulsar Timing Array (NANOGrav 15-year and EPTA) signal. Also, we found that the behaviour of GW spectra satisfies the CMB constraints at the low frequency, $f_*\sim 7.7\times 10^{-17}{\rm Hz}$ corresponding to the pivots scale wave number, $k_*\sim 0.05{\rm Mpc}^{-1}$. Finally, the sharp falling behaviour of the GW spectra within the frequency domain $10^{-7}{\rm Hz}<f<1{\rm Hz}$ validates our theory in the comparatively high-frequency regime as well.

Inflation is a theorized stage of quasi-de Sitter expansion in the very early Universe that was introduced to address the flatness, horizon, and entropy issues [60][61][62][63][64][65][66].During inflation, small quantum mechanical fluctuations generated from linearized tensor and scalar perturbations to the metric are stretched outside the Hubble horizon before reentering much later in the time scale.The inflationary paradigm  is still in excellent shape which is consistent with a large number of cosmological observations [97][98][99][100][101][102][103][104][105][106][107][108][109].Further tests of the inflationary paradigm are one of the prime scientific ingredients of future observational probes [110][111][112].However, the detection of the inflationary SGWB has yet to be made.While this signal has typically been desirable at very low frequencies, f ⩽ O(10 −15 Hz − 10 −16 Hz) [113], based on expectations within the standard framework of inflationary models, such a search need not be restricted to the previously mentioned low-frequency domain, as beyond the canonical single field models of inflation or maybe some interesting version with the canonical one can naturally predict rich phenomenological features at higher frequency domain, including the key frequency range O(10 −9 Hz − 10 −7 Hz) as suggested by very recently observed NANOGrav 15 signal of GW.
The main component and the most important aspects of this essay are summarised below in a brief, point-by-point fashion.This will be beneficial for general readers from a variety of viewpoints: ⋄ In this work we start our discussion with the gauge fixing issue of the Cosmological Perturbation Theory using this we will perform the perturbation of scalar and tensor modes dynamical solution of which we have computed from Mukhanov Sasaki equation along with general quantum initial condition, which we will identify as the Non-Bunch Davies initial condition [90,[114][115][116][117][118][119][120][121][122][123][124][125][126][127][128][129][130][131].The parameter space of such new physics is spanned by SO(1, 4) isometry group of de Sitter space and for this reason, one can expect that the Bunch Davies initial condition will be represented by a point in the larger parameter space of the Non-Bunch Davies initial condition.Further, we will compute the impacts of these perturbations in the inflationary observables, such as in the spectrum, spectral tilt, and tensor-to-scalar ratio in the presence of Non-Bunch Davies initial condition which we believe will show significant deviations compared to the result obtained for the observables with Bunch Davies initial condition.
⋄ From our analysis we will show how one can violate the old consistency relation, r = −8n t with n t < 0 derived in the presence of Bunch Davies initial condition.In our work, we will also point to the preferred range of tensor spectral tilt which is expected to show a blue feature in the tensor power spectrum in the presence of Non-Bunch Davies initial condition.See refs.[132][133][134] to know more about tensor spectral tilt and its physical implications in CMB experiments.
⋄ Further we will compute the new version of the field excursion formula for the inflaton field and will explicitly show in what amount significant deviation can be achieved from the good old Lyth bound derived in the same context from Bunch Davies initial condition.In this computation, our prime objective is to validate the EFT prescription for the canonical single field inflationary paradigm with the help of Non-Bunch Davies initial condition.Though we will do our analysis in a model-independent fashion i.e. we will not specifically talk about any particular structure of the inflationary potential to derive this bound, but we strongly believe that blue tilted tensor power spectrum and to maintain the theoretical constraint from the new consistency relation derived from Non-Bunch Davies initial condition the potential should have a certain structure which will satisfy certain physical properties.
⋄ Next, using the derived new consistency relation between tensor-to-scalar ratio and tensor spectral tilt we will study the behaviour of the Gravitational Waves (GW) abundance (Ω GW h 2 ) with frequency (f ) in the presence of Non-Bunch Davies initial condition.Our claim is that the GW abundance computed for Non-Bunch Davies initial condition can able to justify the various underlying new physical phenomena in the three frequency bands, very low 10 −20 Hz < f < 10 −17 Hz, low 10 −17 Hz < f < 10 −7 Hz and high 10 −7 Hz < f < 1Hz.We will show the status of the Non-Bunch Davies initial condition generated GW in the light of CMB Planck observation, recently observed NANOGrav 15 signal, LISA [135], BBO [136], DECIGO [137], CE [138], ET [139], HLVK and HLV(03) [140][141][142].
The outline of this paper is as follows: In II, we briefly discuss the canonical single-field model of inflation.In III, we discuss in detail the generation of scalar and tensor modes using cosmological perturbation theory.In IV, we move on towards a specific and detailed computation of the scalar power spectrum and later discuss its observational impacts.A similar analysis is followed further in V for the derivation of the tensor power spectrum.In VI, the results derived for the power spectrum in the previous sections are connected with the recent findings of the Pulsar Timing Array (NANOGrav 15-year and EPTA) Data Set.In VII, a new consistency condition is established to explain the blue-tilted tensor power spectrum based on the observed changes in the amplitude and spectral tilt of the tensor perturbations as discussed in previous sections.In VIII, due to the presence of modified Bogoliubov coefficients for the scalar and tensor perturbations, the well-known Lyth bound is also modified and presented with a discussion performed in a model-independent fashion.In IX, we present the numerical analysis required to explain the impact of the blue-tilted power spectrum.After a brief discussion on the generation of stochastic gravitational waves we discuss plots for the tensor-to-scalar ratio based on our new consistency relation and for the abundance of gravitational waves in the presence of Non-Bunch Davies initial condition.Finally, in X we present the conclusions of our findings and discussions.

II. The canonical single field model of inflation: The old wine a new glass
To demonstrate the claim of this work let us first write down the fundamental action of the underlying theory which is described by the following equation: Canonical single field model minimally coupled to gravity: where ϕ is a scalar field that is minimally coupled to gravity and in the present context it describes the single field inflationary paradigm.In the above action the canonical kinetic term of the scalar field is described by the following shorthand notation: and V (ϕ) is the inflationary potential for the scalar field which satisfy the slow-roll conditions during inflation, M pl is the reduced Planck mass which is M pl ∼ 2.43 × 10 18 GeV, and last but not the least, the gravitational interaction is described by the Einstein-Hilbert term which is described in terms of R, which is the Ricci scalar in the following discussion.Most importantly, it is important to mention that all the terms in the above-mentioned action are minimally coupled with the gravitational sector through the term √ −g, which physically represents the determinant of the metric in the corresponding space-time.
In this connection, the best possibility is the spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric which describes our observable Universe, and the corresponding metric of space-time is given by: Qusi de Sitter space-time: where τ is the conformal time coordinate.The scale factor a(τ ) is expressed in terms of the conformal time coordinate and the corresponding expression describes the quasi-de Sitter solution which is extremely useful for describing the inflationary paradigm.Here 'H ′ represents the Hubble parameter in quasi de Sitter space-time, which is not exactly constant, and deviation from the constancy is treated in terms of slow-roll dynamics which we will discuss in brief now.To describe the slow-roll dynamics, let us first write down the classical dynamical equation of motions in spatially flat FLRW space-time, which are known as Friedmann equations and Klein-Gordon equation, and given by the following expressions: Classical field equations: where we have introduced a shorthand notation ′ which represents the conformal time derivative d/dτ .Also, we use the new definition of the Hubble parameter in conformal time, H ≡ a ′ /a = aH which is going to be extremely useful for further understanding the subject material studied in this paper.Now, to initiate and end properly the inflationary paradigm at a specific value in the field space one needs to introduce the following parameters in quasi-de Sitter space-time, known as the slow-roll parameters, Slow-roll parameters: To validate inflation in the slow-roll region, one needs to satisfy the following conditions, On the contrary, slow-roll conditions are violated at the end of inflation when we have either ϵ = 1 or |η| = 1 or both using which one can explicitly compute the field value at the end of inflation.

III. Cosmological Perturbation Theory: Generation of scalar and tensor modes
Let us now have a look at the small perturbations that occur around the spatially flat FLRW background, where the linearized field and the metric perturbations are provided by, Linearized field and the metric perturbations: In this description, we have the following important facts that one needs to remember always henceforth: (a) Here ϕ(τ ) represents the scalar field embedded in a spatially flat background FLRW space-time.
(b) Also, δϕ(x, τ ) quantifies small perturbations in the scalar field which can be able to capture the inhomogeneity and anisotropy.
(c) Here R(x, τ ) represents the scalar comoving curvature perturbation which is often used to describe the scalar fluctuations in inflation: The advantage of this new variable R(x, τ ) is that it is constant outside the horizon and invariant under linearized coordinate transformations.But when working directly in the language of R(x, τ ), taking the φ′ (τ ) → 0 limit can occasionally be misleading because φ′ (τ ) appears in the denominator on the RHS even if the exact de Sitter description should become a good description.
(d) Using two distinct gauges is the most straightforward solution to this problem.
(e) One can work in the gauge where the perturbations are while they are evolving inside the horizon and described by: Gauge A: We refer to this as gauge A. The scalar perturbation in this gauge is given by δϕ(x, τ ) and exhibits a smooth behaviour, which is described by: (f) At this stage one can consider another gauge which is the perfect choice while describing the perturbations outside the horizon and given by: Gauge B: The correlation functions expressed in terms of ζ(x, τ ) are time-independent because in this gauge the scalar perturbation is given by: which is constant outside the horizon.We refer to this as gauge B and in this article, we are going to use this specific gauge to perform the rest of the computation.
(g) There is an equivalent description using which sometimes one describes the Cosmological Perturbation Theory.This is known as Arnowitt-Deser-Misner (ADM) formalism and is extremely useful to describe the Hamiltonian formulation of the General Theory of Relativity, which is commonly used in the context of canonical quantum gravity and numerical relativity.In this description, the corresponding ADM metric is described by the following expression: Arnowitt-Deser-Misner (ADM) metric: where instead of using conformal time coordinate the formalism is completely developed in terms of the physical time coordinate t and for this reason, the perturbed metric is not expressed in terms of conformally flat form.
Using this metric one needs to choose the following gauge fixing condition: where N and N i represent the lapse and shift functions.In this context, the spatial component of the metric g ij is described using the following expression: Here R(x, t) and h ij (x, t) represent the comoving curvature perturbation and tensor perturbation and both of them are written in terms of the physical time coordinate t.Also, the scale factor is given by, a(t) = exp(Ht), where H is not exactly constant for quasi-de Sitter space-time.We will talk about the tensor perturbation in detail in the later points, so up to that point please allow us to explain some of the crucial facts regarding the gauge choices and formalism that one follows within the context of Cosmological Perturbation Theory.
(h) It is further important to note that the gauge fixing conditions written in terms of the lapse and shift functions, within the framework of ADM formalism, do not fix the gauge degrees of freedom completely.After performing the mentioned gauge fixing there is a possibility to have the following spatial reparametrization which is of the following form: and temporal reparametrization which is of the following form: Here ξ i (x) and ξ(x) represent the corresponding parameters of gauge transformations which describe spatial reparametrization and temporal reparametrization respectively.Out of both of these possibilities to write down the Cosmological Perturbation Theory in a correct fashion, one needs to only look into the temporal reparametrization which provides the correct coordinate transformation where one needs to choose the timeindependent parameter ξ(x) suitably.
(i) At the linearized level of the perturbation one can note the following relationship among the parameter ζ(x, τ ) of the gauge B and the parameter δϕ(x, τ ) corresponding to the gauge A, given by: Relationship between gauge A and gauge B: which becomes extremely useful for converting the cosmological observables, particularly the correlation functions and other derived parameters, from one gauge to another just by performing a very simple change of variables.
(j) Finally, in the present context h ij (x, τ ) represent the spin-2 tensor perturbations which describe the primordial gravitational waves (PGW) that satisfy the following properties at the linearized level of perturbation theory: Properties of PGW: In the upcoming section, we are going to study the impacts of scalar and tensor perturbations in detail and their direct connection with the presently observed NANOGrav 15-year Data Set.

IV. Computation of scalar power spectrum with Non-Bunch Davies initial condition
In this section, we are interested in studying the theoretical and observational impacts of the scalar power spectrum.To serve this purpose we begin by expanding the representative canonical scalar field model action up to the second order in the scalar comoving curvature perturbation, R(x, τ ) = ζ(x, τ ), which gives the following simplified form of the action: Second order perturbed action for scalar modes: For further computationally simplification purposes we introduce the rescaled space-time dependent field variable, which is written as: using which the second-order action can be further expressed in terms of the canonically normalized form and the corresponding expression is given by: where we introduce a conformal time-dependent quantity in the quasi-de Sitter space-time: Here the scale factor and slow-roll parameter ϵ are defined using the conformal time as mentioned in the previous part of this paper.Now instead of performing the rest of the analysis in coordinate space, we will perform the further computations in the momentum space as this is the most cosmologically relevant region from the observational perspective.Also, the coordinate space result involves the Infra Red (IR) cut-off in terms of the underlying length scale of the quasi-de Sitter theory, which is of course not good and one cannot able to fix such cut-offs just from the observational point of view.For this reason, we make use of the following Fourier transform ansatz using which the rescaled field variable in the coordinate space can be written in terms of its Fourier space counterpart by the expression: Using this Fourier transform ansatz the previously mentioned second-order action written in terms of the rescaled scalar perturbation can be recast in the following simplified form: This action mimics the role of the action of an oscillator having a conformal time-dependent effective frequency in the present discussion.Variation of the aforementioned action with respect to the modes of the rescaled field variable gives us following the equation of motion written as: Mukhanov Sasaki equation for scalar modes: which in the corresponding literature is commonly referred to as the Mukhanov Sasaki (MS) equation.For further computational purposes, it is important to note that in this computation we have the effective, conformal timedependent, mass which can be expressed in terms of the following simplified language: where we have parameterized the effective mass for the scale mode of perturbation in terms of the conformal timedependent mass parameter function for the scalar modes by the following relation at the leading order of the slow-roll approximation: Here the subscript s corresponds to the contribution computed for the scalar modes only and ϵ, η are the slow-roll parameters that we have explicitly defined in the earlier part of this paper.Using the above-mentioned parameterization and the corresponding slow-roll approximation the most general solution of the second-order MS equation can be written in terms of the following expression: In this solution, H νs (−kτ ) and H νs (−kτ ), represent the Hankel function of the first and second kind which can further be expressed in terms of the Bessel function of the first kind and second kind as: H (1)  νs (−kτ ) = J (1)  νs (−kτ ) + iY (1)  νs (−kτ ), H (2)  νs (−kτ ) = J (1)  νs (−kτ ) − iY (1)  νs (−kτ ), (35) which further implies that: k and β (s) k , are the two Bogoliubov coefficients which are fixed by the proper choice of the quantum initial condition during the slow-roll phase of the inflation.In a deeper sense, the correct choice of the initial condition can able to construct the corresponding quantum mechanical vacuum state.In principle, these Bogoliubov coefficients can be any arbitrary function of the momentum scale which will appear and span within the large classes of SO(1, 4) de-Sitter isometry group.The generator and the conserved charges of this infinite parameter family should contain the mentioned two Bogoliubov coefficients which describe the corresponding general quantum initial condition and the associated states can be fixed in this description very clearly.
In general, one can consider an arbitrary initial condition described in terms of a general quantum vacuum state which is specified by the two Bogoliubov coefficients α where Ĉk represents the corresponding annihilation operators for the scalar modes which become extremely important when we canonically quantize the perturbed modes using Cosmological Perturbation Theory.
In principle, the general quantum state can be expressed in terms of the well-known Euclidean vacuum state, which is commonly referred to as the Bunch Davies quantum vacuum state in the present description using a well-defined Bogoliubov transformation.Here before going to further computational details in this present section, it is important to mention that in the large parameter space of SO(1, 4) de-Sitter isometry group Bunch Davies initial condition is described by a single point, which is actually characterized by the fixing Bogoliubov coefficients as, α (s) k = 1 and β (s) k = 0. However if one is interested in utilizing the effect of a general quantum initial condition one should have an infinite number of possibilities where these coefficients are arbitrary functions of the characteristic momentum scale as appearing all over in the present computation.In this corresponding literature, this possibility is often referred to as Non-Bunch Davies initial condition which we explicitly show to be extremely useful for the present computational purpose.Now let us come back to the previously mentioned Bogoliubov transformation which helps us to express any general Non-Bunch Davies initial quantum states in terms of the Bunch Davies initial quantum state and such transformation is described by the following equation: Bogolibov transformation for scalar modes: BD , (38) where the overall normalization factor of the defined Non-Bunch Davies initial quantum states are described by the following equation: Additionally, it is important for the present description to note that, the initial Bunch Davies quantum state for the scalar modes is often referred to as: Here the careful observation clearly shows that we have omitted the superscript s because within the framework of describing Bunch Davies initial quantum state one cannot able to distinguish between the ground states which describe the generation of scalar and tensor modes i.e.

|1, 0⟩
(s) k,λ and β On the other hand, once we describe the Non-Bunch Davies initial quantum states in this context, where s stands for scalar and t for tensor perturbations, such tags cannot be removed as the corresponding states describing the scalar and tensor sectors are not at all same in general: where the quantum states |α k,λ ⟩ NBD and the corresponding Bogoliubov coefficients α k,λ as appearing in the context of tensor perturbations are explicitly defined in the next section.Here the index λ = +, × corresponds to the two polarizations of the tensor modes.
For further better understanding purpose let us also mention the following constraint condition on the Non-Bunch Davies initial quantum states which needs to be satisfied during performing rest of the computation: Similarly, the annihilation operators of the Non-Bunch Davies initial quantum states and Bunch Davies initial quantum state are related via the following sets of Bogoliubov transformations: Further, it is important to note that the signed coordinate transformation has been recently introduced in refs.[143][144][145] is an interesting and innovative transformation that changes the signature of the Jacobian if considered locally.It resembles the Bogoliubov transformations in a non-trivial curved space background, where it mixes the positive and negative frequency components.Here the non-trivial curved background is important as in the case of the Minkowski flat vacuum there is no mixing of the positive and negative frequency components appearing.As an immediate consequence of this local signed coordinate transformation the Fock space operators of the Non-Bunch Davies initial quantum state are written in terms of its Bunch Davies counter-part through Bogolibov transformations in terms of positive and negative frequency components of the creation and annihilation operators.Initially, the inclusion of the Non-Bunch Davies initial vacuum in the various theoretical frameworks was considered to incorporate the effect of interacting vacuum in different types of computations.Out of various possibilities, α-vacua which describes the squeezed vacuum states is very well studied within the framework of Quantum Field Theory of de Sitter space using which the key concepts of the primordial cosmology are established.However, α-vacua was introduced in a completely phenomenological fashion in the previous literature and strong physical ground was absent to explain the proper theoretical origin of such concepts to describe the Non-Bunch Davies initial vacuum.However, the underlying physical connection between the local signed coordinate transformation and the Bogolibov transformation clearly helps to justify the appearance of the Non-Bunch Davies initial vacuum in the presence of non-trivial curved quasi de Sitter background geometry.It would be interesting if a detailed formal and rigorous derivation of the statement is found and we would certainly keep it in our mind for our future investigations 2 .
It is also important to note that the scalar modes satisfy the following normalization condition: which is referred to as the Klein-Gordon product and can get translated further in terms of the following crucial constraint on the Bogoliubov coefficients describing the Non-Bunch Davies initial condition: Analyzing the derived solutions of the rescaled field and associated canonically conjugate momentum and deriving some physically significant information is challenging.We will therefore take into account the asymptotic limits in the resulting solutions, which will be very beneficial for our research at various cosmic scales.To determine how Hankel functions of the first and second kind of order ν s behave, we use the asymptotic limits −kτ → 0 and −kτ → ∞.After applying the said limits, we obtain the following expressions: Further, imposing the above-mentioned asymptotic limits applicable to the super-horizon (−kτ ≪ 1) and sub-horizon (−kτ ≫ 1) the expressions for rescaled field variable can be obtained by the following expressions: Utilizing this fact the most general solution for rescaled field variable for any general Non-Bunch Davies initial condition can be expressed as: Using the above-mentioned expression for the rescaled scalar modes finally the expression for the comoving curvature perturbation in gauge B is given by the following simplified expression: Solution for scalar modes in gauge B: which represents the most general solution of the MS equation for scalar modes in the presence of any general Non-Bunch Davies initial condition.Now we require the quantization of the comoving scalar curvature perturbation to compute the expression for the power spectrum.For this purpose, we use previously defined Non-Bunch Davies quantum vacuum state.The canonical quantization of the scalar mode and its associated conjugate momenta should satisfy the necessary equaltime commutation relation: using which one can get the following commutation relations between the creation and annihilation operators of the Non-Bunch Davies quantum vacuum state for the scalar modes: Then the two-point correlation function for the comoving curvature perturbation can be expressed as follows, and the associated dimensionless power spectrum can be computed as: Finally, for super-horizon scales (−kτ ≪ 1), the dimensionless power spectrum can be further simplified as: Scalar power spectrum in super-horizon scale: At horizon exit, −kτ = 1 condition is imposed in the above expression which is used for further estimation of the amplitude and connection with the constraints from observations.

V. Computation of tensor power spectrum with Non-Bunch Davies initial condition
In this section, our prime objective will be to study the theoretical and observational impacts of the tensor power spectrum.To serve this purpose, further expanding the representative Einstein-Hilbert classical gravitational action up to the second order gives the following simplified result: Second order perturbed action for tensor modes: Further, we consider the following rescaling in tensor mode perturbation field variable which will be useful for the further computational purpose: where, e (λ) ij (x) represents polarization tensor for the tensor modes for two helicities, λ = +, ×, respectively.Using the above-mentioned field redefinition the second-order perturbed action for the tensor modes can be written as: Now, making use of the following Fourier transform ansatz which helps to convert the above-mentioned second order action for the tensor modes in the momentum space: Using this Fourier transform ansatz the previously mentioned second-order action written in terms of the tensor modes of rescaled perturbation can be recast in the given simplified form: This action mimics the role of the action of an oscillator having a conformal time-dependent effective frequency in the present context of the discussion.Varying the action with respect to the rescaled field variable written in Fourier space the MS equation for the tensor mode can be written as: Mukhanov Sasaki equation for tensor modes: For the further computational purpose, it becomes important to remember that here we have obtained an effective, conformal time-dependent, mass which can further be expressed in the following simplified language: Here the subscript t corresponds to the contribution computed for the tensor modes.Utilizing these facts, the most general solution of the second-order MS equation can be written using the following expression: In this solution, H 3/2 (−kτ ) and H (2) 3/2 (−kτ ), represent the Hankel function of the first and second kind.Here, α k,× , are the two sets of Bogoliubov coefficients describing the two helicity dependent sectors of tensor modes which can be fixed by the proper choice of the quantum initial condition during the slow-roll phase of the inflation.In general, one can consider an arbitrary initial condition described in terms of a general quantum vacuum state which is characterized by the following equation: Definition of Non-Bunch Davies state for tensor modes: where Dk,λ represents the corresponding annihilation operators for the two helicity-dependent sectors of the tensor modes.
Further, it is important to mention that in the large parameter space of SO(1, 4) de Sitter isometry group Bunch Davies initial condition for the tensor modes is described by a single point for each helicity, which is characterized by the fixing Bogoliubov coefficients as, α k,× = (0, 0).However if one is interested to utilize the effect of a general quantum initial condition one should have an infinite number of possibilities where these coefficients are arbitrary functions of the characteristic momentum scale as appearing in both helicity-dependent sectors of the tensor modes and often referred to as Non-Bunch Davies initial condition.Here any general Non-Bunch Davies initial quantum states can be expressed in terms of the well-known Bunch Davies initial quantum state and such Bogoliubov transformation is described by the following equation: Bogoliubov transformation for tensor modes: where the overall normalization factor of the defined Non-Bunch Davies initial quantum states for the two helicitydependent tensor modes is described using the following equation: On top of this, it is important to note that in the present description, the Bunch Davies initial quantum state for the tensor modes is often referred to as, |1, 0⟩ BD = |0⟩.The similar constraints which we have discussed in the case of scalar modes holds good in the case of describing tensor modes with two helicities as well.
Further utilizing the asymptotic properties of the Hankel function in the sub-horizon and super-horizon region, as discussed for scalar modes in the previous section, the most general solution for the rescaled field variable for the tensor perturbation modes for any general Non-Bunch Davies initial condition can be expressed as: These tensor modes satisfy the following necessary normalization condition: which is referred to as the Klein-Gordon product and this further translates into the following crucial constraint on the Bogoliubov coefficients describing the Non-Bunch Davies initial condition: Using the above-mentioned expression for the rescaled field variable finally, the expression for the tensor mode for two helicities is given by the following combined format: Solution for tensor modes: which represents the most general solution of the MS equation for scalar modes in the presence of any general Non-Bunch Davies initial condition.Now we need to quantize the tensor perturbation to compute the expression for the power spectrum.For this purpose, we use previously defined Non-Bunch Davies quantum vacuum state.The canonical quantization between the tensor mode and its associated conjugate momenta has to satisfy the following equal-time commutation relation: using which one can get the following commutation relations between the creation and annihilation operators of the Non-Bunch Davies quantum vacuum state for the tensor modes: Then the two-point correlation function for the tensor perturbation can be expressed as follows, and the associated dimensionless power spectrum can be computed as: Finally, for super-horizon scales (−kτ ≪ 1), the dimensionless power spectrum can be further simplified as: Tensor power spectrum in super-horizon scale: At horizon exit, −kτ = 1 condition is imposed in the above expression which is used for further estimation of the amplitude and connection with the constraints from observations.

VI. Inflationary observables from inflation in presence Non-Bunch Davies initial condition
In this section, our prime objective is to connect the derived results for Non-Bunch Davies initial condition to describe the inflationary paradigm with the recent findings of the NANOGrav 15-year and the EPTA Data Set.Most importantly, in this section, we will provide a quintessential explanation of accommodating the blue-tilted tensor power spectrum within the context of canonical single-field slow-roll models of inflation in a model-independent fashion.We will write down the expressions for each observable and some other important quantities using Non-Bunch Davies initial condition in the super-horizon region which we believe will be extremely helpful to establish our key point in the rest half of this paper: (a) Scalar power spectrum: The previously derived expression for the dimensionless scalar power spectrum in the presence of Non-Bunch Davies initial condition, in the super-horizon region, can be written in terms of a simplified formula as: Scalar power spectrum in super-horizon scale: where at the pivot scale k = k * = 0.05Mpc −1 , the amplitude of the scalar power spectrum is given by the following expression: Scalar power spectrum at the pivot scale: (b) Scalar spectral tilt: The scalar spectral tilt from the aforementioned form of the scalar power spectrum is computed to be as: Scalar spectral tilt at the pivot scale: Here we fix the Non-Bunch Davies initial condition for the scalar counterpart of the Bogoliubov coefficients in such a way that the following constraint condition hast to be obeyed throughout our analysis: which means that the spectral tilt for scalar perturbation is not very sensitive to the choices of the respective Bogoliubov coefficients used to describe them.Here, first of all, it is important to note that both the slow-roll parameters ϵ and η are evaluated at the pivot scale and to avoid any cumbersome notation we will not be using the * notation for the these parameters.We must also mention the important fact that the following approximations are applicable in the slow-roll regime of inflation which helps us to write everything using the effective potential during inflation and its derivatives: Here the ′ is used to describe the derivative with respect to the inflaton field.
(c) Tensor power spectrum: The previously derived expression for the dimensionless tensor power spectrum in the presence of Non-Bunch Davies initial condition in the super-horizon region can be expressed in terms of the following simplified formula: Tensor power spectrum in super-horizon scale: where at the pivot scale k = k * = 0.05Mpc −1 , the amplitude of the tensor power spectrum is given by the following expression: Tensor power spectrum at the pivot scale: (d) Tensor spectral tilt: The tensor spectral tilt from the derived above-mentioned form of the tensor power spectrum is computed as: Tensor spectral tilt at the pivot scale: Here we fix the Non-Bunch Davies initial condition for the tensor counterpart of Bogoliubov coefficients in such a way that the following constraint condition hast to be obeyed throughout our analysis: which means that the spectral tilt for tensor perturbation is extremely sensitive to the choices of respective Bogoliubov coefficients used to describe them.Also, this implies that the tensor spectral tilt is highly blue tilted as during the slow-roll phase of inflation, the slow-roll parameter ϵ is very small for almost all classes of inflationary effective potentials.One needs to always satisfy the above-mentioned constraint on the specific functional form of the Bogoliubov coefficients which describe the tensor perturbation.Now, it is important to provide a warning to the general readers that the tensor spectral tilt described in equation ( 86) represents the first term in the following Taylor series expansion ln ∆ 2 t (k) in terms of ln k, which is given by: which can be further simplified as: In the above computation, we define: where the Bunch-Davies counterpart satisfies the slow-roll limiting condition given by the following expression: Now, if the quantity G is of order unity or larger, that would call for including higher terms in the abovementioned series expansion.The series can only be truncated iff its coefficients n t , α t , β t , etc are small, i.e. in the slow-roll limit for the Bunch Davis vacuum when it takes the values n BD d ln k 2 demand to include higher order terms and strictly pointing towards the violation of standard slow-roll consistency conditions, which is only valid in presence of Bunch Davies vacuum.In this article, we utilize this violation of slow-roll condition in the presence of Non-Bunch Davies vacuum to achieve our final desired result established at the end of the paper for gravitational waves.
Recently in the NANOGrav 15-year and the EPTA Data Set, it is reported that the tensor power spectrum is highly blue-tilted and described using the following suggested fitting formula: NANOGrav 15-year Data Set fitting formula for tensor spectral tilt: where r is the tensor-to-scalar ratio which we will compute from our present prescription.The above-suggested form clearly suggests the existence of the spectral blue tilt for the tensor spectrum.Further, this formula gives hints towards having a violation of the standard slow-roll condition, n t = −r/8, which is only true for Bunch Davies vacuum.In this article, through our established methodology we have suggested that inclusion of the Non-Bunch Davies vacuum is one of the most promising options using which one can violate the standard slowroll consistency condition.We will discuss this fact in a more detailed fashion in the next section where we will justify that our findings strongly support the recently observed result from the NANOGrav 15-year and EPTA Data Set.Additionally, we will provide a new consistency relationship in the next section which violates the standard slow-roll consistency condition in the presence of Non-Bunch Davies vacuum.
(e) Tensor-to-scalar ratio: Further, in the super-horizon region, the tensor-to-scalar ratio in the presence of Non-Bunch Davies initial condition can be expressed as: Tensor-to-scalar ratio in super-horizon scale: where at the pivot scale k = k * = 0.05Mpc −1 , the amplitude of the tensor-to-scalar ratio is given by the following expression: Tensor-to-scalar ratio at the pivot scale: where in presence of Non Bunch Davies initial condition the modification factor F is defined as: We will now mention the immediate important outcomes from the inflationary paradigm in the presence of Non-Bunch Davies initial condition.

VII. New consistency relation from Non-Bunch Davies initial condition
Due to having changes in the amplitude and spectral tilt of the tensor perturbation in the presence of Non-Bunch Davies initial condition, it is expected to have a violation of old consistency relation, r = −8n t , derived for the canonical single field slow-roll models of inflation in presence of Bunch Davies initial condition.Here we found that the new consistency relation between the amplitude and spectral tilt of the tensor perturbation is given by the following simplified formula: New consistency relation: where the modification factors F and G are already defined before in terms of the Bogoliubov coefficients of Non-Bunch Davies vacuum which describe the scalar and tensor perturbations respectively.In the next section, we will suggest the functional forms of both of the modification factors which can be able to explain the blue-titled tensor spectrum through Non-Bunch Davies initial condition in the light of PTA (NANOGrav 15-year and EPTA) Data Set.One important fact of the above relation is in the case of Bunch Davies initial condition where we have F = 1 and G = 0, which gives the old consistency relation and is only able to explain red-tilted tensor perturbation in the power spectrum.The PTA Data Set discarded this possibility and pointed towards having a blue-tilted tensor power spectrum which is almost impossible to explain via Bunch Davies initial condition.Though there exists a very small number of counterexamples using which one can produce a very small blue-titled tensor power spectrum, which is, a quintessential inflationary paradigm [93][94][95][96]146], phantom models of inflation [147,148], alternatives to inflation, such as string gas cosmology [149,150], ekpyrotic [151,152] and bouncing scenarios [153][154][155][156][157][158] and many more in the list [159][160][161][162][163].But the PTA Data Set is suggested to have large spectral blue tilt which cannot be explained by the above-mentioned examples without incorporating the details of reheating phenomenology.Recently in ref [107] the author showed that large blue tilt favours very low reheating temperature, for example, the numerical analysis suggests that to have n t ∼ 2 one should have reheating temperature within the range 10GeV < T reh < 1TeV, which is again almost impossible to generate from the standard known models of canonical single field inflation.However, we strongly believe once the microphysics of reheating will be known, which is till date completely unknown to all of us, some of the ambiguities and problems associated with reheating phenomenology related to the present work can be addressed in a more crystal clear fashion.Some of the efforts have been made in refs.[114,115,164,165] to address this mentioned issue, which we believe will be useful for the purpose of a better understanding of the general readers.In this work instead of talking about such possibilities, we have pointed out that the Non-Bunch Davies initial conditions are one of the best possibilities using which one can reliably explain the origin of the blue tilted tensor power spectrum as observed by the PTA (NANOGrav 15-year and EPTA) Data Set.Not only can Non-Bunch Davies initial conditions explain the blue-tilted tensor spectrum in the low-frequency domain, but also one can able to show using a bit of detailed numerical analysis that it is consistent with the constraints obtained from the high-frequency probes of the gravitation waves spectrum.In the next section, we are going to analyze all of these possibilities together in a detailed fashion.

VIII. Anti Lyth bound: New formula for field excursion using with Non-Bunch Davies initial condition
Since the Non-Bunch Davies initial condition is the key ingredient of this paper then due to having modifications in the presence of different Bogoliubov coefficients which explain the scalar and tensor perturbations correctly the well-known field excursion formula for the inflaton field derived in the presence of the Bunch Davies initial condition, commonly known as the Lyth bound, will also be going to modify.Though we perform our analysis in a completely model-independent fashion, i.e. without using any specific structure of inflationary effective potentials in the slow-roll regime, the result derived in this section will be going to further constrain the behaviour of such effective potentials in the presence of the Non-Bunch Davies initial condition.In this connection, it is important to point out that the Lyth bound is the determining factor that tells whether inflation occurs in the super-Planckian or in the sub-Planckian regime.The old version of the Lyth bound always suggests having a super-Planckian inflationary paradigm, which is against the general notion of the basics of the Effective Field Theory (EFT) prescription.However, enormous efforts have been made in the past to evade such strong outcomes of Lyth bound which renders the EFT prescription trustworthy within the framework of inflationary paradigm governed by Bunch Davies initial condition.See refs.[76,81,82,84,85,146] more details on this aspects.Here we are digging out this good old issue again in this paper again as we believe the Non-Bunch Davies initial condition is one of the promising possibilities using which one can able to evade the good old Lyth bound derived for Bunch Davies initial condition.This means that the existence of Non-Bunch Davies initial condition clearly validates the inflationary paradigm in a region where the basic understandings and general notions of EFT hold good perfectly.Not only validating the EFT framework but also our computation performed in this section is going to be extremely helpful to know how exactly the observed tensor blue tilt by the PTA (NANOGrav 15-year and EPTA) Data Set can be able to constrain the structure of the inflationary potential to validate EFT prescription in the present context of the discussion.Here we start our derivation by the following conversion formula which connects the underlying momentum scale with the inflaton field: which is also used during the computation of running of the scalar and tensor spectral tilt of the primordial spectrum as discussed before in the earlier part of this paper.Using the above-mentioned conversion formula one can further write down the expression for the first slow-roll parameter ϵ V in the following simplified form: Now through the use of the expression mentioned above, in the super-horizon region one can go on to write down the expression for the tensor-to-scalar ratio as a function of momentum scale in the presence of Non-Bunch Davies initial condition: using which one can further compute the following expression for the field excursion formula: where the field excursion in terms of the inflaton field is defined as: where the field value at the end of inflation and at the horizon crossing scale is described by, ϕ end and ϕ * .Now for the further simplification purpose we introduce the number of e-foldings ∆N which in the present context helps us to make the connection with momentum-dependent scale and is described by the following expression: Using the above-mentioned result, the new field excursion formula for the inflaton field can be recast in the following form: New field excursion formula: Anti Lyth bound For the blue-tilted tensor spectrum let us fix, n t ∼ 2 and for the red tilted scalar spectrum, we consider n s ∼ 0.96.Additionally, we fix the number of e-foldings at ∆N ∼ 60.We also use the relation r(k * ) = 16ϵ V F, which further gives the following simplified result in terms of order of magnitude: where in the last step during the simplification we use the constraint on the potential dependent first slow-roll parameter as, ϵ V ∼ 0.0044, which is consistent with the Planck data [166].This implies that blue tensor spectral tilt with Non-Bunch Davies initial condition describes sub-Plancking EFT description of the canonical single field inflation.It is then expected from the above-mentioned estimation that if we push the blue tilt of the tensor power spectrum to the lower values, say n t ∼ 0.22 then we get |∆ϕ| ∼ O(1)M pl where the EFT prescription breaks and the super-Planckian physics start dominating over the sub-Planckian behaviour of the inflationary potential.This means that at, n t ⩽ 0.22, EFT prescription breaks even if we consider a small blue tilt in the mentioned spectrum.This analysis also tells us that if we consider the red tilt in the tensor power spectrum as obtained from Bunch Davies initial condition, then the field excursion formula suggests that EFT description breaks and super-Planckian physics play a very significant role in the present computation.After doing this analysis we have found that the findings of Non-Bunch Davies initial condition are more appropriate than the result obtained from Bunch Davies initial condition as we always want that the EFT prescription holds good perfectly.

IX. Numerical results: Explanation of blue titled tensor spectrum through Non-Bunch Davies initial condition in the light of Pulsar Timing Array Data Set
In this section, our prime objective is to explain the impact of the blue titled tensor spectrum through the Non-Bunch Davies initial condition in the light of the Pulsar Timing Array (NANOGrav 15-year and EPTA) Data Set which we believe is one of the promising possibilities to defend the canonical single-field slow roll inflationary paradigm.To explain the prescribed scenario we emphasize the violation of the consistency conditions as derived in the earlier section.In the upcoming discussion, we suggest some frequency-dependent and tensor spectral tilt-dependent specific features of the Bogoliubov coefficients describing the tensor sector of the perturbation, which will be going to fix some of the unknown quantities as appearing in the expression for the derived new consistency relation in the presence of Non-Bunch Davies initial condition.Before going into further details on the related discussion, below we will first discuss some of the basic details about the generation of stochastic gravitational waves and related constraints using the Pulsar Timing Array (NANOGrav 15-year and EPTA) Data Set.
Let us start with the following conversion formula which helps us to convert the underlying cosmological momentum scale in terms of the GW frequency: where we must point that, 1nHz ≡ 10 −9 Hz.Using this conversion formula the amplitude of the tensor power spectrum can be recast in the following simplified formula: Using this relation in the region k ≫ k eq , the spectral energy density of GW can be expressed as [167]: Spectral density of GW from theory of inflation: where f * ∼ 7.7 × 10 −17 Hz represent the frequency of the pivot scale which is k * ∼ 0.05Mpc −1 .Here η 0 is the conformal time scale at the present day, H 0 is the present epoch Hubble parameter, k eq is the wave number associated with matter radiation equality, which is fixed by k eq ∼ 7.3 × 10 −2 × Ω m h 2 Mpc −1 ∼ O(10 −2 )Mpc −1 corresponding to the frequency scale, f eq ∼ 1.54 × 10 −17 Hz.Here Ω m is the density parameter for matter and h represents the reduced Hubble parameter in this discussion.From Planck 2018 [168] one can take the best-fit values of the cosmological parameters for the further analysis purpose, which are Ω m = 0.315, h = 0.673 for η 0 ∼ 1.4 × 10 4 Mpc, ∆ 2 ζ ∼ 2.2 × 10 −9 , where the pivot scale and matter-radiation equality scale are given by, k * = 0.05Mpc −1 and k eq ∼ 0.01Mpc −1 .Now we will talk about the GW produced by PTA where the spectral density is given by the following expression: Spectral density of GW from PTA portal: where h(f ) is the frequency-dependent GW power spectrum, which appears at the characteristic frequency of PTA, which is f PTA ∼ 10 −9 Hz = 1nHz.For the practical purpose in the mentioned characteristic frequency domain one can use power law parametrization using which the GW power spectrum can be expressed by the following formula: where A G is the amplitude of the PTA-based GW and the corresponding spectral tilt is α, which is related to the spectral tilt γ of the PTA residual cross-power spectral density by, α = 1 2 (3 − γ).Using this specific power law form of the spectrum the spectral density for GW from PTA can be further expressed by the following simplified relationship: Now at the pivot scale k * = 0.05Mpc −1 one can explicitly compute the amplitude of the GW amplitude from PTA, which is given by the following expression: Constraint on GW amplitude from PTA portal: where we use n t = 5 − γ in this computation.This relation suggests that for the allowed upper bound of the tensor-to-scalar ratio r ≲ 0.06, if one can provide a large blue tilt from a fundamental theoretical framework then only with PTA portal in the nHz frequency domain one can able to explain the production of GW amplitude.As we have pointed already out in the previous part of the discussion in this paper, Non-Bunch Davies initial condition is one of the interesting and promising possibilities using which one can able to explain the generation of a large blue tilted tensor spectrum.For this reason, it is highly possible to detect the imprints of such theoretical framework in the GW spectrum in the light of observed by the Pulsar Timing Array (NANOGrav 15-year and EPTA) Data.
In the context of Non-Bunch Davies initial condition, we parameterize the Bogoliubov coefficients associated with scalar and tensor perturbations in such a fashion that the following new consistency relation is maintained in all frequency ranges: where we suggest that the factor F is very small and approximately a constant and G is a function of the tensor spectral tilt n t .We suggest the following functional forms from the numerical best-fitting perspective when considering all  ( In figure (1), we have plotted the behaviour of tensor-to-scalar ratio (r) vs tensor spectral tilt (n t ) plot using new consistency relation for Non-Bunch Davies initial condition.We have also compared our result with the newly proposed fitting formula as suggested in the PTA (NANOGrav 15-year and EPTA) result.This representative plot clearly suggests that our proposed new consistency relation conforms well with the PTA result.Now keeping the above parametrization in mind the spectral energy density of GW can be further written in the following simplified form: Sepctral density of GW amplitude for Non-Bunch Davies initial condition: The suggested feature of the spectral energy density of GW in the three frequency domain for Non-Bunch Davies initial condition with n t > 0 is actually motivated by the scenario of analyzing slow-roll (SRI), ultra slow-roll (USR) and then a slow-roll (SRII) phase one after another to study the generation of large amplitude fluctuation in presence of instantaneous transition necessarily required to produce PBHs.See refs.[169][170][171][172][173][174][175][176][177][178][179][180][181][182][183][184][185][186][187][188] for more details on this issue.However, through deeper analysis, we have understood in ref [172][173][174] that quantum loop effects dominate in between the frequency domain 10 −10 Hz < f < 10 −9 Hz corresponding to the momentum scale, 10 5 Mpc −1 < k < 10 6 Mpc −1 and one cannot able to produce a sufficient number of e-foldings (we have found N ∼ 25) to produce M PBH ∼ M ⊙ , where M ⊙ ∼ 10 31 kg.Renormalization and DRG resummation of the one-loop scalar power spectrum put such strong constraints on the final result.For this reason, one can only able to produce a small mass PBHs (M PBH ∼ 10 2 gm) at very high frequency 10 6 Hz < f < 10 7 Hz corresponding to the momentum scale, 10 21 Mpc −1 < k < 10 22 Mpc −1 .In this case, a sufficient number of e-foldings can be achieved i.e.N ∼ 60 even after applying strong constraints from renormalization and DRG resummation in the one-loop corrected result of the scalar power spectrum.This problem can be evaded in the case of Galileon inflation which we have explicitly shown in ref [174,175,184,186] due to having a strong non-renormalization theorem.In this article, our motivation is not PBHs production but rather violating the well-known slow-roll consistency condition r = −8n t , with red tilt n < 0 derived for canonical single field slow-roll models of inflation with Bunch Davies initial condition.Though Bunch Davies initial condition seems like a natural choice as it asymptotically goes to the Minkowski vacuum, but not capable enough to resolve all the above-mentioned crucial issues.Also, the corresponding red-tilted spectrum could not able to show enhancement to confront the findings of the recently pointed Pulsar Timing Array (NANOGrav 15-year and EPTA) result.Though it is not fully very naturally originated but Non-Bunch Davies initial condition helps us to violate the consistency condition.Instead of inserting a slow-roll (SRI), ultra slow-roll (USR), and then a slow-roll (SRII) phase one equivalently performs the analysis by violating the consistency condition in the presence of Non-Bunch Davies initial condition.Previous analysis also points towards the fact that the transition from SRI to USR and USR to SRII shifts the corresponding vacuum state to a Non-Bunch Davies one if we initially start with the Bunch Davies condition at the starting point of SRI phase.Such a newly created vacua in USR and SRII phases not only deviates from its Bunch Davies counterpart but also makes the corresponding Bogoliubov coefficients of the states momentum and time scale-dependent.In the above-mentioned parametrization of the new consistency relation, all of this information is not only maintained but also helps us to go beyond our general notion, which was established in the ground of very specific choice of Bunch Davies initial condition.
In figure (2) we have depicted the behaviour of the Gravitational Waves (GWs) abundance (Ω GW h 2 ) with frequency (f ) in the presence of Non-Bunch Davies initial condition.By following the new structure of the modification factors in the three frequency bands, very low 10 −20 Hz < f < 10 −17 Hz, low 10 −17 Hz < f < 10 −7 Hz and high 10 −7 Hz < f < 1Hz, the behaviour of the spectral density of the GW and consequently the GW abundance changes its features with respect to the tensor spectral tilt n t which in turn fixes the value of the tensor-to-scalar ratio by making use of the newly proposed consistent relations for Non-Bunch Davies initial condition.Throughout the analysis, we have maintained the positive signature of the tensor spectral tilt n t > 0. The suggested new consistency relation is applicable to the mentioned three frequency regions and the above-mentioned plot clearly shows that in the intermediate low region, 10 −17 Hz < f < 10 −7 Hz, the blue-tilted behaviour of the power spectrum is dominating which finally gives rise to the enhancement in the GW abundance and using such prescription we can easily observe that it is able to explain the signals obtained from the recent Pulsar Timing Array (NANOGrav 15-year and EPTA) results.On the other hand, in the very low-frequency domain, 10 −20 Hz < f < 10 −17 Hz, the spectrum is behaving like a red-tilted though we have maintained the signature of n t > 0. It means that in this region the frequency-dependent behaviour of the spectral density of GW and the corresponding GW abundance is such that the outcome gives a red-tilted spectrum.Careful observation tells us that this is the region where CMB observation plays a very crucial role and one has to respect the constraints from the Planck observation as well in this domain, particularly at the pivot scale f * ∼ 7.7 × 10 −17 Hz which corresponds to the wave number k * ∼ 0.05Mpc −1 .We have obtained the red-tilted behaviour of the tensor spectrum in this particular domain which is also consistent with the Planck observation.Further, it is important to note that in this particular region with the help of Bunch Davies initial condition one can In addition to this, it is important to note that all the benchmark values of the tensor-to-scalar ratio are consistent with the upper bound r < 0.06 as obtained from Planck observation.
able to explain the existence of the red tilted tensor power spectrum which supports CMB observation.That means, in our analysis, the mentioned low-frequency domain is the Bunch Davies limiting condition of the proposed Non-Bunch Davies initial condition where everything seems perfect and the outcomes conform to the Planck observation as well.This is a good news that the proposed Non-Bunch Davies condition captured the information of the good old Bunch Davies initial condition in the very low-frequency domain which validates our analysis from the perspective of the CMB observation as well.On top of that, in the mentioned intermediate low-frequency domain due to having the blue tilted feature in the spectrum it also satisfies the constraint obtained from Pulsar Timing Array (NANOGrav 15-year and EPTA) data with high accuracy.Now we discuss the comparatively high-frequency domain 10 −7 Hz < f < 1Hz, where the behaviour of the spectrum shows rapidly-falling behaviour compared to the behaviour observed from the very low-frequency domain.This particular feature points towards the fact that in this region the spectrum is different but shows a red-tilted feature, which is different from the feature obtained from the very low-frequency domain GW spectrum.This high-frequency domain behaviour of the spectrum is also consistent with other observational constraints obtained from LISA [135], BBO [136], DECIGO [137], CE [138], ET [139], HLVK and HLV(03) [140][141][142].These mentioned probes have not observed any signature from inflation and as the obtained spectrum from the Non-Bunch Davies initial condition falls in a rapid fashion and does not intersect the scanning parameter space of any of them, the high-frequency domain behaviour of the spectrum is also consistent.In figure (2) for the visualization purpose and to show the applicability of our result here we have taken five benchmark points, which are n t = 1.2 for r = 2.3 × 10 −5 , n t = 1.5 for r = 3.6 × 10 −8 , n t = 1.8 for r = 4.3 × 10 −11 , n t = 2.2 for r = 7.3 × 10 −15 , n t = 2.5 for r = 7.9 × 10 −18 .All of these values suggest that blue-tilted tensor spectra obtained from Non-Bunch Davies initial condition is more favourable than the red-tilted spectra as originally derived from the canonical single field model of inflation in the presence of Bunch Davies initial condition.Furthermore, it is important to note that all the benchmark values of the tensor-to-scalar ratio are consistent with the upper bound r < 0.06 as obtained from Planck observation.

X. Conclusion
We finally conclude our discussion with the following key highlighting points from our analysis: ✓ In this work we start our discussion with the gauge fixing issue of the Cosmological Perturbation Theory and mention clearly the usefulness of each of the choices in the context of primordial cosmology.
✓ Next, using a preferred gauge condition we have performed the perturbation of scalar and tensor modes dynamical solution of which we have computed from Mukhanov Sasaki equation along with general quantum initial condition, which we have identified as the Non-Bunch Davies initial condition in this paper.The parameter space of this new initial condition is spanned by SO(1, 4) isometry group of de Sitter space and for this reason one can expect that the Bunch Davies initial condition will represent a point in the large parameter space of the Non-Bunch Davies initial condition.
✓ Further we have computed impacts of these purterbations in the inflationary observables, such as in the spectrum, spectral tilt and tensor-to-scalar ratio in presence of Non-Bunch Davies initial condition.The obtained result show significant modifications/deviations compared to the result obtained for the observables with Bunch Davies initial condition.
✓ From our analysis we expect that the Bogoliubov coefficients representing the scalar and tensor sector of the cosmological perturbations are different.For this reason it will show immediately its impact in the consistency relation between the tensor-to-scalar ratio (r) and tensor spectral tilt (n t ).We have found from our analysis that the old consistency relation, r = −8n t with n t < 0 derived in presence of Bunch Davies initial condition is violated in the present context.We have explicitly derived the expressions for the new consistency relation in the presence of Non-Bunch Davies initial condition where the significant deviation is pointed in the derived result.
✓ Further with the help of newly derived results for inflationary observables derived from Non-Bunch Davies initial condition we have computed the new version of the field excursion formula for the inflaton field and show significant deviation from the good old Lyth bound derived in the same context from Bunch Davies initial condition.This newly derived result in our paper suggests that the field excursion in the red-tilted regime of the tensor spectrum is controlled by the tensor-to-scalar ratio and the total number of e-foldings necessarily required for inflation.On the other hand, the blue-tilted regime of the tensor power spectrum is controlled by the tensor-to-scalar ratio and the values of the tensor spectral tilt and scalar spectral tilt.We have found from our analysis that in the allowed region of tensor-to-scalar ratio, r < 0.06 from CMB observation one is able to achieve sub-Planckian field excursion which is the fundamental requirement to validate EFT prescription within the framework of inflation.Though during the derivation of this bound we have not used any specific structure of the inflationary potential but the present analysis suggests that to get blue tilted tensor power spectrum and to maintain the theoretical constraint from the new consistency relation derived from Non-Bunch Davies initial condition the potential should have a certain structure and physical properties which will respect both the newly derived consistency condition and field excursion formula.
✓ Next, we have suggested certain structure of the modification factors as appearing in the new consistency relation, which are made up of certain combinations of the Bogoliubov coefficients computed from scalar and tensor perturbations.This suggested form is incorporated so that it is able to respect the red tilted behaviour in the frequency ranges, 10 −20 Hz < f < 10 −17 Hz and 10 −7 Hz < f < 1Hz and the blue tilted behaviour in the domain 10 −17 Hz < f < 10 −7 Hz.Though we are currently not interested in PBH formation from the present setup but it is important to note that the present set up somewhat mimics the transient roll inflation where the scenario is describe in terms of SRI, USR and SRII phases and the corresponding transition from one phase to other phases.In near future we will investigate this issue in great detail and establish such underlying connection more clarity.
✓ Finally we have depicted the behaviour of the abundance of Gravitational Waves (GW) (Ω GW h 2 ) with frequency (f ) in the presence of Non-Bunch Davies initial condition.By following the new structure of the modification factors applicable in the three frequency bands, very low 10 −20 Hz < f < 10 −17 Hz, low 10 −17 Hz < f < 10 −7 Hz and high 10 −7 Hz < f < 1Hz the behaviour of the spectral density of the GW and consequently the GW abundance changes its features with respect to the tensor spectral tilt n t which in turn fixes the value of the tensor-to-scalar ratio by making use of the newly proposed consistent relations for Non-Bunch Davies initial condition.We have found that this GW spectra computed from Non-Bunch Davies initial condition is consistent with CMB observation within the frequency range, 10 −20 Hz < f < 10 −17 Hz, where the pivot scale is f * ∼ 7.7×10 −17 Hz.Most importantly, the blue tilted behaviour of this GW spectrum conforms well with Pulsar Timing Array (NANOGrav 15-year and EPTA) signal in the frequency range 10 −9 Hz < f < 10 −7 Hz, where the peak of the spectrum appears at the frequency scale, f ∼ 10 −7 Hz.In the comparatively high-frequency domain 10 −7 Hz < f < 1Hz behaviour of the spectrum is also consistent with other observational constraints obtained from LISA [135], BBO [136], DECIGO [137], CE [138], ET [139], HLVK and HLV(03) [140][141][142].These mentioned probes did not have observed any signature from inflation and as the obtained spectrum from the Non-Bunch Davies initial condition falls in a rapid fashion in this regime and does not intersect the scanning parameter space of any of them, the high-frequency domain behaviour of the spectrum is also consistent.Our analysis suggests blue tilted tensor spectra obtained from Non-Bunch Davies initial condition is more favoured than the red tilted spectra as originally derived from the canonical single field model of inflation in the presence of Bunch Davies initial condition.Finally, we have found that the values of the tensor-to-scalar ratio are consistent with the upper bound r < 0.06 as obtained from Planck observation.

k
as mentioned earlier.Within this description of the Quantum Field Theory of de-Sitter space, a general quantum vacuum state is described by these two numbers as |α (s) k , β (s) k ⟩ NBD and characterized by the following equation: Definition of Non-Buch Davies state for scalar modes:

t , α BD t and β BD t 3 .
The nonnegligible values of the additional terms G, dG d ln k and d 2 G

FIG. 1 .for 10 − 20
FIG.1.Tensor-to-scalar ratio (r) vs tensor spectral tilt (nt) plot using new consistency relation for Non-Bunch Davies initial condition.We have also compared our result with the newly proposed fitting formula as suggested in NANOGrav 15 year result.

16 FIG. 2 .
FIG.2.Representative plot showing the behaviour of the Gravitational Waves (GWs) abundance (Ω GW h 2 ) with frequency (f ) in the presence of Non-Bunch Davies initial condition.This plot clearly shows that the blue tilted behaviour of the spectrum in the frequency range 10 −17 Hz < f < 10 −7 Hz is perfectly consistent with the NANOGrav 15 signal (in black bands) and the EPTA signal (in cyan bands).This plot is also consistent with the CMB observation (Planck) in the very low-frequency range, at the pivot value f * ∼ 7.7 × 10 −17 Hz corresponding to momentum scale k = 0.05Mpc −1 .For the visualization purpose and to show the applicability of our result here we have taken five benchmark points, which are nt = 1.2 for r = 2.3 × 10 −5 , nt = 1.5 for r = 3.6 × 10 −8 , nt = 1.8 for r = 4.3 × 10 −11 , nt = 2.2 for r = 7.3 × 10 −15 , nt = 2.5 for r = 7.9 × 10 −18 .All of these values suggest that blue-tilted tensor spectra obtained from Non-Bunch Davies initial condition is more favoured than the red-tilted spectra as originally derived from the canonical single field model of inflation in the presence of Bunch Davies initial condition.In addition to this, it is important to note that all the benchmark values of the tensor-to-scalar ratio are consistent with the upper bound r < 0.06 as obtained from Planck observation.