Physics Letters B 540 (2002) 1–8 www.elsevier.com/locate/npe
Cardassian expansion: a model in which the universe is flat, matter dominated, and accelerating
Katherine Freese, Matthew Lewis
Michigan Center for Theoretical Physics, University of Michigan, Ann Arbor, MI 48109, USA Received 15 January 2002; received in revised form 30 May 2002; accepted 30 May 2002 Editor: J. Frieman
Abstract A modification to the Friedmann–Robertson–Walker equation is proposed in which the universe is flat, matter dominated, and accelerating. An additional term, which contains only matter or radiation (no vacuum contribution), becomes the dominant driver of expansion at a late epoch of the universe. During the epoch when the new term dominates, the universe accelerates; we call this period of acceleration the Cardassian era. The universe can be flat and yet consist of only matter and radiation, and still be compatible with observations. The energy density required to close the universe is much smaller than in a standard cosmology, so that matter can be sufficient to provide a flat geometry. The new term required may arise, e.g., as a consequence of our observable universe living as a 3-dimensional brane in a higher-dimensional universe. The Cardassian model survives several observational tests, including the cosmic background radiation, the age of the universe, the cluster baryon fraction, and structure formation. 2002 Elsevier Science B.V. All rights reserved.
1. Introduction radiation. Pure matter (or radiation) alone can drive an accelerated expansion if the Friedmann–Robertson– Recent observations of type IA Supernovae [1, Walker (FRW) equation is modified by the addition of 2] as well as concordance with other observations a new term on the right-hand side as follows: (including the microwave background and galaxy H 2 = Aρ + Bρn, (1) power spectra) indicate that the universe is accel- ˙ erating. Many authors have explored a cosmologi- where H = R/R is the Hubble constant (as a function cal constant, a decaying vacuum energy [3,4], and of time), R is the scale factor of the universe, the quintessence [5–7] as possible explanations for such energy density ρ contains only ordinary matter and an acceleration. radiation, and we will take Here we propose an alternative which invokes n<2/3. (2) no vacuum energy whatsoever. In our model the universe is flat and yet consists only of matter and In the usual FRW equation B = 0. To be consistent with the usual FRW result, we take 8π E-mail address: [email protected] A = . (M. Lewis). 2 3mpl 0370-2693/02/$ – see front matter 2002 Elsevier Science B.V. All rights reserved. PII:S0370-2693(02)02122-6 2 K. Freese, M. Lewis / Physics Letters B 540 (2002) 1–8
We note here that the geometry is flat, as required The second term starts to dominate at a redshift zeq n by measurements of the cosmic background radiation when Aρ(zeq) = Bρ (zeq), i.e., when [8], so that there are no curvature terms in the = 1−n + 3(1−n) equation. There is no vacuum term in the equation. B/A ρ0 (1 zeq) . (5) This Letter does not address the cosmological constant From evaluating Eq. (1) today, we have (Λ) problem; we simply set Λ = 0. In this Letter, we first study the phenomenology of 2 n H = Aρ0 + Bρ (6) the ansatz in Eq. (1), and then turn to a discussion 0 0 of the origin of this equation.1 Directions for a future so that search for a fundamental theory will be discussed. = 2 − n−1 A H0 /ρ0 Bρ0 . (7) From Eqs. (5) and (7), we have 2. The role of the Cardassian term2 2 3(1−n) H (1 + zeq) B = 0 . (8) The new term in the equation (the second term n 3(1−n) ρ [1 + (1 + zeq) ] on the right-hand side) is initially negligible. It only 0 comes to dominate recently, at the redshift zeq ∼ We have two parameters in the model: B and n, O(1) indicated by the supernovae observations. Once or, equivalently, zeq and n. Note that B here is the second term dominates, it causes the universe chosen to make the second term kick in at the right to accelerate. We can consider the contribution of time to explain the observations. As yet we have ordinary matter, with no explanation of the coincidence problem; i.e., we have no explanation for the timing of zeq.Suchan −3 ρ = ρ0(R/R0) (3) explanation would arise if we had a reason for the to this new term. Once the new term dominates the required mass scale of B. The parameter B has units 2−4n right-hand side of the equation, we have accelerated of mass . Later, we will discuss the origin of the expansion. When the new term is so large that the Cardassian term in terms of extra dimensions, and ordinary first term can be neglected, we find discuss the origin of the mass scale of B. As discussed below, to match the CMB and supernovae data we take 2 R ∝ t 3n (4) 0.3 zeq 1, but this value can easily be refined to better fit upcoming observations. so that the expansion is superluminal (accelerated) for n<2/3. As examples, for n = 2/3wehaveR ∼ t; for n = 1/3wehaveR ∼ t2;andforn = 1/6wehave 3. What is the current energy density of the R ∼ t4. The case of n = 2/3 produces a term in the universe? FRW equation H 2 ∝ R−2; such a term looks similar to a curvature term but is generated here by matter in a Observations of the cosmic background radiation universe with a flat geometry. Note that for n = 1/3the show that the geometry of the universe is flat with acceleration is constant, for n>1/3 the acceleration is Ω = 1. In the Cardassian model we need to revisit diminishing in time, while for n<1/3 the acceleration 0 the question of what value of energy density today, ρ , is increasing (the cosmic jerk). 0 corresponds to a flat geometry. We will show that the energy density required to close the universe is much 1 As discussed below, we were motivated to study an equation smaller than in a standard cosmology, so that matter of this form by work of Chung and Freese [10] who showed that can be sufficient to provide a flat geometry. terms of the form ρn can generically appear in the FRW equation as The energy density ρ0 that satisfies Eq. (6) is, by a consequence of embedding our observable universe as a brane in extra dimensions. definition, the critical density. From Eqs. (1) and (5), 2 The name Cardassian refers to a humanoid race in Star Trek we can write whose goal is to take over the universe, i.e., accelerated expansion. 2 = + 1−n + 3(1−n) n This race looks foreign to us and yet is made entirely of matter. H A ρ ρ0 (1 zeq) ρ . (9) K. Freese, M. Lewis / Physics Letters B 540 (2002) 1–8 3
Fig. 1. The ratio F(n,zeq) = ρc/ρc,old as given by Eq. (13). The contour labeled 0.3 corresponds to parameters n and zeq roughly consistent with present observations.
= 2 3 Evaluating this equation today with A 8π/(3mpl), above, i.e., we have = × × −29 2 3 ρ0 (1/3, 1/5, 0.15) 1.88 10 h0 gm/cm 2 = 8π + + 3(1−n) H0 ρ0 1 (1 zeq) . (10) 2 for n = (2/3, 1/3, 1/6) and zeq = 1. (16) 3mpl
Defining ρ0 = Ω0ρc we find that the critical density For larger values of zeq, the modification to the ρc has been modified from its usual value, i.e., the value of ρc can be even larger. Note the amusing number has changed. We find result that for zeq = 2andn = 1/12, we have ρc = 2 2 0.046ρc,old so that baryons would close the universe 3H0 mpl (not a universe we advocate). ρ = . (11) c 3(1−n) 8π[1 + (1 + zeq) ] Thus 4. Cluster baryon fraction ρc = ρc,old × F(n), (12) where For the past ten years, a multitude of observations has pointed towards a value of the matter density = + + 3(1−n) −1 F(n) 1 (1 zeq) (13) ρ0 ∼ 0.3ρc,old. The cluster baryon fraction [11,12] as and well as the observed galaxy power spectrum are best fit if the matter density is 0.3 of the old critical density. = × −29 2 −3 ρc,old 1.88 10 h0 gm/cm (14) Recent results from the CMB [8,9] also obtain this value. In the standard cosmology this result implied and h0 is the Hubble constant today in units of 100 km s−1 Mpc−1. In Fig. 1, we have plotted the new critical density ρc as a function of the two 3 An alternate possible definition would be to keep the standard parameters n and zeq. For example, if we take zeq = 1, value of ρc and discuss the contribution to it from the two terms on we find the right-hand side of the modified FRW equation. Then there would be a contribution to Ω from the ρ term and another contribution F = (1/3, 1/5, 0.15) from the ρn term, with the two terms adding to 1. This is the approach taken when one discusses a cosmological constant in lieu for n = (2/3, 1/3, 1/6), respectively. (15) of our second term. However, the situation here is different in that we have only matter in the equation. The disadvantage of this second We see that the value of the critical density can be choice of definitions would be that a value of the energy density = much lower than previously estimated. Since Ω0 1, today equal to ρc according to this second definition would not we have today’s energy density as ρ0 = ρc as given correspond to a flat geometry. 4 K. Freese, M. Lewis / Physics Letters B 540 (2002) 1–8
Table 1 fected. Below we discuss the impact on late structure Values of zeq for various values of n corresponding to a universe formation during the era where the Cardassian term = with ρ0 0.3ρc,old. The age of the universe today t0 corresponding is important. This term accelerates the expansion of to the two parameters n and z is listed in the last column, where eq the universe, and freezes out perturbation growth once H0 is the value of the Hubble constant today it dominates (much like when a curvature term dom- nz H t eq 0 0 inates); this freeze out happens late enough that it is 0.60 1.00 0.73 relatively unimportant. To obtain an idea of the type 0.50 0.76 0.78 0.40 0.60 0.83 of effects that we may find, instead of analyzing the 0.30 0.50 0.87 exact perturbation equations with metric perturbations 0.20 0.42 0.92 included, we will merely modify the time dependence 0.10 0.37 0.95 of the scale factor in the usual Jeans analysis equation. 0.00 0.33 0.99 For now we take the standard equation for perturba- tion growth; as a caveat, we warn that recent structure that matter could not provide the entire closure density. formation may be further modified due to a change in Here, on the other hand, the value of the critical Poisson’s equations as described below. For we now density can be much lower than previously estimated. we take Hence the cluster motivated value for ρ0 is now compatible with a closure density of matter, Ω = 1, 0 δ¨ + 2(R/R)˙ δ˙ = 4πρδ/m2 , (17) all in the form of matter. For example, if n = 0.6 with pl z = 1, or if n = 0.2 with z = 0.4, a critical density eq eq where δ = (ρ −¯ρ)/ρ¯ is the fluid overdensity. Now of matter corresponds to ρ ∼ 0.3ρ , as required by 0 c,old one must substitute Eq. (1) for R/R˙ . In the standard the cluster baryon fraction and other data. In Fig. 1, FRW cosmology with matter domination, R ∼ t2/3, one can see which combination of values of n and z eq and there is one growing solution to δ with δ ∼ R ∼ produce the required factor of 0.3. If we assume that t2/3. This standard result still applies throughout most the value ρ = 0.3ρ is correct, for a given value 0 c,old of the (matter dominated) history of the universe in our of n (that is constant in time) we can compute the new model, so that structure forms in the usual way. value of z for our model from Eq. (13). Table 1 lists eq Modifications set in once the new Cardassian term these values of n and z . Henceforth, we shall focus eq becomes important. When R ∼ tp, Eq. (17) can be on these combinations of parameters. written
2 5. Age of the universe 2p 3p − δ (x) + δ (x) − x 3pδ = 0, (18) x 2 In the Cardassian model, the universe is older due to the presence of the second term. In Table 1, we where x ≡ t/t0 with t0 denoting the time today and present the age of the universe for various values of n superscript prime refers to d/dx. This equation can (under the assumption that ρ0 = 0.3ρc,old). generally be solved in terms of Bessel functions for If one takes t0 > 10 Gyr as the lower bound on constant p (such as is the case once the Cardassian globular cluster ages, then one requires t0H0 > 0.66 term completely overrides the old term). A simple = = for h0 = 0.65. If one requires globular cluster ages example is the case of n 2/3andp 1; in greater than 11 Gyr [23], then t0H0 > 0.73 for h0 = the limit x 3/4, the last term in Eq. (18) can 0.65. All values in Table 1 satisfy these bounds. be dropped and the equation is solved as δ(t) = −1 a1 + a2t . Perturbations cease growing and become frozen in. This result agrees with the expectation 6. Structure formation that in a universe that is expanding more rapidly, the overdensity will grow more slowly with the scale Since the new (Cardassian) term becomes impor- factor. As mentioned at the outset, as long as the tant only at zeq ∼ 1, early structure formation is not af- Cardassian term becomes important only very late K. Freese, M. Lewis / Physics Letters B 540 (2002) 1–8 5 in the history of the universe, much of the structure with F given in Eq. (13). As discussed previously, we see will have already formed and be unaffected. as our standard case we will take F ≡ ρc/ρc,old = Further comments on late structure formation (e.g., 0.3. With this assumption, and by using expression cluster abundances) follow below. Eq. (24) in Eq. (22), we find that d changes by a factor of (1.47, 1.88, 2.04 and 2.23) for n = (0.6, 0.3, 0.2 and 0.1), respectively, compared to the usual 7. Doppler peak in cosmic background radiation (nonCardassian) FRW universe with ρ0 = ρc,old.In addition Here we argue that the location of the first Doppler √ √ 1 1 + R∗ + R∗ + r∗R∗ peak is only mildly affected by the new Cardassian s∗ ∝ √ ln √ , cosmology. We need to calculate the angle subtended F 1 + r∗R∗ by the sound horizon at recombination. In the standard − where r∗ = 0.042(F h2) 1 and R∗ = 30Ω h2 and we FRW cosmology with flat geometry, this value corre- b use h = 0.7andΩ = 0.04. We find that s∗ changes sponds to a spherical harmonic with l = 200. A peak b by a factor of (1.44, 1.62, 1.67, 1.29) for n = (0.6, at this angular scale has indeed been confirmed [8]. In 0.3, 0.2 and 0.1), respectively, compared to the usual the Cardassian cosmology we still have a flat geome- FRW universe with ρ = ρ . The angle subtended try. Hence, we can still write 0 c,old by the sound horizon on the surface of last scattering θ = s∗/d, (19) decreases and the location (l) of the first Doppler peak increases by roughly a factor of where s∗ is the sound horizon at the time of recom- bination tr and d is the distance a light ray trav- (1.02, 1.11, 1.12, 1.13) els from recombination to today. To calculate these n = , , lengths, we use the fact that for a light ray ds2 = 0 = for (0.6, 0.3, 0.2, 0.1) respectively (25) − 2 + 2 2 dt a dx to write compared to the usual FRW universe with ρ0 = ρc,old. t0 This shift still lies within the experimental uncertainty d = dt/a. (20) on measurements of the location of the Doppler peak. We note the following: in the same way that a t r nonzero Λ may make the current CBR observations Following the notation of Peebles [13], we define the compatible with a small but nonzero curvature, indeed redshift dependence of H as a nonzero Cardassian term could also allow for a nonzero curvature in the data. A more accurate study H(z)= H E(z) (21) 0 of the effects of Cardassian expansion on the cosmic so that Eq. (20) can be written background radiation (including the first and higher
zr peaks) is the subject of a future study. 1 dz d = . (22) H0R0 E(z) 0 8. The cutoff energy density Similarly, the sound horizon at recombination is ∞ An alternate way to write Eq. (1) is = s∗ dt/a. (23) 2 n−1 H = Aρ 1 + (ρ/ρcutoff) , (26) zr ≡ = In standard matter dominated FRW cosmology where ρcutoff ρ(zeq) A/B. This notation of- 3/2 fers a new interpretation; it indicates that the sec- with Ωm,0 = 1, E(z) = (1 + z) in Eq. (22) and ond term only becomes important once the energy d = 2/H0R0. For the cosmology of Eq. (1), we have density of the universe drops below ρcutoff,which has a value a few times the critical density. Hence, E(z)2 = F × (1 + z)3 + (1 − F)× (1 + z)3n (24) regions of the universe where the density exceeds 6 K. Freese, M. Lewis / Physics Letters B 540 (2002) 1–8 this cutoff density will not experience effects asso- 10. Best fit of parameters to current data ciated with the Cardassian term. In particular, we can be reassured that the new term will not affect We can find the best fit of the Cardassian parame- gravity on the Earth or the Solar System. The den- ters n and zeq to current CMB and Supernova data. 3 sity of water on the Earth is 1 gm/cm ,whichis The current best fit is obtained for ρ0 = 0.3ρc,old (as 28 orders of magnitude higher than the critical den- we have discussed above) and n<0.4 (equivalently, sity. w<−0.6) [20,21]. In Table 1 one can see the values of zeq compatible with this bound, as well as the resul- tant age of the universe. As an example, for n = 0.2 (equivalently, w =−0.8), we find that z = 0.42. 9. Comparing to quintessence eq Then the position of the first Doppler peak is shifted by a factor of 1.12. The age of the universe is 13 Gyr. We note that, with regard to observational tests, The cutoff energy density is ρcutoff = 2.7ρc,sothat [14–18], one can make a correspondence between the new term is important only for ρ<ρcutoff = 2.7ρc. the Cardassian and quintessence models; we stress, Hence, as mentioned above, the Cardassian term will however, that the two models are entirely different. not affect the physics of the Earth or Solar System in Quintessence requires a dark energy component with a any way. specific equation of state (p = wρ), whereas the only We note the enormous uncertainty in the current ingredients in the Cardassian model are ordinary mat- data; future experiments (such as SNAP [22]) will ter (p = 0) and radiation (p = 1/3). However, as far constrain these parameters further. as any observation that involves only R(t), or equiv- alently H(z), the two models predict the same effects on the observation. Regarding such observations, we 11. Extra dimensions can make the following identifications between the Cardassian and quintessence models: n ⇒ w+1, F ⇒ A Cardassian term may arise as a consequence + Ωm,and1− F ⇒ ΩQ,wherew is the quintessence of embedding our observable universe as a (3 1)- equation of state parameter, Ωm = ρm/ρc,old is the ra- dimensional brane in extra dimensions. Chung and tio of matter density to the (old) critical density in the Freese [10] showed that, in a 5-dimensional universe standard FRW cosmology appropriate to quintessence, with metric = ΩQ ρQ/ρc,old is the ratio of quintessence energy ds2 =−q2(τ, u) dτ 2 + a2(τ, u) dx2 density to the (old) critical density, and F is given by + 2 2 Eq. (13). In this way, the Cardassian model can make b (τ, u) du , (27) contact with quintessence with regard to observational where u is the coordinate of the fifth dimension, one tests. may obtain a modified FRW equation on our observ- All observational constraints on quintessence that able brane with H 2 ∼ ρn for any n (see also [19]). depend only on the scale factor, R(t) (or, equivalently, This result was obtained with 5-dimensional Einstein H(z)) can also be used to constrain the Cardassian equations plus the Israel boundary conditions relating model. However, because the Cardassian model re- the energy–momentum on our brane to the derivatives quires modified Einstein equations (see below), the ofthemetricinthebulk. gravitational Poisson’s equations and consequently We do not yet have a fundamental higher-dimen- late-time structure formation may be changed; e.g., the sional theory, i.e., a higher-dimensional Tµν ,which redshift dependence of cluster abundance should be we believe describes our universe. Once we have this, different in the two models. These effects (and others, we can write down the modified four-dimensional Ein- such as the fact that quintessence clumps) may serve stein’s equations and compute the modified Poisson’s to distinguish the Cardassian and quintessence mod- equations, as would be required, e.g., to fully under- els. The correspondence with quintessence, as well as stand latetime structure formation. discussion of distinguishing tests will be the subject of There is no unique 5-dimensional energy–momen- a future paper. tum tensor Tµν that gives rise to Eq. (1) on our brane. K. Freese, M. Lewis / Physics Letters B 540 (2002) 1–8 7
Hence, in this Letter we construct an example which 10−101, which cancels other small numbers in such a is easy to find but is clearly not our universe, simply way as to again require roughly Eq. (30) to be satis- as a proof that such an example can be written down. fied. The form of Tµν given in Eq. (28) is by no means Following [10] (see Eqs. (24) and (25) there with unique and has been presented merely as an existence F(u)= u), we have constructed an example of a bulk proof; we hope a more elegant Tµν may be found, per- 2 n Tµν for arbitrary n in H ∼ ρ , matter on the brane as haps with a motivation for the required value of B. in Eq. (3), and with q = b in Eq. (27). We display only 0 T0 here (the other components will be published in a future paper): 12. Discussion − 4+n − 2 3 n B n , κ2T 0 =− We have presented Eq. (1) as a modification to 5 0 n2τ 2 the FRW equations in order to suggest an explana- 4 n tion of the recent acceleration of the universe. In the 1 2 1 4 1 2 2 2 × 4 · 81 n B n − 16 n κ n τ + u , 5 nτ Cardassian model, the universe can be flat and yet matter dominated. We have found that the new Car- (28) dassian term can dominate the expansion of the uni- where verse after zeq = O(1) and can drive an acceleration. 2 + n We have found that matter alone can be responsible 2 n − 1 2 1 , = exp −(2/3) n B n κ u (29) for this behavior (but see the comments below). The 5 nτ current value of the energy density of the universe is and the constant κ5 is related to the 5-dimensional then smaller than in the standard model and yet is at Newton’s constant G5 and 5D reduced Planck mass the critical value for a flat geometry. Structure forma- 2 = = −3 M5 by the relation κ5 8πG5 M5 .Thisis tion is unaffected before zeq. The age of the universe merely one (inelegant) example of many bulk Tµν that is somewhat longer. The first Doppler peak of the cos- produce Cardassian expansion. mic background radiation is shifted only slightly and We may now investigate the meaning of the values remains consistent with experimental results. Such a of B(n) required by Eq. (8), where B(n) is the modified FRW equation may result from the existence parameter in front of the new Cardassian term in of extra dimensions. Further work is required to find a Eq. (1). As mentioned previously, the mass scale of B simple fundamental theory responsible for Eq. (1). − n has units of m2 4 . We find that the corresponding Questions of interpretation remain. We have said mass scale is very small for n<1/2, is singular at n = that matter alone is responsible for the accelerated 1/2, and then goes over to a very large value for n> behavior. However, if the Cardassian behavior results 1/2. Specifically, for n = 2/3andzeq = 1, we obtain from integrating out extra dimensions, then one may − − B ∼ 10 52 GeV 2/3 which corresponds to a mass ask what behavior of the radii of the extra dimensions scale of 1078 GeV. In the context of extra dimensions, is required. The Israel conditions connect the energy this large mass scale turns out to cancel against other density on the brane to fields in the bulk. The required large numbers in such a way that it corresponds to behavior of bulk fields is not transparent when one reasonable values of the energy–momentum tensor in writes the modified FRW equation. We have found the bulk. We find that τ is roughly the age of the a large or small mass scale to be required, which universe and we have , ∼ 1forallu.Thenwehave must result from the extra dimensions. In principle − one would like to have a complete 5-dimensional T 0 ∼ 10 5 GeV 5. (30) 0 theory so as to perform post-Newtonian tests on the Although this value is not motivated, it is not unrea- model and also to check other consequences. For sonable. In other words, reasonable bulk values can example, with a 5-dimensional model, one would like generate the required parameters in Eq. (1). Numeri- to compare with limits from fifth force experiments cal values for other components of Tµν are the same and to check that none of the higher-dimensional order of magnitude, with the exception that T04 ∼ 0. fields are overcontributing to the energy density of the For the case of n = 1/3, we obtain a mass scale of universe at any point (the moduli problem). 8 K. Freese, M. Lewis / Physics Letters B 540 (2002) 1–8
One might attempt to use a Cardassian term (the [3] K. Freese, F.C. Adams, J.A. Frieman, E. Mottola, Nucl. Phys. second term in Eq. (1)) to drive an early inflationary B 287 (1987) 797. era in the universe as well. For n<1/2 one could [4] J. Frieman, C. Hill, A. Stebbins, I. Waga, Phys. Rev. Lett. 75 (1995) 2077. have a superluminal expansion during the radiation [5] L. Wang, P. Steinhardt, Astrophys. J. 508 (1998) 483. dominated era. However, once accelerated expansion [6] R. Caldwell, R. Dave, P. Steinhardt, Phys. Rev. Lett. 80 (1998) begins, Eq. (1) without a potential provides no way for 1582. inflation to stop. Hence we have focused on using this [7] G. Huey, L. Wang, R. Dave, R. Caldwell, P. Steinhardt, Phys. new term to generate acceleration today rather than to Rev. D 59 (1999) 063005. [8] C.B. Netterfield et al., astro-ph/0104460; cause an inflationary mechanism early on. However, R. Stompor et al., astro-ph/1015062; it may be possible to combine such a term with a N.W. Halverson et al., astro-ph/0104489. different way to end inflation. [9] C. Pryke et al., Astrophys. J. 568 (2002) 46. [10] D.J. Chung, K. Freese, Phys. Rev. D 61 (2000) 023511. [11] S.D.M. White, J.F. Navarro, A.E. Evrard, C.S. Frenk, Na- ture 366 (1993) 429. Acknowledgements [12] A.E. Evrard, Mon. Not. R. Astron. Soc. 292 (1997) 289. [13] P.J.E. Peebles, Principles of Physical Cosmology, Princeton K.F. thanks Ted Baltz, Daniel Chung, Richard Univ. Press, Princeton, NJ, 1993. [14] M. Davis, J. Newman, Astrophys. J. 513 (1999) L95. Easther, Gus Evrard, Paolo Gondolo, Wayne Hu, [15] D. Huterer, M. Turner, Phys. Rev. D 64 (2001) 123527. Lam Hui, Will Kinney, Risa Wechsler, and especially [16] C. Alcock, B. Paczynski, Nature 281 (1979) 358. Jim Liu for many useful conversations and helpful [17] L. Hui, A. Stebbins, S. Burles, Astrophys. J. 511 (1999) L5. suggestions. [18] P. McDonald, J. Miralda-Escude, Astrophys. J. 518 (1999) 24. We acknowledge support from the Department of [19] P. Binetruy, C. Deffayet, D. Langlois, Nucl. Phys. B 565 (2000) Energy via the University of Michigan. 269. [20] R. Bean, S. Hansen, A. Mechiorri, astro-ph/0201127. [21] S. Hannestad, E. Mortsell, astro-ph/0205096. [22] C. Akerlof et al., Supernova/Acceleration Probe (SNAP): An References Experiment to Measure the Properties of the Accelerating Universe, Resource Book on Dark Energy; C. Akerlof et al., Contributions from the Snowmass 2001, [1] S. Perlmutter et al., Supernova Cosmology Project Collabora- Workshop on the Future of Particle Physics, Snowmass, CO, tion, Astrophys. J. 517 (1999) 565, astro-ph/9812133. 2001. [2] A.G. Riess et al., Supernova Search Project Collaboration, [23] L.M. Krauss, B. Chaboyer, astro-ph/0111597. Astron. J. 116 (1998) 1009, astro-ph/9805201. Physics Letters B 540 (2002) 9–13 www.elsevier.com/locate/npe
Black hole entropy without brick walls
Li Xiang a,b
a CCAST (World Lab.), P.O. Box 8730, Beijing 100080, PR China b Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2735, Beijing 100080, PR China 1 Received 12 April 2002; received in revised form 5 May 2002; accepted 5 May 2002 Editor: J. Frieman
Abstract The properties of the thermal radiation are discussed by using the new equation of state density motivated by the generalized uncertainty relation in the quantum gravity. There is no burst at the last stage of the emission of a Schwarzschild black hole. When the new equation of state density is utilized to investigate the entropy of a scalar field outside the horizon of a static black hole, the divergence appearing in the brick wall model is removed, without any cutoff. The entropy proportional to the horizon area is derived from the contribution of the vicinity of the horizon. 2002 Elsevier Science B.V. All rights reserved.
PACS: 04.70.Dy; 04.62.+v; 97.60.Lf
The title is the same as Ref. [1] where Demers et by introducing some regulators. These fictitious fields al. show that the divergence appearing in the brick are especially designated in the number, statistics and wall model [2] can be absorbed into the renormal- masses. To my surprise, the entropy expressed by the ized Newton’s constant. By using the WKB approx- masses of the regulators can be precisely renormalized imation, ’t Hooft investigates the statistical properties to the Bekenstein–Hawking formula, S = A/(4GR), of a scalar field outside the horizon of a Schwarzschild GR is the renormalized Newton’s constant. However, black hole. The entropy proportional to the horizon it is hard to understand the introduction of the “bare area is obtained, but with a cutoff utilized to remove entropy” in Ref. [1]. The “bare entropy” seems to be the divergence of the density of states. The cutoff is negative and its meaning is unclear.2 Is there a better introduced by hand and looks unnatural. Susskind and method can remove the divergence appearing in the Uglum suggest that the explosive free energy and en- brick wall model? tropy in the model of ’t Hooft are related to the diver- Recently, many efforts have been devoted to the gence of the one-loop effective action of the quantum generalized uncertainty relation field theory in curved space [3]. Their conjecture is confirmed by [1]. The authors of [1] remove the cutoff λ x p h¯ + ( p)2, (1) and regularize the divergent free energy and entropy h¯
E-mail address: [email protected] (L. Xiang). 2 Dr. Fursaev told me, this difficulty can be overcome in the 1 This is the mailing address. Sakharov’s induced gravity [4].
0370-2693/02/$ – see front matter 2002 Elsevier Science B.V. All rights reserved. PII:S0370-2693(02)02123-8 10 L. Xiang / Physics Letters B 540 (2002) 9–13 and its consequences [5–11], especially the effect on then the density of states [10,11]. Here h¯ is the Planck −4 constant, λ is of order of the Planck length. Eq. (1) u = β G(0) + G (0)a √ 4 2 means that there is a minimal length, 2 λ. As well- π − 40π λ = β 4 1 − . (6) known, the number of quantum states in the integrals 15 7β2 d3x d3p is given by In the usual case, above equation does not essentially d3x d3p , (2) change the well-known conclusion for the black body (2πh)¯ 3 radiation because the correction is very slight. For ex- which can be understood as follows: since the uncer- ample, the temperature of the center of the neutron star 9 32 2 tainty relation x p ∼ 2πh¯ , one quantum state cor- is 10 K, but the Planck temperature is 10 K,λ/β ∼ −46 responds to a “cell” of volume (2πh)¯ 3 in the phase- 10 . However, Eq. (6) is no longer valid for the case 2 space. Based on the Liouville theorem, the authors λ/β 1, that is higher than the Planck temperature. of Ref. [11] argue that the number of quantum states We calculate the upper bound of energy density, that is should be changed to the following ∞ 2 − x dx 3 3 4 d x d p u<β 2 , (3) 1 + λx 3 (2πh)¯ 3(1 + λp 2)3 0 β2 −3/2 2 = i = − π λ where p pi p ,i 1, 2, 3. Eq. (3) seriously de- = β 4 · forms the Planckian spectrum of the black√ body ra- 16 β2 diation at the Planck temperature, T = 1/λ (see π − λ = β 1, (7) Ref. [11], Fig. 2). 16λ3/2 Let us discuss the more details than Ref. [11]. This where the inequality is due to ex − 1 >x.This will benefit the following investigation of the black means that when the temperature is higher than the hole entropy. From Eq. (3), we directly write down the Planck temperature the state equation of the thermal density of internal energy of the thermal radiation radiation is essentially different from the well-known ∞ conclusion, u ∼ β−4. This will influence the emission ω3 dω u = of the black hole. According to the Stefan–Boltzmann βω − + 2 3 (e 1)(1 λω ) law, the loss mass rate of a Schwarzschild black hole 0 ∞ reads 3 −4 x dx = β dM −4 1 x 2 3 ∼ β A ∼ , (8) (e − 1)(1 + ax ) 2 0 dt M − = β 4G(a), (4) where M the mass of the hole. At the last stage of emission, M → 0, so the emission rate becomes 2 where a = λ/β ,x = βω. We take the units G = c = divergent. However, from Eq. (7), at the last stage, the h¯ = kB = 1. The above integral cannot be expressed as rate will be changed to a simple formula, but we can investigate its asymptotic dM − behavior in the two different conditions. We first con- ∼ β 1A ∼ M → 0, (9) sider the case a 1. This means that the temperature dt is much less than the Planck temperature. We have here is no burst. ∞ We turn to the problem of black hole entropy. Re- x3 dx π4 G(0) = = , calling the brick wall model, the number of quantum ex − 1 15 states less than energy ω is given by [2,12,13] 0 ∞ L 5 6 3 2 x dx 24π 2ω r dr G (0) =−3 =− , (5) Γ(ω)= , (10) ex − 1 63 3π f 2 0 r0+" L. Xiang / Physics Letters B 540 (2002) 9–13 11 which is for a massless scalar field in a spherical and From Eq. (3), the number of quantum states with static space–time as follows energy less than ω is given by 2 =− 2 + −1 2 + 2 2 + 2 2 2 1 dr dθ dϕdpr dpθ dpϕ ds fdt f dr r dθ r sin θdφ , g(ω) = 3 + 2 3 (11) (2π) (1 λω /f ) 1 dr dθ dϕ where f = f(r). The horizon is located by f(r0) = 0. = (2π)3 (1 + λω2/f )3 " is the cutoff near the horizon. Obviously, the number of states is divergent if " = 0. We carefully check 2 ω2 1 1 1/2 × − p2 − p2 the derivation of Eq. (10) and find that it agrees with f 1/2 f r2 θ r2 sin2 θ ϕ Eq. (2), not (3). The former leads to the following × dpθ dpϕ formula 4πω3 r2 dr 8π3 r2 dr = sin θdθdϕ S = , (12) 3(2π)3 f 2(1 + λω2/f )3 45β3 f 2 2ω3 r2 dr = , (18) which is analogous with√ the usual state equation of the 2 + 2 3 −1 3π f (1 λω /f ) thermal radiation:√ (β f) is the local temperature, 4πr2 dr/ f is the element of the spatial volume of where the integration goes over those values of pθ ,pϕ the spherical shell. The divergent entropy means the for which the argument of the square root is positive = invalidity of the usual state equation near the black (please refer to Refs. [2,12,13]). When λ 0, Eq. (18) = hole horizon. If we take Eq. (3), the situation may be naturally returns to (10). However, in the case λ 0, essentially different. Why not have a try? Substituting Eq. (18) is essentially different from (10): it is conver- the wave function Φ = exp(−iωt)ψ(r,θ,ϕ) into the gent at the horizon without any cutoff! By using the usual method, the free energy is given by equation of massless scalar field 1 √ = 1 − −βω √ ∂ −ggµν ∂ Φ = 0, (13) F(β) dg(ω)ln 1 e −g µ ν β ∞ we obtain =− g(ω)dω βω − ∂2ψ f 2 ∂ψ e 1 + + 0 ∂r2 f r ∂r ∞ 2 r2 dr ω3 dω 1 ω2 1 ∂2 ∂ 1 ∂2 =− . + + + cotθ + ψ 2 βω − + 2 3 2 2 2 2 3π f (e 1)(1 λω /f ) f f r ∂θ ∂θ sin θ ∂ϕ r0 0 (19) = 0. (14) The entropy reads By using the WKB approximation with ∂F S = β2 ψ ∼ exp[iS(r,θ,φ)], ∂β ∞ we have 2β2 r2 dr eβωω4 dω = 2 3π f 2 (eβω − 1)2(1 + λω2/f )3 2 1 ω 1 2 1 2 = − − r0 0 pr 2 pθ 2 pϕ , (15) f f r r2 sin θ ∞ 2β−3 r2 dr x4 dx where = 2 , 3π f 2 (1 − e−x )(ex − 1) 1 + λx 3 ∂S ∂S ∂S β2f p = ,p= ,p= . (16) r0 0 r ∂r θ ∂θ ϕ ∂ϕ (20) where x = βω. Taking into account the following We also obtain the square module of momentum inequalities ω2 2 = i = rr 2 + θθ 2 + ϕϕ 2 = −x x p pi p g pr g p g pϕ . (17) 1 − e > , θ f 1 + x 12 L. Xiang / Physics Letters B 540 (2002) 9–13
ex − 1 >x. (21) near the horizon. This convergency is due to the effect of the generalized uncertainty relation on the quantum We obtain states. This provides an evidence for the idea of Li ∞ 2β−3 r2 dr (x3 + x2)dx and Liu. The more details between the Li–Liu equa- S< tion and the generalized uncertainty relation will be 3π f 2 + λx2 3 1 2 r0 0 β f investigated in the future. − − 2β−3 r2 dr 1 λ 2 π λ 3/2 As pointed by Ref. [6], the generalized uncertainty = + relation (1) may have a dynamical origin since it π 2 2 2 3 f 4 β f 16 β f contains a dimensional coupling constant λ.IfEq.(1) r0 − is indeed due to the string theory, λ should be β λ 3/2 r2 dr = 2 + associated with the stringy scale l2. This implies 2 r dr 1/2 . (22) s 6πλ 24 f that one takes into account the contribution from r0 r0 the stringy excitation when calculating the density of We are only interested in the contribution from the quantum states. The convergent property should be [ + ] vicinity near the horizon, r0,r0 " , which corre- reexamined. There are some evidences (or arguments) sponds to√ a proper distance of order of the minimal for the convergence of the density of states even if length, 2 λ. This is because the entropy closes to the considering the stringy excitation: firstly, the bosonic upper bound only in this vicinity. Furthermore, it is string can be described by a discrete field theory, then just the vicinity neglected by brick wall model. We the number of degrees of freedom of it is smaller have than that of the usual field theory [16]. Secondly, the + + r0 " r0 " entropy of a string is proportional to its mass since √ dr dr 2 λ = √ ≈ √ the degeneracy increases exponentially with the mass − f 2κ(r r0) level. However, the massive string cannot be excited r0 r0 in the low energy effective theory (such as general 2" = , (23) relativity) [17]. Therefore, the contribution from the κ stringy excitation is ignored in the case of the massive where κ is the surface gravity at the horizon of black hole where the semi-classical approximation is black hole and it is identified as κ = 2πβ−1. Thus still valid. As to the black hole at the Planck scale we naturally derive the entropy proportional to the the usual quantum field theory is no longer valid. The horizon area entropies of the black hole and the excited string states β λ−3/2 √ are matched in the correspondence principle [18]. S ∼ r2" + · 2r2 λ 6πλ2 0 24 0 3A = , (24) Acknowledgements 16λπ = 2 where A 4πr0 is the surface area of the black hole. The author thanks X.J. Wang, Z.Q. Bai and H.G. As early as 1992, Li and Liu phenomenally pro- Luo for their zealous help during the research. The posed that the state equations of the thermal radiation author also thank Profs. Y.Z. Zhang, R.G. Cai and near the horizon should be changed to a series of new C.G. Huang for their comments on this research. This formulae rather than Eq. (12), in order to maintain the work is supported by the Post Doctor Foundation of validity of the generalized second law of thermody- China and K.C. Wong Education Foundation, Hong namics [14]. Using the Li–Liu equation, Wang inves- Kong. tigates the entropy of a self-gravitational radiation sys- tem and obtains the Bekenstein’s entropy bound [15]. Here, Parallel to the brick wall model, the scalar field References near the horizon of a static black hole is investigated again, we obtain the entropy proportional to the hori- [1] J.G. Demers, R. Lafrance, R.C. Myers, Phys. Rev. D 52 (1995) zon area. There is no divergence without any cutoff 2245. L. Xiang / Physics Letters B 540 (2002) 9–13 13
[2] G. ’t Hooft, Nucl. Phys. B 256 (1985) 727. [10] S.K. Rama, Phys. Lett. B 519 (2001) 103. [3] L. Susskind, J. Uglum, Phys. Rev. D 50 (1995) 2700. [11] L.N. Chang et al., hep-th/0201017. [4] V.P. Frolov, D.V. Fusaev, A.I. Zelnikov, Nucl. Phys. B 486 [12] T. Padmanabhan, Phys. Lett. B 173 (1986) 43. (1997) 339, hep-th/9607104; [13] J.L. Jing, Int. J. Theor. Phys. 37 (1998) 1441. V.P. Frolov, D.V. Fursaev, Phys. Rev. D 56 (1997) 2212, hep- [14] L.X. Li, L. Liu, Phys. Rev. D 46 (1992) 3296. th/9703178; [15] D.X. Wang, Phys. Rev. D 50 (1994) 7385. V. Frolov, D. Fursaev, J. Gegenberg, G. Kunstatter, Phys. Rev. [16] I. Klebanov, L. Susskind, Nucl. Phys. B 309 (1988) 175. D 60 (1999) 024016. [17] A.W. Peet, Class. Quantum Grav. 15 (1998) 329. [5] A. Kempf, G. Mangano, R.B. Mann, Phys. Rev. D 52 (1995) [18] L. Susskind, hep-th/9309145; 1108. E. Halyo, A. Rajaraman, L. Susskind, Phys. Lett. B 392 (1997) [6] L.J. Garay, Int. J. Mod. Phys. A 10 (1995) 145. 319; [7] H.A. Kastrup, Phys. Lett. B 413 (1997) 267. G.T. Horowitz, J. Polchinski, Phys. Rev. D 55 (1997) 6189. [8] D.V. Ahluwalia, Phys. Lett. A 275 (2000) 31, gr-qc/0002005. [9] R.J. Adler et al., gr-qc/0106080. Physics Letters B 540 (2002) 14–19 www.elsevier.com/locate/npe
Implications of the first neutral current data from SNO for solar neutrino oscillation
Abhijit Bandyopadhyay a, Sandhya Choubey b, Srubabati Goswami c,D.P.Royd,e
a Saha Institute of Nuclear Physics, Bidhannagar, Kolkata 700 064, India b Department of Physics and Astronomy, University of Southampton, Highfield, Southampton S017 1BJ, UK c Harish-Chandra Research Institute, Chhatnag Road, Jhusi, Allahabad 211-019, India d Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400 005, India e Physics Department, University of California, Riverside, CA 92521, USA Received 1 May 2002; received in revised form 13 June 2002; accepted 15 June 2002 Editor: P.V. Landshoff
Abstract We perform model independent and model dependent analyses of solar neutrino data including the neutral current event rate from SNO. The inclusion of the first SNO NC data in the model independent analysis determines the allowed ranges of 8Bflux normalisation and the νe survival probability more precisely than what was possible from the SK and SNO CC combination. We perform global νe–νactive oscillation analyses of solar neutrino data using the NC rate instead of the SSM prediction for the 8B flux, in view of the large uncertainty in the latter. The LMA gives the best solution, while the LOW solution is allowed only at the 3σ level. 2002 Published by Elsevier Science B.V.
= NC CC+NC The neutral current results from the Sudbury Neu- where r σνµ,τ /σνe 0.157 for a threshold en- trino Observatory measures for the first time the to- ergy of 5 MeV (including the detector resolutions tal flux of 8B neutrinos coming from the Sun [1]. In and the radiative corrections to ν–e scattering cross- a recent paper [2] we had examined the role of the sections). All the rates are defined with respect to anticipated NC data from SNO in enhancing our un- the BBP2000 Standard Solar Model (SSM) [4]. We derstanding of the solar neutrino problem. The SNO showed in [2] that because SNO has a greater sensi- NC rate can be expressed in terms of SNO CC and SK tivity to the NC scattering rate as compared to SK, the elastic scattering rates as [3] SNO NC measurement will be more precise and hence incorporation of this can be more predictive than the NC = CC + el − CC SNO CC and SK combination. We took three repre- RSNO RSNO RSK RSNO /r, (1) SNO = ± sentative NC rates, RNC 0.8, 1.0 and 1.2 ( 0.08) and showed that
E-mail addresses: [email protected] ν (A. Bandyopadhyay), [email protected] (1) For a general transition of e into a mixture of (S. Choubey), [email protected] (S. Goswami), active and sterile neutrinos the size of the sterile [email protected] (D.P. Roy). component can be better constrained than before.
0370-2693/02/$ – see front matter 2002 Published by Elsevier Science B.V. PII:S0370-2693(02)02138-X A. Bandyopadhyay et al. / Physics Letters B 540 (2002) 14Ð19 15
(2) For transition to a purely active neutrino the where Pee and Pea denote the probabilities folded with 8B neutrino flux normalisation and the survival the detector response function [5] and averaged over probability Pee are determined more precisely. energy. To extract a model independent bound on Pee (3) We had also performed global two flavour oscilla- one has to ensure an equality of the response functions tion analysis of the solar neutrino data for the νe– which amounts to slight adjustment of the SK thresh- CC νactive case, where instead of RSK and RSNO we old energy and the rate [5,6]. Our approach is slightly el NC CC NC different. We treat P to be effectively energy inde- used the quantities RSK/RSNO and RSNO/RSNO. ee These ratios are independent of the 8B flux nor- pendent. The SK spectrum data indicates a flat proba- malisation and hence of the SSM uncertainty. We bility down to 5 MeV [7]. This is corroborated by SNO showed that use of these ratios can result in drastic [8,9] which now has a threshold of 5 MeV for kinetic reduction of the allowed parameter regions spe- energy of the observed electron. Hence we consider cially in the LOW-QVO area depending on the this assumption as justified and expect the results to be value of the NC rate. insensitive to the differences in the response functions. It should be noted here that in contrast to the SNO CC We now have the actual experimental result events their NC events correspond to a neutrino en- ergy threshold of 2.2 MeV. However it is clear from RNC = 1.01 ± 0.12 (2) Eq. (5) that for a νe to νa transition there is no rea- SNO NC son to expect any energy dependence in RSNO.Onthe while Eq. (1) gives 1.05 ± 0.15. Thus in 306 live days (577 days) the SNO NC measurement has achieved a precision, which is already better than that obtained from the SK and SNO CC combination. This Letter follows closely the analysis that we have done in [2] but incorporating the actual data. In addition we also perform an alternative global analysis for νe–νactive 8 oscillation by letting the B normalisation factor fB NC = vary freely, where the inclusion of RSNO ( fB) in the fit serves to control this parameter. As we shall see below the two methods of global analysis give very similar results. In Section 1 we discuss the constraints on the elec- tron neutrino survival probability, the 8B normalisa- tion factor fB and the fraction of sterile component without assuming any particular model for the proba- bilities. In Section 2 we perform the global analyses assuming two flavour νe–νactive oscillation.
1. Model independent analysis
For the general case of νe transition into a combina- tion of νactive (νa)andνsterile (νs ) states one can write the SK, SNO CC and SNO NC rates as
el = + RSK fBPee fBrPea, (3) CC = Fig. 1. The SNO CC and NC rates shown relative to their SSM RSNO fBPee, (4) predictions. The dashed line is the prediction of the pure νe to νs NC = + transition. The pure sterile solution is seen to be disfavored at 5.3 σ. RSNO fB(Pee Pea), (5) 16 A. Bandyopadhyay et al. / Physics Letters B 540 (2002) 14Ð19
8 Fig. 2. Best fit value of the B neutrino flux fB shown along with its 8 1σ and 2σ limits against the model parameter sin2 α, representing Fig. 3. The 1σ and 2σ contours of solutions to the B neutrino flux fB and the νe survival probability Pee assuming νe to νa transition. νe transition into a mixed state (νa sin α + νs cos α). The dashed line denote the ±2σ limits of the SSM. The 1σ SSM error bar for fB is indicated on the right. other hand for the general case of νe transition into a of fB from this fit with the 2σ upper limit from the combination of νa and νs our approach effectively as- SSM (vertical lines) gives a lower limit of sin2 α> sumes Pes to be energy independent down to 2.2 MeV. 0.45, i.e., the probability of the active component CC NC A comparison of the current values RSNO with RSNO is is > 45%. Note that there is no upper limit on this shown in Fig. 1. It constitutes a 5.3 σ signal for transi- quantity, i.e., the data is perfectly compatible with νe tion to a state containing an active neutrino component transition into purely active neutrinos. or alternatively a 5.3 σ signal against a pure sterile so- Assuming transition into purely active neutrinos lution. (Pea = 1 − Pee) we show in Fig. 3 the 1σ and 2σ Next we consider the general case where νe goes to contours in the fB–Pee plane from the combinations a mixed state = νa sin α + νs cosα. Then one can write SK + SNOCC and SK + SNOCC + SNONC. The 2 Pea = sin α(1 − Pee). Substituting this in Eqs. (3) inclusion of the NC rate narrows down the ranges of and (5) and eliminating Pee using Eq. (4) one gets the fB and Pee. The error in fB after the inclusion of NC 2 following set of equations for fB and sin α [2] data is about half the size of the corresponding error sin2 α f − RCC = Rel − RCC /r, (6) from SSM as is seen from Fig. 3. B SNO SK SNO 2 − CC = NC − CC sin α fB RSNO RSNO RSNO. (7) We treat sin2 α as a model parameter. And for different 2. Model dependent analysis input values of sin2 α we determine the central value 2 and the 1σ and 2σ ranges of fB by taking a weighted In this section we present the results of our χ average of Eqs. (6) and (7). The corresponding curves analysis of solar neutrino rates and SK spectrum data are presented in Fig. 2. Combining the 2σ lower limit in the framework of two flavour oscillation of νe A. Bandyopadhyay et al. / Physics Letters B 540 (2002) 14Ð19 17
Fig. 4. The νe → νa oscillation solutions to the global solar neutrino data using (a) Ga rate, the SK zenith angle energy spectra and the SK and SNO (CC) rates, both normalised to the SNO (NC) rate and (b) total Ga, Cl, SK, SNO (CC) and SNO (NC) rates along with the SK zenith angle 8 8 energy spectra, keeping the B flux normalisation fB free. In both cases we use the SNO (NC) error as the error in the B flux. (c) is similar to (b), but without using the SNO (NC) rate. to an active flavour. We use the standard techniques SK spectrum data. The SNO CC and NC rates have a described in our earlier papers [10,11] excepting for large anticorrelation. We have taken into account this el CC the fact that instead of the quantities RSK and RSNO correlation between the measured SNO rates in our el NC CC NC global analyses. Further details of this fitting method we now fit the ratios RSK/RSNO and RSNO/RSNO. The 8B flux normalisation gets cancelled from these can be found in [2]. In Table 2 we present the best-fit 2 ratios and the analysis becomes independent of the parameters, χmin and goodness of fit (GOF). The best- large (16–20%) SSM uncertainty associated with this. fit comes in the HIGH(LMA) region as before [11,12]. We include in our global analyses the 1496 day SK However as is seen from Fig. 4(a) the incorporation of zenith angle spectra [15]. Since we use both SK rate the NC data narrows down the allowed regions, and in and SK spectrum data we keep a free normalisation particular the LOW region becomes much smaller. 2 factor for the SK spectrum. This amounts to taking We have also performed an alternative χ fit to the the information on total rates from the SK rates data rates of Table 1 [1,7,13,14] along with the 1496 day and the information of the spectral shape from the SK spectra [15], keeping fB as a free parameter. Even 18 A. Bandyopadhyay et al. / Physics Letters B 540 (2002) 14Ð19
Table 1 The observed solar neutrino rates relative to the SSM predictions (BP2000) are shown along with their compositions for different experiments. For the SK experiment the νe contribution to the rate R is shown in parentheses assuming νe → νa transition. In the combined Ga rate we have included the latest data from SAGE and GNO Experiment R Composition Ga 0.553 ± 0.034 pp(55%), Be(25%), B(10%) Cl 0.337 ± 0.030 B(75%), Be(15%) SK 0.465 ± 0.014 (0.363 ± 0.014) B(100%) SNO(CC) 0.349 ± 0.021 B(100%) SNO(NC) 1.008 ± 0.123 B(100%)
Table 2 2 The χmin, the goodness of fit and the best-fit values of the oscillation parameters obtained for the analysis of the global solar neutrino data 2 2 2 2 Data used Nature of solution m in eV tan θχmin Goodness of fit (%) − Ga + LMA 9.66 × 10 5 0.41 35.95 80.08 − SK/NC + LOW 1.04 × 10 7 0.61 46.73 36.09 − CC/NC + VO 4.48×10 10 0.99 54.25 13.84 − − SKspec SMA 6.66 × 10 6 1.35 × 10 3 67.06 1.41 − Cl + Ga+ LMA 6.07 × 10 5 0.41 40.57 65.99 − SK + CC+ LOW 1.02 × 10 7 0.60 50.62 26.14 − NC + SKspec VO 4.43×10 10 1.156.11 12.39 −6 −3 +fB free SMA 5.05 × 10 1.68 × 10 70.97 0.81
though we allow fB to vary freely the NC data serves Letter we have discussed two useful strategies, of 2 to control fB within a range determined by its error. incorporating the NC data in the global χ analysis of As we see from Table 2 and Fig. 4(b) the results of rates and spectrum data, by which one can avoid the this fit are very similar to the previous case. The best fit large 8B flux uncertainty from the SSM. comes from the HIGH(LMA) region, while no allowed • We fit the ratios of the SK elastic and SNO CC region is obtained for the LOW solution at the 99% rates w.r.t. the NC rate, from which the f cancels CL level. Maximal mixing is seen to be disallowed B out. at the 3σ level. To illustrate the impact of the NC • We fit the rates by keeping f as a free parameter, rate on the oscillation solutions we have repeated the B where the inclusion of the SNO NC rate (= f ) free f fit without this rate. The results are shown in B B serves to control this parameter. Fig 4(c). Evidently the NC data plays a pivotal role in constraining the oscillation solutions, particularly in Both the analyses give very similar results. They the LOW/QVO region, which is allowed only at the 3σ clearly favour the HIGH(LMA) solution, while a level. It puts an upper bound on the m2 in the LMA limited region of the LOW solution is also acceptable region and rules out maximal mixing. at the 3σ level. The maximal mixing solution is disfavoured at the 3σ level. As more data accumulate one expects a substantial reduction in the error bar of 3. Summary and conclusions the SNO NC rate, resulting in further tightening of the allowed regions of neutrino mass and mixing. The first SNO NC data constitutes a 5.3 σ signal for transition into a state containing an active neutrino component. The inclusion of this data puts much Note added tighter constraints on fB and Pee from a model independent analysis involving active neutrinos as The paper [16] appeared on the net after completion compared to the SNO CC/SK combination. In this of our work. In the region of overlap our results agree A. Bandyopadhyay et al. / Physics Letters B 540 (2002) 14Ð19 19 with theirs as well as with the updated version of [17]. S. Goswami, D. Majumdar, A. Raychaudhuri, hep- It may be added here that the SNO CC and NC rates ph/9909453; given in Table 1 are obtained assuming undistorted S. Choubey, S. Goswami, N. Gupta, D.P. Roy, Phys. Rev. D 64 (2001) 053002. energy spectra above 5 MeV, which for transitions [11] A. Bandyopadhyay, S. Choubey, S. Goswami, K. Kar, Phys. to active neutrinos has good empirical justification as Lett. B 519 (2001) 83; mentioned above. We thank Prof. Mark Chen of SNO S. Choubey, S. Goswami, K. Kar, H.M. Antia, S.M. Chitre, Collaboration for communication on this point. Phys. Rev. D 64 (2001) 113001; S. Choubey, S. Goswami, D.P. Roy, Phys. Rev. D 65 (2002) 073001; A. Bandyopadhyay, S. Choubey, S. Goswami, K. Kar, Phys. References Rev. D 65 (2002) 073031. [12] G.L. Fogli, E. Lisi, D. Montanino, A. Palazzo, Phys. Rev. D 64 [1] SNO Collaboration, Q.R. Ahmad et al., nucl-ex/0204008, to (2001) 093007; appear in Phys. Rev. Lett. J.N. Bahcall, M.C. Gonzalez-Garcia, C. Pena-Garay, [2] A. Bandyopadhyay, S. Choubey, S. Goswami, D.P. Roy, hep- JHEP 0108 (2001) 014; ph/0203169, to appear in Mod. Phys. Lett. A. P.I. Krastev, A.Yu. Smirnov, hep-ph/0108177; [3] See, e.g., V. Barger, D. Marfatia, K. Whisnant, Phys. Rev. M.V. Garzelli, C. Giunti, JHEP 0112 (2001) 017. Lett. 88 (2002) 011302. [13] J.N. Abduratshitov et al., SAGE Collaboration, astro- [4] J.N. Bahcall, M.H. Pinsonneault, S. Basu, Astrophys. J. 555 ph/0204245; (2001) 990. W. Hampel et al., GALLEX Collaboration, Phys. Lett. B 447 [5] F.L. Villante, G. Fiorentini, E. Lisi, Phys. Rev. D 59 (1999) (1999) 127; 013006. M. Altman et al., GNO Collaboration, Phys. Lett. B 490 (2000) [6] G.L. Fogli, E. Lisi, D. Montanino, A. Palazzo, Phys. Rev. D 64 16; (2001) 093007. E. Belloti, Talk at Gran Sasso National Laboratories, May 17, [7] S. Fukuda et al., Super-Kamiokande Collaboration, Phys. Rev. 2002; Lett. 86 (2001) 5651. T. Kirsten, Talk at Neutrino 2002, Munich. [8] SNO Collaboration, Q.R. Ahmad et al., Phys. Rev. Lett. 87 [14] B. Cleveland et al., Astrophys. J. 496 (1998) 505. (2001) 071301. [15] M.B. Smy, Talk at noon2001, hep-ex/0202020. [9] SNO Collaboration, Q.R. Ahmad et al., nucl-ex/0204009, to [16] V. Barger, D. Marfatia, K. Whisnant, B.P. Wood, hep- appear in Phys. Rev. Lett. ph/0204253. [10] S. Goswami, D. Majumdar, A. Raychaudhuri, Phys. Rev. D 63 [17] P. Creminelli, G. Signorelli, A. Strumia, hep-ph/0102234. (2001) 013003; Physics Letters B 540 (2002) 20–24 www.elsevier.com/locate/npe
Non-standard contributions to ν–e elastic scattering from solar neutrinos observations and LSND measurement
A. Ianni
INFN-Laboratori Nazionali del Gran Sasso, S.S. 17bis Km 18+910, I-67010 Assergi (Aquila), Italy Received 16 May 2002; received in revised form 6 June 2002; accepted 10 June 2002 Editor: L. Rolandi
Abstract We calculate bounds on neutrino magnetic moment and charge radius squared from solar neutrinos observations after Super- Kamiokande and SNO results. We use LSND νe–e elastic scattering measurement to constrain further the charge radius. −10 −10 2 −32 2 With few assumptions, we derive µν 1 × 10 µB (µν 1.2 × 10 µB ), when r =0, and −1.2 × 10 r −32 2 −32 2 −32 2 2.7 × 10 cm (−0.5 × 10 r 3.5 × 10 cm ), when µν = 0, at 90% C.L. for Dirac (Majorana) neutrinos. We −10 also show that µν cannot be larger than 1.2 × 10 µB for any allowed value of the charge radius. Moreover, we show that from our fit a fraction of about 1.3% of νxR (x = µ,τ) could be present in the solar neutrinos flux on Earth. 2002 Published by Elsevier Science B.V.
1. Introduction Dirac neutrino species of the order of [1] −19 mν µν ∼ 3.2 × 10 µB . (1) In the Standard Model neutrinos are massless and 1eV do not have electromagnetic dipole and transition mo- Among the non-standard processes allowed by a non- ments [1], i.e., they do not couple to the electromag- zero µν ,theν–e elastic scattering is of primary netic field. However, there are many possible exten- importance because can be studied extensively in sions of the Standard Model which allow neutrinos different energy ranges. In this exotic scenario the to interact directly to the electromagnetic field [2]. differential cross-section for ν–e elastic scattering is This possibility is of great interest because such a written as [4] new coupling would allow a variety of non-standard 2 dσ(Eν,T) GF me 2 processes [1]. The strong evidence of neutrino flavor = (CV + X + CA) dT 2π transformation reported by the SNO Collaboration [3] shows that neutrino have a non-zero mass and as a con- 2 + − 2 − T sequence they should have a non-zero magnetic mo- (CV CA) 1 Eν ment. In particular, a minimal extension of the Stan- 2 2 meT dard Model predicts a magnetic moment for massive + C − (CV + X) A E2 ν 2 + 2 πα 1 − 1 E-mail address: [email protected] (A. Ianni). µν 2 . (2) me T Eν 0370-2693/02/$ – see front matter 2002 Published by Elsevier Science B.V. PII:S0370-2693(02)02126-3 A. Ianni / Physics Letters B 540 (2002) 20–24 21
2 1 In Eq. (2) CV = 2sin θW + 1/2, CA = 1/2forνe, take into account an energy-dependent deficit. With 2 CV = 2sin θW − 1/2, CA =−1/2forνx (x = µ,τ). this in mind we should consider the possibility of Moreover, dealing with anitneutrinos, we must make both flavor oscillations (MSW) and, as we are deal- the substitution CA →−CA. In Eq. (2) the magnetic ing with a non-zero µν , of resonant spin-flavor pre- moment, µν ,√ is given in units of µB (Bohr magne- cession (RSFP) [10,11,14]. For this latter we distin- 2 παr2 → →¯ ton) and X = ,wherer2 is the mean-square guish νeL νeR,νxR and νeL νxR transitions for 3GF 3 charge radius. With this definition we can see that a Dirac and Majorana neutrinos, respectively, where as − = charge radius of ∼ 6.5 × 10 16 cm2 gives X = 1and above x µ,τ. this corresponds to a significant change in the stan- Both processes (MSW and RSFP) could be used dard cross-section. We notice that Eq. (2) is given as a to fit solar neutrinos data [15–17]. Hence, following function of the neutrino energy, Eν , the recoil electron a reasoning similar to the one presented in [18], kinetic energy, T , and contains two phenomenologi- we call α the survival probability at Earth of νeL’s 2 β cal parameters, namely µν and r . For this reason produced at the core of the Sun and the fraction of a correct approach to study non-standard electromag- νxL’s which have not changed chirality because of an netic interactions should take into account both para- electromagnetic interaction with the Sun’s magnetic 4 meters as they are strongly correlated [5]. Finally, we field. With this model in mind, on Earth, we will also point out that the interaction cross-section propor- measure νeL’s with a reduction factor (with respect to tional to µ2 changes the chirality of the incoming neu- the SSM prediction) α and, as an example, a fraction ν ¯ − − trino. This could play an important role if we study of νxR equal to (1 α)(1 β). electromagnetic processes involving solar neutrinos In the most general scenario both α and β are en- which have been travelling through the magnetic field ergy dependent. In the following, however, as a first inside the Sun before reaching a detector on Earth. approximation we consider α and β to be constant. As 2 a matter of fact in the LMA5 scenario for solar neutri- Existing laboratory limits on µν and r come from ν–e elastic scattering [6–8]. These are of the nos, the survival probability is with good approxima- −10 −9 E order of few ×10 µB for νe and few ×10 µB for tion constant [16] for ν 5 MeV. So, in the energy range of Super-Kamiokande and SNO this is a good νµ. In particular, a study using 825-days data from 2 −10 W(T) 6 Super-Kamiokande gives µν 1.5 × 10 µB at assumption. In this case, if we call, , the de- 90% C.L. [7] On the other hand astrophysical limits tector energy resolution function, the effective cross- −13 −11 ν e are of the order of 10 –10 µB [9–11]. section for solar neutrinos from – elastic scattering, Best limits on r2 are derived using accelerator + and beam dump experiments as the sensitivity on this Eν /(1me/2Eν ) dσ(Eν,T) parameter increases with the neutrino energy. For νe σ(Eν ) = dT W(T) , (3) − × −32 2 × −32 2 dT it is found [12] 3 10 r 4 10 cm , 0 −16 while for νµ [13] the limit is | r | 0.7 × 10 cm.
3 Here we are not considering the possibility to have electron 2. Searching for neutrino non-standard antineutrinos in the final state due to a conversion of the kind νeL →¯νxR →¯νeL as these could be easily identified through the electromagnetic interactions with solar neutrinos inverse β-decay reaction. 4 observations For Dirac neutrinos the final state, νxR, is sterile with respect to weak interactions. 5 In this Letter we analyse data from solar neutri- The Large Mixing Angle scenario is strongly favored by a 2 global analysis of solar neutrinos data [16]. nos observations to search for µν and r . A real- 6 Here, using a gaussian energy resolution function, istic analysis of solar neutrinos observations should Tup − T T − T W(T)= 0.5 Erf √ − Erf low√ , 2σ 2σ 1 2 T T In this Letter we use sin θW = 0.2315. 2 This can be considered a laboratory limit as the solar neutrino where Erf is the error function, Tup and Tlow are the experimental flux is now experimentally established. recoil electron kinetic energy upper and lower limit, respectively. 22 A. Ianni / Physics Letters B 540 (2002) 20–24 can be split, using Eq. (2), into terms proportional to Table 1 2 2 µν and r . Doing this we can write σ from Eq. (3) Parameters for the SNO χSNO function. See text for details as ρx (ρx¯ )ηem η1e η1x (η1x¯ )η2 − . ( . ) . × −4 σ = ασweak + β(1 − α)σweakµ2σ 0.15(0.12) 0.011 0.033 0 011 0 009 2 7 10 e x ν em 2 + r ασ1e + β(1 − α)σ1x + 2 SK + − SK 2 2 r αη1e β(1 α)η1x + r α + β(1 − α) σ2, (4) + r2 2 α + β(1 − α) η2SK − data . (8) and Dirac neutrinos7 and i for Dirac neutrinos, and σ = ασweak + β(1 − α)σweak e x weak 2 SK = + − SK + (1 − α)(1 − β)σ + µ σ ai α β(1 α)ρx x¯ ν em 2 + − − SK + 2 SK + r ασ1e + β(1 − α)σ1x (1 α)(1 β)ρx¯ µν ηem + − − + 2 2 + r2 αη1SK + β(1 − α)η1SK (1 α)(1 β)σ1x¯ r σ2, (5) e x + − − SK weak (1 α)(1 β)η1x¯ for Majorana neutrinos. In Eqs. (4) and (5) σe weak weak 2 2 SK and σx (σx¯ ) are the Standard Model cross- + r η2 − datai (9) sections for ν–e elastic scattering for νe and νµ,τ for Majorana neutrinos. (ν¯µ,τ ), respectively, σem is the cross-section due to the In Eqs. (8) and (9) spin-flip term in Eq. (2). The cross-sections σ1e, σ1x (σ1 ¯ )andσ2 are from the contributions proportional x ρSK ρSK = σ weak/σ weak(σ weak/σ weak), to r2 and r22 for ν (ν¯), respectively. x x¯ x e x¯ e SK = weak Moreover, to derive Eqs. (4) and (5) we have ηem σem/σe , assumed that µν is an effective magnetic moment and SK weak η1 ¯ = σ1e,x,x¯ /σ , is the same for different neutrino flavors, the same e,x,x e SK = weak being true for the charge radius. η2 σ2/σe In this Letter, in particular, we take into account and data is the Super-Kamiokande relative measure- the Super-Kamiokande measurement of the recoil i ment in bin i. In evaluating Eq. (7) we have used electron spectrum [19] for 1258-days data taking and 8Bandhep spectra according to [20], a gaussian en- the SNO measurement of ν–e elastic scattering [3]. √ ergy resolution with σ = 0.47 T/1 MeV and data For the Super-Kamiokande data we calculate the error T from [19]. Cross-sections are folded over the neutri- covariance matrix according to data reported in [19] as nos energy spectra. = 2 + For SNO the χ2 function is written Sij σi δij si sj , (6) where σ 2 is the sum in quadrature of the statistical aSNO − rSNO 2 i χ2 = ES , (10) and uncorrelated errors and si are the correlated ones. SNO SNO σr With Eq. (6) we can define a χ2 function for the Super- SNO SK SNO = SNO SSM Kamionade data as where a is similar to a , rES φES /φ SNO SNO SNO = ± 2 2 2 SK −1 SK and σr is the error on rES .HereφES (2.39 χ α, β, µ , r = a S a , (7) − − − SK ν i ij j 0.27) × 10 6 cm 2 s 1 is the measured elastic scat- ij tering neutrinos flux [3], φSSM = (5.05 ± 1.01) × where 10−6 cm−2 s−1 is the SSM 8B flux [20]. In evaluat- ing Eq. (10) we have used a gaussian energy√ resolution aSK = α + β(1 − α)ρSK + µ2ηSK i x ν em function with σT =−0.0684 + 0.331 T = 0.0425T as in [3] with a lower threshold at 5 MeV for the re- 7 Eq. 4 holds for both νeL → νeR and νeL → νxR with the coil electron kinetic energy. Table 1 reports the values assumptions made in the Letter (see below). of terms used in Eq. (10). A. Ianni / Physics Letters B 540 (2002) 20–24 23
3. The LSND contribution
To improve the sensitivity on r2 we used the measurement of νe–e cross-section performed by the LSND Collaboration [12]. Using for νe a spectrum shape of the form ∝ 2 2 2 − −1 −1 dΓ/dE mµc Eν (mµc 2Eν) MeV s , 2 with mµc muon rest energy, and detector information from [12], we have calculated 167 expected events in the Standard Model scenario against the 191 ± 22 measured [12]. Hence in 1 d.o.f. analysis we derive × −9 µνe 1.1 10 µB at 90% C.L. in agreement with what reported in [12]. However, in a proper 2 d.o.f. 2 2 analysis we will write a χ function similar to χSNO with α = β = 1. Fig. 1. Allowed regions (90–99% C.L.) from Super-Kamiokande, In order to combine results from Super-Kamiok- SNO and LSND data on ν–e elastic scattering electromagnetic phenomenological parameters for Dirac neutrinos. See text for ande, SNO and LSND we add Eqs. (7), (10) and the details. χ2 function for LSND. The SNO measurement of neutral current [3] al- SNO SNO = ± lows us to measure φCC /φNC 0.35 0.02, where SNO SNO φCC and φNC are the measured charge and neu- tral current fluxes, respectively. Therefore, to mini- mize the combined χ2 function we search for a mini- 2 2 mum in µν and r with 0.29 α 0.41 (3σ range) and β in [0.5, 1] for Majorana neutrinos, for in this SNO SNO = case φCC /φNC α. On the other hand, for Dirac SNO SNO = + − neutrinos φCC /φNC α/(α β(1 α)). Hence, for any β in [0.5, 1], we search for a minimum with (β + 0.155)/3.85 α (β + 0.28)/3.125, which cor- responds to a slide in the α–β plane where we can con- strain these parameters on the basis of the 3σ range SNO SNO measured value for the φCC /φNC ratio. Result of minimization for both kind of neutrinos 2 ≈ { = gives χmin 20.24 and best-fit parameters are α = = × −4 2= } Fig. 2. Allowed regions (90–99% C.L.) from Super-Kamiokande, 0.35,β 0.98,µν 3.0 10 , r 1.50 for SNO and LSND data on ν–e elastic scattering electromagnetic Dirac neutrinos and {α = 0.342,β = 0.995,µν = phenomenological parameters for Majorana neutrinos. See text for − 3.4 × 10 5, r2=1.47} for Majorana ones.8 Allowed details. 2 2 regions, in the (r –µν ) plane with α and β equal 2 = to the best fit values, at 90% and 99% C.L. ( χ 10−32 cm2 (−0.5 × 10−32 r2 3.5 × 10−32 cm2) 2 − 2 = χ χmin 4.61, 9.21) are reported in Figs. 1 and 2, when µ = 0 for Dirac (Majorana) neutrinos. In order ν respectively. From these figures we can see that µν to have a better idea of how the fit result combine × −10 × −10 1.0 10 µB (µν 1.2 10 µB ) at 90% C.L. with the experimental data in Fig. 3 we show the 2= − × −32 2 × when r 0, and 1.2 10 r 2.7 measured Super-Kamiokande relative spectrum for solar neutrinos [19] with the predicted spectrum with 8 −10 2 Here as above µν is given in units of 10 µB and r in the best fit parameters reported above. In particular, − units of 10 32 cm2. the solid line is for Dirac neutrinos. 24 A. Ianni / Physics Letters B 540 (2002) 20–24
including for this latter a possible calibration test with a neutrino source [25], could improve further the sensitivity to study lower values of µν ,forthese experiments will probe lower energy regions around 1 MeV. Moreover, a possible combined analysis using CHARM II data [8] could help improving limits on r2.
Acknowledgements
The author thanks J.F. Beacom, D. Montanino and F. Vissani for discussions.
Fig. 3. Measured solar neutrinos spectrum from Su- References per-Kamiokande [19] relative to SSM predictions [20]. Results of fit in case of non-standard electromagnetic contribution to ν–e [1] G. Raffelt, Stars as Laboratories for Fundamental Physics, elastic scattering is shown for Dirac neutrinos (solid line) and Ma- Univ. of Chicago Press, Chicago, IL, 1996. jorana neutrinos (dashed line) using best-fit parameters. See text for [2] R. Mohapatra, P.B. Pal, Massive Neutrino in Physics and details. Astrophysics, World Scientific, Singapore, 1991. [3] SNO Collaboration, Q.R. Ahmad et al., nucl-ex/0204008, see also: http://www.sno.phy.queensu.ca. 4. Conclusions [4] P. Vogel, J. Engel, Phys. Rev. D 39 (1989) 3378. [5] P. Salati, Astrop. Phys. 2 (1994) 269–290. In this Letter we have found bounds for the neu- [6] A. Derbin et al., JEPT Lett. 55 (1993) 206. [7] J. Beacom, P. Vogel, Phys. Rev. Lett. 83 (1999) 5222. trino effective magnetic moment and charge radius ex- [8] CHARM II Collaboration, P. Vilain et al., Phys. Lett. B 345 ploiting a correlated analysis with ν–e elastic scatter- (1995) 115. ing measurements from Super-Kamiokande, SNO and [9] G. Raffelt, Astrophys. J. 365 (1990) 559. LSND. We have shown with few assumptions that the [10] A. Barbieri, R. Mohapatra, Phys. Rev. Lett. 61 (1988) 27; contribution from spin-flipped neutrinos to measured A. Goyal et al., Phys. Lett. B 49 (1995) 312. [11] M.B. Voloshin, Phys. Lett. B 209 (1988) 360. rates in Super-Kamiokande and SNO should be only a [12] LSND Collaboration, hep-ex/0101039. minor component ( 1.3% ) of the signal. This limit is [13] CHARM II Collaboration, P. Vilain et al., Phys. Lett. B 335 calculated without any assumption on the Sun’s mag- (1994) 246. netic field and could be used to develop models of [14] E.Kh. Akhmedov, hep-ph/9705451. analysis where the shape of the field is important to [15] G. Fogli, E. Lisi, D. Montanino, A. Palazzo, Phys. Rev. D 64 (2001) 093007. determine the flavor transition probability. In this kind [16] J.N. Bahcall et al., hep-ph/0204314. of analysis [17,21] one makes the assumption that the [17] J. Pulido, hep-ph/0106201; − neutrino magnetic moment is large enough (10 12– J. Pulido, hep-ph/0112104. −13 10 µB ) to obtain a significant conversion between [18] J. Pulido, A. Mourao, Phys. Rev. D 57 (1998) 1794. neutrino flavors inside the Sun. In this Letter, on the [19] Super-Kamiokande Collaboration, hep-ex/0103032. [20] J.N. Bahcall et al., astro-ph/0010346. contrary, we have used observations to derive a lim- [21] A. Friedland, A. Gruzinov, hep-ph/0202095. its on the magnetic moment and as side result we have [22] A. Pieke, Nucl. Phys. B (Proc. Suppl.) 91 (2001) 99. calculated the fraction of solar neutrinos which may [23] C. Broggini et al., Nucl. Phys. B (Proc. Suppl.) 91 (2001) 105. undergo a spin-flip conversion. [24] BOREXINO Collaboration, Astrop. Phys. 16 (2002) 205. We conclude pointing out that future results from [25] A. Ianni, D. Montanino, G. Scioscia, Eur. Phys. J. C 8 (1999) 609. KamLand [22], MUNU [23] and BOREXINO [24], Physics Letters B 540 (2002) 25–32 www.elsevier.com/locate/npe
+ Measurements of relative branching ratios of c decays into states containing
FOCUS Collaboration J.M. Link a,M.Reyesa,P.M.Yagera,J.C.Anjosb,I.Bediagab,C.Göbelb,J.Magninb, A. Massafferri b,J.M.deMirandab,I.M.Pepeb,A.C.dosReisb,S.Carrilloc, E. Casimiro c,E.Cuautlec, A. Sánchez-Hernández c,C.Uribec,F.Vázquezc, L. Agostino d, L. Cinquini d,J.P.Cumalatd, B. O’Reilly d,J.E.Ramirezd, I. Segoni d, J.N. Butler e, H.W.K. Cheung e, G. Chiodini e, I. Gaines e, P.H. Garbincius e, L.A. Garren e,E.Gottschalke,P.H.Kaspere,A.E.Kreymere,R.Kutschkee, L. Benussi f,S.Biancof,F.L.Fabbrif,A.Zallof,C.Cawlfieldg,D.Y.Kimg,K.S.Parkg, A. Rahimi g,J.Wissg,R.Gardnerh, A. Kryemadhi h,K.H.Changi, Y.S. Chung i, J.S. Kang i,B.R.Koi,J.W.Kwaki,K.B.Leei,K.Choj,H.Parkj, G. Alimonti k, S. Barberis k,A.Ceruttik,M.Boschinik,P.D’Angelok,M.DiCoratok,P.Dinik, L. Edera k,S.Erbak, M. Giammarchi k,P.Inzanik,F.Leverarok, S. Malvezzi k, D. Menasce k, M. Mezzadri k, L. Moroni k,D.Pedrinik, C. Pontoglio k,F.Prelzk, M. Rovere k,S.Salak, T.F. Davenport III l,V.Arenam,G.Bocam, G. Bonomi m, G. Gianini m, G. Liguori m,M.M.Merlom, D. Pantea m,S.P.Rattim, C. Riccardi m, P. Vitulo m, H. Hernandez n,A.M.Lopezn,H.Mendezn,L.Mendezn,E.Montieln, D. Olaya n,A.Parisn, J. Quinones n,C.Riveran,W.Xiongn,Y.Zhangn,J.R.Wilsono, T. Handler p, R. Mitchell p,D.Enghq,M.Hosackq, W.E. Johns q,M.Nehringq, P.D. Sheldon q,K.Stensonq, E.W. Vaandering q,M.Websterq,M.Sheaffr
a University of California, Davis, CA 95616, USA b Centro Brasileiro de Pesquisas Físicas, Rio de Janeiro, RJ, Brazil c CINVESTAV, 07000 Mexico City, DF, Mexico d University of Colorado, Boulder, CO 80309, USA e Fermi National Accelerator Laboratory, Batavia, IL 60510, USA f Laboratori Nazionali di Frascati dell’INFN, Frascati I-00044, Italy g University of Illinois, Urbana-Champaign, IL 61801, USA h Indiana University, Bloomington, IN 47405, USA i Korea University, Seoul 136-701, South Korea j Kyungpook National University, Taegu 702-701, South Korea k INFN and University of Milano, Milano, Italy l University of North Carolina, Asheville, NC 28804, USA m Dipartimento di Fisica Nucleare e Teorica and INFN, Pavia, Italy n University of Puerto Rico, Mayaguez, PR 00681, USA
0370-2693/02/$ – see front matter 2002 Elsevier Science B.V. All rights reserved. PII:S0370-2693(02)02103-2 26 FOCUS Collaboration / Physics Letters B 540 (2002) 25–32
o University of South Carolina, Columbia, SC 29208, USA p University of Tennessee, Knoxville, TN 37996, USA q Vanderbilt University, Nashville, TN 37235, USA r University of Wisconsin, Madison, WI 53706, USA Received 8 June 2002; accepted 10 June 2002 Editor: L. Montanet
Abstract + + ∗0 + + + − We have studied the Cabibbo suppressed decay c → K (892) and the Cabibbo favored decays c → K K , + + + ∗0 + − + + + + − c → φ and c → ( K )K and measured their branching ratios relative to c → π π to be (7.8 ± 1.8 ± 1.3)%, (7.1 ± 1.1 ± 1.1)%, (8.7 ± 1.6 ± 0.6)%and(2.2 ± 0.6 ± 0.6)%, respectively. The first error is statistical and the + − + + + + ∗0 second is systematic. We also report two 90% confidence level limits ( c → K π )/ ( c → K (892)) < 35% + + + − + + + − and ( c → K K )NR/ ( c → π π )<2.8%. 2002 Elsevier Science B.V. All rights reserved.
1. Introduction complete the charged particle tracking and momen- tum measurement system. Three multi-cell, thresh- Past experiments have reported results on non- old Cerenkovˇ counters discriminate between different leptonic branching fractions of the lowest lying particle hypotheses, namely, electrons, pions, kaons + charmed baryon c [1,2]. In this Letter we report on and protons. The FOCUS apparatus also contains one + several c decay channels containing a baryon in hadronic and two electromagnetic calorimeters as well the final state. These measurements may be useful in as two muon detectors. testing theoretical predictions of the contributions to Events are selected using a candidate driven ver- inclusive decay amplitudes. For instance, as pointed texing algorithm where the vector components of the out by Guberina and Stefancic [3], direct measure- reconstructed decay particles define the charm flight + ments of c singly Cabibbo suppressed decay rates direction. This is used as a seed track to find the pro- + can improve our theoretical understanding of the c duction vertex [6]. Using this algorithm we determine lifetime, which can then be compared to recent high the confidence level of the decay and production ver- statistics measurements [4]. tices, and the significance of their separation. For each of the decay modes analyzed, we require the primary vertex to have a confidence level greater than 1% and 2. Reconstruction to contain at least two tracks other than the charm seed track. − → − + → + + → 0 This analysis uses data collected by the FOCUS The nπ , nπ and pπ 1 experiment at Fermilab during the 1996–1997 fixed- decays are reconstructed using a kink algorithm [7] target run and is based on a topological sample of where the properties of the neutral particle in the decay events with a charged Sigma hyperon plus two other are not detected, but rather inferred, with a two-fold + ambiguity in the momentum solution for some decays. charged particles emerging from the c decay vertex. FOCUS is a photo-production experiment equipped Systematic effects due to this ambiguity are reduced + → + + − with very precise vertexing and particle identification by normalizing to the decay mode c π π , detectors. The vertexing system is composed of a sili- where the same effect exists. To aid in fitting the con microstrip detector (TS) [5] interleaved with seg- mass distribution, for channels containing a ,we ments of the BeO target and a second system of twelve microstrip planes (SSD) downstream of the target. Be- 1 Throughout this Letter the charge conjugate state is implied − yond the SSD, five stations of multi-wire proportional unless explicitly stated. Note that the is not the charge conjugate + chambers plus two large aperture dipole magnets partner of the . FOCUS Collaboration / Physics Letters B 540 (2002) 25–32 27 implement a double Gaussian to determine the yield of signal events. By double Gaussian we mean two Gaussian shapes with separate amplitudes, means and widths. In this Letter we will discuss the +π+π−, +K+π−, −K+π+ and +K+K− final states. In + → + + − the c π π mode we let all the parameters float while in the other lower statistics modes we fix some of the parameters to their Monte Carlo values.
+ → + + − 3. c π π normalization mode + → + + − The c π π mode is our highest statis- tics decay containing a + particle. The events are selected requiring a detachment between primary and secondary vertex divided by its error (l/σl ) greater + + − then 5.5. A minimum + momentum cut of 50 GeV/c Fig. 1. π π invariant mass distribution fit with a double c Gaussian for the signal and a linear background. is imposed, as is a minimum secondary vertex confi- dence level of 10%. We also apply a cut on the lifetime + → + ∗0 + → − + + resolution, σt < 120 fs for the run period where we 4. c K (892) and c K π had a silicon detector (TS) in the target region (about decay modes 2/3 of the events) and σt < 150 fs otherwise (NoTS). + Further, we reject events which have a lifetime greater The c is reconstructed in the decay channel + + + − than six times the c lifetime. K π . We fit the invariant mass distribution We identify charged tracks using information from with and without a mass cut (832