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 {\it horizonless} compact objects cannot support spatially regular {\it 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 {\it 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 $\mu$ and proper frequency $\omega$ in the curved background of a spherically symmetric compact reflecting star of mass $M$ and radius $R_{\text{s}}$. It is proved that the regime of existence of these stationary composed star-field configurations is characterized by the simple inequalities $1-2M/R_{\text{s}}<(\omega/\mu)^2<1$. Interestingly, in the regime $M/R_{\text{s}}\ll1$ of weakly self-gravitating stars we derive a remarkably compact {\it analytical} formula for the discrete spectrum $\{\omega(M,R_{\text{s}},\mu)\}^{n=\infty}_{n=1}$ 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.


I. INTRODUCTION
In their physically interesting and mathematically elegant no-hair theorems, Bekenstein [1] and Teitelboim [2] (see also [3][4][5][6][7][8] 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 exterior spacetime regions of asymptotically flat black holes [9][10][11]. 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 [12][13][14] 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 [12][13][14] 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 [15,16] boundary conditions. It is of physical interest to test the general validity of the no-hair property revealed in [12][13][14] 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 [17] 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 [18]. 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 [19] which characterizes the stationary composed reflecting-star-massive-scalar-field configurations in the physical regime M/R s ≪ 1.

II. 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 [20,21] 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 [22][23][24][25][26][27] It is convenient to decompose the scalar field eigenfunction in the form [28] 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 [22-26, 29, 30] 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 [12,31] Ψ(r = R s ) = 0 (7) 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ω 2 − µ 2 )/ µ 2 − ω 2 , see Eqs. (21)-(22) below]. 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)} [19] of resonant oscillation frequencies which characterize the stationary boundstate 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.

III. 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 [32] in terms of which the scalar differential equation (5) takes the form 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-bound-state-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 δ = −[2f (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 the relation (8), one finds that the characteristic resonant oscillation frequencies of the stationary bound-state massive scalar field configurations in the curved spacetime of the spherically symmetric compact reflecting star are restricted to the dimensionless frequency interval

IV. THE CHARACTERISTIC RESONANCE CONDITION OF THE STATIONARY COMPOSED COMPACT-REFLECTING-STAR-LINEARIZED-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 [33] R s ≫ M .
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 and The general mathematical solution of the Whittaker differential equation (21) 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 expression (25) should vanish: One therefore finds that, in the physical regime R s ≫ M of weakly self-gravitating stars [see (17)], the supported bound-state 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 [29]. 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 for the composed spherically-symmetric-reflecting-star-linearized-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.

V. 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 [34] 1 with the asymptotic property ̟ ∞ → 1 [see Eq. (34) below].
In Table I 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 1 − 2M/R s < ̟ < 1 [see Eq. (16)]. In Table II 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 γ.  II: 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 γ ≡ µRs. 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 1 − 2M/Rs < ̟ < 1.
In Table III 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 bound-state 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 conform to the analytically derived bounds 1 − 2M/R s < ̟ < 1.

VI. 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  III: 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 star-field 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 1 − 2M/Rs < ̟ < 1 [see Eq. (16)].
asymptotic regime [35] Substituting the asymptotic relation (see Eq. 13.5.16 of [29]) 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 [36] for the discrete spectrum of resonant oscillation frequencies which characterize the marginally-bound composed spherically-symmetric-compact-reflecting-star-stationary-massive-scalar-field configurations. It is worth emphasizing the fact that, in accord with the exact (numerically computed) data presented in Tables I and  II, 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 .

VII. NUMERICAL CONFIRMATION
It is of physical interest to test the accuracy of the approximated (analytically derived) compact formula (34) for the asymptotic resonance spectrum which characterizes the stationary marginally-bound composed compact-reflectingstar-massive-scalar-field configurations.
In Table IV we present the resonant oscillation frequencies ̟ analytical (n; α, γ, ℓ) of the stationary massive scalar fields as obtained from the approximated (analytical) resonance formula (34) in the asymptotic regime (31) of marginallybound star-field configurations. Also displayed in Table IV 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 [37] between the exact resonant oscillation frequencies [as computed numerically directly from the resonance equation (29)] and the corresponding analytically derived resonant oscillation frequencies (34) which characterize the stationary composed compact-reflecting-star-massive-scalar-field configurations.
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' [12] in the sense that they cannot support static bound-state scalar (spin-0) fields, vector (spin-1) fields, and tensor (spin-2) fields [12][13][14].
One naturally wonders whether the no-scalar-hair behavior observed in [12] 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) It is worth noting that the presence of a positive lower bound on resonant oscillation frequencies of the composed starfield configurations agrees with the no-scalar-hair theorem presented in [12], 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-star-massive-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 starfield 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 I, II, and III).
(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)] ; n ≫ 1 (37) 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 marginallybound composed compact-reflecting-star-linearized-massive-scalar-field configurations, as deduced from the analytically derived resonance formula (37), agree remarkably well [37] with the corresponding exact resonant oscillation frequencies of the stationary massive scalar fields as obtained numerically from the resonance condition (36) (see Table IV).
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 [38]. As we have explicitly proved in this paper, using these two simplifying assumptions, one can study analytically the physical and mathematical properties of the stationary composed compact-reflecting-star-linearized-massive-scalar-field system. As a next step, it would be interesting to use numerical techniques in order to extend our results to the physically important regime of non-linear stationary bound-state scalar field configurations supported by strongly self-gravitating horizonless compact stars.