1. Introduction
The rapidly expanding phase of the universe can be modelled by de Sitter (dS) space, and the simplest form of matter by a real scalar. It is believed that basic effects that have left an imprint on the Cosmic Microwave Background (CMB) were of a thermal nature. Therefore, a real scalar field
in the expanding Poincare patch of dS space [
1,
2,
3] represents a simple model that can explain observed features of the CMB as formulated in the context of thermal quantum field theory [
4]. The action is
which we quantize while taking into account finite temperature effects. Here,
m is the mass of the scalar field,
is the non-minimal coupling of the field to gravity, and
is the scalar curvature.
Consider a
dimensional FRW spacetime with the metric
with
being the conformal time and
the scale factor. The expanding Poincare patch of de Sitter space corresponds to
, where
is the Hubble constant with
. The expanding Poincare patch of dS space is parameterized by
. The scalar field mode in d-dimensional momentum space
in this background yields the classical Klein–Gordon equation of motion (with
being the four-momentum and the prime being derivative with respect to
):
where
and the time-dependent mass is provided by
. The dS mass parameter is
,
, and
H is the inverse curvature parameter of dS space satisfying
. The solutions to Equation (
3) are linear combinations of the Hankel function
and its complex conjugate of weight
with
and as we will revisit later on,
is a part of the classical scaling dimensions of the bulk and boundary operators.
Focusing on a given vacuum at a specific time, we can use Taylor expansion in the usual way for the scalar field and calculate the Hamiltonian [
2]:
where
are the annihilation and creation operators that satisfy the commutation relations
In addition,
and
are defined as
where
are the mode functions that pair with the ladder operators
. The frequency
is defined below, and carries the time dependence. One of the main key points of a QFT in a curved spacetime is that there is not a single choice for a vacuum state, while different choices lead to different ladder operators
and mode functions that are connected by the Bogolyubov Transformation (BT) [
2]:
with
being the Bogolyubov coefficients. Using the commutation relations defined above, it can be shown in a straightforward manner that the Hamiltonian (
5) is hermitian for a given ground state defined at a given time
. An analysis concerning the full time range that spans the entire dS space, instead of only its expanding patch, as we consider here, is more delicate, and hermiticity may be lost.
The work presented in this paper follows the related work of [
5], where a possible connection between a bulk dS theory and a boundary Ising model was examined along with its implications for the value of the cosmological spectral index
. In particular, experiments [
6,
7] have found that
deviates slightly from unity, which indicates that the CMB is nearly scale-invariant. Although it is true that other endeavors to explain this deviation exist, e.g., [
8,
9,
10], the idea of using the critical exponent of a boundary Ising field in order to predict cosmological observables is new. Furthermore, it is important to note that the calculations and results in the current paper are model-independent, meaning that we only take for granted the experimental value of
and the assumption that the inflation era of the universe can be explained by the expanding Poincare patch of dS spacetime.
In
Section 2, we compute thermal propagators in dS spacetime by generalizing methods first produced in flat spacetime (as discussed in [
4]), i.e., the Schwinger–Keldysh (SK) path integral and the Thermofield Dynamics (TFD). To the best of our knowledge, while there are cases in which the SK path integral has been used in dS before, e.g., [
11], the use of the TFD formalism in a general FRW spacetime is novel. Moreover, a plain though non-trivial connection between the two formalisms is presented, which can hold for other spacetime choices other than dS, enabling the ambiguities rising in the usual SK construction to be avoided. In
Section 3, we make use of the dS thermal propagators in order to incorporate the thermal corrections that arise when proceeding to calculate the scalar spectral index
. In contrast to previous works [
12], where the spectral index to the leading order was found to be equal to unity, we argue here that thermal effects actually slightly break scale invariance, resulting in
. This leads to a parametric freedom which can be fixed by an RG flow argument that has its origins in the
Ising model. In this way, it is possible to match the deviation of
from unity to the experimental data. Finally, using the newly formed thermal dS propagator, it is possible to extract more cosmological observables determined by the scalar metric fluctuations, namely, the running of the spectral index
and the non-Gaussianity parameter
.
2. Propagators and Temperature
Quantization of this system results in the notion of a time-dependent vacuum state and a doubled Hilbert space. Regarding the vacua, we are concerned with the so called “in” vacuum defined at
and the “out” vacuum defined at the boundary (i.e., the horizon) of the expanding patch at
. These are empty vacua from the perspective of corresponding local (in conformal time) observers. The
is chosen to be the maximally symmetric Bunch–Davies (BD) vacuum [
13,
14]. The two vacua are related via the BT
, where
is a label of the vacuum and
is the field operator with mode function
. Note that the field is the same in both vacua, with the mode functions and the creation and annihilation operators inside it being the vacuum dependent quantities. Common notation is
and
.
The doubled Hilbert space can be understood in the context of the SK path integral as being related to a + (or forward) branch and a − (or backward) branch in conformal time evolution. The field propagator
in such a basis has a
matrix structure and is (
(
), denoting time (anti-time) ordering, while
is a generic vacuum):
and
where
,
and
=
,
. Note that because the vacuum state
is time-dependent, for now we construct the above formulas without choosing a specific vacuum. In addition, there is no need for a time ordered product in Equation (
9) because the two fields commute, as they are defined in different parts of the SK contour.
The above matrix elements satisfy the relation
Hidden in these expressions is the
shift, implementing the projection on the vacuum at
. It can be chosen such that in the flat limit the propagator becomes diagonal with
. The above construction of the propagator at zero temperature in dS spacetime has been recently studied in [
11].
The thermal generalization of the propagator components in Equations (
9) and (
10) is our next goal. If the Hamiltonian of the system were time-independent, it would be possible to simply follow the process described in
Appendix A and show that the propagator satisfies the KMS condition [
15,
16], which ensures that it is a good thermal propagator. Here, however, we are dealing with a time-dependent Hamiltonian, and this is not straightforward. Instead, we use the method introduced in [
17], which takes advantage of the SK contour by adding an extra “thermal” leg to it. In particular, if
is the forward branch, where time evolution follows the path
, then
is the backward branch, where
; we attach an extra part to the contour
, where
and
is the inverse temperature parameter
Figure 1:
Furthermore, we introduce the propagators
where
is the field operator on
and
; in addition, we demand that the following junction conditions for
:
must be satisfied at time instance
, where the
and
contours meet. Moreover, the conditions
need to be satisfied at
, where
and
meet. Finally, in order for the SK analogue of the KMS condition to hold, we need to sew together
and
, which results in the conditions
that ensure the consistency of the deformed contour and yield a good thermal propagator.
The above conditions introduce corrections of a thermal nature into the propagators in Equations (
9) and (
10), which we compute by making two assumptions. Because the chosen contour allows for an imaginary time flow, we assume that there is no inflation in that direction. This means that the mode functions on the
leg of the contour can be taken as having the form of a plane wave. In addition, at
we assume the BD vacuum such that the mode functions are expressed in terms of the Hankel functions of
order. According to these assumptions, the solution to the conditions results in the in–in thermal propagator components [
17]:
with
being the Bose–Einstein distribution parameter
We can conveniently express this propagator collectively in matrix notation as
with
and the parametrization
and
. The flat limit of this propagator is diagonal, and its
component is such that the
shift of the zero-temperature propagator denominator becomes
in the thermal state. It is easy to see that this thermal propagator satisfies a condition such as Equation (
11).
Here, we are actually interested in the out–out thermal propagator. We first derive the result using a novel shortcut, then show that it yields the correct result. The shortcut uses the TFD formalism, where the doubled Hilbert space is seen as the tensor product of the Hilbert spaces of positive and negative momenta, respectively
and
. The corresponding fields residing in these Hilbert spaces are
and
. The validity of this strategy is based on the fact that the SK structure can be read as a TFD structure, in which case the passage to finite temperature is via the transformation
and [
18]
That this is an allowed operation on dS propagators is supported by the fact that a transformation by the matrix
is a BT with coefficients
and
. Hence, we essentially calculate the thermal corrections that the BT has on the propagator via the TFD formalism. The result of the rotation provides us with the out–out thermal propagator
It is immediately noticeable that the two expressions in Equations (
18) and (
20) disagree in the thermal correction, as the latter has an extra term along
. While this might seem troublesome at first, they both contain the same physical information. Taking advantage of the trivial identity
the propagators in Equations (
18) and (
20) are seen to be equal for
. Note that the above identity does hold in the
and
parameterization, where it reads
. Therefore, we have proved that the known form of the dS thermal propagator of [
17] can be equivalently obtained via TFD rotation of the zero-temperature SK propagator of the half-thermal parameter. Of course, the equivalence of the two expressions reflects the universal nature of the dS temperature as measured at an arbitrary time instance by the in and out observers. The advantage of the TFD rotation operation is that it is very simple and can be easily generalized to any background. Thus, we use this point of view in what follows.
The result of all allowed thermal transformations of
are correlators of the form
The doublet field, now in the language of TFD, is
, while
is a thermal index associated with any combination of thermal transformations of the form Equation (
19). The label (not index)
I on the field is a reminder of the vacuum state to which the mode functions belong. The two types of thermal transformations that are relevant to us here are the insertion of an explicit density matrix, resulting in a transformation by a unitary operator
U, as
, where the eigenvalue of
U is
and the Gibbons–Hawking (GH) effect [
19] (for which we momentarily use the parameter
to distinguish it from
) is expressed as
with
. The only temperature that dS space can sustain is the GH temperature, which means that
. As such, it is sufficient to know the form of the thermal dS–scalar propagator for some generic temperature and then set
.
3. The Spectral Index with Thermal Effects
The propagators in Equations (
18) and (
20) determine several important observables. At equal space-time points and at the time of horizon exit, defined as
and concentrating on horizon exiting modes specified by
, they determine various cosmological indices derived from the following thermal scalar power spectrum [
12] (here,
is the
matrix with unit elements):
in terms of a single parameter (when the temperature takes its natural value
):
where
is the frequency defined under Equation (
3). This parameter can be traded for the weight of the Hankel function, as determined by the Klein–Gordon equation, in Equation (
4). Of special importance in
is the choice of
, or
, which is known to generate a scale invariant CMB spectrum. This corresponds to
and decaying modes at the time of exit.
A particularly useful point of view [
20] is to recognize the system at
as being related to a UV Conformal Field Theory (CFT) labeled by the weight
and associated with the Gaussian fixed point of the
real scalar theory, which flows towards an interacting IR fixed point and the corresponding
at
. It is clear that in the present context exact scale invariance is realized in the
vacuum, with the deviations generated by a spontaneous shift in
M that, according to Equation (
20), should have a finite temperature origin. Deviations can be encoded in general by a shift of the weight
, which can be interpreted as a shift in the scaling dimension of a dS scalar field
There is a corresponding shadow partner solution to this with
. In this paper, we are concerned with
.
In order to understand
(which turns out to be a non-trivial zero), we first point out that the
(
) state is a BT of the Bunch–Davies vacuum. The mode functions before and after the transformation solve the same Bessel equation with frequency
. Upon a time-dependent BT however, the frequency that an observer sees for a time other than his own is [
21]
As a result, the horizon exit parameter is transformed as
where we have defined the dimensionless temperature parameter
, which takes values in
. The transformed state in general has a reduced isometry with respect to the Bunch–Davies state. This can be seen by the fact that the BT introduces a non-zero mass term
in the Lagrangian with exit parameter
and that the late time equations of motion
have no non-trivial solution with
and a mass term that is non-zero and finite.
The two limiting values of
x are interesting. Its natural value
, where
, provides
for
. This is a special case in which we recover a dS solution of maximal isometry that corresponds to
. As in the BD vacuum, no modes are seen to exit the horizon, this time due to their ultra-short wavelength. In the limit
, on the other hand, the out observer sees modes of any wavelength as exiting modes, as in this limit the time of exit approaches the horizon. This means that if the observer calls their frequencies
, then their horizon exit parameter is forced to be
(in this limit,
x becomes an odd multiple of
). This suggests the construction of a trajectory from
to
along which the value of some yet to be defined thermal effect is kept non-zero and constant, starting from a position somewhat shifted away from the scale invariant limit
. Deviations from exact dS isometry due to finite temperature effects can be encoded in the shift of the spectral index of scalar curvature fluctuations [
22]:
where
is the thermal scalar power spectrum defined in Equation (
23).
In the previous section, we have shown that the SK and TFD formalisms result in equivalent propagators. Consequently, from Equation (
23), they both determine the same thermal deviation
of
away from unity. Observe that in
, where
and
,
vanishes and we see a scale-invariant spectrum. Moving somewhat away from it, to
(here, it is implicitly assumed that moving away from
is a result of spontaneous breaking of scale invariance, which is expected to lower the temperature), the state is
and
becomes a one-parameter expression of
. We can fix this freedom by determining the value
by interpreting its deviation from unity as an anomalous dimension in the dual field theory in the spirit of the dS/CFT correspondence. Then, we can reach
along a trajectory that keeps this value constant for all temperatures.
In [
5], it was proposed that within the dual field theory on the horizon, the anomalous dimension that shifts the spectral index is the critical exponent
, for which the non-perturbative value is around
. Thus, near the horizon we have
while the experimentally measured value is equal to [
6]
A known fact about the dual field theory of dS is that it is expected to be non-unitary [
8]. Thus, the argument could be made that the Ising model (or its large
N relatives), which is unitary, is not a good candidate as the complement to dS theory. The suggestion made here is that the complementary field theory to the AdS version of the model under discussion, which could be in the universality class of the Ising model (or its large
N relatives), is an analytic continuation of dS. This is supplemented by the fact that Equation (
31) is invariant under an analytic continuation; see [
5] for a more detailed justification.
This is a constraining statement that leaves us with no free parameters. In [
5], it was shown that the quantity by which
shifts is the operator anomalous dimension of the trace of the Ising stress energy tensor
, which is an exact zero, realized as the cancellation
= 0, where
is the “total” operator anomalous dimension and
is the field anomalous dimension (or the so-called wave function renormalization). Therefore, it is in this sense that
is a non-trivial zero, being related to the vanishing anomalous dimension of a special operator that is the energy–momentum tensor. In [
5], it was further demonstrated that it is the total anomalous dimension
that ends up shifting the spectral index
. Of course, outside of the fixed point where the Ising field is massive,
M deviates from zero in the bulk, and the solution to Equation (
28) is not dS. It is important, however, to understand that the main effect on
comes from the critical value
of
and that the deviation from the critical value is small as long as the system sits in the vicinity of the fixed point. For this reason, the leading order results are independent of the source of the breaking. In a sense, the only assumption here is that there is a mechanism of spontaneous breaking of scale invariance. From the point of view of the boundary, this could (for example) be justified as some sort of Coleman–Weinberg mechanism.
4. Line of Constant Physics and Other Observables
A line of constant physics (LCP) is a set of points on the phase space upon which the value of a physical quantity remains fixed. What we demonstrate now is that, in the bulk, there is an LCP labelled by the fixed value
along which the system is heated up from zero temperature, where
and
, up to the dS temperature. A few points on this line and a picture of the LCP can be seen in
Figure 2.
We stress that for a given x the corresponding value of is fixed by the label of the LCP. Thus, near the endpoint of the LCP, where , the value is a fixed output. It is important to emphasize that the LCP is really meaningful up to just outside of its two limiting points. Up to around , it is characterized by a non-zero ; however, at exactly this becomes equal to zero, as the trace of the boundary stress–energy tensor to which the bulk scalar couples vanishes. Analogously, the interpretation of each point on it as a dS space of the same is possible everywhere except at , where the intrinsic temperature must become abruptly unobservable.
Because there are no free parameters, several other observables that are determined by
are expected to be fixed. For example, we can define the moment
and compute it using
. The result, evaluated under the same conditions as
, is
which, substituting
and
, provides us with
for the running of the index. The constraints in [
6] are
Finally, the universal contribution to the non-Gaussianity parameter [
23] can be expressed in terms of
and its derivatives in the in-vacuum as [
24]
with
,
and
. (The parameter
is defined by a more general expression [
8]. The expression we use is provided in [
24] for the special case of a single scalar field.) It is computed as being
which for
and
provides us with
while one of the experimental results for
in a certain analysis is [
7]