Interacting fields in real-time AdS/CFT

We compute time-ordered 2- and 3-pt correlation functions of CFT scalar operators between generic in/out states. The calculation is holographically carried out by considering a non backreacting AdS scalar field with a $\lambda \phi^3$ self-interaction term on a combination of Euclidean and Lorentzian AdS sections following the Skenderis-van Rees prescription. We show that, although working in an essentially different set up, the final result for the 3-pt correlators agree with those of Rastelli et al. for Euclidean AdS. By analyzing the inner product between the in/out excited states in the large $N$ approximation, we argue that a cubic bulk interaction deforms the excited states from coherent into \emph{squeezed}. Finally, a diagrammatic interpretation of the results suggests some general properties for the $n$-point correlation functions between excited states.


Introduction
At its inception, large N two-and three-pt correlators for local CFT operators in the vacuum were computed in the AdS/CFT framework using a prescription that requires to solve for classical dual fields on Euclidean AdS in the presence of boundary sources [1,2,3,4]. Even though these results can be extrapolated to the Lorentzian AdS spacetime by analytic extension, some important ingredients of the real time description are left out from the beginning. For instance, amplitudes in arbitrary excited states 〈out|O (x 1 ) . . . O (x n )|in〉 and the analysis of intrinsic real time phenomena such as response functions. See [5,6,7] for previous related work. Skenderis and van Rees (SvR) have proposed a setup that extends the Gubser-Klebanov-Polyakov-Witten (GKPW) prescription [2,3] to Lorentzian signature and allows the calculation of n-point correlation functions of local CFT operators in real time. Among its virtues one finds the natural emergence of the causal Feynman propagator [8,9].
For the case of vacuum to vacuum scattering amplitudes, the SvR prescription requires to consider a Lorentzian AdS cilinder M L smoothly glued to two halves of Euclidean AdS M ± along the past/future spacelike surfaces Σ ± that limit the Lorentzian region, as shown in Fig.1(a). The outcome of their construction is that time ordered correlators for CFT operators O , inserted at the timelike boundary ∂M L , are computed in the large N limit, through the formula (1.1) On the right hand side S φ L ; φ ± = 0 stands for the on-shell action for Φ, the bulk field dual to the CFT operator O , which is solved on a multi-pieces geometry with: (i) vanishing Dirichlet conditions φ ± = 0 on the Euclidean asymptotic boundaries ∂M ± and (ii) arbitrary non-trivial φ L boundary condition on ∂M L , the latter being the external source from the CFT perspective. A picture of the whole manifold M , and Φ as a single field, is presented in Fig. 1(b).
In [8] it was also conjectured that in order to describe non-vacuum states, one should consider turning on non-trivial boundary conditions φ ± = 0 on ∂M ± . A detailed proof of this claim was given in [10] where it was explicitly analyzed for the case a free scalar field 1 , and it was also shown that states constructed this way were precisely given by where P stands for Euclidean time ordering. Consistency with the alternative holographic prescription called BDHM [14], then permitted infer that (1.2) are coherent states [10]. The prescription for the generating functional given by (1.1) thus generalizes, for arbitrary in/out states φ ± , as Here the right hand side denotes the bulk field action evaluated on the classical solution for Φ satisfying non-trivial boundary data φ ± on the Euclidean boundaries of M as shown in Fig. 1(a). As customary, time ordered n-pt O correlators are obtained by taking n derivatives with respect to φ L and then setting φ L → 0. The glued manifold M can be thought of as dual to a (complex) time evolution as shown in Fig. 2.(a). The purpose of the present work is to study the formalism presented above in the presence of bulk interactions 2 . To get some intuition we recall that the (free) field theory solution to the Lorentzian bulk equation of motion reads Φ 0 (z, t , x) = dt dx K (z, t , x;t ,x) φ L (t ,x) + Φ 0(N ) (z, t , x) . (1.4) The first term consists of the familiar boundary-bulk propagator K carrying the information contained in the source φ L into the bulk, while the second term, Φ 0(N ) , encodes the excited states structure into the normalizable mode 3 . In fact one can verify that whenever Φ 0(N ) = 0 the expectation value for O (t , x), obtained 1 For previous work regarding excited states in Lorentzian signature see [11,12,13] 2 Bulk interactions have been already been studied, for example in [15,16], but not in the context of semiclassical excited states as far as the authors are aware. 3 We would like to remind the reader that the normalizable modes are uniquely defined once the propagator ambiguities present in the Lorentzian context are fixed. Figure 1: Bulk dual to the In-Out formalism as a gluing of Lorentzian and Euclidean AdS manifolds.
by the standard prescription [2,3,8], becomes non-trivial This well known result [10,12,17] is interpreted as a manifestation of the CFT being shifted from its symmetric vacuum state. Bulk interactions of the form λΦ m will induce corrections to (1.4) so that, to first order in λ, the on-shell configuration Φ will now contain a polynomial in φ L of degree m, thus modifying every n-point function n ≤ m. In the following we will compute in the SvR setup time ordered 2-and 3-pt amplitudes of conformal (scalar) local operators for states (1.2), and explicitly verify that adding interactions to the gaussian bulk action result in the states losing their coherent character 4 . In particular, we show that to first order, the cubic self-interacting field turns (1.2) into squeezed states. The paper is organized as follows: in Sec. 2 we specify the scalar field model in which we are going to study eq. (1.3), this is, a cubic self-interacting real massive scalar field in AdS, dual to a scalar local operator in the CFT. In Sec.3 we construct the free field solution to the field in the presence of boundary sources, this is the necessary building block for computing the first order λ-correction to the on-shell action. In Sec. 4.1 we analyze the standard boundary term of the on-shell action arising from the quadratic piece, and in Sec. 4.2 the (bulk) self-interacting contribution. An analysis of the results is made in Sec. 5. Finally, Sec. 6 summarizes the results and suggest prospects for future work. We relegate to the Appendixes many explicit computations and technical details.

Interacting Scalar Field Theory on AdS
Let us consider the simplest example of interacting fields on a global AdS spacetime background: a real massive scalar field with a cubic self interaction, which should be enough to see corrections to the 1-and 2-pt functions and also to the inner product between the (1.2) states. The generalization to multiple scalar field becomes straightforward after considering this minimal example. The action we are going to work with in the In-Out formalism ( Fig. 2(a)) is, where the integration runs over the d + 1 dimension manifold M ≡ M + M L M − obtained by appropriately gluing Euclidean and Lorentzian AdS regions as shown in Fig.2(b) (see [8,10] for details). The equation of motion following from (2.1) is which we solve perturbatively in λ expanding the field as Here the free field solution Φ 0 meets the {φ ± , φ L } boundary conditions on ∂M , whereas every other Φ i have Dirichlet conditions over the asymptotic boundary.
The on shell action results where γ is the induced metric on the boundary ∂M and n µ is the outgoing unit normal vector. When writing the third line we have used the equation of motion (2.2). These expressions are rather formal though, as an appropriate prescription is required for imposing the asymptotic boundary conditions on Φ to avoid divergences. We will choose to work with the so called -prescription [2,20] which consists in regularizing the problem by setting the boundary conditions at a fixed radial distance z = 1 in Poincare coordinates, with the → 0 limit taken at the end of the computations. Inserting (2.3) into (2.5) one finds where the first term is the contribution of the free field action, the second term 5 will be shown to be zero in App. A and the third term will give rise to first order corrections in correlation functions. Being the second term in (2.6) absent, Φ 0 is the only relevant piece of Φ necessary to compute the first order corrections in λ to correlation functions. We will devote the next section to build such a solution for the In-Out path. A somewhat technical but relevant comment regarding the structure of this work is related to the treatment of the expression coming from the free on shell action, containing only boundary terms, and the one coming from the interaction terms, containing bulk integrals. As noticed in [4], there is a slight difference between the Asymptotic prescription as performed in [3] and the -prescription as defined in [2] when regularizing the divergences appearing from the asymptotic boundary of AdS. The two techniques give different normalizations for the CFT 2-pt function (which gets its contribution from the boundary term), but give identical results for any other higher point function, which arise from bulk terms. While the first prescription leads to easier computations, the latter is more natural in the sense that it automatically meets the 5 The contribution ∂M γ Φ 1 n µ ∂ µ Φ 0 is absent in (2.6) since, as stated above, Φ 1 = 0 on ∂M .
Ward identities between two and higher order point functions. With this in mind, we are going to follow the -prescription when treating the free contribution of the on-shell action and follow the asymptotic prescription in the interacting terms.

Free Field Solution Φ 0 on In-Out path
In this section we build the leading order solution Φ 0 in M by solving the first equation in (2.4). To this end we will first find the most general solutions on M L and M ± and afterwords impose continuity of the fields and their conjugated momenta across the Σ ± gluing surfaces. This procedure will determine uniquely the solution.

Lorentzian Section
The Klein-Gordon equation of motion in M L is with the Poincaré AdS d +1 metric given by Inserting the ansatz Φ L 0 = e −i ωt +i kx f ωk (z) in (3.1) one finds In the -prescription equation (3.1) is supplemented by a boundary condition at z = given by where ∆ ≡ d /2 + ν and ν = d 2 /4 + m 2 . In the present work we consider the case where ν ∈ N 6 . The most general solution to (3.1) satisfying (3.4) can be written as (3.5) The first line, which we will refer as NN-solution, fulfills (3.4) provided we take the boundary-bulk propagator K to be given by [2,20] K (z, t , Here J ν is the Bessel function of the first kind and K ν the modified Bessel function of the second kind. A small imaginary piece −i 0 + is added to q to properly define the momentum integrals appearing in (3.6) as shown in Fig. 3(a). As discussed in [9,12], this choice leads to the Feynman propagator. The second line in (3.5), or (Normalizable) N-modes, specify the form of Φ 0(N ) that we mentioned in (1.4). They involve a combination of Bessel functions and correspond to solutions which by construction vanish at z = . We have written the second independent solution to (3.1) (second term in the last parentheses of (3.5)) as a K ν Bessel function with imaginary argument. Normalizability demands the momentum domain for N-modes to be timelike, i.e. (ω 2 − k 2 ) ≥ 0. We have explicitly separated the positive and negative frequencies in (3.5), and in the following, ± denote that the integration in the first differential variable runs over R ± . Every other integration variable should be taken over R. We will show below that the coefficients L ± ωk in (3.5), arbitrary in principle, will get determined after imposing continuity conditions for Φ 0 across Σ ± (see [8,10]).

Euclidean Sections
The metric and the equations of motion for the field in the Euclidean section can be obtained by Wick rotating (3.2) and (3.3). The crucial feature for M ± is that τ runs over the half line R ± , allowing for the existence of Normalizable modes. The general solutions to the Klein-Gordon equation in M ± are therefore with Euclidean boundary-bulk propagator The main difference with the Lorentzian form (3.6) is that no cuts appear along the integration contour, as can be seen from Fig. 3(b). Notice that the N-modes involve real exponentials in τ with ω > 0, thus normalizability admits the e ∓ωτ behavior for M ± . The coefficients E ± ωk are the euclidean counterparts of L ± ωk and are fixed by the continuity conditions across Σ ± .

Gluing the Euclidean and Lorentzian sections
The continuity conditions across Σ ± , uniquely fix the L ± ωk , E ± ωk coefficients in (3.5) and (3.7) to be where φ L ωk and φ ± ωk are the Fourier components of the sources defined in (B.5) and (B.9). See App. Bfor a derivation of these expressions. Similar relations were obtained for the In-In formalism in [18].

On-shell action perturbation theory 4.1 Free field boundary term
In this section we generalize to Poincare coordinates and d dimensions the results obtained in [10]. In particular, we compute for the excited states (1.2) the inner product and the 1-and 2-point functions of local operators of conformal dimension ∆. To this end we need to evaluate the on-shell action.
We start with the Gaussian contribution S 0 in (2.6). Inserting the solutions (3.5) and (3.7), with L ± ωk , E ± ωk coefficients given by (3.9), back into S 0 we find a sum over the three sections displayed in Fig.2 We have explicitly replaced the volume element γ = −d , inserted the boundary conditions (3.4) and omitted the arguments of the functions to shorten the expression. Only the asymptotic (z = ) boundary contribution appear in the expression since the gluing procedure guarantees that the on shell action pieces arising from the Σ ± surfaces cancel each other [9]. The asymptotic analysis of the boundary terms in (4.1) shows that both the NN-solution and N-modes contribute to the physical observables. We work out explicitly the Lorentzian piece, the Euclidean cases being analogous. The NN-solution in (3.5), for z ≈ behaves as As it is well known [2,3], the above expansion shows two distinct features: an analytic (leading) piece given by a polynomial P ν−1 q 2 of order (ν − 1) and a non-analytic (subleading) piece given, for integer ν, by the ln q terms. The former give rise to contact terms in S 0 and will be dropped in the following. This is formalized in App. A. The latter give rise to the familiar CFT propagator when transformed back to configuration space. As for the N-modes, second line in (3.5), an -expansion gives Notice that a leading ∆ behavior was factored out, both in (4.2) and (4.3), which compensates the divergent −∆ factors appearing in (4.1). Higher orders in (4.2) and (4.3) are thus unimportant as they will give vanishing contributions to (4.1) in the → 0 limit. Notice that the contribution from (4.3) does not yield contact terms in spite of being analytic in (ω 2 − k 2 ) since the N-modes are solely integrated over a timelike domain (see (3.5)). See App. C for detailed calculations. We now have all the ingredients to obtain the on-shell action S 0 . Inserting (4.2) and (4.3) in (4.1) leads to momentum integrals that we carry out in detail in App. C. The expressions adopt a compact form if we define the distance between points on the contour in Fig.2(a) as where the complex time variable η parametrizes the In-Out path shown in Fig.2(a) 8 . For example, we shown 7 The on-shell action, both the free and interacting parts, should be understood as a functional only of the boundary sources , for the ease of notation in what follows we will suppress these dependences. 8 One may alternatively parametrize the time path as in the App. C that the term bilinear in the Lorentzian sources come up Feynman regulated with For the Lorentzian section M L the free on shell action becomes and similar expressions are obtained for the Euclidean sections M ± Notice the appearance of crossed terms between the Lorentzian and Euclidean sources in these expressions which add up in (4.1). Their consequences will be explored below. In particular, using (1.3) we will evaluate the inner product between excited states and compute the 1-and 2-pt correlation functions.
Inner product: The inner product between excited states (1.2) can be computed by collapsing the Lorentzian section (∆T = (T + − T − ) → 0) in the absence of Lorentzian sources [10]. This amounts to consider only the first terms in S + 0 and S − 0 . Defining the inner product can be rearranged to give [10] ln 〈φ + |φ − 〉| λ=0 = lim in agreement with the well known expression in [3].

1-pt correlation function:
The 1-point function arises from the linear terms in φ L in S 0 , Notice that this expression corresponds to a propagation of the boundary conditions φ ± on the Euclidean sections to the Lorentzian section. When performing the integral recall thatx µ = (T ± ∓ iτ) forτ ≷ 0.

Connected 2-pt function:
The second term in (4.1) is the relevant one for computing the 2-pt connected correlator. The result is Eqs. (4.6)-(4.8) review and generalize the results in [10], where they have been thoroughly analyzed.

Self-interaction bulk contribution
In this section we compute the first order correction in λ to the on shell action. We already mentioned that the second term in (2.6) vanishes (see App. A), thus the term we are going to work with is where we have used the complex time variable η introduced in (4.4). As mentioned in Sec 2, we will work out the bulk integrals involving the non-linear interaction in the Asymptotic prescription. Appealing to results in [4], we will find a closed analytic expression for the contributions. In passing we will give a diagrammatic interpretation of the results.
To compute (4.9) we need the expression for Φ 0 in the Asymptotic prescription. This can be obtained from expressions (3.5),(3.7),(3.9) by taking → 0. For the NN-piece, the first line in (3.5), the limit leaves a momentum integral which can be computed analytically, giving [3,9] notice that the boundary-bulk propagator comes Feynman regulated as mentioned in the paragraph below (4.4). We now show that the N-modes containing the information of the excited states (second line in (3.5)) can be written as a convolution involving the boundary sources φ ± and a generalized boundary-bulk propagator. Taking the → 0 limit on the N-modes piece, with L ± ωk given by (3.9), one again finds a momentum integral which can be explicitly carried out (see App. C) with φ E defined in (4.5) and the distance |x µ −x µ | 2 in (4.4). Analogous expressions are obtained for the fields on M ± . Recalling the complex time variable η introduced in (4.4) motivates the definition of a single source which allows to write Φ 0 over the mixed signature manifold M as a single complex integral which puts together the contributions from the Lorentzian and Euclidean sources. Summarizing, the final expression for Φ 0 in the Asymptotic prescription taking into account all sources is (4.13) Inserting the expression above in (4.9) leads to an integral over M which is known to give [4], where the arguments x µ i of the sources φ defined in (4.12) lie on ∂M = ∂ z M − ∂ z M L ∂ z M + , and the expression for R is given by Example of a first order correction diagram contributing to the 1-pt function. It is obtained from (4.14) by stripping away one Lorentzian source and removing its accompanying integral. The final result for the first order correction is given by (4.18).
Since (4.14) contains powers of φ L up to cubic order, we get corrections to the inner product and connected 1-, 2-and 3-pt correlators. The leading order on shell action for our system is thus Inner product λ-correction: From the on-shell action we can obtain the inner product between the excited states. This is found by turning off the Lorentzian sources and setting the (real time) interval of evolution to zero [10] ln〈φ (4.5). In analogy with the non-interacting case [10], the left hand side, which reduces to the computation of the generating functional for a source φ E on the boundary of the Euclidean AdS, matches [4].
The excited states given by (1.2) are not normalized. Defining |φ ± 〉 N ≡ N ± |φ ± 〉, the appropriate normalizing factors follow from N 〈φ ± |φ ± 〉 N = |N ± | 2 e i S 0 (1 + i λS 1 where we have expanded to first order in λ. Recall that upon taking ∆T → 0 limit, i S 0 and i S 1 become real, see (4.6),(4.16). After some manipulations, the inner product of normalized states becomes ) (4.17) which naturally gives one when (φ + (τ, x)) ≡ φ + (−τ, x) = φ − (τ, x). This is the conjugation rule given by [21]. [10], is induced by the inner product on the space The factor in parentheses modifying the Gaussian in (4.17) indicates that, in the presence of interactions, the excited states deviate from being strictly coherent.

1-pt function λ-correction:
The first order correction to the 1-pt function (4.7) is obtained by taking the derivative of (4.14) with respect to φ L . It yields In the ∆T → 0 limit, i.e. when no time evolution takes place, the expression above becomes the matrix element of O between the states |φ ± 〉. The correction is again in line with the deformation of coherence, as the matrix element are no longer linear in the Euclidean sources [10,12].

2-pt function λ-correction:
Taking the second derivative of (4.14) with respect to φ L yields the first order correction, Notice that the correction, depending on the arbitrary profile φ E , does not correspond to an anomalous dimension as other type of corrections coming from bulk interactions, cf. [16].

3-pt function:
A third φ L differentiation of (4.14) leads to the 3-pt function of the scalar operator O , This is the time-ordered Lorentzian extension of the result in [4]. One notices no dependence on φ E in (4.20). This happens as a consequence of the cubic nature of the bulk interaction vertex form, and is similar to the result found in [10] for the 2-pt function in a free bulk theory. However, in the present case, the φ E sources independence in (4.20) is only valid up to linear order in λ. This will become clear below from a diagrammatic point of view.

Diagrammatic Interpretation of the results:
The result (4.13) and the bulk interaction vertex (4.9) allow a diagrammatic interpretation of the corrections found above (see Figs.4 and 5). The boundary-bulk propagator, given by the parentheses in (4.13), carrying the information of φ(η, x) from ∂M to a point inside M , is represented as a wavy line. The cubic interaction (4.9) maps to three wavy lines meeting at a point. Higher order λ corrections involve bulk-bulk propagators joining bulk vertices which can also be represented by wavy lines. The colored wavy lines (red/blue) emphasize dependence on the initial/final states. By construction (see (1 It is easy to understand why no Euclidean sources appear in the 3-pt function correction (4.20) as no vertices legs are available to reach the Euclidean regions from a single bulk vertex point. This also gives a alternative way for understanding the result in [10] showing that (4.8) was independent of the euclidean sources at 0-th order in λ.

Squeezed states
Generically the AdS/CFT correspondence in the large N limit amounts to consider classical bulk dynamics at leading order. This manifests in the full fledged quantum theory in terms of coherent states as discussed in [10]. In the present section we will consider the consequences of turning on bulk interactions of the form (4.9) on the quantum states and show that the coherent states (1.2) deform into squeezed ones.
To this end, we take profit of the BDHM prescription [14,22]. The equivalence between the boundary sources technique (GKPW) developed in [2,3,8] and the Fock space construction (BDHM) [14,22] was shown in [23] . The BDHM framework consists in canonically quantizing AdS bulk fields and obtaining dual CFT operator through the relation 10 The normalization factor N ν ≡ 2ν is required for a precise match between the prescriptions [10,23,25].
From the action (2.1), one builds an interacting quantized field Φ as a perturbative λ-expansion. The free field Φ 0 is a linear combination of the orthonormal set of eigenfunctions of the KG equation with ladder operators as coefficients, i.e. [12] Φ 0 (z, t , x) ≡ π From this last expression and (5.1) we obtain the 0-th order CFT operator, 3) The first order quantum correction to the quantum field which follows from (2.4) is where G(z, t , x;z,t ,x) is the (standard) bulk-bulk AdS propagator [43] and :Â : denotes normal ordering of the operatorÂ [26]. Using the relation [20] lim z→0 the first λ correction to the CFT operator becomes To first order we therefore get showing that the correction O 1 leads to quadratic terms in the ladder operators. The final step is to notice that the path ordered exponential in (1.2) can be thought, for computational purposes, as a time evolution operator with Hamiltonian d x O (τ, x)φ − (τ, x), evolving in Euclidean time the ground state to our initial excited state (1.2). From (5.5) we see that the "time evolution" is quadratic in the ladder operators. We now make use of the results in [27] where, by use of disentangling theorems, it was shown that a generic quadratic operator with the parameters α ωk (t ) and ζ ωk;ωk (t ) determined by the coefficients in (5.6) [27]. In short, to linear order in λ, the states (1.2) can be dissentangled as squeezed states. In the present case we get The "time" dependence in (5.7) is understood to be evaluated at Euclidean time τ = 0 corresponding to an evolution along the whole M − manifold.

Multiple scalar fields and entanglement
Through this work we have consistently worked with a single scalar field. Nevertheless, the formalism can straightforwardly be generalized to the case of many fields by using (4.13) 11 . We show below that generically, interactions among bulk fields lead to entangled CFT states.
As a simple example, let us consider three scalar fields Φ I , I = 1, 2, 3 with an interaction term given by Therefore, the states space can be expressed as The BDHM prescription [14,22] gives us the operator By expressing Φ I ;0 (z, t , x) in terms of the ladder operators a I , the λ correction O J K ;1 becomes a linear combination of the operators a † J a † K , a J a K , a † J a K (5.9) So for instance, the excited states associated only to sources φ 1 = 0 read as thus the second term in the exponent, proportional to the coupling λ 1J K , involves products a † J a † K that typically characterize maximally entangled states in the space H J ⊗ H K . As an aside, notice that thermal states in the TFD formalism also have this form [28,29,30,31].
This result suggests a possible connection with some recent ideas and results pointing that (emergent) classical geometry, and field configurations on it, should be intimately related to entangled states in the boundary quantum field theory [32,33,34,35,36]. More recently these ideas have been developed in the context of MERA networks [37,38,39].

Conclusions
In a previous article [10] it was argued that the states (1.2), prescribed by the SvR set up for a free bulk theory are coherent states. This can be understood by the fact that the large-N approximation describes the semiclassical regime of the bulk theory, such that only semiclassical states should make sense. Moreover, these states are strictly coherent in the bulk representation of the Hilbert space, since the N → ∞ limit of supergravity reduces to a free theory, in the sense that, with a suitable normalization [18], only two point correlation functions remain non-vanishing. In the present work, motivated by previous literature [4,14], we have considered the simplest toy model for an interacting bulk theory 12 . We have analyzed the nature of these excited states and explicitly shown in Sec. 5 that the bulk interaction slightly deviates the in/out states from coherent to squeezed states.
After setting notations in Secs. 2 and 3 we computed in Sec 4.1 the 0-th order contributions to the inner product, 1-and 2-pt function. In Sec. 4.2 we analyzed the first order corrections arising from the self interacting (bulk) contribution to the on-shell action, in particular we gave a (Witten) diagrammatic representation for the contributions to correlation functions in the SvR framework in Figs. 4 and 5, matching the intuition one has from classical perturbation theory. It is important to stress that Euclidean and Lorentzian sources stand in different footing. Euclidean sources are never turned off and prescribe the in/out states, concomitantly fixing the normalizable modes propagating into the Lorentzian bulk. On the other hand, Lorentzian sources serve as tools to obtain the correlation functions and are set to zero at the end of computations.
From the results in Sec 4.2, we can infer some general properties of the CFT n-pt correlation functions arising from a self interacting λΦ m bulk dual: up to linear order in λ, the correction to the n-pt correlation function will arise from the on-shell action term involving m − n Euclidean sources. These can be pictured as diagrams analogous to Figs. 4 and 5. Therefore, up to leading order in λ we can conclude that n-pt functions with n < m will depend on the excited states profile, while the (connected) m-pt function will be non-zero but independent of the excited states. It is worth noticing that the squeezed character of states (1.2) only follows for a cubic interaction.
In a future work, we aim to study a more realistic set up based on supergravity models that captures the space-time back-reaction and how the coupling with the spin 2 field h αβ would affect the form of excited states (1.2). Typically, the bulk interaction between many SUGRA fields should produce entangled states as pointed out in Sec 5.

Acknowledgments
Special thanks are due to Martin Kruczenski for useful comments and discussions. PJM also would like to thank Ariana Christodoulou for helpful correspondence. This work was supported by projects PICT 2012-0417 ANPCyT, PIP0595/13 CONICET and X648 UNLP.

A Holographic Renormalization in -prescription
In this appendix we pursue two objectives: first, we show explicitly the vanishing of the second term in (2.6), which leads to the first order correction to the on-shell action arising solely from (4.9); secondly, we construct the adequate counter-terms needed to take care of the contact terms in (4.2). The Holographic Renormalization method in the regularized space (i.e. -prescription), which we will review in what follows, will take care of both issues. We start by ignoring the N modes in the solutions, thus considering only the NN solutions involving the bulk-boundary propagator K defined in (3.6). We then show that considering the N modes leaves the results unaffected. For the sake of self-consistency, we present some formulae that will become useful in what follows. The first order solution to (2.2) disregarding the N modes is with Φ 0 given only by the first line in by (3.5) and G the Feynman bulk-bulk Green function in the regularized space, i.e. 13 − m 2 G (z, t , x;z,t ,x) = 1 The relevant property of G we will use in what follows is the relation between the bulk-bulk and boundarybulk propagators z∂ z (G (z, t , x;z,t ,x))| z= = − ∆ K (z,t ,x; t , x) .
The on-shell action (4.1) involves radial derivatives of the boundary-bulk propagator K , which lead to divergences (see (4.2)) upon taking the → 0 limit. The expansion that will become useful below is where ψ(x) is the DiGamma function. The divergent terms in the on-shell action will mainly come form the second line above, containing integer powers of q 2 that lead to contact terms in configuration space and an additional logarithmic divergence arise from the first line for k = 0 that leads to a matter conformal anomaly [44]. We build below appropriate counter-terms for each of these divergent terms. Inserting (A.1) in (2.5) we find to first order in λ where the outward normal is n µ ∂ µ = −z∂ z . The cancellation between the second and third term in (A.4) follows by using (A.2) in the second term in (A.4), giving [20] λ which exactly cancels the third term in (A.4). One may worry that in the presence of N-modes, the integral in parentheses no longer gives Φ 0 , since the second line in (3.5) is missing. Nevertheless, the cancellation persists upon considering the S ± contributions in (2.1). A nice outcome of the SvR construction is that one can write, on the complete manifold M , the field solution to (2.2) as (4.13) packaging both the NN-and N-modes information in terms of a boundary-bulk propagator. Therefore (A.5) remains valid as long as we integrate over the whole manifold M . We now devote ourselves to build the counter-terms needed to obtain a finite on shell action. Thedivergences in (A.4) arise from the first term 14 . Each of these divergences take the form of boundary terms, and as such, can be subtracted without altering the equations of motion. We will work out the ν = 1 example, the general ν ∈ N case follows the same procedure. For completeness we mention that the treatment of the ν = 0 case, corresponding to the Breitenlohner-Freedman lower mass bound [19], differs slightly from the general integer case. The reason being a logarithmic decay in the field when approaching the boundary. The boundary condition (3.4) modifies to Φ( , t , x) = ∆ ln( ) φ(t , x) and an interesting outcome to point out is that the coefficient in (4.8) changes to Γ(∆)/(2π ∆ ), which can be readily seen to be necessary since ν/Γ(ν) would make the result trivial otherwise [4].
Counter-terms for ν = 1 By using the expansion (A.3) for ν = 1, the leading -terms for the normal derivative in the first term on (A.4) turn into where we used the relation (3.4),˜ ≡ 1 2 e γ , γ is the Euler-Gamma number and γ ≡ 2 η i j ∂ i ∂ j , with i , j = 1, . . . , d is the induced D'Alambertian on the boundary. The first term of (A.4) can therefore be written as The first line of (A.7) shows that the divergent terms of the on shell action (in the → 0 limit) can be written as local functions of the boundary values Φ 0 | = d −∆ φ L , therefore they can be removed by adding identical terms with opposite signs. The second line in (A.7) is the relevant term in (4.2) and gives rise to the expected 2-pt function for a conformal operator with conformal weight ∆, as we show in App. C. We would like to mention that in the -prescription we follow in the present work, the relationship between the boundary condition for the field at the radial boundary and the CFT source is simply Φ 0 | = d −∆ φ L as compared to Φ 0 = z d −∆ φ L + . . . for z 1 in the Asymptotic prescription. The simple boundary condition (3.4) avoids the iterative procedure one finds in determining the counter-terms in [44]. The second term in the first line of (A.7) involving ln gives rise to the matter conformal anomaly of the dual CFT [20,44]. Notice that this term appears as a consequence of ν being a positive integer. The counter-term action to be added, for ν = 1, can be written as

B Solutions near Σ ±
In this appendix we derive the equations (3.9) from the continuity conditions (3.8). We start by reminding the reader that the continuity condition following from the path integral formulation is imposed on the field and its conjugated momentum. Starting from the definition, and using the complex time variable η defined in (4.4), it follows from the action (2.1) that which is the second equation in (3.8).
The continuity conditions (3.8) on Σ ± give linear relations between L ± and E ± which determine them uniquely in terms of the Euclidean sources φ ± . To compute the relations we can safely take the → 0 limit from the outset. This is allowed because the sources φ L,± turn off at the vicinity of Σ ± . The absence of sources guarantee that the field configuration can be expanded in terms of N-modes, for which the regularization is superfluous. From a mathematical point of view, this manifests in (4.3) where one sees that the leading ∆ in front only requires the 0 information of the coefficients L ± , see discussion below (4.3).
We will transform to momentum space via and make use of the following properties [45,46] with q defined in (3.6). We will now find the behavior of (B.1) near the surfaces Σ ± for the Lorentzian and Euclidean solutions (3.5) and (3.7) in the → 0 limit.
Lorentz Section: We will analyze the NN-and N-pieces in Φ separately. The momentum components of the NN-piece in (3.5) read Θ ω 2 − k 2 (ω 2 − k 2 ) Theω integral in (B.4) is solved by the Residue theorem, with the integration contour closing in the upper/lower half plane depending on Σ ± . As an example, taking into account that the CFT sources φ L (t ,x) lie before Σ + , for t close to Σ + it follows that (t −t ) ∼ (T + −t ) > 0, therefore theω-path must be closed from below. For Σ − the opposite contour is needed. Summarizing, the NN contribution to (B.1) gives (B.4) = Θ(ω 2 − k 2 ) 2πω i χ ωk φ L ±ωk e ∓i ωt , t near Σ ± ; χ ωk ≡ π ω 2 − k 2 where we dropped the i 0 + and the absolute value since ω > 0 for N-modes. Notice that only the time-like Fourier components of φ L excite the normalizable modes (see [12]). The momentum components of the N-modes are where we have used (B.2). Notice that one of the Θ functions can be thrown away as it becomes redundant.

Euclidean Section
We perform similar calculation for the Euclidean solution (3.7). The NN piece gives where we defined the Euclidean Fourier transform of the source as

C Useful mathematical results
For completeness, we devote this Appendix to carry out some sample integrals that lead to the results (4.6)-(4.8).