Hidden freedom in the mode expansion on static spacetimes

We review the construction of ground states focusing on a real scalar field whose dynamics is ruled by the Klein-Gordon equation on a large class of static spacetimes. As in the analysis of the classical equations of motion, when enough isometries are present, via a mode expansion the construction of two-point correlation functions boils down to solving a second order, ordinary differential equation on an interval of the real line. Using the language of Sturm-Liouville theory, most compelling is the scenario when one endpoint of such interval is classified as a limit circle, as it often happens when one is working on globally hyperbolic spacetimes with a timelike boundary. In this case, beyond initial data, one needs to specify a boundary condition both to have a well-defined classical dynamics and to select a corresponding ground state. Here, we take into account boundary conditions of Robin type by using well-known results from Sturm-Liouville theory, but we go beyond the existing literature by exploring an unnoticed freedom that emerges from the intrinsic arbitrariness of secondary solutions at a limit circle endpoint. Accordingly, we show that infinitely many one-parameter families of sensible dynamics are admissible. In other words, we emphasize that physical constraints guaranteeing the construction of full-fledged ground states do not, in general, fix one such state unambiguously. In addition, we provide, in full detail, an example on $(1 + 1)$-half Minkowski spacetime to spell out the rationale in a specific scenario where analytic formulae can be obtained.


Introduction
Quantum field theory on curved spacetimes has lead to significant improvements in our understanding of different physical phenomena ranging from particle production in cosmology to Hawking radiation in black hole physics.In the analysis of the vast majority of the available models, the first step consists in constructing full-fledged quantum states for free fields.Under the mild assumption that the correlation functions are Gaussian, it reduces to the identification of an on-shell two-point correlation function that has to abide to physically motivated constraints.The prime example in this direction is the Hadamard condition, which ensures not only that the quantum fluctuations of all observables are finite, but also that Wick ordered fields can be constructed following a locally covariant scheme.In turn, it entails control of the underlying renormalization group and of interactions that are studied at a perturbative level.Yet, in many concrete scenarios one is limited to abstractly argue the existence of such distinguished two-point functions and an explicit construction is, at best, elusive.
Major improvements occur when one concentrates on static spacetimes M R×Σ, regardless of whether they are globally hyperbolic or not.Let us consider, for simplicity, a free, scalar field Ψ : M → R that abides to the Klein-Gordon equation.By calling t the time coordinate along R, the latter simplifies to where K is an elliptic, second order partial differential operator.The key rationale consists of reading K as a symmetric operator on the Hilbert space H := L 2 (Σ, dµ Σ ) of square-integrable functions with respect to the measure dµ Σ induced by the Lorentzian metric tensor of M on Σ.This leads to two notable advantages, one at a classical and one at a quantum level, as described next.At a classical level, solutions of Equation ( 1) can be constructed as follows.Since K is a real and symmetric operator, it admits a non-necessarily unique self-adjoint extension K. Assuming, for convenience, that K has positive spectrum, and given initial data (Ψ 0 , Ψ0 ) ∈ C ∞ 0 (Σ) × C ∞ 0 (Σ) ∩ D( K) × D( K), for each t ∈ R we have where each term is well-defined using spectral calculus.Moreover, there exists a unique Ψ ∈ C ∞ (M) such that, Ψ| Σ t = Ψ t and ∇ n Ψ| Σ t = Ψt , where Σ t ≡ {t} × Σ, t ∈ R, while n is the unit vector field normal to Σ t .If (M, g) is a globally hyperbolic spacetime without boundary, than there exists a unique choice for K and the dynamics is therefore unambiguously determined [1].On the contrary, if K has more than one self-adjoint extension, then multiple, physically inequivalent scenarios do exist.Markedly, the latter is not a remote possibility: it occurs for example when (M, g) is a globally hyperbolic spacetime with a timelike boundary, see [2].This class of backgrounds encompasses several physically interesting scenarios such as AdS spacetime, which was first analyzed in the language above by Ishibashi and Wald in [3].It must be stressed that, using the language of boundary triples [4], the infinite set of different choices of self-adjoint extensions for K can be put in correspondence with the choice of a boundary condition for Equation (1).
At a quantum level, the assumption that (M, g) is static guarantees a considerable advantage: the existence of a ground state.Under the same premises of the previous paragraphs, the associated two-point function ψ 2 (t, x, t , x ) can be constructed directly as the integral kernel of the operator [5] Observe that there exists a different ground state for each self-adjoint extension K, in accordance with the fact that each K characterizes a different physical system.Furthermore, if the underlying background is a globally hyperbolic spacetime with or without timelike boundary, then ground states are of Hadamard form as a consequence of the results of [6,7].While at this stage the analysis of a scalar quantum field theory on a static spacetime seems a rather well-understood problem, the drawback lies in two crucial details.On the one hand, when non unique, an explicit construction and characterization of all self-adjoint extensions of K in Equation ( 1) is a daunting task.On the other hand, the quantitative evaluation of physical observables in concrete scenarios, such as on black hole spacetimes, requires a deeper and more hands on knowledge of the two-point function, far beyond the spectral level as per Equation (3).To bypass this conceptual hurdle, it is customary to consider static backgrounds with a high degree of symmetry.Beyond reasons of mathematical simplicity, this class includes many physically relevant backgrounds, such as cosmic strings, black holes and asymptotically AdS spacetimes.
In this paper, we consider the class of n-dimensional static spacetimes M that are isometric either to where I ⊆ R, and Σ n−2 j are Cauchy-complete, connected, (n − 2)-dimensional Riemannian manifolds of constant sectional curvature j.The line element associated to the metric tensor on M reads where f and h are suitable positive functions.Barring some technical aspects that will be specified in the next sections, we emphasize that a large class of spacetimes is characterized by the line-element above, including black hole backgrounds ranging from the three-dimensional static BTZ spacetime to the n-dimensional Schwarzschild or Schwarzschild-AdS spacetime.
On top of these manifolds, we consider a real, scalar field Ψ whose dynamics is ruled by the Klein-Gordon equation, which, as before, can be written as per Equation (1).With the construction of a quantum field theoretical framework in mind, we are interested in obtaining distinguished two-point functions that correspond to full-fledged ground states.Although the procedure outlined above is applicable, especially when considering scenarios where boundary conditions needs to be imposed, it is common to follow a more computationally oriented approach that exploits the underlying symmetries.In the following, we sketch the steps usually followed in the literature in these scenarios.More details will be given in the next sections of this work.
Consider a solution of the Klein-Gordon equation on (M, g) assuming that it admits a mode expansion: where Y j (ϕ 1 , ..., ϕ n−2 ) are the eigenfunctions of the Laplace operator on Σ n−2 j whose corresponding eigenvalue is denoted by η j , while ω ∈ R plays the standard rôle of frequency.
The only unknown function R ωη j can be shown to satisfy an eigenvalue equation AR ωη j = λR ωη j where A is a second order differential operator in the radial coordinate r whose domain is the interval I, here taken for definiteness as (a, b).Most notably A can be written in the form of a, possibly singular, Sturm-Liouville operator, see [8].
Similarly to the rôle played by K in Equation (1), one reads A as a symmetric operator on a space of square-integrable functions over the interval I = (a, b).[9], three options are possible.

Following von Neumann's theory of deficiency indices
A can admit just one self-adjoint extension, a one-parameter or a two-parameters family of self-adjoint extensions.In most applications, the last option does not occur.
At this stage it is necessary to pause the description of the procedure to construct a two-point function and draw the attention to self-adjoint extensions.More precisely, from the viewpoint of the differential equation AR ωη j = λR ωη j , the existence of such extensions can be inferred by looking at the behavior of solutions close to a and b, the endpoints of the interval I. Henceforth, for definiteness, we focus on a.The general theory of Sturm-Liouville operators guarantees that, for any λ ∈ C, there is always a distinguished function u, called principal solution.In Section 3 we dwell on the technical details of this concept.For now it suffices to say that u tends to zero as r → a + faster than any other solution that is linearly independent from it, and it is square-integrable in any neighborhood of the endpoint a. Notwithstanding, the existence of another solution, called secondary solution, that is linearly independent from u and square-integrable in any neighborhood of a depends on the differential problem at hand.If such a solution does not exist at both endpoints, then it happens that there is a unique self-adjoint extension for A, see [8].More interesting is the scenario for which that is not the case for some value of λ ∈ C, on account of the fact that if a secondary solution exists, then it is highly non unique.
The intrinsic arbitrariness of the secondary solution lies at the core of this work.Suppose there exists a secondary solution at a, but only the principal one at b. From the viewpoint of the Sturm-Liouville operator, the choice of a specific secondary solution at a is irrelevant as its rôle lies only in establishing a one-to-one correspondence between self-adjoint extensions of A and boundary conditions of Robin type, assigned at the endpoint a.These exhaust all possibilities at the level of the ordinary differential equation and therefore one can read the choice of two different secondary solutions as two different, albeit equivalent, ways to span the same space of boundary conditions and, consequently, of solutions of the underlying ordinary differential equation.In most of the recent literature in quantum field theory on curved spacetimes, a consequential dogma has always been to consider the secondary solution as an innocuous and physically insignificant abstraction.Our main goal is to argue that this by far not the case since, tracking back the problem at a fully covariant level, making different choices of the secondary solution, while keeping Robin boundary conditions, impacts significantly the analysis of scalar field Ψ.More precisely one is able to codify at a covariant level a much larger class of admissible boundary conditions than just a one-parameter family as the standard analysis might suggest.
The physical relevance of choosing a secondary solution is limpid when one constructs the two-point correlation function of the underlying ground state.Let us thus focus once more on the procedure we started sketching at Page 4 and let us state the subsequent steps, as follows.
Since the two-point correlation function ψ 2 must obey the Klein-Gordon equation it is suitable to use the same mode expansion as for the construction of the solutions Ψ.This, together with the ansatz that only positive frequencies contribution are of relevance, identifies a full-fledged ground state.On account of the large isometry group of the background, ψ 2 is completely determined up to a kernel along the radial direction.
As explained in Section 5, such kernel can be constructed using an algorithmic scheme once a self-adjoint extension for the operator A has been chosen in terms of a Robin boundary conditions imposed at the level of principal and secondary solutions.
We emphasize that the procedurehas been extensively applied on static spacetimes with a timelike boundary in the last years.To mention a few, exhaustive works based on Von Neumann deficiency index theory are [11,12], on AdS spacetimes, and [13], on a static BTZ black hole.Beyond these, analyses based on a mode expansion and on Robin boundary conditions and used in the construction of physically-sensible two-point functions within quantum field theory on asymptotically AdS spacetimes can be found in [14,15,16,17,18], and, more recently, also in [19,20,21,22,23].
The scope of this work is to highlight the fundamental rôle played, in this whole procedure, by the secondary solutions of the underlying Sturm-Liouville problem.To strengthen this statement we consider a simple, yet illustrative example, namely the two-dimensional half-Minkowksi spacetime R × R + .The advantage is that, in this scenario, we can address the problem using explicit, analytic formulae that allow to make clear that even minor adjustments to the choice of secondary solution yield, at the fully covariant level, boundary conditions that have a completely different physical interpretation.In a nutshell, we aim to convey that the choice of secondary solution is of physical consequence.
This paper is organized as follows.In Section 2 we show that the Klein-Gordon equation untangles into a Sturm-Liouville problem for the radial part of the Klein-Gordon operator on static spacetimes with maximally symmetric sections.Subsequently, in Section 3 we provide straightforward generalizations of main results from singular Sturm-Liouville theory that allow us to obtain all self-adjoint representations of the latter.Markedly, singular endpoints give rise to an ambiguity in the definition itself of generalized Robin boundary conditions.We clear up this ambiguity in Section 4 by defining generalized (γ, v)−Robin boundary conditions and explaining its connection with the regular case.In Section 4.2 we show how a boundary condition on the radial part translates to a boundary condition on the full solution Ψ. Insofar as solutions of the Klein-Gordon equation characterize a classical dynamics, it is meaningful to include a discussion on the canonical quantization procedure.Hence, in Section 5 we explain the connection between the imposition of the canonical commutation relations and the spectral resolution of the identity given by the before-mentioned Sturm-Liouville problem.We illustrate the main points of this work in a detailed example given in Section 6.Most importantly, this example clarifies in which sense the generalized (γ, v)−Robin boundary conditions imposed on the radial part may render time-dependent boundary conditions on Ψ. Final remarks are given in Section 7.

The Klein-Gordon equation
In this initial section we introduce both the geometric data of the spacetimes we are interested in and the Klein-Gordon equation.In addition we show that, under the specific assumptions on the background metric, the Klein-Gordon equation can be reduced to a Sturm-Liouville problem.
In this work, for n > 2, Σ n−2 j denotes a Cauchy-complete, connected, (n − 2)dimensional Riemannian manifold of constant sectional curvature j, parametrized by (ϕ 1 , ..., ϕ n−2 ), and whose standard metric has an associated line element dΣ n−2 j (ϕ 1 , . . ., ϕ n−2 ).Unless stated otherwise, we shall assume j has been normalized so that j ∈ {−1, 0, +1}.The symbol M refers to an n-dimensional, static spacetime isometric to the warped geometry R × I × Σ n−2 j , where I ⊆ R, while the line element of the underlying metric read in global Schwarzschild-like coordinates (t, r, ϕ 1 , ..., ϕ n−2 ): For convenience, we shall call t ∈ R, r ∈ I ⊆ R, and (ϕ 1 , ..., ϕ n−2 ), respectively, the time, the radial and the angular coordinates.Equation ( 4) is completely specified aside from the two functions f, h.For simplicity, we assumed them to be elements of C ∞ (I; (0, ∞)), although in many instances throughout this work less regularity would suffice.Note that we also allow for the case n = 2, in which M is isometric to R × I endowed with the line element Remark 1.We consider I to be an open interval, say I = (a, b), which might suggest that we are discarding scenarios of notable interest such as globally hyperbolic manifolds with timelike boundary, e.g. the universal cover of AdS n .In these cases, the counterpart of I would include one or both endpoints in the domain of the coordinate r.Yet, if one is interested in the analysis of boundary conditions and their effects, it suffices to focus the attention on the interior of the underlying manifold.Therefore, our analyses can be straightforwardly applied to such cases as well.
On M, we consider a free, scalar field with mass m 0 ≥ 0, Ψ : M → R whose dynamics is ruled by the Klein-Gordon equation: where ξ ∈ R, while R is the scalar curvature built out of the metric as per Equation (4).In the case in hand, the D'Alembert wave operator, denoted by , reads where ∆ Σ n−2 j is the Laplace operator on Σ n−2 j .In the following, we construct the solutions of the Klein-Gordon equation.Although Equation ( 5) can be recast in the form of Equation ( 1), we take a route alined with the proceduredescribed in the Introduction.Assuming that the regularity of Ψ is such that we can work at the level of modes, we consider the ansatz where Y j (ϕ 1 , ..., ϕ n−2 ) are the eigenfunctions of ∆ Σ n−2 j with corresponding eigenvalues denoted by η j .Observe that ∆ Σ n−2 j has a continuous spectrum if j ∈ {−1, 0}, while a discrete one if j = 1.
Equation ( 6) in combination with Equation ( 5) yields that R ωη j obeys to a secondorder, ordinary differential equation, dubbed radial equation: We can rewrite Equation ( 7) as where If λ = ω 2 , the functions p, q and µ are given by In the remaining cases the corresponding expressions can be derived directly from Equation ( 7), but we omit listing them as they are a rather straightforward modification of Equation ( 10a) and (10b).To conclude this section we observe that Equation ( 7) identifies a, possibly singular, Sturm-Liouville problem, following the standard nomenclature of ordinary differential equations, see e.g.[8].

Self-adjoint extensions
Envisioning the construction of ground states for the Klein-Gordon field Ψ, it is essential to obtain the advanced and retarded fundamental solutions associated to the operator P as in Equation (5).To this end we bear in mind a procedure that has been considered in several examples in the literature [16,17,18,19,20,21], mainly when the underlying spacetime possesses a conformal, timelike boundary.The starting point is A, as per Equation ( 8), which shall be read as an operator on the Hilbert space L 2 (I, µ(r)dr).
The operator A, and consequently also L, is manifestly symmetric when taken with the dense domain C ∞ 0 (I).Herein, we scrutinize whether L admits self-adjoint extensions and, if so, how many of them.This is a mathematical question that can be answered by combining tools of Sturm-Liouville theory with the theory of unbounded operators on Hilbert spaces, see e.g.[9].Accordingly, in the following we recall the main results known in the literature that are of relevance to our investigation as well as necessary to make this work self-contained.All definitions and lemmas introduced here culminate in Theorem 3.5, which constitutes the resolution to the question in hand.
First, let us pose the question more precisely.Consider the Sturm-Liouville problem, as in Equation (8), where as of now we omit the subscripts ω, η j from the radial function for decluttering.
In addition, letting L 1 loc refer to locally integrable functions, we assume that 1/p, q, µ ∈ L 1 loc (I), µ > 0, and As alluded to in the previous paragraphs, L as per Equation ( 11) identifies either a minimal or a maximal operator respectively indicated by L min and L max , with corresponding domains where the closure is taken with respect to the graph topology, and AC loc (I) denotes the set of functions that are absolutely continuous on all compact intervals of I. Specifically, our quest is to find self-adjoint extensions L S.A. of L min , whose domain shall be denoted by D S.A. (L S.A. ).If their spectrum satisfies σ(L S.A. ) ⊆ [0, ∞), then we say L S.A. is a positive, self-adjoint extension of L min .
In order to address the quest stated above, we introduce some additional tools tailoring the analysis of [8] to the case of interest, i.e.Equation (11) together with the assumptions of the previous sections.
For y, z ∈ AC loc (I) we denote the Lagrange sesquilinear form and the Wronskian, respectively, by ii) singular if it is not regular; iii) limit circle if all solutions of Equation ( 11) lie in L 2 (I c , µ(r)dr), ∀c ∈ I; iv) limit point if it is not limit circle.
The following definition is especially relevant for our analysis since it differentiates among the solutions of Equation ( 11) depending on their behaviour close to an endpoint.
Definition 3.2.Let y be a non-vanishing solution of Equation (11) in I c , ∀c ∈ I. Then we say y is a i) principal solution at e if, for any other solution z of Equation ( 11) that is linearlyindependent from y, y(x) z(x) ii) secondary (or non-principal) solution at e if it is not a principal solution.
Definitions 3.1 and 3.2 relate by the fact that at a limit point only the principal solution belongs to L 2 (I c , µ(r)dr), ∀c ∈ I. Manifestly, the classification given by Definition 3.2 is of relevance only when at least one of the endpoints is a limit circle.Thus, for the remainder of this section, we assume that on I = (a, b), a is a limit circle while b is a limit point.In addition, we denote by u and v, respectively, the principal and secondary solutions at the limit circle endpoint and we set [u, v](c) = Λ ∈ C for all c ∈ I.We observe that although a and b might be singular endpoints, for any y, z ∈ D max (L max ), since the following limits exist, it makes sense to define Next, we report three results concerning the interplay between the Lagrange sesquilinear form and Equation (11).Their proofs, omitted here, are a direct adaptation to the case in hand respectively of Lemmas 10.
Lemma 3.3.Let L be as per Equation (11).Then for any λ ∈ R and α, β ∈ C, there exists f ∈ D max (L max ) such that Moreover, if a is a regular endpoint, then there exists g ∈ D max (L max ) such that Remark 2. It is interesting to notice that if Λ = 1, Equation ( 16) follows directly from Equation ( 15) by setting g = f .This is not the case if Λ = 1 and this plays a significant part in the discussion of generalized versus regular boundary conditions in the next sections.
The following result concerns properties of the self-adjoint extensions of L min , whereas their existence is a direct consequence of Von Neumann lemma [9,Thm.5.43] since the differential operator L in Equation (10a) has real coefficients.For its proof we refer to [8, Th. 10.4.1], and references therein.Conversely, for any g ∈ D max (L max ) abiding to the conditions in item (1), there exists a self-adjoint extension of L min whose domain D S.A. (LS.A.) is defined as per item (2).
Our quest reaches a finale with the following paramount result, which is specially tailored to befit singular Sturm-Liouville problems.In particular, it is instrumental to relating the existence of multiple self-adjoint extensions to the choice of specific boundary conditions.We include a detailed proof due to its relevance and in light of the fact that it is not exactly the well-known result as per [8, Thm.10.4.5], but rather a slight generalization of it.Namely, the principal and secondary solutions are not necessarily normalized to [u, v] = 1.Theorem 3.5.Let L be as in Equation (11).As per Definition 3.2, let u and v be principal and secondary real-valued solutions at r = a such that [u, v] = Λ.Then for any identifies the domain of a self-adjoint extension of L min .Moreover, all self-adjoint extensions of L min are of this form.
Proof.We divide the analysis in two separate parts: proof of the first statement, and proof of the "moreover" statement.
Let (B 1 , B 2 ) ∈ R 2 \ {(0, 0)}.To prove that Equation ( 17) identifies the the domain of a self-adjoint extension of L min , we shall use Lemma 3.4.In other words we set g = B 1 u+B 2 v Λ and, using Lemma 3.1, we can make sure that Observe that item (2) of Lemma 3.4 is automatically fulfilled by Equation ( 17).
Consider now L S.A. a self-adjoint extension of L min .By Lemma 3.2, for y, g ∈ D max (L max ), it holds where C 1 := [g, u](a) and C 2 := [g, v](a).On account of Lemma 3.4 there exists g / ∈ D min (L min ) such that [g, g](a) = 0.By setting y = g in the equation displayed above and assuming both u and v to be real-valued, it descends In addition, still Lemma 3.4 guarantees that f ∈ D S.A. (L S.A. ) if and only if [f, g](a) = 0.This reduces to Equation ( 17) setting therein Remark 3. Note that the proof of Theorem 3.5 assuming [u, v] = Λ is analogous to that on [8, Thm.10.4.5] that assumes [u, v] = 1.It is worth mentioning that this normalization does not select a secondary solution.It is easy to see that this is the case if we take into account that [u, u] = 0 for real-valued u.In turn, the reality of both u and v is an essential aspect of the validity of the proof.In addition, a consequence of this restriction is that Equation ( 17) can be equivalently written in terms of the Wronskians instead of the Langrage sesquilinear form, i.e.

Generalized (γ, v)-Robin boundary conditions
In this section we take a closer look at the boundary condition stated in Equation ( 17) and we reiterate two important facts: • although Equation ( 17) depends on the choice of the pair (B 1 , B 2 ) = (0, 0), it is always possible to rescale y ∈ D max (L max ) so to fix one of the parameters to 1, i.e. we can consider only pairs of the form (1, B 2 B 1 ); • the characterization of D S.A. (L S.A. ) also depends on the chosen secondary solution.
At the mere level of the Sturm-Liouville problem under consideration, the freedom in the choice of secondary solution is inconsequential if one is interested in characterizing all self-adjoint extensions of the corresponding operator.Notwithstanding, it plays a distinguished, physically relevant rôle when we turn back to analyzing the dynamics of the Klein-Gordon field ruled by Equation (5).The following definitions aim at highlighting this freedom in the overall process and the difference that occurs when considering a singular rather than a regular Sturm-Liouville problem.
In particular, we say that y abides to a (i) generalized Dirichlet boundary condition at a if it satisfies a generalized (0, v)-Robin boundary condition: lim It is worth stressing that Definition 4.1 is applicable only to regular Sturm-Liouville problems since it implicitly requires differentiability of the solution at the endpoint a.In addition, consistently with what one could a priori expect, the "Dirichlet boundary condition" is actually independent of the choice of the secondary solution.In the following, we elucidate more in detail the connection between the two definitions above.

Reduction to the regular case
Consider the setting of Theorem 3.5, and assume that r = a is a regular endpoint as per Definition 3.1.For real-valued u and v such that {u, v} = Λ, a solution y = cos(γ)u + sin(γ)v ∈ L 2 (I, µ(r)dr) satisfies a generalized (γ, v)-Robin boundary condition, as per Equation (18).Hence That is, at a regular endpoint, for a given u there is a choice of secondary solution v for which the generalized (γ, v)-Robin boundary condition, as per Definition 4.2, yields a regular γ-Robin boundary condition, as per Definition 4.1.Conversely, if we do not choose it in such a way and if γ = 0, then a generalized boundary condition does not necessarily reduce to a regular one.

Generalized (γ, v)-Robin boundary conditions and the Klein-Gordon equation
In view of the foregoing discussion, it is natural to wonder what is the consequence of choosing a specific generalized (γ, v)-Robin boundary condition at the level of the fully covariant Klein-Gordon equation.This question becomes especially relevant when we are working on a globally hyperbolic spacetime with timelike boundary [2], such as the Poincaré patch of an n-dimensional anti-de Sitter spacetime (PAdS n ).In this case, it is known that the dynamics is completely specified when, and only when, initial data are supplemented with a boundary condition assigned at conformal infinity.Then, one would slavishly follow the analysis outlined in the previous sections.For definiteness, let us assume we are working under conditions for which Definition 4.2 is meaningful.Explicitly, take Ψ to be a solution of the Klein-Gordon equation ( 5) written as the mode expansion given in Equation ( 6).In addition, let u and v be realvalued principal and secondary solutions at an endpoint a for the radial equation, such that the radial mode satisfies a generalized (γ, v)-Robin boundary condition.One can infer that the latter translates to a boundary condition on Ψ: where ) is the integral over the spectrum of the Laplace operator on . Details regarding the latter are left to the reader since they play no rôle in our discussion, yet we observe that if j = 1 this integral reduces to a sum of (hyper-)spherical harmonics, whereas if j ∈ {0, 1} it is nothing but an ordinary Lesbegue integral.Similarly, σ(A) dλ is formally the integral over the spectrum of the self-adjoint extension A with respect to the associated spectral measure.For all practical purposes in many instances σ(A) = (0, ∞) and the integral reduces to a standard Lesbegue integration on the half real line.
With the discussion of Section 4 in mind, we see that Equation ( 22) reduces to a regular Robin boundary condition only under special conditions.In addition, we highlight that generalized (γ, v)-Robin boundary conditions on R(r) translate at the level of the Klein-Gordon equation, to a wide variety of boundary conditions, including time-dependent ones.Although at this stage, this statement might be elusive and hidden in the meanders of Equation ( 22), it is manifest in the concrete example thoroughly discussed in Section 6.2.

The quantum dynamics
The analysis of the classical solutions to the Klein-Gordon equation is just the starting point to obtain a full-fledged, covariant quantization framework.In this paper we shall not give all the details of the latter, see [16,17,18,19,20,21], rather we focus on the construction of ground states admitting generalized (γ, v)-Robin boundary conditions.Definition 5.1.Let M and P be as defined in Section 2. A two-point function of a quantum state is a bidistribution ψ 2 ∈ D (M × M) such that (i) it solves the Klein-Gordon equation in both entries: (ii) it satisfies the canonical commutation relations: where E is the advanced minus retarded fundamental solution associated to P ; (iii) it is positive: In turn, the bi-distribution E ∈ D (M × M) is a solution of the initial value problem: where Σ t is any constant-time hypersurface, while δ Σ is the Dirac delta thereon.It is important to stress that E is a priori not unique, depending both on the underlying geometry and on the parameters ξ and m 0 of the Klein-Gordon equation, see Equation (5).The details for its construction using tools of Sturm-Liouville and spectral theories can be found in [10] and references therein.
As detailed in [10, Ch.2] and hinted at in the Introduction, among the plethora of two-point functions on a spacetime admitting Schwarzschild-like coordinates, as per Equation ( 4), one can always distinguish the ones that characterize ground states.They are of the form where Θ denotes the Heaviside step function and, for compactness, we have introduced the notation θ = (ϕ 1 , ..., ϕ n−2 ).Using Equation (23a) in combination with the canonical commutation relations in Definition 5.1 and with the completeness of the eigenefunctions of the Laplace operator ∆ Σ n−2 j , it turns out that the unknown ψ 2 (r, r ) can be obtained by the spectral resolution of the Green function G(r, r ) associated to Equation (7), see [24,Ch.7].Namely, promoting λ to a complex variable, one make use of the chain of identities where C ∞ is an infinitely large circle in the λ-plane with a counter-clockwise orientation.In the next section we give a neat and tangible example that unveils how the choice of different secondary solutions, even with the same value of γ, yields well-defined but inequivalent ground states.It corroborates our statement that, in a system where states can be constructed following a mode decomposition, the choice of a secondary solution for the radial equation remains free even after imposing all the physical constraints necessary to guarantee a sensible framework.

An illustrative example: the wave equation on R × R +
In this section, we outline a simple yet most illustrative example aimed at highlighting the relevance of the generalized (γ, v)-Robin boundary conditions: a massless, real, scalar field on the 2-dimensional half-Minkowski spacetime R × R + .Although the endpoints to be considered in the corresponding Sturm-Liouville problem are regular, since the underlying manifold is globally hyperbolic with a timelike boundary, it is conceivable to impose thereon generalized (γ, v)-Robin boundary conditions.
As functions of x and considering ω as a possibly complex parameter rather than a Fourier variable conjugated to time t, it follows that y 1 (x) and y 2 (x) are not squareintegrable at x = ∞ unless Im(ω) > 0 and Im(ω) < 0, respectively.Therefore, as per Definition 3.1, this endpoint is a limit point and the most general square-integrable solution therein can be written concisely as Both y 1 (x) and y 2 (x) are square-integrable in a neighborhood of x = 0. Still according to Definition 3.1, x = 0 is a limit circle.

Generalized versus regular Robin boundary conditions
Since the limits of y 1 , y 2 and of their derivatives exist as x → 0, we can cast the generalized (γ, v κ )-Robin boundary condition above as a regular γ-Robin boundary condition: We find that which highlights the difference between a generalized and a regular Robin boundary condition.Markedly, in the regular scenario we can choose a secondary solution based on the property of obtaining a frequency-independent parameter β κ ; in this case, β 1 .Yet, frequency-dependent boundary conditions are also of physical relevance and hence there is no a priori reason to discard them.

Time-dependence of the boundary conditions
Analogously to the discussion in Section 4.2, given a radial solution of the wave function that satisfies a generalized (γ, v κ )-Robin boundary condition, it is legitimate to wonder which is the corresponding boundary condition satisfied by the solution of the wave equation on R × R + .For ψ κ given by Equation ( 27), let us consider a general solution: It can be written as where, with a slight abuse of notation, we have denoted the Fourier transform with a hat.We find We can now read which boundary condition is satisfied by Ψ κ (t, x) at x = 0 for each κ.
Formally, working at the level of distributions, we look for operators such that The above are satisfied if we take c κ and ζ κ as Accordingly, the solution Ψ κ (t, x) satisfies, for κ ∈ {2, 3}: That is, the solutions Ψ 2 (t, x) and Ψ 3 (t, x) do not satisfy a regular, time-independent, γ-Robin boundary condition at the boundary as it does Ψ 1 (t, x).

The Green functions
To construct the ground state for a Klein-Gordon field admitting generalized (γ, v κ )-Robin boundary conditions on the two-dimensional half-Minkowski spacetime, we follow the rationale outlined in Section 5. On account of Equation ( 25), the only unknown is ψ 2 (x, x ), which in turn can be constructed from the Green function G(x, x ) of the radial equation.Note that, in the case in hand, the rôle of the radial coordinate r is played by the cartesian coordinate x.
Let the Wronskian between the principal and the secondary solutions be given by {u, v κ } =: Λ κ .Then the one between the two general square-integrable solutions given by Equation (26) and Equation (27) reads Therefore, the Green functions G κ (x, x ) for the different choices of secondary solutions, which satisfy where (x < , x > ) = (x, x ) if x < x and (x < , x > ) = (x , x), otherwise.Observe that the dependence on ω is implicit in the solutions ψ κ and ψ ∞ .Moreover, for γ ∈ R, it holds true that For notational convenience, but with no loss of generality, let us consider x < x and γ ∈ [0, π 2 ] in the remainder of this section.Once more we fix s := sign(Im(ω)), Im(ω) = 0. Explicitly, the Green functions can be written as .

The resolution of the identity
For all three values of κ, the suitable contour to be considered for the integration of G κ is the "pac-man" in the ω 2 -complex plane, which is tantamount to integrating 2ωG κ in one semi-disk in the upper or lower ω-complex plane, as illustrated in Figure 1.Note that although G 1 and G 3 diverge in the limit ω → 0, ωG κ has no poles in the ω-complex plane.In addition, since Jordan's lemma holds true, it follows that Next, let us outline the computation of the integral given in Equation (31), to highlight the fact that all three possible choices of a secondary solution does yield the resolution of the identity, as per Equation (25).For G 1 , this computation is standard, [24,Ex.7.3.2,Pg.454].Yet, we do provide a step-by-step solution for all κ in a supplementary notebook available online [25].

The two point functions
Directly from the spectral resolution, as stated in Section 5, we obtain the spatial part of the two-point function in each case: With the three choices of secondary solutions we obtain three one-parameter families of two-point functions for three different ground states in the 2-dimensional half-Minkowski spacetime: Ψ κ (t, x, t , x ) =  33) for different values of κ.One may ponder on the significance of such discrepancy, since it could be the case that different integrands yield equal integrals (as it happens for the resolution of the identity, for example).Yet, as it turns out, the term ψ κ (x, x) has a physical interpretation: it characterizes the probability of de-excitations of a two-level system with energy gap Ω, at a fixed spatial position x and interacting for an infinite time with the quantum field in the ground state specified by Ψ κ (t, x, t , x ).Such physical observable has been extensively used in the last years to probe a wide range of characteristics of the underlying quantum field theoretical framework, see [10,26] and the references therein.

Conclusion
In this work we have highlighted the existence of a hidden freedom in the standard procedure of constructing ground states for a real scalar field in a large class of static backgrounds.In particular we have observed that, when working at the level of the so-called radial equation, boundary conditions of Robin type can be imposed by using an arbitrary secondary solution.While this choice appears to be moot at the level of the underlying ordinary differential equation, it bears notable consequences at a fully covariant level.Interestingly, we have argued and shown via the concrete example of the two-dimensional half-Minkowski spacetime that, by exploring such freedom, one can account for a large class of boundary conditions that are structurally quite different from regular Robin boundary conditions, possibly including time-dependent ones.
From a structural viewpoint, the choice of secondary solution does not alter the effectiveness of the methods used until now for the construction of ground on static spacetimes.Nevertheless, it does open the possibility of studying a much larger class of boundary conditions and of investigating the physical consequence of the various different choices.To conclude we emphasize that a rationale similar to the one considered in this paper can be taken also in the investigation of boundary condition of Wentzell type, see e.g.[19,27].Yet, a full-fledged analysis of this scenario would require a lengthy discussion that is worth leaving to a future work.

Definition 4 . 1 .
(Regular γ-Robin boundary condition) Let L be as per Equation(11).Given any self-adjoint realization L S.A. , we say that y ∈ D S.A. (L S.A. ) satisfies a regular γ-Robin boundary condition at a if lim r→a {cos(γ)y + sin(γ)y } = 0 for γ ∈ [0, π), where the prime indicates the derivative along the r-direction.In particular, we say that y abides to a (i) regular Dirichlet boundary condition at a if it satisfies a regular 0-Robin boundary condition: lim r→a y = 0 and lim r→a y = c ∈ R. (ii) regular Neumann boundary condition at a if it satisfies a regular π 2 -Robin boundary condition: lim r→a y = 0 and lim r→a y = c ∈ R. Definition 4.2.(Generalized (γ, v)-Robin boundary condition) Let L be as per Equation (11) and let u be the principal solution at a and v any secondary solution at a, real-valued and such that {u, v} = Λ.Given any self-adjoint realization L S.A. , we say that y ∈ D S.A. (L S.A. ) satisfies a generalized (γ, v)-Robin boundary condition at a if lim r→a {cos(γ){y, u} + sin(γ){y, v}} = 0 for γ ∈ [0, π).

Remark 4 .
Each generalized (γ, v)-Robin boundary condition on the radial mode R(r) yields a different Green function G(r, r ) for the radial equation and, consequently, a different two-point function-a different ground state.