Stationary bound-state massive scalar field configurations supported by spherically symmetric compact reflecting stars

It has recently been demonstrated that asymptotically flat neutral reflecting stars are characterized by an intriguing no-hair property. In particular, it has been proved that these horizonless compact objects cannot support spatially regular static matter configurations made of scalar (spin-0) fields, vector (spin-1) fields and tensor (spin-2) fields. In the present paper we shall explicitly prove that spherically symmetric compact reflecting stars can support stationary (rather than static) bound-state massive scalar fields in their exterior spacetime regions. To this end, we solve analytically the Klein–Gordon wave equation for a linearized scalar field of mass μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu $$\end{document} and proper frequency ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} in the curved background of a spherically symmetric compact reflecting star of mass M and radius Rs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{\text {s}}$$\end{document}. It is proved that the regime of existence of these stationary composed star–field configurations is characterized by the simple inequalities 1-2M/Rs<(ω/μ)2<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1-2M/R_{\text {s}}<(\omega /\mu )^2<1$$\end{document}. Interestingly, in the regime M/Rs≪1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M/R_{\text {s}}\ll 1$$\end{document} of weakly self-gravitating stars we derive a remarkably compact analytical equation for the discrete spectrum {ω(M,Rs,μ)}n=1n=∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\omega (M,R_{\text {s}},\mu )\}^{n=\infty }_{n=1}$$\end{document} of resonant oscillation frequencies which characterize the stationary composed compact-reflecting-star–linearized-massive-scalar-field configurations. Finally, we verify the accuracy of the analytically derived resonance formula of the composed star–field configurations with direct numerical computations.


Introduction
In their physically interesting and mathematically elegant nohair theorems, Bekenstein [1] and Teitelboim [2] (see also [3][4][5][6][7][8][9] and references therein) have explicitly proved that static non-linear matter configurations which are made of spatially regular massive scalar fields cannot be supported in the extea e-mail: shaharhod@gmail.com rior spacetime regions of asymptotically flat black holes 1 (see  for the case of stationary black-hole spacetimes). This intriguing physical property of the static sector of the non-linearly coupled Einstein-scalar field equations is often attributed to the fact that classical black-hole spacetimes are characterized by one-way membranes (event horizons) that irreversibly absorb matter and radiation fields.
Interestingly, it has recently been demonstrated that the presence of an absorbing (attractive) horizon is not a necessary condition for the validity of the no-hair property for compact objects in general relativity. Specifically, it has been explicitly proved in [40][41][42] that asymptotically flat horizonless compact stars with reflecting (that is, repulsive rather than attractive) boundary conditions are also characterized by an intriguing no-hair property. In particular, the no-hair theorems presented in [40][41][42][43] have revealed the interesting fact that spatially regular static configurations made of scalar (spin-0) fields, vector (spin-1) fields and tensor (spin-2) fields cannot be supported by spherically symmetric compact stars with reflecting 2,3 boundary conditions. It is of physical interest to test the general validity of the no-hair property revealed in [40][41][42] for horizonless compact reflecting objects in general relativity. In particular, in the present paper we raise the following physically intriguing question: Can compact reflecting stars 4 support stationary (rather than static) scalar fields in their exterior spacetime regions?
In order to address this physically interesting question, we shall analyze below the Klein-Gordon wave equation for a stationary linearized scalar field of mass μ and proper frequency ω in the background of a spherically symmetric compact reflecting star of mass M and radius R s . We shall first prove that the stationary (that is, non-decaying in time) massive scalar field configurations, if they exist, must be characterized by the compact relation 1 − 2M/R s < (ω/μ) 2 < 1. 5 Interestingly, we shall then prove explicitly that spherically symmetric compact reflecting stars can support stationary bound-state configurations made of spatially regular linearized massive scalar fields. In particular, we shall derive a remarkably compact analytical formula [see Eq. (34) below] for the discrete resonant oscillation spectrum {ω(M, R s , μ; n)} n=∞ n=1 6 which characterizes the stationary composed reflecting-star-massive-scalar-field configurations in the physical regime M/R s 1.
Footnote 3 continued existence of spatially regular static scalar hair outside the compact reflecting stars. The no-hair theorem for the compact reflecting objects can easily be extended to the case of Neumann boundary conditions with dψ/dr = 0 at the surface r = R s of the star, in which case the scalar Klein-Gordon differential equation is violated at the surface itself. Interestingly, the regime of validity of the no-hair property for spherically symmetric asymptotically flat compact reflecting objects can further be extended to the more generic case of Robin boundary conditions with ψ + b · dψ/dr = 0 and b < 0 at the surface r = R s of the star. In this more generic case, the radial scalar eigenfunction is characterized by the functional relation ψ · dψ/dr >0 at r = R s , a property which, together with the characteristic asymptotic behavior ψ(r → ∞) → 0 of finite mass matter configurations [see Eqs. (7) and (8) of [40]], imply that the scalar eigenfunction must have an extremum point, r = r peak ≥ R s , outside the reflecting star. As explicitly proved in [40], the Klein-Gordon equation for the static scalar field is violated at this extremum point, a fact which rules out the existence of static spatially regular scalar hairy configurations outside the compact reflecting objects. 4 The term 'reflecting star' is used here to describe an object with a well defined compact surface on which the scalar field vanishes [see Eq. (7) below]. 5 It is worth noting that, as expected, this result agrees with the noscalar-hair theorem presented in [40], according to which static (ω = 0) field configurations cannot be supported by compact reflecting stars. 6 Here the resonance parameter n is an integer which characterizes the stationary bound-state massive scalar field configurations in the spacetime of the spherically symmetric compact reflecting star [see Eq. below].

Description of the system
We consider a physical system which is composed of a spherically symmetric compact reflecting star which is linearly coupled to a stationary bound-state massive scalar field configuration. The curved spacetime outside the compact star of mass M is characterized by the spherically symmetric line element [44] 7 where R s is the radius of the star and The dynamics of the linearized massive scalar field (t, r, θ, φ) in the curved spacetime of the compact star is governed by the familiar Klein-Gordon wave equation [45][46][47][48][49][50][51][52] It is convenient to decompose the scalar field eigenfunction in the form 9 (t, r, θ, φ) = lm e imφ S lm (θ )R lm (r ; ω)e −iωt dω. (4) Substituting the scalar eigenfunction (4) and the metric components of the curved line element (1) into the Klein-Gordon wave equation (3), one finds that the spatial behavior of the radial scalar function R lm (r ) is determined by the ordinary differential equation [45][46][47][48][49][50][51][52][53][54] d dr where K l = l(l + 1) (with l ≥ |m|) is the familiar angular eigenvalue of the spatially regular angular scalar eigenfunction S lm (θ ). Clearly, the radial scalar equation (5) is invariant under the symmetry transformation ω → −ω. We shall henceforth assume, without loss of generality, that We are interested in stationary (that is, non-decaying in time) bound-state configurations of the massive scalar fields which are supported by the central compact reflecting star. These field configurations are characterized by the following two boundary conditions: (1) an inner boundary condition of a vanishing scalar eigenfunction [40,55,56] (r = R s ) = 0 at the reflecting surface of the spherically symmetric compact star, and (2) the asymptotic large-r boundary condition of an exponentially decaying (normalizable) scalar eigenfunction at spatial infinity [here κ ≡ M(2ω The ordinary differential equation (5) for the radial scalar eigenfunction, supplemented by the boundary conditions (7) and (8) at the surface of the spherically symmetric compact reflecting star and at spatial infinity, determines the discrete spectrum {ω(M, R s , μ, l; n)} (see footnote 6) of resonant oscillation frequencies which characterize the stationary bound-state massive scalar field configurations in the curved spacetime of the central supporting star. Below we shall use analytical techniques in order to determine this characteristic star-field resonance spectrum.

Lower bound on the resonant oscillation frequencies of the stationary bound-state massive scalar field configurations
In the present section we shall use simple analytical techniques in order to prove that the stationary composed compact-reflecting-star-massive-scalar-field configurations, if they exist, are characterized by resonant oscillation frequencies which, for given values of the physical parameters {M, R s , μ}, cannot be arbitrarily small. To this end, it proves useful to define the new radial function 10 in terms of which the scalar differential equation (5) takes the form 10 For brevity, we shall henceforth omit the harmonic indices {l, m} of the stationary linearized scalar field modes.
Taking cognizance of Eqs. (7), (8), and (9), one deduces that the radial eigenfunction (r ), which characterizes the spatial behavior of the stationary bound-state massive scalar field configurations in the curved spacetime (1) of the spherically symmetric compact reflecting star, must have (at least) one extremum point in the interval In particular, at this extremum point the scalar eigenfunction (r ) is characterized by the simple relations Substituting the characteristic functional relations (12) into the scalar differential equation (10), one deduces that the composed compact-reflecting-star-stationary-boundstate-massive-scalar-field configurations, if they exist, are characterized by the (rather cumbersome) relation at the extremum point r = r ext . Taking cognizance of Eq.
(2), one finds that the inequality (13) yields the upper bound on the resonant oscillation frequencies of the stationary bound-state massive scalar field configurations. The r.h.s. of (14) is maximized for the simple choice δ = −[2 f (r ext )] −1 , in which case one obtains from (14) the simple lower bound Finally, using the inequality r ext > R s in (15) and taking cognizance of equation (8), one finds that the characteristic resonant oscillation frequencies of the stationary boundstate massive scalar field configurations in the curved spacetime of the spherically symmetric compact reflecting star are restricted to the dimensionless frequency interval (see foot-

The characteristic resonance condition of the stationary composed compact-reflecting-starlinearized-massive-scalar-field configurations
Interestingly, in the present section we shall explicitly show that the ordinary differential equation (5), which determines the spatial behavior of the stationary bound-state massive scalar fields in the spacetime of the compact reflecting star, is amenable to an analytical treatment in the large-radii regime 11 In particular, we shall now derive the fundamental resonance condition which determines the discrete spectrum {ω(M, R s , μ, l; n)} of resonant oscillation frequencies that characterize the composed reflecting-star-stationary-scalarfield configurations in the physical regime (17). Substituting the composed radial function into the scalar radial equation (5) and neglecting terms of order O(M/r 3 ) in the large-r regime (17), one obtains the characteristic Schrödinger-like differential equation which determines the spatial behavior of the massive scalar field configurations. It is convenient to use the dimensionless radial coordinate in terms of which the scalar radial equation (19) takes the compact form where the dimensionless physical parameters {κ, β}, which characterize the composed star-field system, are given by 11 Taking cognizance of the analytically derived inequalities (16)], one realizes that the assumed large-radii regime (17) corresponds to stationary bound-state massive scalar field configurations with ω μ. and The general mathematical solution of the Whittaker differential equation (21) can be expressed in terms of the familiar confluent hypergeometric functions (see Eqs. 13.1.32 and 13.1.33 of [53]): where {A, B} are normalization constants.
The asymptotic spatial behavior of the radial function (24) Taking cognizance of the boundary condition (8), which characterizes the asymptotic spatial behavior of the stationary bound-state (normalizable) massive scalar field configurations, one deduces that the coefficient of the exponentially diverging term in the asymptotic radial equation (25) should vanish: One therefore finds that, in the physical regime R s M of weakly self-gravitating stars [see (17)], the supported boundstate configurations of the stationary massive scalar fields are characterized by the radial eigenfunction where U (a, b, z) is the confluent hypergeometric function of the second kind [53]. Defining the dimensionless physical parameters and taking cognizance of the inner boundary condition (7) which characterizes the behavior of the massive scalar fields at the reflecting surface of the central compact star, one obtains from the radial solution (27) the non-linear resonance condition Table 1 Resonant oscillation frequencies of the composed stationary compact-reflecting-star-linearized-massive-scalar-field configurations with γ = 10 and l = 0. We present the smallest resonant oscillation frequency min (α) of the stationary bound-state massive scalar fields for various values of the dimensionless star-field mass parameter α ≡ Mμ.
One finds that the characteristic resonant frequency min (α) of the composed star-field system is a monotonically decreasing function of the dimensionless mass parameter α.  Table 2 Resonant oscillation frequencies of the composed stationary compact-reflecting-star-linearized-massive-scalar-field configurations with α = 0.1 and l = 0. We display the smallest resonant oscillation frequency min (γ ) of the stationary bound-state massive scalar fields for various values of the dimensionless mass-radius parameter γ ≡ μR s .
One finds that the characteristic resonant frequency min (γ ) of the composed star-field system is a monotonically increasing function of the dimensionless mass-radius parameter γ . As expected, the resonant oscillation frequencies conform to the analytically derived bounds Dimensionless for the composed spherically-symmetric-reflecting-starlinearized-massive-scalar-field configurations. Interestingly, as we shall explicitly show below, the mathematically compact resonance equation (29), which is valid in the large-radii regime (17), determines the discrete family of resonant oscillation frequencies { (α, γ, ; n)} which characterize the stationary bound-state massive scalar field configurations in the curved spacetime of the spherically symmetric compact reflecting star.

The characteristic resonance spectrum of the stationary composed compact-reflecting-star-massive-scalar-field configurations
The analytically derived resonance condition (29), which determines the resonant oscillation spectrum of the stationary bound-state linearized massive scalar fields in the background of the compact reflecting star, can easily be solved numerically. Interestingly, one finds that a composed compact-reflecting-star-massive-scalar-field system with a given set of the dimensionless physical parameters {α, γ , } is characterized by a discrete resonant oscillation spectrum of the form 12 12 As noted above, the radial scalar equation (5) is invariant under the symmetry transformation ω → −ω. We shall assume that ω>0 without loss of generality.
In Table 1 we display, for various values of the dimensionless star-field mass parameter α ≡ Mμ, the value of the smallest resonant oscillation frequency min (α). Interestingly, one finds that the resonant oscillation frequency min (α), which characterizes the composed stationary compact-reflecting-star-massive-scalar-field configurations in the large-radii regime (17), is a monotonically decreasing function of the dimensionless physical parameter α. As a consistency check, we note that the numerically computed resonant oscillation frequencies min (α) of the composed star-field system conform to the analytically derived bounds (16)]. In Table 2 we present, for various values of the dimensionless mass-radius parameter γ ≡ μR s of the composed star-field system, the value of the smallest resonant oscillation frequency min (γ ) which characterizes the stationary bound-state field configurations. One finds that the characteristic resonant oscillation frequency min (γ ) of the stationary composed compact-reflecting-star-massive-scalar-field configurations is a monotonically increasing function of the dimensionless mass-radius parameter γ .
In Table 3 we display, for various values of the dimensionless angular harmonic index l of the stationary scalar field modes, the value of the smallest resonant oscillation frequency min (l). Interestingly, one finds that the resonant oscillation frequency min (l) of the stationary boundstate massive scalar field configurations is a monotonically increasing function of the spherical harmonic index l. Again, one finds that the numerically computed resonant oscillation frequencies min (l) of the composed star-field system Table 3 Resonant oscillation frequencies of the composed stationary compact-reflecting-star-linearized-massive-scalar-field configurations with α = 1 and γ = 10. We present the smallest resonant oscillation frequency min (l) which characterizes the stationary bound-state starfield configurations for various values of the dimensionless harmonic index l. One finds that the characteristic resonant frequency min (l) of the composed star-field system is a monotonically increasing function of the angular harmonic index l. As a consistency check, we note that the numerically computed resonant oscillation frequencies conform to the analytically derived bounds conform to the analytically derived bounds √ 1 − 2M/R s < < 1.

Analytical treatment for marginally-bound stationary scalar field resonances
In the present section we shall explicitly show that the resonance equation (29), which determines the discrete family of resonant oscillation frequencies { (α, γ, ; n)} that characterize the stationary bound-state linearized massive scalar fields in the curved spacetime of the spherically symmetric compact reflecting star, can be solved analytically in the asymptotic regime 13 Substituting the asymptotic equation (see Eq. 13.5.16 of of the confluent hypergeometric function into the analytically derived resonance condition (29), one finds the characteristic equation for the marginally-bound ( 2 → 1 − ) stationary scalar field configurations. The resonance condition (33) yields the compact expression 14 13 Note that the asymptotic frequency regime (31) corresponds to the physically interesting case of marginally-bound stationary composed compact-reflecting-star-massive-scalar-field configurations. 14 Here the integer n 1 is the resonance parameter which characterizes the stationary composed star-field configurations.
for the discrete spectrum of resonant oscillation frequencies which characterize the marginally-bound composed spherically-symmetric-compact-reflecting-star-stationarymassive-scalar-field configurations.
It is worth emphasizing the fact that, in accord with the exact (numerically computed) data presented in Tables 1 and  2, the analytically derived expression (34) for the asymptotic resonance spectrum implies that the resonant oscillation frequencies which characterize the stationary composed compact-reflecting-star-massive-scalar-field configurations are a monotonically decreasing function of the dimensionless star-field mass parameter α ≡ Mμ and a monotonically increasing function of the dimensionless mass-radius parameter γ ≡ μR s .

Numerical confirmation
It is of physical interest to test the accuracy of the approximated (analytically derived) compact equation (34) for the asymptotic resonance spectrum which characterizes the stationary marginally-bound composed compact-reflectingstar-massive-scalar-field configurations.
In Table 4 we present the resonant oscillation frequencies analytical (n; α, γ , ) of the stationary massive scalar fields as obtained from the approximated (analytical) resonance equation (34) in the asymptotic regime (31) of marginallybound star-field configurations. Also displayed in Table  4 are the corresponding resonant oscillation frequencies numerical (n; α, γ , ) of the stationary bound-state massive scalar field configurations as obtained by solving numerically the analytically derived resonance equation (29). In the physical regime (31) of marginally-bound star-field configurations, one finds a remarkably good agreement 15 between the Table 4 Resonant oscillation frequencies of the composed stationary compact-reflecting-star-linearized-massive-scalar-field configurations with α = 1, γ = 10, and l = 0. We display the resonant oscillation frequencies analytical (n; α, γ , ) of the stationary scalar fields as obtained from the analytically derived resonance equation (34) in the asymptotic regime (31) of marginally-bound star-field configurations. We also present the corresponding numerically computed resonant oscillation frequencies numerical (n; α, γ , ) of the stationary composed star-field configurations. In the physical regime (31) of marginally-bound massive scalar field configurations, one finds a remarkably good agreement (see footnote 15) between the exact resonant field frequencies [as obtained numerically from the characteristic resonance condition (29)] and the corresponding analytically derived resonant oscillation frequencies (34) which characterize the stationary bound-state star-field configurations  (29)] and the corresponding analytically derived resonant oscillation frequencies (34) which characterize the stationary composed compact-reflecting-star-massive-scalar-field configurations.

Summary and discussion
The elegant no-scalar-hair theorems of [1,2] (see also [3][4][5][6][7][8][9] and the references therein) have revealed the interesting fact that asymptotically flat classical black-hole spacetimes with one-way membranes (absorbing event horizons) cannot support static configurations made of massive scalar fields. Intriguingly, it has recently been demonstrated explicitly that the presence of an absorbing boundary (an event horizon) is not a necessary condition for the validity of the no-hair property for compact gravitating objects in general relativity. In particular, it has been proved that asymptotically flat horizonless neutral reflecting (rather than absorbing) stars are 'bald' [40] in the sense that they cannot support static bound-state scalar (spin-0) fields, vector (spin-1) fields, and tensor (spin-2) fields [40][41][42]. One naturally wonders whether the no-scalar-hair behavior observed in [40] is a generic property of compact reflecting stars? In particular, in the present paper we have raised the following physically intriguing question: Can horizonless compact objects with reflecting boundary conditions (as opposed to the absorbing boundary conditions which characterize the horizons of classical black-hole spacetimes) support stationary (rather than static) massive scalar field configurations?
In order to address this physically interesting question, in the present paper we have studied analytically the characteristic Klein-Gordon wave equation for stationary linearized massive scalar field configurations in the curved spacetime of a spherically symmetric compact reflecting star. The main results derived in this paper and their physical implications are as follows: 1. It has been explicitly proved that the stationary composed compact-reflecting-star-linearized-massive-scalar-field configurations, if they exist, are restricted to the dimensionless frequency regime [see Eq. (16)] (see footnote 12) It is worth noting that the presence of a positive lower bound on resonant oscillation frequencies of the composed star-field configurations agrees with the no-scalarhair theorem presented in [40], according to which spherically symmetric compact reflecting stars cannot support spatially regular static (ω = 0) scalar field configurations. 2. We have then proved that spherically symmetric compact reflecting stars can support stationary (rather than static) bound-state linearized massive scalar fields. In particular, it has been shown that the resonant oscillation modes of the composed compact-reflecting-starmassive-scalar-field system depend on three dimensionless physical parameters: the star-field mass parameter α ≡ Mμ, the mass-radius parameter γ ≡ μR s , and the spherical harmonic index l of the field mode. Solving analytically the Klein-Gordon wave equation for the stationary massive scalar field modes in the large-radii regime M/R s 1, we have explicitly proved that the resonance equation [see Eqs. (28) and (29)] determines the discrete spectrum { (α, γ, ; n)} n=∞ n=1 of dimensionless resonant oscillation frequencies which characterize the stationary bound-state scalar field modes in the spacetime of the spherically symmetric compact reflecting star.
3. Solving numerically the analytically derived resonance condition (36) for the composed star-field configurations, we have demonstrated that the smallest possible resonant oscillation frequency min [see Eq. (30)] which characterizes the stationary bound-state scalar field configurations is a monotonically decreasing function of the dimensionless star-field mass parameter Mμ and a monotonically increasing function of the dimensionless mass-radius parameter μR s and the dimensionless angular harmonic index l (see Tables 1, 2, 3). 4. It has been explicitly proved that the large-radii resonance equation (36) for the composed star-field configurations is amenable to an analytical treatment in the physically interesting regime (ω/μ) 2 → 1 − of marginally-bound stationary massive scalar field configurations. In particular, we have derived the remarkably simple dimensionless analytical formula [see Eqs. (28) and (34)] for the discrete resonance spectrum which characterizes the stationary marginally-bound linearized massive scalar field configurations in the spacetime of the spherically symmetric horizonless reflecting star. 5. It has been explicitly demonstrated that the characteristic resonant frequencies of the stationary marginally-bound composed compact-reflecting-star-linearized-massivescalar-field configurations, as deduced from the analytically derived resonance equation (37), agree remarkably well (see footnote 15) with the corresponding exact resonant oscillation frequencies of the stationary massive scalar fields as obtained numerically from the resonance condition (36) (see Table 4).
Finally, we would like to emphasize again that in the present exploration of the physical and mathematical properties of the composed star-field system we have made two simplifying assumptions: (1) The scalar field was treated at the linear level, and (2) the central supporting star was treated as a weakly self-gravitating object. 16 As we have explicitly proved in this paper, using these two simplifying assumptions, the physical and mathematical properties of the stationary composed compact-reflecting-star-linearized-massivescalar-field system can be studied analytically. As a next step, it would be interesting to use numerical techniques in order