Thermal representations in group field theory: squeezed vacua and quantum gravity condensates

We apply the formalism of thermofield dynamics to group field theory quantum gravity and construct thermal representations associated with generalised equilibrium Gibbs states using Bogoliubov transformations. The newly constructed class of thermal vacua are entangled, two-mode squeezed, thermofield double states. The corresponding finite temperature representations are inequivalent to the standard zero temperature one based on a degenerate vacuum. An interesting class of states, coherent thermal states, are defined and understood as thermal quantum gravity condensates.


I. INTRODUCTION
An increasing number of studies are hinting toward intimate links between entanglement and geometry. Spacetime is thought to be highly entangled, with quantum correlations in the underlying quantum gravitational system being crucial. Particularly in discrete quantum gravity approaches, a geometric phase of the universe is expected to emerge from a quantum, pre-geometric one via a phase transition. This geometric phase must then also be highly entangled. Squeezing techniques have been used often in such contexts to define suitably entangled non-perturbative vacua. Also, the emergent phase must encode a sufficient notion of semi-classicality. This notion is commonly introduced via coherent states. Moreover, fluctuations in relevant observables are expected to be important in the physics of a quantum spacetime and thus must be included in the description of the system. In this work we show how to construct quantum gravitational phases in group field theory (for now, kinematically) using thermofield dynamics, which display all these features, namely entanglement, coherence and statistical fluctuations in given observables. These phases can be of deep relevance in the study of semi-classical and effective continuum descriptions of quantum spacetime, especially in the context of cosmology and black holes.
Group field theories (GFT) [1][2][3] are statistical field theories of combinatorial and algebraic quanta of geometry [4][5][6], formally defined by a partition function They are strictly related to various other approaches like loop quantum gravity [7][8][9][10], spin foams [11,12], causal dynamical triangulations [13], tensor models [14] and lattice quantum gravity [15]. Like in usual field theories, the * mehdi.assanioussi@desy.de † isha.kotecha@aei.mpg.de kinematics is specified by a choice of the fields ϕ κ , each defined in general over a domain space of direct products of Lie groups, and taking values in some target vector space. The dynamics is specified by propagators and interaction vertices encoded in a function S, which can be understood as a Landau-Ginzburg free energy function in the present statistical context [6]; or as a Euclidean action from the point of view of standard quantum field theories [1][2][3]16]. However unlike in usual field theories, S is non-local in general with respect to the base manifold. This non-locality is essential and encodes the nontrivial combinatorial nature of the fundamental degrees of freedom and their dynamics. Moreover, the base manifold is not spacetime, but carries algebraic information associated with discrete geometric and matter degrees of freedom. Such a complete absence of any continuum spacetime structures a priori is a manifestation of background independence in group field theory, like in various other non-perturbative approaches to quantum gravity. The partition function Z GFT perturbatively generates Feynman diagrams that are labelled 2-complexes (dual to labelled stranded diagrams), with boundary states given by labelled graphs [1][2][3]16]. For the choice of models closer to loop quantum gravity and spin foam setups, the boundary states are abstract spin networks (but organised in a second quantised Hilbert space, that of a field theory [4]) and bulk processes are spin foams, both of which in turn are dual to polyhedral complexes when restricting to loopless combinatorics [17,18]. Thus, a group field theory generates discrete quantum spacetimes made of fundamental polyhedral quanta 1 .
We can describe the same structures from a many-body perspective [19], and treat as more fundamental an in-teracting system of many such quanta. This viewpoint enables us to import formal techniques from standard many-body physics for macroscopic systems, by treating a quantum polyhedron or an open spin network node as a single atom or particle of interest [5,6,20,21]. This further leads to a modelling of an extended region of discrete quantum space (a spin network) as a multi-particle state with a large number of such quanta, and a region of dynamical quantum spacetime (a spin foam) as an interaction process. This is the perspective that we employ here to successfully use techniques from thermofield dynamics, even when working with a radically different kind of system, one that is background independent, and devoid of any standard notion of space, time and other associated geometric structures and standard matter couplings.
One of its main advantages lies in the fact that its formulation parallels that of zero temperature quantum field theory, thus powerful tools of the temperature-independent setup can be translated to the thermal case, including perturbative Feynman diagrammatic techniques, symmetry breaking analyses, and for what concerns us here, Fock space techniques. TFD has been applied in various fields such as superconductivity, quantum optics, and string theory.
The core idea of the TFD formalism is to represent statistical ensemble averages as temperature-dependent vacuum expectation values. That is, given a system described by a Hilbert space h and an algebra of observables A h on h, one looks for a vector state |Ω ρ (thermal vacuum) in a Hilbert space H, corresponding to a density operator ρ on h, such that the following condition holds for all observables A of the system, where the subscript h or H denotes a suitable representation of a classical observable in the respective Hilbert spaces. In this work, we are mainly concerned with statistical equilibrium, thus with density operators of the Gibbs form e −βO h . Note that we will not make use of hats to denote operators in order to declutter the overall notation in the following. A vector state satisfying equation (1.1) can only be defined in an extended Hilbert space, by adding the socalled tilde degrees of freedom to the original ones [22]. Importantly, this doubling, or in general an enlargement of the space of the relevant degrees of freedom, is a characteristic feature of finite temperature description of physical systems. This was also discovered in algebraic quantum field theory for equilibrium statistical mechanics [27]. Moreover, equation (1.1) is strongly reminiscent of the construction of a GNS representation [28] induced by an algebraic statistical state, with the vector state |Ω ρ being the cyclic vacuum of a new thermal representation. These intuitions have in fact led to tangible relations between the two formalisms, with the tilde degrees of freedom of TFD being understood as those of the conjugate representation of Tomita-Takesaki theory in the algebraic framework [23,26,[29][30][31].
These structures are also encountered commonly in quantum information theory [32], which in turn is utilised heavily in various areas of modern theoretical physics, like holography. Specifically, constructing a state |Ω ρ is nothing but a purification of ρ. A prime example of a vector state is the thermofield double state, which is the purification of a Gibbs density operator. These states are used extensively in studies probing connections between geometry, entanglement, and more recently, complexity [33,34]. For instance in AdS/CFT, this state is important because a bulk with an eternal AdS black hole is dual to a thermofield double in the boundary quantum theory [35].
Our goal is to construct entangled, thermofield double states |Ω ρ that encode statistical fluctuations with respect to an observable of interest; and to subsequently introduce a family of coherent configurations of the quanta of geometry built upon the thermal vacuum. Then, such thermal quantum gravity condensates, say with respect to the spatial volume observable, can be particularly valuable in cosmology [36][37][38][39][40].
The article is organised as follows. We begin with an overview of the relevant essentials of the formalism of thermofield dynamics in section II, and use them in the rest of the paper to present a systematic extension of group field theories to construct finite temperature phases associated with generalised equilibrium Gibbs states. In section III we give a brief overview of the original setup of a bosonic group field theory coupled to a scalar matter field, which provides the basis for a zero temperature phase as described in section IV. In section V, we briefly recall the subtleties surrounding the foundational issue of defining statistical equilibrium in background independent systems to consider generalised Gibbs states, while focusing on extensive observable as generators. We construct their corresponding thermal vacua (thermofield double states), along with the inequivalent quantum gravity phases that they generate in section VI. We conclude with an outlook on future prospects for further applications in group field theory and other quantum gravity approaches in section VIII.

II. THERMOFIELD DYNAMICS
Below we review the basics of the TFD formalism [22][23][24][25][26] using a simple example of an oscillator. We will extend it to a field theory setup directly for GFTs in the subsequent sections.
Consider a single bosonic oscillator, described by ladder operators a, a † satisfying the commutation algebra, with the a-particle vacuum being specified by a |0 = 0. A Fock space h is generated by actions of polynomial functions of the ladder operators on |0 . Thermal effects are then encoded in density operators defined on h. In particular, an equilibrium state at inverse temperature β is a Gibbs state, where H is a Hamiltonian operator on h, possibly of grand-canonical type. But we must note that in the context of background independent systems with an absence of a notion of absolute energy, one may work with other observables like the spatial volume, as will be discussed briefly in section V [6,20]. This system is extended by including tilde degrees of freedom, with a Hilbert spaceh generated by a vacuum |0 and the ladder operatorsã,ã † satisfying also the bosonic algebra, All tilde and non-tilde degrees of freedom commute with each other, that is The Hilbert spaceh is conjugate to the original one via the action of the modular conjugation operator of KMS theory [29]; or equivalently via the tilde conjugation rules of TFD: , for all nontilde and tilde operators defined on h andh respectively, and z i ∈ C. The zero temperature (or its inverse, β = ∞) phase of the system is described by the enlarged Hilbert space, built from the Fock vacuum, by actions of the ladder operators a, a † ,ã,ã † , such that Temperature is introduced via thermal Bogoliubov transformations of the algebra generators, along with expressions for their adjoints b † andb † . The Bogoliubov transformations are canonical, thus leaving the algebra unchanged, so that the β-ladder operators also satisfy bosonic commutations relations, along with their adjoints. The temperature-dependent annihilation operators now specify a thermal vacuum, which is cyclic for the thermal Hilbert space H β . It is an example of a two-mode squeezed state, which can be restored to its full generality most directly by adding a net phase difference between the cosh and sinh terms in the transformations (2.9) and (2.10).
The Hilbert space H β can most directly be organised as a Fock space with respect to the β-dependent ladder operators that create and annihilate b-quanta over the thermal background |0 β . Like in any other Fock space construction, one can define useful classes of states in H β . We will return to this point in section VII where we define one such interesting class of states in the quantum gravity system, namely the coherent thermal states.
The thermal Bogoliubov transformations (2.9) are parametrized by θ(β), which must thus encode complete information about the corresponding statistical state. In the present case, it must be uniquely associated with the Gibbs state ρ β , which has a well known characteristic Bose number distribution, using a Hamiltonian of the form H = ωa † a, and the number operator N = a † a of the physical, non-tilde system. This is how θ is usually determined in TFD, using equation (1.1) for the number operator. Then, using inverse Bogoliubov transformations, the right hand side of equation (1.1) for the present case gives This specifies θ via the equation, The non-tilde degrees of freedom can be understood as being physically relevant in the sense that they describe the subsystem of interest which is accessible to the observer. In other words, it is the subsystem under study in any given situation. Then, the physically relevant observables belong to the algebra restricted to the non-tilde degrees of freedom, which can be retrieved from the full description by partially tracing away the complement, here the tilde degrees of freedom. Thus, one is usually interested in observable averages of the form 0 β | O(a, a † ) |0 β , for operators O that are in general polynomial functions of the generators of the physical non-tilde algebra. Notable geometric examples include relativistic quantum field vacua in Minkowski and Schwarzschild spacetimes, where the non-tilde algebra has support on spacetime regions exterior to the respective horizons. Here, the tilde subsystems are CPT conjugates (modulo a rotation) of the non-tilde ones, and belong to the interior of the horizons [41][42][43]. Another physical interpretation of the tilde subsystem that is more common in condensed matter theory, is that of a thermal reservoir [24,44].
In the present discrete quantum gravity context, for now we retain the elementary, quantum informationtheoretic interpretation of the non-tilde and tilde degrees of freedom, simply as describing a given subsystem and its complement respectively, without assigning any further geometric meanings.
The two sets of ladder operators are related to each other explicitly via equations (2.9) and (2.10). This suggests that their respective vacua are also related by an associated transformation. Indeed they are, via the following unitary transformation 2 (2.17) In terms of this thermal squeezing operator, we then have It is clear then that such unitary operators map the different β-representations into each other. In finite quantum systems this is simply a manifestation of von Neumann's uniqueness theorem. However, when extending the above setup to quantum field theory, one would expect that representations at different temperatures are inequivalent. This is certainly the case for physical systems in general. It is also true in the present quantum gravity case, as we show in section VI. Lastly, we note that in the field theory extension, equations (2.9)-(2.13) hold mode-wise and still remain well-defined. Together, they describe the β-phase of the system. However, the operator U (θ) is no longer well-defined in general (before any cut-offs). Thus, without any suitable regularization, equations (2.17)-(2.20) technically do not hold in full field theory.

III. BOSONIC GROUP FIELD THEORY
The most commonly studied models in GFT are scalar field theories defined over multiple copies of the local gauge group of gravity, which is SL(2, C) for 4-dim Lorentzian models and Spin(4) for Euclidean ones. The subgroup SU (2) is chosen often for boundary configurations related to loop quantum gravity. For models of pure geometry (that is, before any matter coupling), the group field ϕ is then commonly taken to be a complex-valued 2 The form of this unitary operator, along with equation (2.18), shows that the thermal vacuum |0 β is a two-mode squeezed state in which aã-pairs have condensed [25]. scalar field on SU (2) d , where the number of copies of the group denotes the dimension of the emergent spacetime. This choice of ϕ corresponds to taking an atom of space to be a geometric polyhedron with d faces (under an additional closure condition, see below), or equivalently, an open d-valent node with each of its incident half-link labelled with an SU (2) element.
In this work, we consider the scalar field of the theory to be defined over the manifold SU (2) d × R. From a technical point of view, this allows a more general treatment in the formalism, where the base manifold has a non-compact direction. From a physical point of view, this can be viewed as coupling gravity to a scalar matter field, which can be used to define a relational frame of reference [37,45], crucial for our subsequent application to thermal condensate cosmology [36]. Thus, the group field under consideration here is (3.1) In the following, we will use the notations g ≡ g i (i = 1, . . . , d) interchangeably as per convenience. The field ϕ is further chosen to be invariant 3 under a diagonal right action of SU (2) on SU (2) d , that is This is the geometric condition of closure, which in the corresponding quantised theory leads to the understanding of a quantum of this field as a quantised polyhedron, in turn dual to a gauge-invariant spin network node [17,46,47]. This invariance effectively reduces the geometric part of the domain space to SU (2) d /SU (2). In the following however, we will continue to use a redundant parametrisation for convenience and consider the gauge-invariant functions to be defined on full SU (2) d while explicitly satisfying equation (3.2). The corresponding quantum operator theory is based on a commutations relations algebra for bosonic 4 tetrahedra, with [ϕ, ϕ] = [ϕ † , ϕ † ] = 0, and I being a delta distribution on SU (2) d compatible with gauge invariance, given by We note that this quantum theory is in fact based on an abstract Weyl algebra for GFT [48]. Our constructions here to access an inequivalent thermal phase using TFD may in principle also be formulated more rigorously using modular structures of Tomita-Takesaki theory and their relation to TFD [29][30][31]. In the present work however, we are content with a more physical approach to the problem, directly in line with usual TFD studies.

IV. DEGENERATE VACUUM AND ZERO TEMPERATURE PHASE
The zero temperature phase of the system is based on an enlargement of the Fock representation of the above bosonic (Weyl) algebra, along the lines presented in section II but generalized here to a field theory.
The Hilbert space for a single gauge-invariant quantum is the state space of geometries of a quantum polyhedron with an additional real degree of freedom, where the quotient by SU (2) ensures gauge invariance.
In order to work with formally well-defined quantities, we smear the operator-valued distributions ϕ( g, φ) with a suitable basis of functions in H. 5 For gauge-invariant functions on SU (2) d , a useful basis is given by a set of gauge-invariant Wigner functions where j i ∈ N/2 labels irreducible representations of SU (2), m i , n i ∈ (−j i , ..., +j i ) are matrix indices in representation j i , D ji mini are complex-valued Wigner matrix coefficients (multiplied by a factor of √ 2j i + 1 for normalisation) in representation j i , and C j nι are intertwiner basis elements indexed by ι arising due to the closure condition in (3.2). Orthonormality and completeness are respectively given by, Similarly for the matter part, let us consider a basis of complex-valued smooth functions T α (φ) in L 2 (φ), labelled by a discrete index α, satisfying orthonormality 5 For more functional rigour, the test functions may be defined on the dense subspace of smooth functions C ∞ (.) ⊂ L 2 (.). But this technicality is often overlooked in practice, as is also done here, thus taking the L 2 space to be the space of single particle wave functions and the space of test functions. and completeness, We thus have a complete orthonormal basis on H consisting of functions f Jα of the tensor product form, Then a suitable set of mode ladder operators can be defined by smearing the operators ϕ, ϕ † with this basis, , .

(4.9)
This essentially decomposes the operators ϕ, ϕ † in terms of the modes f Jα , which can be seen directly by inverting the above two equations. The algebra relations are, and [a, a] = [a † , a † ] = 0. The vacuum is specified by, which is a degenerate state with no discrete geometric or matter data, dubbed often as a no-space state. It generates the symmetric Fock space sym H ⊗N (4.12) by cyclic action of the generators {a Jα , a † Jα , 1} of this representation 6 . For instance, a single particle (N = 1), single mode state is, while a generic single particle state with a wavefunction (4.14) The zero temperature phase is then given by extending the above with the conjugate representation spaceH F , as discussed in section II. This gives the zero temperature (β = ∞) description in terms of a Hilbert space, which is a Fock space on the cyclic vacuum with ladder operators {a, a † ,ã,ã † } J,α that satisfy, The vacuum satisfies The action of all polynomial functions of non-tilde and tilde ladder operators on |0 ∞ generates H ∞ , all in complete analogy with standard H F , including the construction of multi-particle states, coherent states, squeezed states, and so on.

V. GENERALISED GIBBS STATES
The familiar way to include thermal effects is with statistical states, as density operators in a given representation or in general as algebraic mixed states. As we discussed in section II, an equivalent way is with their corresponding vector states (thermal vacua) in an enlarged representation of the system, in which the additional degrees of freedom are known to be integral for encoding finite temperature effects.
Here we are interested in Gibbs density operators for describing equilibrium phases of the quantum gravity system. But defining statistical equilibrium in background independent settings, like the present one, is a subtle issue (see [6], and references therein, for detailed discussions). For the purposes of this work however, what concerns us is the so-called thermodynamical characterisation, based on the maximum entropy principle, for defining generalised Gibbs states of the form where β ℓ are generalized inverse temperatures conjugate to a given set of observables O ℓ [5,6,20]. This state is a result of maximising the information entropy, − ln ρ ρ , under the set of constraints O ℓ ρ = constant and 1 ρ = 1. Particularly, observables O need not be dynamical energies. This is especially valuable in background independent systems where the notion of energy is ambiguous at best, or not defined at all in more radical setups such as the current quantum gravitational one. In fact, these observables can be geometric operators such as area and volume, which could be particularly relevant in quantum gravity. For instance, in the context of applying these constructions to homogeneous cosmology, O could be chosen to be the spatial volume, thus giving a volume Gibbs state [20]. Such a state would then encode thermal fluctuations in spatial volume of the underlying discrete quantum space, which may be expected to be important in cosmological dynamics, especially at very early times [36].
Below we consider the more general case with selfadjoint and semi-bounded operators for group field theories coupled with a scalar matter, so that for such generalised Gibbs states, the above machinery of TFD can be applied to construct various different phases characterised by the corresponding thermofield doubles. Naturally, this leaves open the possibilities of different observables, and which ones are relevant in any given situation is also an important part of the broader problem of investigating the statistical mechanics of quantum gravity for an emergent, thermodynamical spacetime with features compatible with semi-classical studies [6,20].
Given a self-adjoint and semi-bounded O, and using the maximum entropy procedure recalled above, under the constraints O ρ = constant, N ρ = constant and 1 ρ = 1, a generalised Gibbs state can be defined as, where β ∈ R and µ ∈ R are chosen such that the combination β(O − µN ) is a positive operator, thus ensuring proper normalisation of the state. A particularly interesting class of observables where the partition function can be evaluated rather straightforwardly is for positive, extensive operators. That is, we can consider the class of Gibbs states characterised by a self-adjoint, positive and extensive operator on the original system H F given by, with λ J,α ∈ R + , ∀J, α. Extensive operators in GFT are a second quantisation of those loop quantum gravity operators which are diagonal in some intertwiner basis, such as the spatial volume operator [4,20,49].
Extensive operators are compatible with the total number operator, that is, [P, N ] = 0. They are thus diagonal in an occupation number basis of H F consisting of multi-particle, multi-mode states of the general form, J2α2 ... |0 (5.5) which are orthonormal, The partition function for (5.2) with operators P can then be evaluated 7 in the basis of states (5.5) to give, where µ < min(λ Jα ) for all J and α. This is a grandcanonical state of a-particles, which essentially describes a gas of these atoms of space with a changing total number in the given (non-tilde) subsystem. Technically, the number operator in (5.2) simply implements a constant shift by µ in the observable spectrum, which for now can be neglected by fixing µ to an arbitrary value, or by simply replacing λ Jα − µ → λ Jα . It is evident that by construction, the parameter β controls the strength of statistical fluctuations in P, regardless of any other interpretations. Note however, that one can reasonably inquire about its geometric meaning, especially if the operator P, which is its thermodynamic conjugate, has a clear geometric interpretation. It would thus be interesting to investigate this aspect in a concrete example where the choice of the observable is adapted to a physical context, like cosmology (see for instance [36]).
Finally, the ensemble average for number density will be useful for the construction of the associated thermal representation as outlined in section II. For the state above, it is given by the characteristic Bose distribution, Tr HF (ρ β,µ a † Jα a Jα ) = 1 e β(λJα−µ) − 1 . (5.8)

VI. THERMAL SQUEEZED VACUUM AND FINITE TEMPERATURE PHASE
Now that we have chosen a suitable class of thermal states of equation (5.2) characterised with operators P, we can proceed to define its associated thermal phase generated by a thermal vacuum and β-dependent ladder operators {b Jα , b † Jα ,b Jα ,b † Jα } β , along the lines detailed in section II. 7 Given a suitable set of properly normalised modes f J α as done in IV above, the calculation for this partition function is along the lines of that showed for a volume Gibbs state in [20], to which we refer for details.
Thermal Bogoliubov transformations 8 give the new ladder operators, mode-wise, along with their adjoints b † Jα andb † Jα . Inverse transformations are, and their respective adjoints. The β-dependent annihilation operators specify the thermal vacuum via the parameters θ Jα can be determined from Lastly, the β-phase that we have constructed here, being described kinematically by This can be seen directly from the transformation equations between the two vacua: Clearly, the pre-factor of the product of inverse cosh functions vanishes in general, without any cut-offs in the modes. This means that the overlap between the two vacua is zero, and the two representations built upon them are inequivalent. In other words, this transformation in field theory is ill-defined in general due to an infinite number of degrees of freedom, giving rise to inequivalent representations describing distinct phases of the system.

VII. COHERENT THERMAL STATES
We now define a family of coherent states in the thermal representation, called coherent thermal states. We understand them as defining thermal quantum gravity condensates, expected to be relevant in the studies of semi-classical and continuum approximations in discrete quantum gravity models based on polyhedral quanta of geometry. Indeed, unlike the purely thermal state |0 β , coherent states can encode a notion of semi-classicality, with which one can attempt to extract effective dynamics from an underlying quantum gravity model. For instance in group field theory, it has been shown that a coherent condensate phase of the a-quanta can support FLRW cosmological dynamics, and thus represent a viable choice of quantum gravitational phases that are relevant in the cosmological sector [37][38][39][40]48]. But a coherent state over the degenerate vacuum (4.11) is unentangled, and we expect a geometric phase of the universe to be highly entangled. Moreover these states cannot in themselves encode statistical fluctuations in different observable quantities. Therefore, from a physical point of view, the construction of coherent thermal states is amply justified.
Coherent thermal states [26,[50][51][52] are a coherent configuration of quanta over the thermal vacuum, implemented by displacing |0 β with displacement operators of the form, for σ ∈ H. To recall, the usual coherent states |σ ∈ H F of a-particles are, while |σ ∈H F forã-particles are, The tilde in the ket notation |σ simply means that the state is an element of the conjugate Hilbert space, and Dã is a displacement operator of the same form as (7.1) but for tilde ladder operators. The most useful property of these states is that they are eigenstates of their respective annihilation operators, a Jα |σ = σ Jα |σ (7.4) a Jα |σ = σ Jα |σ (7.5) which is at the heart of their extensive use as robust, most classical-like, quantum states.

VIII. SUMMARY AND OUTLOOK
In this article we have presented an implementation of the thermofield dynamics formalism in the context of group field theory quantum gravity. It opens the door to using such techniques in discrete quantum gravity, thus facilitating exploration of the phase structure of quantum gravity models characterised by generalised thermodynamic parameters β ℓ , and complementing renormalization investigations in group field theory [53][54][55] and possibly other related approaches [56][57][58][59].
Here, we have constructed finite temperature, equilibrium phases associated with a class of generalised Gibbs states in group field theory, along with identification of their non-perturbative thermal vacua. The vacua are squeezed states encoding entanglement of quantum geometric data, which in turn is expected to be a characteristic property of a physical quantum description of spacetime in general. We have introduced coherent thermal states which, in addition to carrying statistical fluctuations in a given set of observables, are also condensates of quantum geometry.
In group field theory, zero temperature coherent states have been used to obtain an effective description of flat, homogeneous and isotropic cosmology (flat FLRW), where certain quantum corrections arise naturally and generate a dynamical modification with respect to classical general relativity, preventing the occurrence of a big bang singularity along with cyclic solutions in general [37][38][39][40]. Encouraged by these results, the introduction of coherent thermal states may bring further insight and progress in the understanding and analysis of semiclassical and effective continuum dynamics of GFT models. The particular case of flat FLRW is studied in an upcoming paper [36].
It would also be interesting to understand better the connection of our thermal vacua with other works in loop quantum gravity concerning kinematical entanglement between intertwiners (related further to discrete vector geometries) [60][61][62], especially since our squeezed thermal vacua essentially encode entanglement between gauge-invariant spin network nodes, that is intertwiners, but now at a field theory level.
The thermal vacua can be extended even further to more general two-mode squeezed vacua between tilde and non-tilde quanta, between different tilde quanta, or even between different modes of the tilde quanta. For instance, condensates of correlated quanta, like dipole condensates [39], may be directly constructed and studied in this setup. Considering correlations between different modes of the quanta, which encode quantum geometric data, might also make comparisons with studies in loop quantum gravity mentioned above [60][61][62] more direct.
This setup could also prove useful for the study of quantum black holes. In group field theory for example, black holes have been modelled as generalised conden-sates [63], which must also possess related thermal properties. Suitably modified thermal coherent states may provide just the right type of technical structure to explore their statistical and thermal aspects further.
Finally, by providing quite a straightforward handle on collective, quasi-particle modes in discrete quantum gravity, while still allowing for access to different inequivalent representations, this framework may bring closer the studies of microscopic theories of quantum gravity and analogue gravity models [64].