Boundary conditions and vacuum fluctuations in AdS4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm {AdS}}_4$$\end{document}

Initial conditions given on a spacelike, static slice of a non-globally hyperbolic spacetime may not define the fates of classical and quantum fields uniquely. Such lack of global hyperbolicity is a well-known property of the anti-de Sitter solution and led many authors to question how is it possible to develop a quantum field theory on this spacetime. Wald and Ishibashi took a step towards the healing of that causal issue when considering the propagation of scalar fields on AdS. They proposed a systematic procedure to obtain a physically consistent dynamical evolution. Their prescription relies on determining the self-adjoint extensions of the spatial component of the differential wave operator. Such a requirement leads to the imposition of a specific set of boundary conditions at infinity. We employ their scheme in the particular case of the four-dimensional AdS spacetime and compute the expectation values of the field squared and the energy-momentum tensor, which will then bear the effects of those boundary conditions. We are not aware of any laws of nature constraining us to prescribe the same boundary conditions to all modes of the wave equation. Thus, we formulate a physical setup in which one of those modes satisfy a Robin boundary condition, while all others satisfy the Dirichlet condition. Due to our unusual settings, the resulting contributions to the fluctuations of the expectation values will not respect AdS invariance. As a consequence, a back-reaction procedure would yield a non-maximally symmetric spacetime. Furthermore, we verify the violation of weak energy condition as a direct consequence of our prescription for dynamics.


Introduction
One of the most remarkable outcomes of string theory was the proposition of the AdS/CFT correspondence [1]. It is conjectured that a theory of quantum gravity B Vitor S. Barroso vitor.barrososilveira@nottingham.ac.uk Extended author information available on the last page of the article on n-dimensional AdS displays an underlying equivalent conformal quantum field theory without gravity, taking place at the (n − 1)-dimensional conformal boundary of AdS. Accordingly, applications to high energy and condensed matter physics appeared within the efforts to test the limits of this new conjecture, placing the anti-de Sitter spacetime under the scientific spotlight.
Although most of the developments in AdS rely on string theory techniques, on a recent work [2], the authors have focused on studying semiclassical properties of the spacetime. Using the mathematical apparatus of quantum field theory (QFT) in curved spaces, they have found the fluctuations of the expectation values of the energymomentum tensor and the field squared in AdS n . However, they did not discuss in depth the implications of the causal structure of the spacetime, i.e., the effects of non-globally hyperbolicity.
Since AdS has a conformal boundary, we may not be able to determine much about the history of a physical quantity without specifying its behavior at infinity. Such a circumstance poses a fundamental issue on the quantization procedure, namely the solutions of the wave equation will not be uniquely defined by initial conditions in AdS, i.e., the Cauchy problem is not well-posed. Thus, unless we give extra information at the conformal boundary, the lack of predictability makes it impracticable to build a quantized field whose dynamical evolution comprises the entire history of the spacetime.
Avis et al. [3] were the first ones to address the causal pathology of AdS when solving field equations. They developed QFT on AdS 4 by regulating information leaving or entering the spacetime by hand. Their approach proposes the imposition of boundary conditions at the spatial infinity in order to control whether information flows through (or is reflected by) the conformal boundary. Even though Avis et al. provide us with physically consistent solutions to the wave equation, works by Wald [4] and Ishibashi [5,6] reveal that a broader category of boundary conditions might be employed to obtain a physical dynamical evolution.
In [5], the authors present a prescription for dynamics of fields in general nonglobally hyperbolic spacetimes based on the grounds of physical consistency. In order to fulfill some reasonable physical requirements (to be explained later), they argue that the spatial component of the differential wave operator must be self-adjoint. Besides, in [6], they show that the prescription for dynamics in AdS translates into specifying boundary conditions at the conformal boundary. While Kent and Winstanley [2], impose the Dirichlet boundary condition at infinity, they are neglecting an entire set of non-equivalent dynamical outcomes. According to Ishibashi and Wald [6], those outcomes would correspond to the various boundary conditions that one could have specified at infinity.
In this paper, we study physical effects that may arise due to non-Dirichlet boundary conditions at the conformal boundary. We investigate those effects by computing the vacuum fluctuations of the expectation values of the quadratic field and the energymomentum tensor for conformally coupled scalar fields in AdS 4 . Also, we will keep Ref. [2] as a basis for our results and shall return to it for further comparison.
An arbitrary Robin boundary condition (which is neither Dirichlet nor Neumann) introduces an energy scale. We expect that this should be responsible for the break of AdS invariance. In two recent papers [7,8], the authors showed that this is indeed the case, at least for AdS 2 in Poncaré coordinates, where there is only one boundary condition parametrized by β ∈ R. In particular, it was shown that the Green's function does not respect the spacetime symmetries since it does not depend exclusively on the geodesic distance.
Is it possible to obtain a similar effect in global AdS n ? The answer is yes, as we will see later in this paper. However, as can be seen in Ref. [9], the expansion of the Green's functions in terms of the normal modes is highly nontrivial for a general Robin boundary condition and the search of a closed form seem to be hopeless. Nevertheless, it is still possible to break AdS invariance by choosing a nontrivial boundary condition. Since a scalar field in global AdS n has a simple decomposition in terms of angular momentum modes and each mode has its own dynamics, it is possible, as a first attempt, to change the boundary condition in a finite number of these l modes. Clearly, each l mode with a nontrivial boundary condition will have its frequencies spectrum modified. The resulting system becomes numerically tractable, but is this choice physical? Since this question does not have an immediate answer, we will try to set up at least the compatibility of our model with basic concepts from quantum field theory.
In this way, we consider ou model as a first step in showing how non trivial boundary conditions may bring highly nontrivial effects to quantum field theory. Besides, by considering backreaction, we expect to see non trivial effects even in the classical theory (future research should consider these effects more carefully).
We have organized this article as follows. In Sect. 2, we briefly review some of the fundamental aspects of the anti-de Sitter solution. Then, in Sect. 3, we display the systematic procedure that describes the dynamics of scalar fields in non-globally hyperbolic spacetimes-such as AdS-first presented by Wald and Ishibashi. With that scheme in hands, we show the implications their prescription has on scalar fields propagating on AdS, in Sect. 4. Our next step is to build the proper Green's functions in Sect. 5, and employ them in the computations of the renormalized quantities of interest, namely the fluctuations of the expectation values of the field squared and the energymomentum tensor, both shown in Sect. 6. Finally, we discuss our results in Sect. 7.

Anti-de Sitter spacetime
Surfaces of constant negative curvature are well-known in geometry and comprise the set of hyperbolic spaces. In the context of General Relativity, the equivalent to those spaces is the n-dimensional anti-de Sitter space, which appears as a solution to Einstein equations when choosing a negative cosmological constant (Λ < 0) in the absence of matter and energy. Setting Λ := − (n−1)(n−2) 2H 2 , we may write the Einstein equations as The outcome is an n-dimensional maximally symmetric pseudo-Riemmanian metric defined over a Lorentzian manifold with constant negative curvature, i.e., the AdS n spacetime. In a suitable set of parametrized coordinates {x μ }, 1 the line element for the induced metric g μν on AdS n is where dΩ 2 n−2 is the line element on a unit (n − 2)-sphere.

Topology
We may understand AdS n as an isometric embedding of a single sheeted ndimensional hyperboloid in an (n + 1)-dimensional flat space provided with metric diag(−1, 1, . . . 1, −1). Timelike curves in AdS are transverse sections of the hyperboloid, and they are always closed. The periodicity of the timelike coordinate, τ , suggests that given a point in spacetime, we can return to it by only traveling along a timelike geodesic of length 2π in τ . Accordingly, the topology of AdS n becomes apparent, namely S 1 × R n−1 , which is compatible with the existence of closed timelike curves. Thus, unphysical events can take place in the spacetime, such as a particle returning to the same position through a periodic motion in time.

Causal structure
Wald remarks in [10] that observers following closed timelike geodesics would have no difficulty altering past events hence breaking causality. In an attempt to solve this primary issue, we can 'unwrap' the hyperboloid along the timelike direction, and patch together unwrapped hyperboloids one after the other. In other words, we construct a spacetime spatially identical to AdS but extended in time, i.e., the temporal coordinate no longer ranges from −π to π but from −∞ to ∞. We refer to such procedure as the universal covering of AdS, and the resulting spacetime as CAdS.
Even though the unwrapping of AdS prevents the existence of closed timelike curves, another fundamental causality issue remains, namely the lack of predictability associated with fields propagating on the spacetime. Indeed, no Cauchy hypersurfaces exist in AdS (and CAdS) hence portraying it as a non-globally hyperbolic spacetime. The Cauchy problem will not be well-posed, yielding non-unique dynamics for a given set of initial conditions. We can understand this scenario as a result of information leaking through the spatial infinity of the spacetime, i.e., flowing in (out) from (through) the boundary. In order to solve such a pathological behavior, we shall discuss in the next sections how to adequately address causality issues associated with field equations in non-globally hyperbolic spacetimes.

Scalar fields in non-globally hyperbolic static spacetimes
An extensive literature (see, for instance, [11] and references therein) provides a complete guide on QFT in curved spaces, and conduct us through a generalized quantization procedure based on that of QFT in Minkowski spacetime. Nevertheless, several researchers developed most of it in a category of spacetimes whose causal structure is thoroughly well-defined, namely globally hyperbolic spacetimes. Indeed, as we discussed previously if a spacetime does not feature global hyperbolicity, then basic field equations might not have causal solutions, which jeopardizes the quantization of fields. On what follows, we use works by Wald [4] and Ishibashi [5,6] to prescribe the appropriate dynamics of scalar fields in non-globally hyperbolic spacetimes.
Let us consider a static spacetime (M, g μν ), which admits the following decomposition of its metric [12] In Eq. (3), h i j is the metric induced on a spacelike hypersurface Σ orthogonal to a given timelike Killing field τ μ of the metric, and we define V 2 = −τ μ τ μ . In this particular case, Klein-Gordon equation, reduces to in which A := −V D i (V D i φ) + m 2 V 2 + ξ RV 2 is the spatial component of the wave operator, and D i is the covariant derivative in a spatial slice of Σ.
Since Σ is a subset of the spacetime, singular or boundary points cannot be part of Σ. However, it is clear that these points may be in causal contact with some event in M. Wald's strategy [4] relies on the fact that, since we do not know what happens at the boundary of Σ, a conveinent choice for the domain of the operador A is D 0 (A) = C ∞ 0 (Σ), i.e., our "initial" domain is insensitive to the boundary points. By choosing the Hilbert space ) is a symmetric operator; (iii) as a consequence there is at least one positive self-adjoint extension for (A, D 0 (A)). 2 An extensive literature on Functional Analysis (e.g., see [13,14]) discusses the properties of such operators and present a systematic procedure for obtaining their self-adjoint extensions, accredited to Weyl and von Neumann. Let us first define the adjoint operator (A † , D(A † )) by: 2 We should emphasize that the choice D 0 (A) = C ∞ 0 (Σ) makes perfect sense since the classical theory does not say which effects may arise from boundary points. However this is not the only possible option. In fact there are infinitely many dense domains D(A) where A is symmetric which would also fit to our purposes.

D(A
where ·, · is the usual inner product in H. It can be easily checked (simply using integration by parts) that D 0 (A) D(A † ) so that (A, D 0 (A)) is clearly not selfadjoint. However, as mentioned above, (A, D 0 (A)) may be extended to a self-adjoint operator.
In order to find these extensions, let us define the deficiency subspaces of A, denoted N ± ⊂ H, by and the deficiency indices as n ± = dim(N ± ). For n-th order ordinary differential equations, there are three cases to be considered: 1. If n + = n − , then A has no self-adjoint extension.
2. If n + = n − = 0, then A is essentially self-adjoint, and we obtain it by taking the closure,Ā, of A. 3. If n + = n − = n ≥ 1, then infinitely many self-adjoint extensions of A exist.
They are in one-to-one correspondence to the isometries between N + and N − parametrized by an n × n unitary matrix U .
For partial differential equations in which the Hilbert space can be split in the form H = n i H n i (as in the case to be considered in this paper) we must consider each subspace n i individualy.
Certainly, the third case is more complex than the others, and we must follow a method for obtaining the self-adjoint extensions (see [14] for a proper description of it). They are given by A E , with E being a parameter labeling the extension, defined by and for all Φ ∈ D(A E ). This procedure can always be followed to find whether an operator has self-adjoint extensions and identify them, in case they exist. Wald [4] proved that the one parameter family which is a solution of the equation reproduces the solution of Eq. (4) determined by ordinary Cauchy evolution [with initial data in D( It is straightforward to notice that for each extension A E there will be an associated dynamical evolution of Eq. (10). Consequently, the dynamics of the field is not uniquely determined by initial conditions. We identify those non-equivalent solutions as a result of various boundary conditions that one can impose at a region in space, such as a singularity or a boundary [4]. Ishibashi and Wald [5], argue that Eq. (10) is the only one that prescribes a physically sensible dynamics of scalar fields in non-globally hyperbolic static spacetimes. By comparison with the globally hyperbolic case, they establish a set of conditions that determine whether a time evolution is consistent or not, namely: 1. Solutions of the wave equation must be causal; 2. The prescription for dynamics must be invariant under time translation and reflection; 3. There exists a conserved energy functional also respecting time translation and reflection invariance, in agreement with the globally-hyperbolic case; 4. Solutions satisfy a convergence condition, as proposed in [4].
We should emphasize that Wald's prescription gives us a reasonable evolution for the classical field. By reasonable we mean respecting the four requirements listed above. In particular, the conserved energy functional is consistent with the unitarity of the field in the corresponding Hilbert space. This is, of course, a necessary requirement in the construction of a "closed" quantum theory-not loosing or gaining information from spatial infinty.
Nevertheless, there are other ways to evolve and quantize the classical field. In Ref. [3], for example, the authors studied the quantization of a classical scalar field respecting the so called "transparent" boundary condition. In this scheme the AdS spacetime becomes transparent when mapped into the Einsteins static universe. In this way, the energy is not time independent (there is an exergy flux due to transparency). Besides, the resulting quantization is not "closed" in the sense that the one particle state becomes a mixture of particle states of the closed AdS (obtained through the reflective Dirichlet and Neumann boundary conditons). We therefore choose Wald's prescrition since it leads to a very concise quantization. We must, however, stress out that other works address quantum field theory in spacetimes with timelike boundaries differently, e.g., see [15][16][17].

Boundary conditions at infinity of anti-de Sitter
Let us now consider Klein-Gordon equation (5) in AdS n , as follows where m ξ is the effective mass of the field defined by m 2 ξ = m 2 − ξ n(n − 1)H −2 , and is the Laplace-Beltrami operator on the unit (n − 2)-sphere whose eigenfunctions are Generalized Spherical Harmonic functions, Y l (θ j , φ), with eigenvalues l(l +n−3). We may recall that a static slice of AdS n can be decomposed into a real interval [0, π/2), labeled by the radial coordinate ρ, and an (n − 2)-dimensional unit sphere S n−2 , parametrized by the angular coordinates θ j and ϕ. It is also worth pointing out that, as the spacetime is static, there exists a timelike Killing field ∂ τ , whose eigenfunctions e −iωτ with positive energy, ω > 0, can be used to expand the solution φ. Thus, φ will be an eigenfunction of the quadratic operator ∂ 2 τ with eigenvalue −ω 2 . With those considerations in hand, let us write the solution as Under the transformationf and omitting temporal and angular dependence, Eq. (12) reduces to which is a differential operator whose domain is C ∞ 0 (0, π/2) defined over a Hilbert space H = L 2 ([0, π/2], dρ), and the coefficients of the equation are defined as and From Eq. (19), it is straightforward to check that The coefficient ν is taken to be the positive square root of ν 2 and will depend on the mass and coupling factor of the field. In such conditions, there are four relevant cases to be analyzed, namely (i) ν 2 ≥ 1: in this case, the effective mass of the field satisfies the relation H 2 m 2 ξ ≥ −(n + 1)(n − 3)/4, which comprise the minimally coupled, massless scalar field for n ≥ 3.
In [6], the authors examine the positivity of the operator A in terms of ν. They demonstrate that, in all cases in which ν 2 ≥ 0-i.e., in (i), (ii) and (iii)-A is a positive operator. Meanwhile, in case (iv), the operator is unbounded bellow. Consequently, A has no positive, self-adjoint extensions in case (iv). On the other hand, at least one self-adjoint extension to A exists-that is, the Friedrichs extension [13]-in all other cases: (i), (ii) and (iii).
According to Eq. (7), to construct the deficiency subspaces N ± , we must take ω 2 = ±i, so ω ∈ C. In the case (i), i.e., ν ≥ 1, Eq. (17) has no solution in the Hilbert space so that the deficiency indices are null. In this case and the operator admits a unique self-adjoint extension for each value of l. In other words, the repulsive effective potential in A, i.e., (cos ρ) −2 , prevents the fields from reaching spatial infinity. Hence, they vanish there, and no additional boundary conditions are required. Conversely, in cases (ii) and (iii), the square integrable solution to Eq. (16) with ω = ±i is given by The other linear independent solution is not square-integrable, i.e. it does not belong to the Hilbert space, so we neglect it here. The deficiency indices in these cases are n + = n − = 1, so infinitely many positive self-adjoint extensions of A exist. Now, the effective potential is not as strong as in case (i); hence we may associate the extensions to boundary conditions prescribed at infinity. A one-parameter family of self-adjoint extensions, A β , of A exists for 0 ≤ ν 2 < 1 [cases (ii) and (iii)]. Equation (8) provides us with the appropriate domain of A β . Since the domain of A consists of functions in C ∞ 0 , all additional information needed to prescribe a physically consistent dynamical evolution must come from the asymptotic behavior of f + and U f + , for all isometries U .
Let U β denote the isometries between N + and N − , given by for β ∈ (−π, π]. Let us consider the function whose behavior near infinity (ρ = π/2) dictates the boundary conditions satisfied by all solutions φ t of the form (10). For 0 < ν < 1, the asymptotic behavior at ρ = π/2 is where the coefficients of the leading terms, a ν and b ν , are functions of ν, σ , the spacetime dimension n and the parameter β. The leading powers in ρ of f + are from which we can see that the asymptotic boundary condition depends on the ratio a ν /b ν , which may take any real value. For ν = 0, we have and an analogous procedure reveals that the asymptotic boundary condition depends on a 0 /b 0 also in this case. However, the function (sin ρ) −σ −1/2 · (cos ρ) −1/2 · f β and its first derivative in ρ both scale with a 0 when approaching infinity ρ = π/2. Setting a 0 = 0, we recover Dirichlet and Neumann boundary condition imposed simultaneously, which is precisely Friedrichs extension.
On what follows, we shall denote the ratio a ν /b ν by α ν , hence all self-adjoint extensions of the operator will be parametrized by α instead of β, although α ≡ α(β). From Eq. (25), we can check that 4 which we identify as generalized Robin boundary conditions for 0 < ν < 1. One recovers generalized Dirichlet or Neumann boundary conditions by setting α ν equals to 0 and ±∞, respectively. In the particular case ν = 1/2, Eq. (25) reduces to an even simpler form of the boundary conditions given by 5 which is the usual Robin boundary condition, hence mixing Dirichlet (α = 0) and Neumann (α = ±∞) conditions. Even though the extensions A α are now parametrized by a real parameter α ν , not all of them are positive. Except for ν 2 ≥ 1, whose unique self-adjoint extension is already positive, the remaining cases satisfy the positivity conditions shown in [6]: For where γ is the Euler gamma and ψ is the digamma function.
It is worth pointing out that Eqs. (27) and (28) must be satisfied mode by mode, i.e., for each spherical label l-and for each σ , indirectly (see Eq. 19)-, the conditions are satisfied by f β,ω,l . Accordingly, there are infinitely many parameters α ν,l associated to each f β,ω,l , and they all satisfy different positivity conditions, given in Eqs. (29) and (30).

Green's functions in AdS
In [20], Allen and Jacobson show that, in a maximally symmetric spacetime, two-point functions such as G F (x, x ) = −i ψ|T {φ(x)φ(x )}|ψ , where |ψ is a maximally symmetric state, may be written in terms of the geodetic interval s(x, x ), 6 i.e., Their proposition simplifies the computations considerably since the wave equation becomes an ODE of the variable s. They also require that the Green's function falls off as fast as possible at spatial infinity, which in AdS translates into: G F → 0 as s → ∞. In other words, they are choosing Dirichlet boundary condition for the field φ. Kent and Winstanley [2], exploit this simplicity to find the fluctuation of the field squared and the energy-momentum tensor in all spacetime dimensions of AdS. They 5 In case ν = 1/2, we drop the index of α 1/2 and replace it simply by α. 6 In AdS, s is constructed so that it goes to zero as x → x and goes to infinity as we approach the boundary. also verify that their results are compatible with the ones of Burgess and Lütken, whose approach in [21] was to perform a summation of modes of the wave solutions.
We are not aware of any law of nature that restricts the boundary conditions of all modes to Dirichlet ones. Indeed, Ishibashi and Wald showed in [5] that there is an entire category of boundary conditions that prescribe a physically consistent dynamical evolution. Additionally, there is no guarantee that all modes must satisfy the same boundary condition.
Let us then consider a setup in which one of the modes of the wave equation, u ω α ,l α , is chosen so that its radial component f ω α ,l α (ρ) satisfies a generalized Robin boundary condition with parameter α. Meanwhile, the components f ω,l (ρ) of all other modes u ω,l (x) (l = l α ) satisfy Dirichlet boundary conditions. The Green's function in this case is given by mode sum (from now on, we consider τ > τ ) where N ω,l are normalization constants. We may complete the last term in the summation for all Dirichlet modes by adding them to and subtracting them off Eq. (32), i.e., Let us denote the last term in Eq. (33) by G may not be a maximally symmetric function. It seems reasonable for us to write that Equation (34) illustrates the break of AdS invariance of the Green's function, as it may not depend on the geodetic interval s entirely anymore. We attribute the break on the maximal symmetry of G F to the imposition of different boundary conditions for each angular mode.
Given the decomposition (34), one question arises. 7 Is it compatible with quantum field theory by means of the canonical commutation relations where x a = (ρ, θ, ϕ)? The answer is yes, and to prove this, consider the Green's function , x ). Without loss of generality, let us consider a nontrivial boundary condition in the l = 0 mode only. This will also be our choicefor simplicity-later in the paper.
Within our notation, it is given by The first term in the above equation is the zero-th mode corresponding to an arbitrary Robin boundary condition, the second one is zero-th mode corresponding to Dirichlet boundary condition and the third one is the complete Green's function corresponding Dirichlet boundary condition (the one respecting the spacetime symmetries). The first and second term have the form Now, given the normalization of u ω (ρ) and the fact that they form a complete set of a self-adjoint operator (corresponding to l = 0) we have In this way since G Dirichlet is obviously consistent with the commutation relations. This proves that our complete solution is also consistent with the commutation relations.
Another important question is whether the Green's function has the Hadamard form. First of all, since we are dealing with a non-globally hyperbolic spacetime, we are interested in the local Hadamard form only. Global Hadamard form is not equivalent in this case since we can face reflecting singularities of the Green's function [16]. Besides, by choosing the positive self-adjoint extensions and constructing the Green's function with the positive frequency modes with respect to the Killing field ∂/∂τ , we are constructing a Green's function of a ground state. Using Appendix A of Ref. [22] (where it is proved that such construction always leads to a ground state with the correct Hadamard form), we can state that our Green's function has the expected Hadamard form.

Renormalized quantities for a conformal massless scalar field in AdS 4
In order to shed light on what we have discussed so far, we shall specialize to four spacetime dimensions, AdS 4 . For simplicity on the computation of quantities of interest, let us restrict ourselves to a conformally invariant, massless scalar field, φ, i.e., m = 0 and ξ = 1 6 . In this case, from Eq. (18), we get ν = 1/2, and from Eq. (19), we find that σ = (2l + 1)/2. Equation (16) becomes and its solutions are where C 1 and C 2 are constants to be determined, and P and Q are the associated Legendre functions of the First and Second kinds, respectively. Square integrability requires f to fall off at the origin ρ = 0, hence C 1 → 0. 8 A complete set of eigenfunctions is then for normalization constants N ω,l to be determined. As discussed in Sect. 4, boundary conditions at infinity are necessary to prescribe the dynamical evolution of the field in AdS n . In case ν = 1/2, Robin boundary conditions (28) are the appropriate ones. We aim to provide an example of the setups discussed in the last section. For that, we will consider that all non-spherically symmetric modes respect Dirichlet boundary conditions. However, the l = 0 mode will be chosen to satisfy Robin condition with a parameter α. As discussed above, the vacuum will not be AdS invariant in this case. However, since the non-trivial boundary condition is on l = 0 mode, we still preserve spherical symmetry. Formulas 14.5.3 and 14.5.4 in Ref. [23] allow us to describe the behavior of f ω,l and its derivative at the boundary, as follows For l > 0, all modes satisfy f ω,l (ρ → π/2) = 0 (Dirichlet boundary condition), thus the quantized frequencies are given by ω = 0, ±2, ±4, . . . , ±l, ±(l + 2), ±(l + 4), ±(l + 6), . . . .
These frequencies come from the relations (clearly, for l = 0, the quantized frequencies for Dirichlet boundary condition are given by the set 2Z). For l = 0, we calculate the ratio between derivative (44) and function (43) to use it in (28), i.e., Positivity condition (29) requires that In our analysis, we consider α ≥ 0, which includes Dirichlet, α = 0, and Neumann, α → ∞, cases. Equation (47) imposes a quantization condition for the frequencies ω in terms of the parameter α. Except for α = 0 and α = ∞, it cannot be solved analytically for an arbitrary value α. One can readily verify that, in the Neumann case (α → ∞), the frequencies are odd integers. Meanwhile, for Dirichlet, they are even integers, which is consistent with Eq. (45).
In our procedure, we employed the software Mathematica [24] to solve Eq. (47) numerically in a determined range of ω for several values of α. As shown in Fig. 1, the solutions of (47) are given by the intersection points between the two functions. We can see that ω values for arbitrary α always lie between an odd number and its next even integer, which are precisely the frequencies for Neumann and Dirichlet conditions, respectively. Thus, given a Neumann frequency, ω N ,r = 2r − 1, and a Dirichlet one, ω D,r = 2r , for r > 0, we may denote an α frequency between them as ω α,r , even though it is not an integer number.

Quadratic field fluctuations 2
Before computing the Green's function, it is useful to write solution f ω α,r ,0 in a more convenient form and normalize it accordingly. Using Ref. [23], we find 9 Now, we recall our discussion from last section to construct the appropriate Green's function. We can decompose our Green's functions in two parts, i.e., 2 ω α,r π − sin(ω α,r π) sin(ω α,r ρ) sin(ω α,r ρ )e −iω α,r (τ −τ ) and where ω ≥ 0 means that we are summing over the positive frequencies in Eq. (45). Our 'Dirichlet' Green's function (51) s(x, x )). As Kent and Winstanley show in [2], approaching the coincidence limit s → 0, the function G (D) F diverges according to the Hadamard form. Thus, point-splitting renormalization can be employed to compute finite quantities. Furthermore, they obtain the Hadamard forms in AdS for any spacetime dimension through a systematic method, based on [25].
In the particular case of AdS 4 , for a conformally invariant field, the Green's function G (D) F has the Hadamard form given by After renormalization, it may be written as [2] G (D) F ren We may find the expectation value of the quadratic field fluctuations as follows which is naturally in accordance with the results in Ref. [2]. Analogously, the effect of our Green's function (50) on φ 2 appears when taking the coincidence limit x → x. However, calculating G (α) F analytically is impossible, since the summation is taken over numerical values of frequencies. Hence, we adopt a numerical approach to find our results.
We expect G F . On the other hand, we cannot perform the infinite sum in (50) numerically, so a residual divergent behavior might appear. Through our computations, we noted it was convenient to take the coincidence limit in the radial coordinate first, i.e., ρ → ρ, and then in the time coordinate. Thus, our final step would be to take the limit of τ → τ . It is more convenient though, to analytically extend the function on the complex plane and take the limit through the imaginary axis, i.e., τ → τ + i , hence τ − τ → −i . Finally, by multiplying G (α) F by i, we will have an entirely real-valued function that, in the limit → 0, yields directly the quadratic fluctuations of the field, and it is much simpler for us to handle it numerically.
Before implementing the numerical routine, we considered the only case that can be treated analytically, which is the Neumann condition, α → ∞. In this situation, the frequencies are ω ∞,r = 2r − 1, for r > 0, and the Green's function reduces to the following summation which we calculated using Mathematica [24], resulting It is straightforward to find the expectation value φ 2 (N ) by simply taking → 0, i.e., First of all, notice that φ 2 (N ) (ρ) fails to be a constant so that the spacetime is not homegenous anymore-at least for the quantum field. Therefore, there is a preferred frame-the one in which Neumman boundary condition has been chosen for tha particular mode l = 0. What if we had chosen Neumann boundary condition in every l mode, instead? Since this particular boundary condition does not introduces an energy scale, we expect that the relevant physical quantities will respect the symmetries of the spacetime and φ 2 (ρ) will be constant in this case (see Ref. [3], where it was shown that the Green's function respects the spacetime symmetries for reflective Neumann boundary condition). The function φ 2 (N ) is finite because both terms inside the sum in Eq. (55) diverge with same strength. Naturally, their subtraction eliminates the infinities. In particular, the last term in Eq. (55), the Dirichlet counterpart of G (α) F , denoted G (α,D) F , appears for all values of α and dictates the divergent behavior at → 0. We find its form by calculating the infinite summation and expanding it in powers of , i.e., Our numerical approach to find the expectation value φ 2 (α) proceeded as follows: 1. Given a value for α, solve Eq. (47) to find the frequencies ω α,r up to r max = 5000; 2. Given a value of ρ between 0 and π/2, compute numerically the truncated summation for 50 values of equally spaced in the range 0.002 to 0.1. 10 3. Fit the function f (α) ρ [ ] using a model that reproduces the divergent behavior in Eq. (58) followed by a Taylor expansion up to order 2 , i.e., As G (α) F is a finite quantity, we expect the divergent behavior of f (α) ρ [ ] to be extremely attenuated. We have found coefficients b ranging between 10 −9 and 10 −12 , recovering the expected almost-finite behavior. The coefficients c and d were effective on reducing the residuals of the fit. Finally, a gives the approximated finite numerical value of φ 2 (α) at the point ρ. 4. Repeat steps 2 and 3 for as many values of ρ between 0 and π/2 as desired. 5. Repeat the entire procedure for another value of α. We followed the scheme described above for 14 values for the parameter α. We chose 80 equally spaced points in the range (0, π/2) to obtain a good resolution of the behavior of φ 2 (α) (ρ). Our results are plotted in Fig. 2. The curve corresponding to α = 1000 reproduces almost perfectly the analytic Neumann result (57). Accordingly, as we approach the other extreme, α = 0-corresponding to Dirichlet conditionswe can see the curves getting closer to zero. Consistently, if α = 0, then G (α) F indeed vanishes, as one can see from Eq. (50).
From Fig. 2, we notice that the correction to φ 2 arising from the nontrivial boundary condition becomes more relevant close to ρ = 0. This may seem strange, since the boundary condition is necessary at the conformal boundary. However, since any observer (close to ρ = 0 or not) is infinitely far away from the boundary, this does not seem to be a contradiction.

Energy-momentum tensor fluctuations T (˛)
In [2], the authors obtain the renormalized energy-momentum tensor T μ ν ren in AdS n . They use the formula from Ref. [25] T μν ren = − where and Θ μν is a purely geometric tensor constructed to be conserved. Kent and Winstanley find that the non-geometrical component of the tensor is proportional to the metric tensor, which is completely consistent with the maximal symmetry of AdS. In our particular case of a conformally invariant field in four spacetime dimensions, we have and the geometric tensor Θ μν is identically zero. We may obtain this renormalized expectation value from Green's function G (D) F ren , hence is associated with Dirichlet conditions in all modes of the wave equation.
Here, we want the contributions to the energy-momentum tensor coming from G (α) F . Our approach will be analogous to that of the Green's functions: we decompose T μ ν ren into two parts, one carrying the boundary condition, denoted T μ ν (α) ren with 16 components T μ ν , and another one reproducing the Dirichlet results as in Eq. (64). In our case, Eqs. (62) and (63) may be written as and from which it follows that [G](ρ) ≡ φ 2 (α) (ρ). According to formula (61), we have here ren , we readily verify that the only non-vanishing components are diagonal terms and the term T τρ (= T ρτ ). Let us recall the temporal inversion (τ → −τ ) symmetry of AdS, denoted I, given in four dimensions by the transformation matrix I μ μ = diag(−1, 1, 1, 1). As none of our quantities depend explicitly on τ , we expect this discrete symmetry to be preserved. In particular, we expect T τ x j = T −τ x j = T τ x j , for x j = (ρ, θ, ϕ). On the other hand, T μ ν (α) ren transforms as a tensor, so we have That yields T τρ = −T τρ , which then implies T τρ = T ρτ ≡ 0. At this point, we have a diagonal tensor, whose remaining components may be calculated using Eq. (67). Our computational efforts were not successful when trying to compute the numerical expressions directly. However, we came up with a solution based on some properties that T which is our first constrain on the remaining diagonal components. We may use the symmetries of AdS as well. Although our Green's function breaks AdS invariance of the radial coordinate ρ, all other symmetries should remain valid. In AdS 4 there exist 10 Killing fields corresponding to the following isometries: one temporal translation, three rotations, four boosts and four spatial translations. From which, we only expect the first two to be preserved after imposing Robin boundary conditions in only one of the modes. The temporal Killing field, t = ∂ τ , yields a conservation equation along with its flow, given by the Lie derivative of the tensor with respect to t, i.e., which shows that all components of T μ ν (α) ren are independent of τ . Additionally, we have the generators of spherical symmetry, given by the following Killing fields Since a combination of them is still a Killing field, we may use χ 2 and χ 3 to obtain χ 4 = ∂ θ . We can use χ 1 and χ 4 to find other two conservation equations similar to that of t, as follows and q(ρ) = 2 (sin(2ρ)F (ρ) + (7 − cos(2ρ))F(ρ)) (5 − cos(2ρ)) cot ρ .
One can verify that solves the equation. In our case, we have which vanishes at ρ = 0 and diverges at ρ = π/2. Naturally, the inverse function exp − dρ p vanishes at the boundary, but diverges at ρ = 0 with strength 1/ρ. As it is physically reasonable to ask for a finite T τ τ at ρ = 0, we set C to zero. Finally, we compute T τ τ using the following expression For the Neumann case, we used our previous analytic results and found F to be F(ρ) = cot 2 ρ 48π 2 csc 2 ρ + 2 + 2(csc 4 ρ − 1) log(sec ρ) .
Applying it in Eq. (91), and then using (84) and (83), we find Now, we have a result to compare our numerical ones with. To compute the function F numerically, we used our previous results of φ 2 (α) (= [G]), but we also need [G] τ τ . According to (63), we find it by taking the second derivative of G (α) F (τ, τ , ρ, ρ) with respect to τ and, then, taking the coincidence limit. In the convention we adopted, ∂ τ τ = −∂ . Indeed, we find it to be Our numerical procedure to find the expectation value of the energy-momentum tensor fluctuations was: 1. Given a value for α, use the frequencies ω α,r found before; 2. Given a value of ρ between 0 and π/2, compute numerically the truncated summation for 50 values of equally spaced in the range 0.002 to 0.1.

Fit the function F (α)
ρ [ ] using a model that reproduces the divergent behavior followed by a Taylor expansion up to order 2 , i.e., As expected, the divergent behavior of F (α) ρ [ ] is extremely attenuated, and the coefficient b is negligible compared to the others. Again, the coefficients c and d were effective on reducing the residuals of the fit. Finally, a gives the finite approximated numerical value of [G] τ τ at the point ρ. 4. Repeat steps 2 and 3 for as many values of ρ between 0 and π/2 as desired to obtain the complete [G] τ τ (ρ). 5. Use our previous results for [G] together with [G] τ τ in Eq. (83) to find a numerical interpolation of F(ρ), denoted F[ρ]. 6. Given a value of ρ between 0 and π/2, use F[ρ] in Eq. (91) and perform a numerical integration to obtain an approximate value of T τ τ at that specific ρ. 7. Repeat step 6 for several values of ρ to find a complete numerical function T τ τ . With that in hands, compute T ρ ρ and T θ θ using Eqs. (84) and (83). 8. Repeat the entire procedure for a different value of α.
Similarly to our results for the expectation value of the field squared, we followed the numerical procedure for 14 values of α. We have found all components of T μ ν (α) ren . In Fig. 3, we can see T τ τ for several values of α, it is clear that the form of the function follows the analytic result for the Neumann condition (plotted in gray). For all the remaining components of the tensor, the results were the same: by increasing α, they approach the expected analytic behavior of Eq. (93) (Neumann condition).

Discussion and further remarks
Avis, Isham, and Storey took a first-step, in Ref. [3], towards the development of a quantum field theory in anti-de Sitter spacetime. They acknowledged that the conformal infinity poses a serious causality issue to the wave equation but solve it by regulating the information flow through the boundary 'by hand.' They imposed the Fig. 3 Contribution to the expectation value T τ τ due to Robin boundary conditions at infinity for the spherically symmetric mode. H is set to one so-called 'transparent' and 'reflective' boundary conditions at infinity in analogy to a box in Minkowski spacetime. In this way, they quantized the fields in the Einstein Static Universe and restricted it to the AdS later.
Conversely, in this article, we considered the developments made by Ishibashi and Wald in [5], where they propose a physically consistent prescription for the dynamical evolution of fields. In the particular case that we have considered, they show that the imposition of mixed boundary conditions at the spatial infinity is sufficient to determine the evolution of quantum fields uniquely.
In the setup studied by Kent and Winstanley [2], all angular modes of the wave equation satisfy the same Dirichlet boundary condition at infinity. Their results are consistent with the maximal symmetry of AdS. Hence, the expectation values of field-dependent quantities fluctuate in the same way throughout spacetime, i.e., they are coordinate-independent. In the light of Wald's and Ishibashi's developments, we presented a setup here that puts up to question how necessary it is to impose the same boundary conditions to all modes of the wave equation. Indeed, we are not aware of any requirement of nature that precludes us from considering various setups in terms of boundary conditions. Our analysis indicated a violation of AdS invariance in the Green's functions, which carried out implications on the related quantities: the quadratic fluctuations of the field and the energy-momentum tensor. Both of them are now dependent on the radial coordinate for any values of the parameter α, as shown in Figs. 2 and 3. At this stage, any attempt of obtaining a back-reacted metric using Einstein's semi-classical equations, G μν = 8π G T μν ren ≡ 8π G T μν (D) Fig. 4 Energy densitiy − T τ τ ren of a massless scalar field conformally coupled to AdS 4 for several Robin boundary conditions. H is set to one would not yield a maximally symmetric metric anymore, but a spherically symmetric one. In these conditions, the coordinate system used to define the angular modes of the wave equation will be privileged. In particular, in this system, the energy density reaches its minimum at the origin ρ = 0, as shown in Fig. 4.
In Fig. 4, we can see a clear violation of the weak energy condition in most of the spacetime, except close to the boundary, where the Dirichlet contribution, − T τ τ (D) ren , pushes the energy density back up over zero. Even though such violation is no stranger to us-as can be observed in the Casimir effect-it appeared as a consequence of the contribution from the Robin boundary condition exclusively. Indeed, the Dirichlet term, − T τ τ (D) ren , of the energy density is positive throughout the entire spacetime. Thus, it is safe to assert that the violation of the weak energy condition is a direct consequence of the imposition of non-Dirichlet boundary conditions at infinity.