Arxiv:Astro-Ph/0411548V1 18 Nov 2004 Yterlaeo H Rvttoa Nryadproduce and Sur- Into Space

arXiv:astro-ph/0411548v1 18 Nov 2004 yterlaeo h rvttoa nryadproduce and sur- into space. energy heated propagate rounding gravitational freely interior, that the Usual radiation bubble of electromagnetic the release quantization. into the redshift by falls the picture matter for this baryonic into account well can the fit of candidate, and one matter than axions, dark by more Hypothetical leading scale collapse time years. against a million on supported hundreds decays is and pressure bubble quantum The quasars. select red- (distance we nearby different If most with with following. associated the groups QSOs concluded two and properties into shift (QSOs) quasars of hc bec fayproiiywscamd[] The (op- faint [4]. intrinsically luminosity claimed of tical was existence periodicity suggest any observations of absence which ab-odtn oos(NB)[] nalsc mod- such all pseudo In with [7]. connection (PNGBs) in bosons theory Nambu-Goldstone filed quantum in rived idcitrcinptnil assaot10 about masses potential, interaction 6]. [5, riodic nucleus the galactic in the or from galaxies ejection active of nearby process from such ejected the that are of indicate quasars % Observations few a population. in QSO only known constitute quantization and redshift samples) distances QSO of large distant from disappearance detected yields not (which are intrinsically objects Being such faint, redshift. noncosmological quantized n ai 10 radii and where eidcitrcinpotential interaction periodic otnab aaisweeteqatzto a con- (0 was with distant quantization and associated 3] the QSOs [2, where of firmed galaxies studies nearby later division most a by Such pe- supported observed. no Meanwhile, is is galaxies redshift distant intrinsic below. with in 3 riodicity associated val- Fig. samples certain in QSO to shown in close as is (quantized), redshift ues their then galaxies Mpc), nti etrw td asv elsaa field scalar real massive study we Letter this In division noted has [1] Karlsson observations, on Based eew hwta ube fdr atrwt pe- with matter dark of bubbles that show we Here nrnial an usr:eiec for evidence quasars: faint Intrinsically V 0 > 3 osrie yatohscladcsooia arguments. cosmological and astrophysical by constrained htitisclyfitpitlk usr soitdwi quant observed associated 10 quasars the masses point-like explain faint results st intrinsically yields Our that and years. collapse million against dreds v bubble any the dar have supports leading can which the pressure redshift of gravitational one wit equal as possesses axions considered from currently made are particles bubbles redshift describing quantized equations noncosmological Einstein with quasars) known of nteUies n il h xo mass axion the i yield unambiguously and quasars Universe such the of in Properties space. rounding .Ti oeta sqiegnrladde- and general quite is potential This 0. eateto hsc,IsiuefrQatmSuis Texa Studies, Quantum for Institute Physics, of Department − rwn muto bevtosidct rsneo intrin of presence indicate observations of amount Growing V 10 L ( ϕ 4 10 = R = ) ⊙ 8 -10 a xli h nrnial faint intrinsically the explain can V 5 0 − 9 [1 M . 10 01 − ⊙ 7 cos( z < n ai 10 radii and L ⊙ S ugopwith subgroup QSO ) ϕ/f gal < )] , 0 . )glxe for galaxies 3) 3 -10 d 8 Dtd aur 5 2019) 25, January (Dated: ∼ < 4 ntl .Svidzinsky A. Anatoly − R 50 ⊙ 10 ϕ hyaebr natv aaisadeetdit sur- into ejected and galaxies active in born are They . − 9 with m M 100 (1) ≈ ⊙ e,wihfisi h pnainms window mass axion open the in fits which meV, 1 meV ciigtesl-rvttn field self-gravitating de- equations the Einstein and scribing Klein-Gordon static the units stant, [9] where ihbudr conditions boundary with l,tekyigeinsaetesae fgoa symme- global of scales the are breaking ingredients try key the els, omlgclagmnscntanisms ob nthe in be and to astrophysical of mass exists, its range axion constrain the arguments If mat- cosmological dark is [8]. well-motivated is PNGB candidate and ter hypothetical couplings radiation feeble and light matter extraordinarily a to possess of which axion examples the the of One ymtybekn cl olei window a in lie to scale symmetry-breaking × where of nrdc iesols oriae n en h unit and the time define distance, and of coordinates dimensionless introduce 0 ecnie peial ymti ytmwt metric with system symmetric spherically consider We t hnab aaisaeaincbblswith bubbles axionic are galaxies nearby th . 62 and xo akmte nteUniverse the in matter dark axion /m lebtenzr n nnt.Quantum infinity. and zero between alue g m atrcniae h ubeinterior bubble The candidate. matter k c tssal ntesaemr hnhun- then more scale the on stable ates eew n naayia ouinof solution analytical an find we Here . h ailmti,and metric, radial the , dct rsneo xo akmatter dark axion of presence ndicate eidcitrcinptnil Such potential. interaction periodic h g stesedo light, of speed the is zto fqaa esitadsuggest and redshift quasar of ization ϕ ial an usrsbru afw% few (a subgroup quasar faint sically = e)=2 = (eV) r r 2 & nvriy X77843-4242 TX University, A&M s N ′ m 0 g with √ = ′ ′ ds = = g f V g 10 = 2 2 mc 0 0 = (0) ¯ h n xlctsmer raig( breaking symmetry explicit and /f r g N 2 = 1 + 2 r t , × − en h icmeeta ais We radius. circumferential the being stepril as ndimensionless In mass. particle the is 1 − g N − 10 − g 6 2 r ϕ 0 ( 2 2 r − 9 − ∞ g dt = rg as 2 − 1 2 = ) 3 2 + mc 6 + V + ¯ h × × r g 2 N r G 2 10 N ϕ , 10 + ϕ dr ( ϕ ∞ stegaiainlcon- gravitational the is ′ h as,aefunctions are lapse, the , − 12 2 ′ ϕ 2 g 2 3 ϕ 2 + 1 = ) 2 ′′ + 0 − e [8]. GeV Vadteglobal the and eV n h ercare metric the and − = r 2 g 2 g 2 d ∂V , ∂ϕ √ 2 V Ω V 4 1 2 πG , f 0 = , ≈ , , , 10 V 7 0 GeV ) 1 (3) (6) (4) (5) (2) / 4 . 2

g′(0) = N ′(0) = ϕ′(0) = 0, V (ϕ( ))=0, ∞ where prime denotes ∂/∂r, 1 1 m V = [1 cos(αϕ)], α = = pl (7) α2 − √4πGf √4πf is the dimensionless potential and the coupling parameter respectively, m = ¯hc/G = 1.2 1019 GeV is the pl × Planck mass. The interactionp potential has degenerate minima at ϕ = 2πn/α, where n is an integer number. Here we show that in the limit of strong nonlinearity, α 1, Eqs. (4)-(6) have an approximate static solution≫ that describes a spherical bubble with surface width much smaller then its radius R. The bubble surface is an interface between two degenerate vacuum states with ϕ =2πn/α (r < R) and ϕ = 0 (r > R). Outside the bubble Eqs. (4)-(6) lead to the known FIG. 1: Scalar field ϕ as a function of distance r to the bubble Schwazschild solution center for bubbles with equal radius and different quantum numbers n = 1, 2, 3. The unit of length is ¯h/mc. Note, 1 2M g2 = ,N 2 =1 , (8) we plot the field ϕ only in the vicinity of the bubble surface 1 2M/r − r where it undergoes variation. − 2 where M is the bubble mass in units of mpl/m. Let us assume that R α 1. Then, near the surface where Rm is a position of the mth kink, m =1, 2,...,n. ≫ ≫ one can omit terms with 1/r in Eqs. (4)-(6) and take When the coordinate r passes through the point Rm r R, we obtain the scalar field ϕ(r) changes from 2π(n m + 1)/α to ≈ 2π(n m)/α (see Fig. 1). Eq. (13) yields− the following ′′ − ′ ϕ ∂V expression for g as a function of ϕ inside the bubble: 2Rϕ V + =0, (9) − g2 − ∂ϕ 1 4R =1 [2m 1 + cos(αϕ/2)]. (15) g − α2 − ′ NR ′ N = ϕ 2 2g2V , (10) Outside the bubble ϕ = 0, m = n and 1/g = 1 2 − 8nR/α2. The solution is valid if 1/g > 0, that− is 2 R < Rmax = α /8n. Match of the inner solution (13) ′ gR ′ with the Schwazschild solution (8) determines the mass- g = ϕ 2 +2g2V . (11) 2 radius relation Eqs. (9)-(11) can be solved analytically. Their first inte- M =4πnuR2 8π2n2u2R3, (16) gral is − where u is the surface energy density given by an integral N = const, ϕ′2 =2g2V, g′ = Rgϕ′2. (12) over one potential period u = √2V dϕ/4π. For the cosine potential (7) u =2/πα2. R We assume ϕ(0) = 2πn/α, where n = 1, 2, 3,... is the Redshift of the bubble interior z = 1/N 1 can be number of kinks at the bubble surface, and ϕ(r) mono- found by matching the inner N = const and− the outer tonically decreases with r. Eqs. (12) yield (8) solutions:

ϕ(0) √ 1 1 1 √ ′ 2V z = 1= 1. (17) =1 R 2V dϕ, ϕ = ϕ(0) . 1 2M/R − 1 4πnuR − g − Zϕ −1 R √2V dϕ − − − ϕ p R (13) The internal redshift monotonically increases from zero For V (ϕ) given by Eq. (7) the ﬁnal solution is to inﬁnity when the bubble radius R changes from zero to Rmax. Fig. 2 shows the redshift of space as a function of 4R 4R the distance r to the bubble center. The redshift is con- ln sin(αϕ/2) + 1 (2m 1) arctanh[cos(αϕ/2)] α2 | | − α2 − stant in the bubble interior and monotonically decreases outside the bubble. Let us make rescaling M M/4πu, R R/4πu, then = sign[sin(αϕ/2)](r R ), → → − m Eqs. (16), (17) yield 1 M = nR2 n2R3/2, z = 1. (18) ϕ [2π(n m + 1)/α, 2π(n m)/α], (14) − 1 nR − ∈ − − − 3

n R, in α2¯h/mc z 1 0.0329 0.357 2 0.0241 0.629 3 0.0204 0.96 4 0.0182 1.40 5 0.0168 2.06 6 0.0159 3.24 7 0.0153 6.11 8 0.0151 26.6

TABLE I: Redshift of the bubble interior z and its radius R 2 2 for M = 0.00752α mpl/m and diﬀerent kink numbers n.

FIG. 2: Redshift z of space as a function of distance r to the bubble center for bubbles shown in Fig. 1.

For a given M the redshift depends on the integer number n, which implies the redshift is quantized. In early samples of QSOs associated with nearby spiral galaxies, Karlsson showed that the redshift distribution has a periodicity log(1 + z ) log(1 + z )=0.089, n+1 − n where n = 0, 1, 2,... and z0 = 0.061 [10]. It has been later confirmed by other groups [2, 11]. In a recent paper, Burbidge and Napier [3] tested for the occurrence of this periodicity in new QSO samples and found it to be present at a high confidence level. The peaks were found at z 0.30, 0.60, 0.96, 1.41 and 1.96 in agreement with Karlsson’s≈ empirical formula. The formula also includes the peak at z0 =0.061, however, this peak does not occur FIG. 3: Histogram of the redshift distribution of QSOs close for quasars, but for morphologically related objects. to bright nearby active spiral galaxies or multiple QSOs with The redshift periodicity is observed only in QSO sam- small angular separation from Ref. [12]. The solid lines rep- ples satisfying certain selection criteria, in particular, the resent position of the peaks from Table 1. galaxies which are assumed to be paired to the QSOs must be most nearby spirals [1, 12]. This implies that redshift quantization is a property of intrinsically faint with nearby galaxies are dark matter bubbles composed QSOs which are not detected from large distances. of scalar particles with periodic interaction potential. It is naturally to assume that QSOs born in the same One should mention an alternative possibility of quasar type of galaxies have approximately equal masses because evolution. Bubbles can be originally born with the same their formation mechanism must be similar. Such phe- mass and number of kinks n = 5 that corresponds to nomenon is well known for type Ia supernovae or neu- the 5th peak. During evolution the kinks tunnel to the tron stars: practically all measured neutron star masses bubble center and quasars sequentially decay into states cluster around the value of 1.4M⊙ with only a few per- with smaller n but the same mass. cent deviation [13]. If dark matter bubbles are born with For axions with m = 0.1 3 meV and f = 2 109 equal masses then, according to Eqs. (18), their red- 6 1010 GeV Eq. (7) yields−α = 5.6 107 1.×7 10−9. shift must be quantized. For M = 0.0601 (in dimension Hence,× an axion bubble with the internal× redshift− × z = 2 2 7 units M =0.00752α mpl/m) Eqs. (18) have solutions for 0.36 and n = 1 would have the mass M = 3 10 9 2 4 × − n =1, 2,..., 8, they are given in Table 1. 10 M⊙ and the radius R =3 10 10 R⊙. Such radius In Fig. 3 we plot the most recent histogram of the range agrees with the size of the× emission− region expected redshift distribution from Ref. [12] in which five peaks for the intrinsically faint QSOs. Indeed, for Seyfert 1 are clearly seen. The solid lines show the redshifts from galaxies the size of the broad-line region is R 10 100 our Table 1, they match well the observed peaks. The light days. The luminosity L of the QSOs is∼ 5-6 orders− agreement is remarkable because the theory has only one smaller. Based on the empirical relation for Seyfert 1 free parameter, the bubble mass M. Such coincidence galaxies R L1/2 [14], we obtain for the quasars R 4 ∝ ∼ strongly suggests that the point-like quasars associated 10 R⊙. 4

Now we discuss the bubble life time. Under the in- energy mc2 by a moving bubble surface is governed by the 2 fluence of surface tension and gravitational attraction an Boltzmann factor exp( mc /Teff ), where Teff = a/c is initially static bubble starts to collapse. In the thin-wall the effective temperature− and a is the acceleration of the approximation the initial acceleration is given by [15] surface [19]. For the bubble a c2/R(t) and the Boltz- mann factor reduces to exp( R≈(t)/l) where l =¯h/mc is 2N 3 N 2M the surface width. Hence, emission− of scalar particles is R¨ = , (19) − R − (1 + N)R2 exponentially suppressed apart from small regions where R(t) < l. As a result, during one period of oscillation, where N = 1 2M/R. For one kink thin-wall con- tc R/c∼ , the energy loss is ∆E (l/R)E, which yields − ∼ ∼ tracting bubblep the conserved mass is [15] the bubble life time

2 4πuR 2 2 3 R R2 M = 8π u R , (20) t t = . (22) 1 (dR/dτ)2 − l c cl − ∼ 2 p For an axionic bubble with R > 10 R⊙ and l < 0.7 cm where τ is the interior coordinate time. Based on Eqs. Eq. (22) yields t> 108 yrs which is the time we need to (19), (20) one can expect a continuous contraction of the account for the phenomenon of quasars. bubble to the origin on an astronomically short time scale Properties of the intrinsically faint point-like QSOs, R/c 1yr. However, so far we treated the scalar field ≪ combined with equations for the bubble mass M = as classical. Quantum corrections suppress the collapse 2 2 11 0.00752α mpl/m =2.94m(eV) 10 M⊙ and the radius and result in appearance of long-lived bubbles, stable on 2 × 6 R = 0.0329α ¯h/mc = 2.73m(eV) 10 R⊙, allow us to a scale more then hundreds million years. To include determine the axion mass m. The quasar× luminosity sug- quantum effects it has been suggested that the expres- 3 gests that the bubble radius is larger then 10 R⊙ which sion (20) be interpreted as the canonical hamiltonian of 8 yields m > 0.4 meV and M > 10 M⊙. From the other the bubble at the quantum level [16, 17, 18]. The bub- hand, the quasar ejection from active galaxies implies ble wave function Ψ(R) satisfies the following stationary that the bubble mass M must be much smaller then the quantum mechanical equation in one dimension (¯h = 1) 9 galactic mass. It is reasonable to constrain M < 10 M⊙ [18]: 4 which leads to m < 3 meV and R < 10 R⊙. We con- 2 clude, the axion mass is m =0.4 3 meV. This value fits 2 ∂ E +8π2u2R3 + 16π2u2R4 Ψ(R)=0. (21) in the open window for the axion− mass constrained by as- ∂R2 − trophysical and cosmological arguments [8], which unam- This equation possesses stationary solutions that are not biguously points towards the axionic nature of dark mat- possible in the classical model. Bubbles of non-negligible ter composing the intrinsically faint point-like quasars. redshift correspond to highly excited stationary states Current cavity search experiments in Livermore [20] and for which the energy spectrum can be treated as quasi- Kyoto University [8] are looking for the axion in the mass range 1 10 µeV which deviates by two orders of mag- continuous. At the quantum level the collapse is pre- − vented by quantum pressure that balances the surface nitude from our result. Probably now, when the axion tension and gravitational attraction producing station- mass is established from quasar observations, the axion ary configurations. has a better chance to be discovered. Let us estimate the decay time of an excited station- We mention that data on central “black hole” masses ary state of the quantum bubble. The decay occurs by in small companion galaxies allow us to determine the means of scalar particle emission. We estimate the de- axion mass more accurately and yield m = 1.0 1.9 cay time using the Bohr correspondence principle as the meV. We will discuss this elsewhere. Moreover, observa-− time of energy loss by the classical bubble with the ra- tions show that apart from the intrinsically faint point- dius R(t) oscillating between the turning points R(t)= R like objects considered here there is a subgroup of bright and R(t) = 0, where R is determined by Eq. (16). In quasars which probably also possess noncosmological red- the quantum picture, however, there are no such oscil- shift. Tachyons, another dark matter candidate, can ex- lations. The probability of creation a particle with the plain their nature. We discuss this in a detail paper [21].

[1] K.G. Karlsson, A&A 239, 50 (1990). 391, 833 (2002); A&A 418, 877 (2004). [2] H. Arp et al., A&A 239, 33 (1990). [6] H. Arp, Quasars, redshifts and controversies, Interstellar [3] G. Burbidge and W.M. Napier, AJ 121, 21 (2001). Media, Berkeley, 1987; Seeing red: redshifts, cosmology [4] E. Hawkins, S.J. Maddox and M.R. Merriﬁeld, MNRAS and academic science, Apeiron, Montreal, 1998. 336, L13 (2002). [7] C.T. Hill and G.G. Ross, Nucl. Phys. B311, 253 (1988). [5] G.R. Burbidge et al., ApJS 74, 675 (1990); E.M. Bur- [8] R. Bradley et al., Rev. Mod. Phys. 75, 777 (2003). bidge et al., ApJ 591, 690 (2003); H. Arp et al., A&A [9] E. Seidel and W. M. Suen, Phys. Rev. D 42, 384 (1990). 5

[10] K.G. Karlsson, A&A 13, 333 (1971); 58, 237 (1977). [16] V.A. Berezin et al., Phys. Lett. B 212, 415 (1988). [11] J.M. Barnothy and M.F. Barnothy, PASP 88, 837 (1976). [17] A. Aurilia and E. Spallucci, Phys. Lett. B 251, 39 (1990). [12] W.M. Napier and G. Burbidge, MNRAS 342, 601 (2003). [18] A. Aurilia et al., Phys. Lett. B 262, 222 (1991). [13] N.K. Glendenning, “Compact Stars: Nuclear Physics, [19] A. Gorsky and K. Selivanov, Phys. Rev. D 62, 071702 Particle Physics, and General Relativity”, Springer Ver- (2000). lag; New York, 2nd edition, (2000). [20] S.J. Asztalos et al., Phys. Rev. D 69, 011101(R) (2004). [14] T.G. Wang and X.G. Zhang, MNRAS 340, 793 (2003). [21] A.A. Svidzinsky, astro-ph/0409064. [15] S. K. Blau et al., Phys. Rev. D 35, 1747 (1987).