Next Article in Journal
Feature-Based Classification Neural Network for Kepler Light Curves from Quarter 1
Previous Article in Journal
Space–Time Structure of Particle Emission and Femtoscopy Scales in Ultrarelativistic Heavy-Ion Collisions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Effects in Ising Cosmology

1
Department of Physics, School of Applied Mathematical and Physical Sciences, National Technical University of Athens, Zografou Campus, GR-15780 Athens, Greece
2
Institute of Theoretical Physics, Faculty of Physics, University of Warsaw, Pasteura 5, PL 02-093 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Universe 2023, 9(10), 434; https://doi.org/10.3390/universe9100434
Submission received: 27 August 2023 / Revised: 20 September 2023 / Accepted: 26 September 2023 / Published: 28 September 2023
(This article belongs to the Section Cosmology)

Abstract

:
We consider a real scalar field in de Sitter background and compute its thermal propagators. We propose that in a dS/CFT context, nontrivial thermal effects as seen by an ‘out’ observer can be encoded in the anomalous dimensions of the  d = 3  Ising model. One of these anomalous dimensions, the critical exponent  η , completely fixes a number of cosmological observables, which we compute.

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
S = d 4 x g 1 2 g μ ν μ ϕ ν ϕ 1 2 ( m 2 + ξ R ) ϕ 2 ,
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  R  is the scalar curvature.
Consider a  d + 1  dimensional FRW spacetime with the metric
d s 2 = a 2 d τ 2 d x 2 ,
with  τ  being the conformal time and  a ( τ )  the scale factor. The expanding Poincare patch of de Sitter space corresponds to  a = 1 H τ , where  H = a a 2  is the Hubble constant with  a = d a d τ . The expanding Poincare patch of dS space is parameterized by  τ ( , 0 ] . The scalar field mode in d-dimensional momentum space  ϕ | k | = χ | k | a  in this background yields the classical Klein–Gordon equation of motion (with  k = ( k 0 , k )  being the four-momentum and the prime being derivative with respect to  τ ):
χ k + ω | k | 2 χ k = 0 ,
where  ω | k | 2 = | k | 2 + m dS 2  and the time-dependent mass is provided by  m dS 2 = 1 τ 2 ( M 2 d 2 1 4 ) . The dS mass parameter is  M 2 = μ H 2 + 12 ξ μ H 2 = m 2 H 2 , and H is the inverse curvature parameter of dS space satisfying  R = 12 H 2 . The solutions to Equation (3) are linear combinations of the Hankel function  H ν cl ( τ , | k | )  and its complex conjugate of weight  ν cl  with
ν cl = d 2 1 4 M 2 d 2
and as we will revisit later on,  ν cl  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]:
H = 1 4 d 3 k Ω | k | 2 α k α k + δ ( 3 ) ( 0 ) + Λ | k | α k α k + Λ | k | * α k α k
where  α k , α k  are the annihilation and creation operators that satisfy the commutation relations
α k , α q = δ ( k q ) , α k , α q = α k , α q = 0 .
In addition,  Ω | k |  and  Λ | k |  are defined as
Ω | k | = | u | k | | 2 + ω | k | 2 | u | k | | 2 , Λ | k | = u | k | 2 + ω | k | 2 u | k | 2
where  u | k | , u | k | *  are the mode functions that pair with the ladder operators  α k , α | k | . The frequency  ω | k | 2  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  β k , β k  and mode functions that are connected by the Bogolyubov Transformation (BT) [2]:
β k = c | k | α k + d | k | * α k , β k = c | k | * α k + d | k | α k
with  c | k | , d | k |  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  n S . In particular, experiments [6,7] have found that  n S  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  n S  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  n S . 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  n S 1 . This leads to a parametric freedom which can be fixed by an RG flow argument that has its origins in the  d = 3  Ising model. In this way, it is possible to match the deviation of  n S  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  n S ( 1 )  and the non-Gaussianity parameter  f N L .

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  τ = 0 . These are empty vacua from the perspective of corresponding local (in conformal time) observers. The  | in  is chosen to be the maximally symmetric Bunch–Davies (BD) vacuum [13,14]. The two vacua are related via the BT  J | Φ I = I | Φ J , where  I , J = in , out  is a label of the vacuum and  Φ I  is the field operator with mode function  χ | k | I . 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  χ | k | in = u | k |  and  χ | k | out = v | k | .
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  D  in such a basis has a  2 × 2  matrix structure and is ( T  ( T * ), denoting time (anti-time) ordering, while  0 |  is a generic vacuum):
0 | Φ + ( τ 2 ) Φ ( τ 1 ) | 0 = D + ( τ 1 ; τ 2 ) 0 | Φ ( τ 1 ) Φ + ( τ 2 ) | 0 = D + ( τ 1 ; τ 2 )
and
0 | T [ Φ + ( τ 1 ) Φ + ( τ 2 ) ] | 0 = D + + ( τ 1 ; τ 2 ) 0 | T * [ Φ ( τ 1 ) Φ ( τ 2 ) ] | 0 = D ( τ 1 ; τ 2 )
where  D + ( τ 1 ; τ 2 ) = D + * ( τ 1 ; τ 2 ) D ( τ 1 ; τ 2 ) = D + + * ( τ 1 ; τ 2 )  and  D + ( τ 1 ; τ 2 )  =  χ | k | ( τ 1 ) χ | k | * ( τ 2 ) D + + ( τ 1 ; τ 2 ) = θ ( τ 1 τ 2 ) D + ( τ 1 ; τ 2 ) + θ ( τ 2 τ 1 ) D + ( τ 1 ; τ 2 ) . Note that because the vacuum state  | 0  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
D + + + D D + D + = 0 .
Hidden in these expressions is the  i ε  shift, implementing the projection on the vacuum at  τ = . It can be chosen such that in the flat limit the propagator becomes diagonal with  D + + = i k 2 m 2 + i ε . 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  C +  is the forward branch, where time evolution follows the path  τ in τ out , then  C  is the backward branch, where  τ out τ in ; we attach an extra part to the contour  C 3 , where  τ in τ in i β 2  and  β = 1 / T  is the inverse temperature parameter Figure 1:
Furthermore, we introduce the propagators
0 | T [ Φ 3 ( τ 1 ) Φ 3 ( τ 2 ) ] | 0 = D 33 ( τ 1 ; τ 2 ) 0 | Φ + ( τ 1 ) Φ 3 ( τ 2 ) | 0 = D 3 + ( τ 1 ; τ 2 ) 0 | Φ ( τ 1 ) Φ 3 ( τ 2 ) | 0 = D 3 ( τ 1 ; τ 2 )
where  Φ 3 ( τ )  is the field operator on  C 3  and  τ 1 , τ 2 C ; in addition, we demand that the following junction conditions for  a { + , , 3 } :
D a + ( τ 1 ; τ 2 ) | τ 2 = τ out = D a ( τ 1 ; τ 2 ) | τ 2 = τ out τ 2 D a + ( τ 1 ; τ 2 ) | τ 2 = τ out = τ 2 D a ( τ 1 ; τ 2 ) | τ 2 = τ out
must be satisfied at time instance  τ = τ out , where the  C +  and  C  contours meet. Moreover, the conditions
D a ( τ 1 ; τ 2 ) | τ 2 = τ in = D a 3 ( τ 1 ; τ 2 ) | τ 2 = τ in τ 2 D a ( τ 1 ; τ 2 ) | τ 2 = τ in = τ 2 D a 3 ( τ 1 ; τ 2 ) | τ 2 = τ in
need to be satisfied at  τ = τ in , where  C  and  C 3  meet. Finally, in order for the SK analogue of the KMS condition to hold, we need to sew together  C +  and  C 3 , which results in the conditions
D a + ( τ 1 ; τ 2 ) | τ 2 = τ in = D a 3 ( τ 1 ; τ 2 ) | τ 2 = τ in i β / 2 τ 2 D a + ( τ 1 ; τ 2 ) | τ 2 = τ in = τ 2 D a 3 ( τ 1 ; τ 2 ) | τ 2 = τ in i β / 2
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  C 3  leg of the contour can be taken as having the form of a plane wave. In addition, at  τ = τ i n  we assume the BD vacuum such that the mode functions are expressed in terms of the Hankel functions of  ν cl = 3 / 2  order. According to these assumptions, the solution to the conditions results in the in–in thermal propagator components [17]:
D + + β / 2 = D + + + n B ( β / 2 ) D + + + D D β / 2 = D + n B ( β / 2 ) D + + + D D + β / 2 = D + + n B ( β / 2 ) D + + + D D + β / 2 = D + + n B ( β / 2 ) D + + + D
with  n B  being the Bose–Einstein distribution parameter
n B ( β ) = e β ω | k | 1 e β ω | k | .
We can conveniently express this propagator collectively in matrix notation as
D β / 2 = D + s 2 ( β / 2 ) D + + + D + + * 1 1 1 1
with
D = D + + D + D + D ,
and the parametrization  s ( β / 2 ) sinh θ | k | ( β / 2 ) = n B ( β / 2 )  and  c ( β / 2 ) cosh θ | k | ( β / 2 ) . The flat limit of this propagator is diagonal, and its  + +  component is such that the  i ε  shift of the zero-temperature propagator denominator becomes  i E = i ε coth ( β ω | k | / 2 )  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  H  and  H ˜ . 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  D β = U β D U β T  and [18]
U β cosh θ | k | ( β ) sinh θ | k | ( β ) sinh θ | k | ( β ) cosh θ | k | ( β ) .
That this is an allowed operation on dS propagators is supported by the fact that a transformation by the matrix  U β  is a BT with coefficients  sinh θ | k | ( β ) = e β 2 ω | k | 1 e β ω | k |  and  cosh θ | k | ( β ) = 1 1 e β ω | k | . 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
D β = D + ( s 2 ( β ) + s ( β ) c ( β ) ) D + + + D + + * 1 1 1 1 .
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  sinh θ | k | cosh θ | k | . While this might seem troublesome at first, they both contain the same physical information. Taking advantage of the trivial identity
e β ω | k | 1 e β ω | k | + e β 2 ω | k | 1 e β ω | k | = e β 2 ω | k | 1 e β 2 ω | k | ,
the propagators in Equations (18) and (20) are seen to be equal for  β = β . Note that the above identity does hold in the  sinh θ | k |  and  cosh θ | k |  parameterization, where it reads  s 2 ( β ) + s ( β ) c ( β ) = s 2 ( β / 2 ) . 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  D  are correlators of the form
D J , γ I = J ; γ | T [ Φ I ( Φ I ) T ] | J ; γ .
The doublet field, now in the language of TFD, is  ( Φ I ) T = ( Φ I , Φ ˜ I ) , 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  | I ; β = U | I , where the eigenvalue of U is  U β ( θ )  and the Gibbons–Hawking (GH) effect [19] (for which we momentarily use the parameter  δ  to distinguish it from  β ) is expressed as  | I = | J ; δ  with  I J . The only temperature that dS space can sustain is the GH temperature, which means that  1 / β dS = T dS = H / 2 π = 1 / δ . As such, it is sufficient to know the form of the thermal dS–scalar propagator for some generic temperature and then set  β = β dS .

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  | τ | H = 1  and concentrating on horizon exiting modes specified by  | k τ | 1 , they determine various cosmological indices derived from the following thermal scalar power spectrum [12] (here,  1  is the  2 × 2  matrix with unit elements):
P S , β 1 = D β 1 | τ 1 = τ 2 ,
in terms of a single parameter (when the temperature takes its natural value  T = T dS ):
κ ω | k | | τ | | | k τ | = 1 = | k | 2 + m dS 2 | τ | 2 | | k τ | = 1 = 5 d 2 4 + M 2
where  ω | k |  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  d = 3  is the choice of  M = 0 , or  κ = i , which is known to generate a scale invariant CMB spectrum. This corresponds to  ν cl = 3 2  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  ν cl  and associated with the Gaussian fixed point of the  d = 3  real scalar theory, which flows towards an interacting IR fixed point and the corresponding  CFT  at  τ = 0 . It is clear that in the present context exact scale invariance is realized in the  | in  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  ν cl ν = ν cl + ν q , which can be interpreted as a shift in the scaling dimension of a dS scalar field
Δ = d 2 ν = d 2 ν cl ν q = Δ cl , ν q .
There is a corresponding shadow partner solution to this with  Δ + = d 2 + ν . In this paper, we are concerned with  ( Δ , Δ + ) cl = ( 0 , 3 ) .
In order to understand  ν q  (which turns out to be a non-trivial zero), we first point out that the  | out ; β  ( β > β dS ) state is a BT of the Bunch–Davies vacuum. The mode functions before and after the transformation solve the same Bessel equation with frequency  ω | k | . Upon a time-dependent BT however, the frequency that an observer sees for a time other than his own is [21]
Ω | k | = ω | k | ( | c | 2 + | s | 2 ) .
As a result, the horizon exit parameter is transformed as
κ Λ = κ 1 + 2 e 2 x κ 1 e 2 x κ = κ coth ( x κ ) ,
where we have defined the dimensionless temperature parameter  x = π H 2 π T , 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  ( μ H 2 + ξ R H 2 ) a 2 H 2 ϕ 2  in the Lagrangian with exit parameter  Λ 2 = | k τ | 2 + a 2 μ H 2 + ( ξ 1 6 ) R H 2 )  and that the late time equations of motion
ϕ + 2 a H ϕ + μ H 2 + ξ R H 2 a 2 H 2 ϕ = 0 , H = 1 2 a ϕ 2
have no non-trivial solution with  H = const .  and a mass term that is non-zero and finite.
The two limiting values of x are interesting. Its natural value  x = π , where  T = T dS , provides  Λ =  for  κ = i . This is a special case in which we recover a dS solution of maximal isometry that corresponds to  | out ; β dS . As in the BD vacuum, no modes are seen to exit the horizon, this time due to their ultra-short wavelength. In the limit  x , 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  Ω | k | , then their horizon exit parameter is forced to be  Λ 0 lim τ 0 ( Ω | k | τ ) 0  (in this limit, x becomes an odd multiple of  π / 2 ). This suggests the construction of a trajectory from  ( Λ , x ) ( , π )  to  ( 0 , )  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]:
n S , β = 1 + d ln | k | 3 P S , β d ln | k |
where  P S , β  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
δ n S n S , β 1 = 2 x Λ e x Λ 1 e 2 x Λ ,
of  n S  away from unity. Observe that in  | out ; β dS , where  x = π  and  Λ = δ n S  vanishes and we see a scale-invariant spectrum. Moving somewhat away from it, to  x π  (here, it is implicitly assumed that moving away from  T dS  is a result of spontaneous breaking of scale invariance, which is expected to lower the temperature), the state is  | out ; β  and  δ n S  becomes a one-parameter expression of  Λ . We can fix this freedom by determining the value  n S , β  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  x =  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  0.036 . Thus, near the horizon we have
n S 1 η = 0.964
while the experimentally measured value is equal to [6]
n S , exp = 0.9649 ± 0.0042 .
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  Δ + , cl  shifts is the operator anomalous dimension of the trace of the Ising stress energy tensor  Θ , which is an exact zero, realized as the cancellation  Γ Θ = γ O 2 γ σ  = 0, where  γ O  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  ν q  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  γ O  that ends up shifting the spectral index  n S . 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  n S  comes from the critical value  η  of  2 γ σ  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  δ n S = η  along which the system is heated up from zero temperature, where  Λ 0 = 0  and  x = , 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  x π , the value  Λ π 1.5117  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  x π , it is characterized by a non-zero  δ n S ; however, at exactly  x = π  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  T dS  is possible everywhere except at  x = , where the intrinsic temperature must become abruptly unobservable.
Because there are no free parameters, several other observables that are determined by  P S , β  are expected to be fixed. For example, we can define the moment
n S , β ( 1 ) = d n S , β d ln | k |
and compute it using  n S ( 1 ) = 0 . The result, evaluated under the same conditions as  n S , β , is
n S , β ( 1 ) = δ n S 2 1 Λ 2 x Λ 1 + 2 e 2 x Λ 1 e 2 x Λ ,
which, substituting  x π  and  Λ = Λ π 1.5117 , provides us with
n S , β ( 1 ) = 0.0186
for the running of the index. The constraints in [6] are
n S , exp ( 1 ) = 0.013 ± 0.012 .
Finally, the universal contribution to the non-Gaussianity parameter [23] can be expressed in terms of  N = t i t f d t H  and its derivatives in the in-vacuum as [24]
f N L un = 5 6 N ρ ρ N ρ 2 ,
with  N ρ = N ρ N ρ ρ = 2 N ρ 2  and  ρ P S , β . (The parameter  f N L  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
f N L un = 5 x ( 1 + Λ 2 ) 2 1 + x Λ cot ( x Λ 2 ) + 2 Λ 3 sinh ( x Λ ) 6 Λ 2 x ( 1 + Λ 2 ) + Λ sinh ( x Λ ) ,
which for  x π  and  Λ = Λ π 1.5117  provides us with
f N L un = 1.7138 ,
while one of the experimental results for  f N L un  in a certain analysis is [7]
f N L , exp un = 1.7 ± 5.2 .

5. Conclusions

In this paper, we have considered a thermal scalar in de Sitter background. Starting from the Bunch–Davies  | in  vacuum, the Bogolyubov transformation places us in the interior of the finite temperature phase diagram in a thermal state  | out ; β . This state can be connected through holography to the vicinity of an interacting IR fixed point in the universality class of the 3D Ising model. The system in this state is rather special, in the sense that the boundary operator that couples to the scalar curvature perturbations in the bulk has a classical scaling dimension. The critical exponent  η  is the order parameter for breaking the scale-invariant spectrum of curvature fluctuations, and a simple argument from the dS/CFT correspondence fixes the parametric freedom in the dS scalar theory, yielding the prediction  n S = 0.964 . In the same context, we compute a number of additional cosmological observables, such as the first moment of the scalar spectral index and the non-Gaussianity bispectrum parameter  f N L , then evaluate them numerically. Our predicted values of  n S n S , β ( 1 )  and  f N L  are well within current experimental bounds [6,7].

Author Contributions

Conceptualization, N.I., A.K. and F.K.; formal analysis, N.I., A.K. and F.K.; Writing—original draft, N.I., A.K. and F.K. All authors have read and agreed to the published version of the manuscript.

Funding

The research of F.K. leading to these results has received funding from the Norwegian Financial Mechanism for the years 2014–2021 under grant number DEC-2019/34/H/ST2/00707.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank I. Dalianis for discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this Appendix, we discuss the real time construction in the Hamiltonian formulation. First, we provide a shortcut derivation of the thermal propagator Equation (18) that starts from flat space and the definitions
D + β ( τ 1 , τ 2 ) = W 2 ( τ 1 , τ 2 ) + W 1 ( τ 1 , τ 2 ) D + β ( τ 1 , τ 2 ) = W 1 ( τ 2 , τ 1 ) + W 2 ( τ 2 , τ 1 )
with the Wightman functions defined as
W 1 ( τ 1 , τ 2 ) Tr { a ( τ 1 ) a ( τ 2 ) ρ } Tr { ρ } = n B e i ω ( τ 1 τ 2 ) W 2 ( τ 1 , τ 2 ) Tr { a ( τ 1 ) a ( τ 2 ) ρ } Tr { ρ } = ( 1 + n B ) e i ω ( τ 1 τ 2 )
where  ρ = e β H  is the thermal density matrix,  H  is the (harmonic oscillator) Hamiltonian, and the second equalities show the result of the trace computations. Now, because the time-dependent part of the mode function in flat space is  u ( τ ) = e i ω τ , we can write  e i ω ( τ 1 τ 2 ) = u ( τ 1 ) u * ( τ 2 )  and pass to dS space via the substitution  u ( τ ) χ | k | ( τ ) . Then, we have
D + β ( τ 1 , τ 2 ) = χ | k | * ( τ 1 ) χ | k | ( τ 2 ) + n B ( β ) χ | k | ( τ 1 ) χ | k | * ( τ 2 ) + χ | k | * ( τ 1 ) χ | k | ( τ 2 ) D + β ( τ 1 , τ 2 ) = χ | k | ( τ 1 ) χ | k | * ( τ 2 ) + n B ( β ) χ | k | ( τ 1 ) χ | k | * ( τ 2 ) + χ | k | * ( τ 1 ) χ | k | ( τ 2 )
and by imposing  D + + β ( τ 1 ; τ 2 ) = θ ( τ 1 τ 2 ) D + β ( τ 1 ; τ 2 ) + θ ( τ 2 τ 1 ) D + β ( τ 1 ; τ 2 ) D β ( τ 1 ; τ 2 ) = D + + β * ( τ 1 ; τ 2 )  and applying for  β / 2 , we arrive again at Equation (18).
A thermal propagator has to satisfy a variant of the KMS condition. The KMS condition originates from the definition
ϕ ( t 1 , x 1 ) ϕ ( t 2 , x 2 ) β = Tr ϕ ( t 1 , x 1 ) ϕ ( t 2 , x 2 ) ρ Tr { ρ }
which in principle leads to the thermally corrected dS propagator. However, the direct computation of the trace is not obvious for the time-dependent Hamiltonians.
The condition takes a simple form though near  τ , which can be shown explicitly. In Thermofield Dynamics, the form of the KMS condition depends on a gauge parameterized by a real number, say,  α . It is a well known fact that TFD propagators satisfy such a condition in any of these  α -gauges [25,26]. The condition holds due to the relation [18]
a k | 0 ; β = e α β ω | k | a ˜ k + | 0 ; β , 0 ; β | a k + = 0 ; β | a ˜ k e ( 1 α ) β ω | k |
between the standard annihilation operator acting on the vacuum of the Hilbert space  H  and the tilded creation operator acting on the vacuum of  H ˜ . The thermal vacuum  | 0 ; β  is defined by the action of a unitary operator on the tensor product of the vacuum states in  H × H ˜ . Using the above relations, it can be straightforwardly shown that the Wightman function between the fields  Φ , Φ ˜  for  α = 1 / 2  satisfies the following condition (with the mode functions of the scalar field near  τ  reducing to plane waves):
0 ; β | Φ ( τ 1 , x ) Φ ˜ ( τ 2 , y ) | 0 ; β = 0 ; β | Φ ( τ 2 + i β 2 , y ) Φ ˜ ( τ 1 i β 2 , x ) | 0 ; β ,
which is the KMS condition in the  α = 1 / 2  gauge. This is relevant for our case, as the transformation matrix Equation (19) is in this gauge [25]. A different gauge choice is to take  α = 1 , where
a k | 0 ; β = e β ω | k | a ˜ k + | 0 ; β , 0 ; β | a k + = 0 ; β | a ˜ k .
Then, the same Wightman function as above needs to satisfy
0 ; β | Φ ( τ 1 , x ) Φ ˜ ( τ 2 , y ) | 0 ; β = 0 ; β | Φ ( τ 2 , y ) Φ ˜ ( τ 1 i β , x ) | 0 ; β ,
that is, the KMS condition in the α = 1 gauge. Note that this is a relation in which the usual form of the KMS condition of thermal field theory can be recognised. These two gauges can, however, be readily seen to correspond to the equivalent Wightman functions by  τ 1 , 2 τ 1 , 2 i B 2 , as they can be related by a shift in the imaginary time. From this shift freedom, it can additionally be seen that the diagonal elements of the propagator do not satisfy any non-trivial constraint. In conclusion, to the extent that Equation (A4) applies to dS space and the trace is computable, the thermal propagator that it defines satisfies a KMS condition.

References

  1. Birrell, N.D.; Davies, P.C.W. Quantum Fields in Curved Space; Cambridge Monographs on Mathematical Physics; Cambridge University Press: Cambridge, UK, 1982. [Google Scholar]
  2. Mukhanov, V.F.; Winitzki, S. Introduction to Quantum Effects in Gravity; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  3. Akhmedov, E.T. Lecture notes on interacting quantum fields in de Sitter space. Int. J. Mod. Phys. D 2014, 23, 1430001. [Google Scholar] [CrossRef]
  4. Das, A. Finite Temperature Field Theory; World Scientific: Singapore, 1997. [Google Scholar]
  5. Irges, N.; Kalogirou, A.; Koutroulis, F. Ising Cosmology. Eur. Phys. J. C 2023, 83, 431. [Google Scholar] [CrossRef]
  6. Akrami, Y. et al. [Planck Collaboration] Planck 2018 results. X. Constraints on inflation. Astron. Astrophys. 2020, 641, 1807.06211. [Google Scholar]
  7. Akrami, Y. et al. [Planck Collaboration] Planck 2018 results. IX. Constraints on primordial non-Gaussianity. Astron. Astrophys. 2020, 641, 1905.05697. [Google Scholar]
  8. Maldacena, J.M. Non-Gaussian features of primordial fluctuations in single field inflationary models. J. High Energy Phys. 2003, 5, 13. [Google Scholar] [CrossRef]
  9. Larsen, F.; van der Schaar, J.P.; Leigh, R.G. De Sitter holography and the cosmic microwave background. J. High Energy Phys. 2002, 4, 047. [Google Scholar] [CrossRef]
  10. Larsen, F.; McNees, R. Inflation and de Sitter holography. J. High Energy Phys. 2003, 7, 051. [Google Scholar] [CrossRef]
  11. Chen, X.; Wang, Y.; Xianyu, Z. Schwinger-Keldysh Diagrammatics for Primordial Perturbations. J. Cosmol. Astropart. Phys. 2017, 12, 006. [Google Scholar] [CrossRef]
  12. Anderson, P.R.; Molina-Paris, C.; Mottola, E. Short distance and initial state effects in inflation: Stress tensor and decoherence. Phys. Rev. D 2005, 72, 043515. [Google Scholar] [CrossRef]
  13. Chernikov, N.A.; Tagirov, E.A. Quantum theory of scalar field in de Sitter space-time. Ann. L’Inst. Henri Poincaré Sect. A Phys. Thér. 1968, 9, 109–141. [Google Scholar]
  14. Allen, B. Vacuum States in de Sitter Space. Phys. Rev. D 1985, 32, 3136. [Google Scholar] [CrossRef]
  15. Kubo, R. Statistical-Mechanical Theory of Irreversible Processes. I. General Theory and Simple Applications to Magnetic and Conduction Problems. J. Phys. Soc. Jpn. 1957, 12, 570–586. [Google Scholar] [CrossRef]
  16. Martin, P.C.; Schwinger, J. Theory of Many-Particle Systems. I. Phys. Rev. 1959, 115, 1342–1373. [Google Scholar] [CrossRef]
  17. Semenoff, G.; Weiss, N. Feynman Rules for Finite Temperature Green’s Functions in an Expanding Universe. Phys. Rev. D 1985, 31, 689. [Google Scholar] [CrossRef] [PubMed]
  18. Takahashi, Y.; Umezawa, H. Thermo field dynamics. Int. J. Mod. Phys. B 1996, 10, 1755–1805. [Google Scholar] [CrossRef]
  19. Gibbons, G.W.; Hawking, S.W. Cosmological Event Horizons, Thermodynamics, and Particle Creation. Phys. Rev. D 1977, 15, 2738–2751. [Google Scholar] [CrossRef]
  20. Antoniadis, I.; Mazur, P.O.; Mottola, E. Conformal Invariance, Dark Energy, and CMB Non- Gaussianity. J. Cosmol. Astropart. Phys. 2012, 9, 024. [Google Scholar] [CrossRef]
  21. Garbrecht, B.; Prokopec, T.; Schmidt, M.G. Particle number in kinetic theory. Eur. Phys. J. C 2004, 38, 135–143. [Google Scholar] [CrossRef]
  22. Basset, B.; Tsujikawa, S.; Wands, D. Inflation dynamics and reheating. Rev. Mod. Phys. 2006, 78, 537–589. [Google Scholar] [CrossRef]
  23. Creminelli, P.; Zaldarriaga, M. Single field consistency relation for the 3-point function. J. Cosmol. Astropart. Phys. 2004, 10, 006. [Google Scholar] [CrossRef]
  24. Kehagias, A.; Riotto, A. The Four-point Correlator in Multifield Inflation, the Operator Product Expansion and the Symmetries of de Sitter. Nucl. Phys. B 2013, 868, 577–595. [Google Scholar] [CrossRef]
  25. Henning, P.A. Thermo Field Dynamics for Quantum Fields with Continuous Mass Spectrum. Phys. Reps. 1995, 253, 235–280. [Google Scholar] [CrossRef]
  26. Ojima, I. Gauge Fields at Finite Temperatures: Thermo Field Dynamics, KMS Condition and their Extension to Gauge Theories. Ann. Phys. 1981, 137, 1–32. [Google Scholar] [CrossRef]
Figure 1. The Schwinger-Keldysh (closed path) contour realized as the coexistence of a forward (+) and a backward (−) branch in conformal time evolution. The thermal generalization of this contour is constructed by adding an extra “thermal” leg to it or, in other words, by attaching the extra part C 3 at the contour such that τ in τ in i β 2 .
Figure 1. The Schwinger-Keldysh (closed path) contour realized as the coexistence of a forward (+) and a backward (−) branch in conformal time evolution. The thermal generalization of this contour is constructed by adding an extra “thermal” leg to it or, in other words, by attaching the extra part C 3 at the contour such that τ in τ in i β 2 .
Universe 09 00434 g001
Figure 2. (Left) A few points of the nearly conformal LCP defined by  δ n S = η  and (Right) the Bogolyubov transformation  | in | out ; β  and the LCP on the complex plane, where  κ = Λ + i Im κ .
Figure 2. (Left) A few points of the nearly conformal LCP defined by  δ n S = η  and (Right) the Bogolyubov transformation  | in | out ; β  and the LCP on the complex plane, where  κ = Λ + i Im κ .
Universe 09 00434 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Irges, N.; Kalogirou, A.; Koutroulis, F. Thermal Effects in Ising Cosmology. Universe 2023, 9, 434. https://doi.org/10.3390/universe9100434

AMA Style

Irges N, Kalogirou A, Koutroulis F. Thermal Effects in Ising Cosmology. Universe. 2023; 9(10):434. https://doi.org/10.3390/universe9100434

Chicago/Turabian Style

Irges, Nikos, Antonis Kalogirou, and Fotis Koutroulis. 2023. "Thermal Effects in Ising Cosmology" Universe 9, no. 10: 434. https://doi.org/10.3390/universe9100434

APA Style

Irges, N., Kalogirou, A., & Koutroulis, F. (2023). Thermal Effects in Ising Cosmology. Universe, 9(10), 434. https://doi.org/10.3390/universe9100434

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop