Stephen Wolframstephenwolfram.com
Publications by Stephen Wolfram * Articles * Cosmology * Abundances of New Stable Particles Produced in the Early Universe (1979)
Abundances of New Stable Particles Produced in the Early Universe (1979)


Main Text

The standard model of the early universe has recently been used to place constraints on the masses and lifetimes of possible nearly-stable heavy neutrino-like particles predicted by various gauge models of weak interactions [1]. Several models of this kind imply the existence of absolutely-stable charged and/or strongly-interacting particles more massive than the proton (e.g. [2]). In this note, I show that rather large numbers of such particles would have been produced in the early universe, so that experimental limits on their terrestrial abundances may place stringent bounds on their masses.

Any new stable charged particles with masses below about should already have been seen in interactions. The next generation of accelerators (PETRA, PEP) could extend this limit to masses up to . Attempts to produce pairs of new stable hadrons in 400 GeV proton interactions have probed up to masses [2,3], but the production cross-sections for heavy hadrons near threshold are not known with sufficient accuracy for definite conclusions to be drawn [4].

The number density () of any species of stable particles spread uniformly throughout a homogenous universe should obey the rate equation [1,5]

where is the expansion scale factor for the universe and is the product of the low-energy annihilation cross-section and relative velocity for the particles, averaged over their energy distribution at time . is their number density in thermal equilibrium. The first term in eq. (1) accounts for the dilution in due to the expansion of the universe, while the second term arises from the annihilation and production of particles in interactions. Let

where is the equilibrium temperature, and in the upper (lower) sign is for fermions (bosons). Then, ignoring the curvature of the universe, which has no effect at the times we consider, eq. (1) becomes (1)

If the temperature of the universe was arbitrarily high at early times (2) (and the cross-sections for particle interactions do not decrease too rapidly at very high energies), then all particle species should then have been in thermal equilibrium, so that the boundary condition in eq. (1) was or . The solutions of eq. (3) for various values of subject to this boundary condition are shown in fig. 1. As the universe cooled, the equilibrium number density of particle species fell dramatically around . The more strongly-interacting (higher ) the particles were, the longer they will have remained in thermal equilibrium, and thus the lower their final number density will have been.



[ Figure 1 ] Solutions to the differential eq. (3) for various values of the parameter (in units of ). The curves give the number densities (divided by ) of particle species with a variety of low-energy annihilation cross-sections as a function of the average temperature of the universe . The equilibrium number density, , is also given.

The parameter appearing in eq. (3) is the effective number particle species in thermal equilibrium at temperature . It determines the energy density and hence the expansion rate of the universe. Ultrarelativistic fermion (boson) spin states contribute 7/16(1/2) to . (The observed spectrum of particles suggests that for ; for ; for (according to QCD quarks and gluons should contribute to as if they were free for GeV); and for .)

The present number density of a particle species is given approximately by , where is the temperature which the microwave background radiation would now have if it had frozen out of thermal equilibrium at the same time and temperature as the particle species under consideration. The difference between and the present temperature of the actual microwave background radiation arises from the heating of the universe by the annihilation of other species. Specific entropy conservation gives

Eq. (3) may be solved approximately by assuming for , and neglecting compared to for . This gives

which is the correct solution to eq. (3) within about a factor of 20 for the cases considered below.

To obtain estimates of for particular types of particles, one must estimate . Charged stable heavy leptons with should annihilate primarily into two photons, and through a virtual photon to hadrons and lighter leptons, giving

This cross-section, together with the form for discussed above, may now be used to solve eq. (3) and to obtain an estimate for the present abundances of any charged stable heavy leptons. (The exact results are well-approximated by eq. (4).) One finds that for , corresponding to an abundance of about one new stable charged heavy lepton in nucleons. For , the estimated number density rises roughly linearly with , except for slight decreases due to increases in the annihilation cross-section associated with the opening of new channels. The abundances of any produced in the early universe should therefore be rather large, and hence easily amenable to experimental investigation.

To estimate the present abundances of any stable heavy hadrons (H) (containing heavy quarks Q), one must assume a form for the low-energy annihilation cross-section. An upper bound on is probably provided by the low-energy limit of for protons [6] . If , then the universe at the freezing temperature for H should have consisted of almost free quarks and gluons, so that a better estimate of annihilation may be given by the rate for electromagnetic annihilation (5) and for (obtained from the first term of eq. (5) by replacing by ). The first estimate for leads to for , decreasing (roughly as ) to for . The second estimate for suggests for , increasing roughly as , and perhaps reaching for . Since it seems most unlikely that the annihilation cross-section is smaller than its value according to the first estimate, any stable heavy hadrons (with masses below about ) should exist in concentrations above one in about nucleons.

These estimates for heavy hadron abundances may be applied to protons. They give a result too small. The discrepancy is due to the assumption of homogeneity made in eq. (1); in fact, there must either be a net excess of baryons over antibaryons in the universe, or protons and antiprotons must have become spatially separated (presumably at ) thereby preventing their annihilation [7]. Similar phenomena may have occured for other stable particles. (An indication that they were not important comes from the result that the present chemical potential for all species of neutrinos is below (3) , while for [8].) Inhomogeneity can serve only to increase , so that our estimates should be considered in fact as lower bounds on .

The observed average mass density in the present universe is around . The requirement that yet unobserved new stable particles produced in the early universe should not contribute a larger mass density than is observed yields (from eq.(4)) , which is irrelevant for all species of particles except those undergoing only weak interactions [1].

After their production in the early universe, stable heavy particles will presumably have followed the gravitational clumping of ordinary matter. Their number densities should not, however, usually have become sufficiently high for much annihilation to occur. Any produced should have been combined into tightly-bound systems, while should occur in or, in the absence of many , composites. The fact that the lightest strange and charmed baryons do not undergo strong decay indicates that the lightest baryon carrying a new absolutely-conserved quantum number should not be able to decay into a meson carrying the same quantum number and should therefore be stable. These new stable baryons and mesons should be bound into ordinary nuclei. Any L and H produced in the early universe should therefore occur in terrestrial material.

Another source of heavy stable particles is pair production by the interaction of cosmic ray particles with the earth's atmosphere. Assuming that all will eventually get into water, this gives [4] (4). The cosmic-ray-induced heavy hadron abundance should be about. These abundances are insignificant compared to those expected from heavy particle production in the early universe.

There have been a number of searches for heavy integer-charged stable particles, mostly in sea water. The best published experiment [9] found no such particles in nucleons, for almost all masses between 6 and . When combined with the abundances expected form the early universe, this result suggests that no stable integer-charged particles exist with . The most sensitive search yet made is presently being performed [10] using a mass spectrometer to scan the equivalent of kg of sea water. This experiment should detect concentrations down to one particle in nucleons, for . Modern nuclear physics accelerator techniques, if applied to the same sample, should allow the sensitivity of new particles per nucleon to be reached [11]. Even if no heavy stable particles were produced in the early universe, a null result in this experiment would show that their abundance was in many cases below that expected just from their production in cosmic ray interactions. The conclusions that no such particles exist (with masses less than several hundred ) would then surely be inescapable, placing an important constraint on present and future models in particle physics.

I am grateful to N. Isgur and H.J. Rose for discussions.

previous  l  next