Issue 
A&A
Volume 551, March 2013



Article Number  A68  
Number of page(s)  13  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201219985  
Published online  22 February 2013 
Mass limits for dark clusters of degenerate fermions
^{1}
Departamento de Física TeóricaUniversidad de Zaragoza,
50009
Zaragoza,
Spain
email:
membrado@unizar.es
^{2}
BIFI, Instituto de Biofísica y Física de Sistemas Complejos,
Universidad de Zaragoza, 50009
Zaragoza,
Spain
email:
amalio@unizar.es
Received:
11
July
2012
Accepted:
29
December
2012
We calculate the range of possible masses for dark spheres of bound fully degenerate fermions as a function of the fermion mass. The cosmological constant is included in our calculations. We deduce that the minimum fermion mass that is able to give rise to degenerate fermion clusters is ~0.02 g^{−1/4} eV, where g is the spin degeneracy parameter. We show that degenerate fermions of mc^{2} ≈ (15−30) g^{−1/4} eV can build bound degenerate dark haloes that could reproduce the values of the rotation velocities of galaxies. The masses and radii derived for these degenerate dark haloes of typical galaxies, without considering any cosmological information, agree with the Jeans masses and radii of a cosmological background of ~(20−30) g^{−1/4} eV degenerate fermions at redshift z ≈ 50. However, degenerate fermion objects of 10^{15} M_{⊙} composed of these particles are too small to constitute the dark halo of galaxy clusters. We also derive degeneracy conditions for hot and cold dark matter fermions.
Key words: dark matter / dark energy
© ESO, 2013
1. Introduction
It is accepted that most matter in the universe is dark. This belief comes from the fact that more mass than visible matter is required to explain the observed structure. At the end of the 1980s, it was suggested that baryonic dark matter could also be clumped in dark clusters (for a review, see for example, Carr 1994). The search for baryonic dark objects started at the beginning of the 1990s (see, for example, Moscoso 1993; Kerins & Carr 1994; Wasserman & Salpeter 1994). However, baryonic dark matter does not seem to be sufficient to describe the rotation curves observed in galaxies (see, for example, Caldwell 1995). Moreover, baryondominated models are ruled out by the observations of anisotropies of the cosmic microwave background (CMB; see, for example: Bond 1988; Strukov et al. 1987). Indeed, cosmological evidence favors nonbaryonic dark matter (e.g., weakly interacting massive particles; see, for example, Trimble 1987), which reproduces rotational curves. The most natural nonbaryonic dark matter candidate is the massive neutrino, but there are other alternatives such as the axion and various supersymmetric particles. The main argument in favor of nonbaryonic dark matter comes from the primordial or big bang nucleosynthesis (BBN): it estimates a baryonic energy density much lower than present estimates of the energy density of matter. In addition, dominant nonbaryonic dark matter models are consistent with other cosmological observations (the angular power spectrum of CMB, the largescale structure of the universe (LSS), lensing clusters), which also impose constraints on them (see, for example, Peacock 1999). In this work we will deal with fully degenerate fermionic dark matter. Semidegenerate massive neutrinos, with different degrees of degeneracy, have been extensively addressed; see, for example, Freese et al. (1983); Kang & Steigman (1992); Lesgourgues & Pastor (1999); Lesgourgues & Liddle (2001); Dolgov et al. (2002b); and Dolgov (2002a).
The evolution of the nonbaryonic dark matter distribution is collisionless. According to the Liouville theorem, the phasespace density of a collisionless system is preserved in its evolution. Therefore, the maximum phasespace density is conserved. Tremaine & Gunn (1979) applied this result to leptons. They used the condition that the maximum phasespace density must not increase and that the energy density of the cosmological background of leptons must be lower than that of matter. They assumed that initially the distribution function is that of FermiDirac, very far from degeneracy, i.e., neglecting the chemical potential, and they took a Maxwellian velocity dispersion for the final relaxed cluster. Tremaine and Gunn concluded that massive galactic halos cannot be composed of stable neutral leptons of mass ≤1 MeV, which are relativistic when they decouple from radiation. Assuming a typical isothermal massive galactic halo, the phasespace constraint for neutrinos is mc^{2} > 33 eV (see, for example, Peacock 1999), close to the result required for a criticaldensity universe. Thus, nondegenerate isothermal neutrinos cannot be the components of these dark haloes. However, neutrinos could constitute dark haloes in clusters of galaxies; for them, the constraint gives mc^{2} > 1.5 eV. Leptons with masses >1 MeV are nonrelativistic at the decoupling, and they would not be ruled out by the phasespace density arguments.
In any case, the Tremaine and Gunn conclusion cannot be claimed for degenerate fermions such as those treated here. It must be taken into account that in a fully degenerate state, all quantum phasespace cells, from momentum 0 up to the Fermi momentum, are occupied. Therefore, the maximum phasespace density of fully degenerate fermions in bound clusters is equal to g/h^{3}, where g is the spin degeneracy factor. Because in the evolution the maximum phasespace density cannot increase, the maximum phasespace density of the cosmological background of such fermions should also have been g/h^{3}. Hence, no constraint such as that derived by Tremaine and Gunn can be derived for the mass of the fermion. Nevertheless, in this work we look for restrictions on the mass of degenerate fermions by imposing the degeneracy condition on fermions that are relativistic when they decouple and on fermions that are nonrelativistic at the decoupling.
There is strong evidence that the expansion of our flat universe (Efstathiou et al. 1999) is speeding up (see, for example: Perlmutter et al. 1997, 1999; Schmidt et al. 1998), which cannot be explained only in terms of matter and radiation. Thus, it is necessary to assume that another component is involved. This new component has been named “dark energy”, and its pressure (from a hydrodynamical point of view) must be a negative function of its density. The first models for dark energy were hydrodynamical (Steinhardt 1996; Turner & White 1997), but dynamical scalar field models were soon proposed (see, for example, Coble et al. 1997; Caldwell et al. 1998). When the dark energy density is assumed to be a constant and pressure is taken as minus density, a cosmological constant model is considered.
Pacheco et al. (1986) found that a cosmological constant imposes a minimum mass for the existence of Newtonian selfgravitating degenerate fermion systems. Though bound states were allowed by the repulsive effect of the cosmological constant, clusters could also suffer gravitational instability due to an excessive mass (see, for example, Shapiro & Teukolski 1983). In this work, we study both the lower bound for the cluster mass imposed by the repulsive effect of the cosmological constant and the upper bound imposed by the gravitational instability.
Finally, we calculate the Jeans mass and Jeans length of a cosmological background of Newtonian degenerate fermions as a function of their redshift of formation.
The paper is organized as follows. In Sect. 2, Einstein equations are written for a fluid composed of a pure degenerate Fermi gas where the cosmological constant is present; the boundary conditions to solve these equations are also shown. In Sect. 3, we assume a weak field approximation to deal with Newtonian clusters. Section 4 shows the process we have followed to impose the stability conditions that limit the existence of these clusters. Numerical results are presented in Sect. 5. Section 6 deals with the degeneracy of a cosmological background of fermions as a function of the fermion mass. The Jeans mass for Newtonian degenerate fermions is calculated in Sect. 7. Finally, Sect. 8 contains our conclusions.
2. Einstein equations
We consider a static spherically symmetric cluster whose components are perfect fluids fulfilling different equations of state. The energystressmomentum tensor T^{μν} is then given by (see, for example, Glendenning 1997) (1)In Eq. (1), P and ρ are the pressure and the energy density of the fluid; g_{μν} is the metric which, using polar coordinates (x^{0} = ct, x^{1} = r, x^{2} = θ, x^{3} = φ), is given by (2)(ds^{2} = B(r)(cdt)^{2} − A(r)(dr)^{2} − r^{2} [(dθ)^{2} + sin^{2}(dφ)^{2}] is therefore the line element); U^{0} = (dx^{0}/ds) = B^{−1/2} and U = (dx/ds) = U^{0}(v/c) = 0 are the components of the fourvelocity tensor (they come from assuming a static cluster, so v = (dx/dt) = 0, and from U^{μ}U_{μ} = 1). B(r), A(r), ρ(r) and P(r) must be solutions of Einstein equations.
We consider a cluster composed of degenerate fermions of mass m in a cosmological background of dark energy treated by the cosmological constant model. Then, the total energy density and pressure are: where F and Λ represent fermions and dark energy, respectively.
In a degenerate state, fermions have the lowest energy without violating the Pauli principle. In this state, all quantum phasespace cells, from zero momentum up to Fermi momentum p_{F}(r) (maximum momentum that a fermion can have at a distance r from the cluster center), are occupied by g fermions, where g is the degeneracy per kinematic state. Hence, the distribution function, f(p), called Fermi distribution function (see, for example, Kippenhahn & Weigert 1994), is (5)Thus, the number of fermions per unit volume whose absolute value of momentum lies between p and p + dp is 4πf(p)p^{2}dp. Eq. (5) can be compared with the FermiDirac distribution function (see Landau & Lifshitz 2002), (6)In Eq. (6), μ and T are the chemical potential and the temperature of fermions, and E is the kinetic energy related with the momentum, p, by (7)(Equation (6) can also be expressed as f_{FD} = (g/h^{3}) [exp(U/K_{B}T − μ^{∗}/K_{B}T) + 1] ^{1}, where U = E + mc^{2} is the total energy and μ^{∗} = μ + mc^{2}; see, for example, Raychaudhuri et al. 1992). As can be seen, when μ ≫ K_{B}T, Eq. (6) coincides with Eq. (5) expressed as a function of the kinetic energy, E, (8)In Eq. (8), E_{F} is the Fermi kinetic energy, so (9)The number density of particles, n_{F}, is then (10)and as a function of the dimensionless Fermi momentum, x = p_{F}/mc, (11)(the end of a gravitationally bound system at r = R is then characterized by the fact that no cell is occupied, so p_{F}(R) = 0, and therefore n_{F}(R) = 0). Energy density, ρ_{F}, and pressure, P_{F}, of degenerate fermions are given by where U(p) = (p^{2}c^{2} + m^{2}c^{4})^{1/2} and v(p) = pc/(p^{2} + m^{2}c^{2})^{1/2} are the total energy and velocity of a fermion with momentum p, and For dark energy, we assume a cosmological constant model. Because the metric tensor g^{μν} fulfills , we can write a modified set of field equations, G^{μν} + Λg^{μν} = −(8πG/c^{4})T^{μν} from the Einstein equations (G^{μν} being Einstein tensor) that are also consistent with the conservation laws, . Λ is the cosmological constant and was introduced by Albert Einstein. Therefore, the existence of an energymomentum tensor of the vacuum can be considered, . From a hydrodynamical point of view, will agree with Eq. (1) if we take Defining Ω_{Λ0} = ρ_{Λ}/ρ_{c0}, where is the critical energy density of the universe at present (see, for example, Padmanabhan, 1995) and H_{0} is the Hubble function at present, ; i.e. (18)where (19)Values for H_{0} and Ω_{Λ0} can be taken from Spergel et al. (2003).
In the absence of other forces, except gravitation, the nontrivial components of the field equation, , read as (see, for example, Weinberg 1972)
where R_{μν} is the Ricci tensor and T_{μν} is that of Eq. (1), with ρ and P given by Eqs. (3) and (4). In Eqs. (20)−(22), a prime represents d/dr. A fourth nonindependent equation to (20)−(22) comes from the conservation of the energymomentum tensor, ; is identical to null, and gives the equation of the hydrostatic equilibrium of the system, i.e., (23)And a fifth equation is also obtained from Eqs. (20)−(22), as follows: (24)Now, let us define a length scale, a, given by Inserting (25) and (26) in (20)−(24), and using Eqs. (3), (4), (12), (13), (16) and (17), we have In Eqs. (27)−(31), a dot means d/dy. The value of a given in Eq. (26) was chosen so that the coefficient on the right hand term in Eqs. (27)−(29) is 3/2.
From Eqs. (27)−(31) coupled differential equations for x and A can be obtained. Using Ḃ (from Eq. (30)), Ȧ (from Eq. (31)) and (from the differentiation of Eq. (30) with respect to y) in Eq. (27), and after some manipulations, we obtain (32)where (33)The differentiation of Eq. (31) with respect to the dimensionless radius y gives (34)where (35)If the cosmological constant terms were not considered in Eqs. (32) and (34), we would have the wellknown TolmanOppenheimerVolkoff equations and universal solutions for x and A, i.e., solutions that are m independent. However, as cosmological constant terms (a^{2}Λ) depend on the fermion mass, it will be necessary to solve Eqs. (32) and (34) for each m.
From Eqs. (32) and (34), the behavior of x and A for y → 0 can be deduced. Up to second order in y, we obtain According to Eqs. (36) and (37), the coupled differential equations (32) and (34) must be integrated with the four boundary conditions: (38)The mass of the clusters inside a distance r = ay is (39)It can also be calculated by integrating Eq. (31), (40)from 0 up to some y. Thus, (41)When A is close to unity, it is better to use Eq. (39) to calculate M(y). The radius of the cluster will be the distance, R_{c} = ay_{c}, that, as has been said above, fulfills x(y_{c}) = 0; so, the cluster mass will be M_{c} = M(y_{c}).
The number of fermions inside a sphere of dimensionless radius y is determined by (42)where n_{F} is given by Eq. (10).
From Eq. (41), A(y) can be expressed as (43)With respect to B, the integration of Eq. (30), (44)from a point y up to y_{c} (i.e., from B(y) up to B(y_{c}) and from x(y) up to x(y_{c}) = 0) leads to (45)B(y_{c})can be calculated from Eqs. (27) and (29) for y ≥ y_{c} where u = v = 0. Adding the product of B and Eq. (27) to the product of A and Eq. (29), we have (46)i.e., d(ln(AB))/dy = 0, and therefore, (47)At y = ∞, the metric given by Eq. (2) has to be the Minkowski metric (that of special relativity), so A(y = ∞) = B(y = ∞) = 1; therefore (48)And from Eq. (43) evaluated at y_{c}, (49)
3. Weak field approximation
When fermions are nonrelativistic (x ≪ 1), the weak field approximation can be considered; i.e., So, taking into account Eqs. (43), (45) and (49), the leading terms in γ and β are Using Eqs. (50) and (51) in Eq. (29), and taking into account that we have (56)Equation (56) is the classical potential equation of gravity with (57)φbeing the gravitational potential; i.e., (58)To derive a differential equation for x from Eq. (58), we can replace β by φ (see (57)) in Eq. (53), i.e., (59)and insert this equation in (58). Thus, we have (60)From Eq. (60), the behavior of x in the neighborhood of the cluster center, up to second order in y, is given by (61)Evidently, Eq. (36) converts into Eq. (61) when x_{0} ≪ 1 (taking into account that ). Thus, from Eq. (61), Eq. (60) must be integrated with the boundary conditions (62)Using (54) in (39), (39) in (52), and (52) in (50), A can be expressed as (63)so, (64)The integral in Eqs. (63) and (64) can be expressed as a function of x and y as follows: integrating Eq. (56), (65)from 0 up to some y, we have (66)and as from Eq. (53), , (67)Hence Instead of solving Eq. (60), we can also make the following changes of variable Equating Eqs. (59) and (70), and using (71), we obtain that (72)And using Eqs. (70)−(72) and (25) in (58), we have (73)Equation (73) depends on two length scales (a, given by Eq. (26), and the free parameter b), on the mass of the cluster (M_{c}), and on the cosmological constant (Λ). However, by taking (74)we can make Eq. (73) depend on just one parameter. Thus, Eq. (73) reads (75)where (76)Finally, using the value of a, given in Eq. (26), in Eqs. (74) and (76), b and α are i.e., In the neighborhood of the cluster center, the behavior of χ up to second order in z is (81)Comparing (61) and (81), and taking into account (72) and (74), it can be seen that (82)Equation (75) has to be integrated for each value of α, imposing the boundary conditions (83)which derives from Eq. (81). For each α, χ_{1} must be that which causes the mass condition to be fulfilled; i.e., (84)where R_{c} = bz_{c} with χ(z_{c}) = 0. So using Eqs. (11), (72) and (79) in Eq. (84), (85)must be fulfilled.
Finally, when the cosmological constant is neglected, χ_{1} = 0.264 and z_{c} = 5.10, so the cluster radius is given by (86)
4. Minimum and maximum mass allowed for degenerate fermion clusters
As pointed out in Sect. 2, because the Λ terms in Eqs. (32) and (34) depend on fermion mass, it is necessary to solve these equations for each m. Thus, the discussion presented in the following paragraph must be made for each fermion mass.
The repulsive effect of the cosmological constant is already present in the neighborhood of the cluster center, as can be seen from Eqs. (36) and (37). According to Eq. (36), models with , fulfilling (87)show dx/dy_{y = 0} > 0 and are therefore unstable (see, for example, Jeans 1919; Shapiro & Teukolski 1983). As we take greater x_{0}’s, with , the instability condition dx/dy > 0 appears at greater dimensionless radii y’s (by solving Eqs. (32) and (34)). This will continue to happen until some x_{0} = x_{0,m}, for which the instability condition appears exactly at y = y_{c,m} with x(y_{c,m}) = 0. The mass inside the sphere of dimensionless radius y_{c,m}, calculated from Eq. (39), is denoted by M_{m}. Clusters with x_{0} > x_{0,m} do not fulfill dx/dy > 0 at any y ≤ y_{c}, with x(y_{c}) = 0; therefore, M_{m} is a lower limit for the mass of bound degenerate fermion clusters. The minimum number of fermions, N_{m}, that build a bound fermion cluster can then be calculated from Eq. (42) taking y = y_{c,m}.
Pacheco et al. (1986) found, using Eq. (75) for nonrelativistic clusters, that α > α^{∗} = 6.29 × 10^{3} does not allow bound states. Indeed, for α > α^{∗}, dx/dz > 0 at some z, so such a cluster will be unstable. Therefore, for each fermion mass, m, α^{∗} fixes the lower mass, M_{m}, for bound nonrelativistic clusters. Taking α = α^{∗} in Eq. (75) and integrating that equation, the condition (85) is fulfilled for χ_{1} = 0.180 and z_{c} = 7.71. Thus, minimum mass, (from Eq. (80), taking α = α^{∗} = 6.29 × 10^{3}), minimum fermion number, , and maximum radius, (using (79) in (71), with z = z_{c} = 7.71 and (88)), are Therefore, (91)Equation (60) could also be used to determine M_{m}; however, it would be necessary to calculate x_{0,m} for each fermion mass, so the calculation time would be longer.
Fig. 1 Mass limits, M_{M} (upper continuous line) and M_{m} (lower continuous lines), versus fermion mass, m, for dark clusters of degenerate fermions. is assumed. For details see the text. 

Open with DEXTER 
So far, we have dealt with lower bounds in the cluster mass. Now, we study upper bounds caused by the gravitational instability. There are x_{0} > x_{0,m} models that are not stable. An upper limit for x_{0}, and therefore for the mass of clusters, comes from the condition (92)Equation (92) indicates that stability changes at x_{0,M}. Our x_{0,m}<x_{0} models are stable for x_{0} < x_{0,M} fulfilling dM_{c}/dx_{0} > 0, and unstable for x_{0} > x_{0,M} fulfilling dM_{c}/dx_{0} < 0. Therefore, M_{M} = M_{c}(x_{0,M}) is the OppenheimerVolkoff limit for the mass of degenerate fermion clusters. When the repulsive effects of the cosmological constant are negligible, , and are given by As a brief summary of this section, we emphasize that Eqs. (88)−(90) are nonrelativistic results including repulsion effects, whilst Eqs. (93)−(95) are pure relativistic results that do not take repulsion into account.
A dimensional analysis of the lower bounds and upper bounds for fermions and bosons is described in Membrado & Pacheco (2012).
5. Numerical results and discussion
Numerically exact results for M_{m} and M_{M} versus fermion mass, m c^{2}, up to 10 g^{−1/4} GeV, are shown in the lower and upper curves of Fig. 1.
If relativistic effects are not taken into account in M_{m} (from Eq. (88), ) and if repulsive effects in M_{M} are not considered (from Eq. (93), ), then, for eV. However, for this fermion mass, relativistic and repulsive effects can not be neglected. Our exact calculation, for , gives m_{∗}c^{2} = 2.38 × 10^{2} g^{−1/4}eV, M_{∗} = 1.66 × 10^{21}M_{⊙} (see the point in Fig. 1 where the two continuous lines intersect) and R_{∗} = 1.26 × 10^{3}Mpc. m_{∗} is, therefore, the minimim fermion mass that can build bound degenerate fermion clusters.
We can now compare values of M_{m} obtained by neglecting relativistic effects, (see Eq. (88)), with those shown in Fig. 1. For fermion masses 1, 0.1, 0.04 and 0.02 g^{−1/4} eV, is about 0.1%, 0.6%, 10% and 72%. Thus, we can say that relativistic effects on M_{m} appear for mc^{2} < 0.1 g^{−1/4} eV. With respect to M_{M}, the repulsive effects of the cosmological constant begin to be appreciable for mc^{2} < 0.1 g^{−1/4} eV; for m_{∗}c^{2}, M_{M} = M_{∗} = 1.66 × 10^{21} M_{⊙} and , so .
Fig. 2 Cluster radii, R, versus fermion mass, m, for several cluster masses, M. is assumed. Horizontal lines correspond to 1, 10, 30 kpc and 1 Mpc. For details see the text. 

Open with DEXTER 
Numerical results for the radii of clusters, R, versus fermion mass, m, for several cluster masses, M, are shown in Fig. 2. In this figure, each line represents radii allowed for degenerate fermion clusters of the same mass, M. The maximum radius shown in each line corresponds to the limit imposed by the cosmological constant. The minimum radius in each line is the OppenheimerVolkoff limit.
We consider the mass of a typical galaxy dark halo (G). This mass would be on the order of ~10^{11} M_{⊙}. According to Eq. (88) the minimum fermion mass that creates a degenerate cluster with is ; from Eq. (90) we know that its radius is . Thus, the contributions to the rotational velocity would be about , too small to reproduce the rotational curve of a typical galaxy. The stability limit for is reached with m_{M,G}c^{2} = 2.97 × 10^{3} g^{−1/4} eV (see Eq. (93)); in that case, from Eq. (95), , so v_{h} = 10^{5} km s^{1}. Therefore, that degenerate cluster has too small a radius and too great a v_{h} to be the dark halo of a typical galaxy. However, as we show below, there are fermion masses between m_{m,G} and m_{M,G} that could give rise to dark clusters whose v_{h} could lead to the observed values of the rotational velocities of typical galaxies.
We now look at the mass of a dark halo of a typical galaxy cluster (GC), i.e., ~10^{15} M_{⊙}. If in Eq. (88) we take , we conclude that the minimum fermion mass that can build a degenerate fermion object with such a mass is ; and from Eq. (90), its radius is . v_{h} is then ≈650 km s^{1}, which is too small to reproduce the velocities of galaxies in clusters. It can also be seen that the maximum fermion mass that creates a 10^{15} M_{⊙} degenerate cluster is, from Eq. (93), m_{M,GC}c^{2} = 2.97 × 10^{1} g^{−1/4} eV; in that case, from Eq. (95), , so v_{h} ≈ 10^{5} km s^{1}. Therefore, the radius is too small and v_{h} too high with respect to the size of a typical galaxy cluster and the velocities of their galaxies. We can also consider the case of 1 Mpc, which has about the visible size of a typical galaxy cluster. Figure 2 shows that a degenerate cluster of 10^{15} M_{⊙} and 1 Mpc (so, v_{h} = 2 × 10^{3} km s^{1}) can be built using a 2 g^{−1/4} eV fermion (see Eq. (86)). However, this particle cannot build clusters with masses smaller than 5.8 × 10^{13} M_{⊙} (see Eq. (88)), so it cannot be the component of the dark halos of galaxies.
It can be seen from Eq. (86) that fermions with masses between 15 g^{−1/4} and 30 g^{−1/4} eV could create bound degenerate dark haloes that would contribute to the rotational velocity with values, v_{h} = (GM_{c}/R_{c})^{1/2}, between 150 and 300 km s^{1}; these haloes would have radii between 15 and 50 kpc and masses between 8 × 10^{10} M_{⊙} and 7 × 10^{11} M_{⊙}. However, they could not compose the dark halo in a galaxy cluster, unless they were composed of dark galaxies. We consider some examples. For 20 g^{−1/4} eV, M_{c} = (2−6) × 10^{11} M_{⊙} degenerate dark clusters show R_{c} = (38−27) kpc and v_{h} = (150−313) km s^{1}; and those with M_{c} = (10^{12}−10^{15}) M_{⊙} present R_{c} = (22−2.2) kpc and v_{h} = (4.4 × 10^{2}−4.4 × 10^{4}) km s^{1}, and therefore are incompatible with the observations of the largest galaxies and clusters of galaxies. For 25 g^{−1/4} eV, M_{c} = (1−4) × 10^{11} M_{⊙} clusters are built with R_{c} = (26−17) kpc, showing v_{h} = (129−318) km s^{1}; while for M_{c} = (10^{12}−10^{15}) M_{⊙}, R_{c} = (12−1.3) kpc and v_{h} = (5.9 × 10^{2}−5.9 × 10^{4}) km s^{1}, so they are again incompatible with observations.
For fermion masses smaller than 15 g^{−1/4} eV, only clusters with masses greater than 10^{12} M_{⊙} are able to give the kind of v_{h} that is needed to reproduce the values of the rotational curves of typical galaxies; so, they cannot be the dark haloes of galaxies. The same conclusion is obtained for fermions with masses greater than 30 g^{−1/4} eV, but now because degenerate clusters show radii smaller than those of typical galaxies.
We now look at the problem of dwarf galaxies. Dwarf galaxies are the most abundant type of galaxy in the Universe. The largest dwarf irregulars have stellar masses of ~10^{9} M_{⊙} and are rotationally supported, while dwarf spheroidal galaxies have stellar masses from 10^{7} M_{⊙} to below 10^{6} M_{⊙} with no rotational support (the lightest known is Segue 1 with a mass of ~6 × 10^{5} M_{⊙}; see, for example, Simon et al. 2011).
Highresolution simulations of the standard ΛCDM (cosmological constant and cold dark matter) cosmology (see, for example: Klypin et al. 1999; Moore et al. 1999; Diemand et al. 2007; Springel et al. 2008) were applied to the Local Group; they showed a total number of dark matter subhaloes much greater than that of dwarf galaxies predicted using luminosity functions corrected by completeness and bias (Tollerud et al. 2008; Koposov et al. 2008) (one order of magnitud higher than the number of observed dwarf galaxies). It has been attempted to solve this problem found in ΛCDM simulations by suggesting that the supernova feedback could lead to the loss of baryonic components in haloes with masses smaller than a critical value, and that the effect of the photoionisation could prevent star formation in the smallest haloes (see, for example: Larson 1974; Dekel & Silk 1986; Ferrara & Tolstoy 2000; Efstathiou 1992; Somerville 2002; Hoeft et al. 2006; Simon & Geha 2007). These mechanisms would make the dark matter masslight ratio very high in haloes on the order of or smaller than ~10^{9} M_{⊙} (Sawala et al. 2010).
Another problem found in simulations of dark halos using the standard ΛCDM cosmology comes from dark halos of ~10^{10} M_{⊙}. ΛCDM simulations of ~(1−2) × 10^{10} M_{⊙} dark halos show the formation of dwarf galaxies with stellar masses between 5 × 10^{7} M_{⊙} and 2 × 10^{8} M_{⊙} (Pelupessy et al. 2004; Stinson et al. 2007; Governato et al. 2010; Sawala et al. 2011). Guo et al. (2010) compared the abundance of dark haloes in the Millennium and MillenniumII Simulations to the observed abundances of galaxies as a function of the stellar mass obtained from SDSS DR7 by Li & White (2009). They found that the abundance of ~10^{10} M_{⊙} dark haloes obtained from ΛCDM simulations, and whose stellar mass is ~5 × 10^{7} M_{⊙}, is about the abundance of observed galaxies with stellar mass of ~10^{6} M_{⊙}. Therefore, ΛCDM simulations predict galaxies with about fifty times more stellar mass than observed. In any case, if dark halos of 10^{10} M_{⊙} contain dwarf galaxies with stellar masses of 10^{8} M_{⊙}, a ΛCDM universe would predict about four times more of those galaxies than observed.
Thus, ΛCDM simulations present serious problems for describing dwarf galaxies. In this respect, it should be said that warm dark matter models (ΛWDM) predict about three times fewer dark haloes of 10^{10} M_{⊙} than ΛCDM models for m_{WDM} = 1 keV (Zavala et al. 2009).
We now look at our results for small degenerate dark haloes. We take as an example the case of fermions of 20 g^{−1/4} eV. According to Eq. (88), bound dark haloes with masses smaller than cannot exist due to the repulsive effect of the cosmological constant. These small degenerate dark haloes would be very extensive (a degenerate dark halo of M_{c} = 10^{10} M_{⊙} would have a radius R_{c} ≈ 100 kpc). However, these extensive degenerate dark haloes should be truncated by the tidal effects of larger neighboring galaxies (below we show as an example the case of the Large Magellanic Cloud). With respect to their contributions to the rotational velocities, v_{h} = (GM_{c}/R_{c})^{1/2}, we see that they are small; thus, for a dark cluster of 10^{10} M_{⊙}, v_{h}(100 kpc) ≈ 20 km s^{1}. We may also wonder about the mass and velocity contribution at 5 and 10 kpc, which could be the size of dwarf galaxies contained in the dark haloes. For a 10^{10} M_{⊙} cluster, M(5 kpc) ≈ 6.7 × 10^{6} M_{⊙} and v_{h}(5 kpc) ≈ 2.4 km s^{1}; M(10 kpc) ≈ 5.4 × 10^{7} M_{⊙} and v_{h}(10 kpc) ≈ 4.8 km s^{1}.
We now consider the case of the Large Magellanic Cloud (LMC). This is a nearby barred spiral galaxy, a satellite of the Milky Way (MW). It is located at D ≈ 48.5 kpc from the center of MW. Up to 8.9 kpc its stellar disk mass is ~2.7 × 10^{9} M_{⊙} (de Vaucouleurs & Freeman 1972; Bothum & Thompson 1988; Zaritsky 1999; Bell & de Jong 2001) and the neutral gas mass is ~0.5 × 10^{9} M_{⊙} (Kim et al. 1998); these masses contribute to the circular velocity of 39.3 km s^{1}. The circular velocity derived by Van der Marel et al. (2002) at 8.9 kpc is 64.8 ± 15.9 km s^{1}; consequently, a dark halo must exist in this galaxy.
The LMC suffers the tidal effects exerted by the MW. The gravitational interaction between LMC and MW is a very complex problem. But, for the sake of simplicity and merely to illustrate the problem, we have made very considerable simplifications. Therefore, the results described in the following must be considered as very rough estimates.
For the Milky Way we take, as an example, a degenerate dark halo composed of fermions of 23 g^{−1/4} eV with a mass of M_{MW} = 2.5 × 10^{11} M_{⊙}, and therefore a radius (see Eq. (86)) of R_{MW} = 24.4 kpc. At the Sun position, R_{⊙} ≈ 8.5 kpc, the contribution to the rotational velocity is v_{h} ≈ 158 km s^{1}. Using the galactic bulgedisk model developed by Klypin et al. (2002), and taking for the bulge and disk masses, M_{B} = 7 × 10^{9} M_{⊙} and M_{D} = 3.5 × 10^{10} M_{⊙}, the rotational velocity at 8.5 kpc is , v_{d} and v_{b} being the contributions of the bulge and the disk (we have taken M_{D}/M_{B} = 5, as happens with the luminosity ratio L_{D}/L_{B} ≈ 5 (see Kent et al. 1991); the kband luminosity of the Galaxy measured by Kent et al. (1991) is L_{K} = 6.7 × 10^{10} L_{⊙}, and measured by Drimmel & Spergel (2001), L_{K} = 8.9 × 10^{10} L_{⊙}). At 20 kpc, v_{h} ≈ 225 km s^{1}, which together with the contributions of the disk and bulge make v_{r} ≈ 250 km s^{1} (the rotation curve derived by Honma & Sofue (1996) based on the geometry of the HI disk shows v_{r} ≈ 265 ± 30 km s^{1}).
For the LMC, we could propose a dark halo composed of 23 g^{−1/4} eV degenerate fermions which, if isolated, would have M_{c} = 8 × 10^{10} M_{⊙} (according to Eq. (86), its radius would be R_{c} = 35.7 kpc). However, this halo should be truncated by the gravitational effects of the MW. We can derive a rough estimate of a tidal radius, r_{T} and a mass of the degenerate dark halo, M_{LMC} ≈ M(r_{T}), as those that fulfill . This equation is the balance of strengths per unit mass in the direction that joins the centers of the LMC and the MW. The first term on the left hand denotes the rotation of the LMC around the MW () and the other terms are the gravitational strength per unit mass due to the LMC and the MW at a distance r_{T} from the center of the LMC. For this calculation it is necessary to know M(r). This can be derived by solving Eq. (75), neglecting the cosmological constant, and with the boundary conditions given by (83) and (85). Thus, we have obtained that R_{LMC} = r_{T} ≈ 14.5 kpc and M_{LMC} ≈ 2.3 × 10^{10} M_{⊙}. Assuming that the inner structure of the degenerate dark cluster is not changed by the truncation, the halo mass and the contribution to the circular velocity of LMC at r = 8.9 kpc are M(8.9 kpc) ≈ 6.5 × 10^{9} M_{⊙} and v_{h}(9 kpc) ≈ 56 km s^{1}. Accordingly, the total circular velocity would be v_{r} = 69 km s^{1}, within the range given by Van der Marel et al. (2002). Smaller masses than M_{c} lead to smaller r_{T}, v_{h} and v_{r}, while greater masses than M_{c} lead to greater r_{T}, v_{h} and v_{r}.
With respect to dwarf spheroidal galaxies, we can consider the ultrafaint Milky Way satellite galaxy Segue 1 (Simon et al. 2011; Martinez et al. 2011). This is located at a distance of 23 kpc from the MW center. The stars in this galaxy show a velocity dispersion σ ≈ 3.7 km s^{1}. This velocity dispersion could imply a dark halo mass of M_{1/2} ≈ 5.8 × 10^{5} M_{⊙} within a sphere with a radius r_{1/2} ≈ 36 pc that encloses half the stellar luminosity of the galaxy. As the Vband luminosity of Segue 1 is 240 L_{⊙}, the Vband masstolight ratio within the halflight radius is ≈3400 M_{⊙}/L_{⊙}. Assuming for the dark halo of Segue 1 a mass of 6 × 10^{5} M_{⊙} and using the same formula as we used above to derive the tidal radius of LMC, we obtain that the tidal radius of Segue 1 is ≈240 pc; this is much greater than r_{1/2}. This means that Segue 1 should not suffer tidal effects from the MW.
We assume that the dark halo of Segue 1 is composed of degenerate fermions like those proposed here to describe the dark halo of the MW. In that case, the condition that within a sphere of 36 pc there is a mass of 5.8 × 10^{5} M_{⊙} is fulfilled by a degenerate dark halo with a total mass of ≈3 × 10^{12} M_{⊙} and a radius of ≈10 Kpc. This has no sense because the MW would suffer strong tidal effects. Therefore, the dark halo of Segue 1 cannot be composed of degenerate fermions like those proposed in this paper to describe the dark halo of typical galaxies. However, if the halo were composed of degenerate fermions of ≈1.3 keV, a cluster with total mass M_{1/2} and radius r_{1/2} could be reproduced (see Eq. (86)). The problem with this assumption is that, as can be seen from Eq. (125) shown below, such fermions are not degenerate; so Eq. (86) cannot be used. Therefore, such a nondegenerate fermion cluster with mass M_{1/2} should have a radius greater than r_{1/2}.
Hence, it is apparent that to adequately model the dark halo of Segue 1 it is necessary to have more compact selfgravitating spheres than those provided by degenerate fermions. In this sense, spheres formed by fermions that are not completely degenerate are of no help because they are obviously less compact than their degenerate counterparts. Selfgravitating boson spheres, however, might be an answer. One kind of these objects is the selfgravitating collisionless Newtonian boson cluster in which all bosons occupy the lowest lying oneparticle Hartree orbital (Membrado et al. 1989a,b; Membrado & Aguerri 1996). For these condensates, the radius R_{BC}, which covers 99% of the assembly, is , where M_{BC} is the mass of the boson cluster and m_{B} the boson mass. An upper limit for the mass of the boson that could reproduce the dark halo of Segue 1 can be obtained assuming M_{BC} = M_{1/2} and R_{BC} = r_{1/2}. Thus, m_{B}c^{2} < 2 × 10^{20} eV. It should be said that invoking the uncertainly principle, it may be deduced that the linear momentum of each boson in the cluster is on the order of p_{B} ~ ħ/R_{BC} and hence that the deBroglie wavelength (nonrelativistic particles) for the boson is derived as R_{BC}. Then, the phasespace condition for degeneracy in this condensate of bosons is trivially fulfilled: within a deBroglie cell there are M_{BC}/m_{B} particles.
Scalar fields fulfilling different potentials have been proposed to model dark matter (see, for example, Lee & Koh 1996; Sahni & Wang 2000; Matos & UreñaLópez 2000, 2001). The mass of this scalar dark matter, m_{B}, ranges from 10^{26} eV up to 10^{23} eV. The dark clusters composed of these fields show a minimal cutoff radius, greater than their Compton legth Λ_{C} = ħ/m_{B}c. Indeed, the enormous Compton length of these particles would prevent structure formation on small subgalactic scales. In this sense, collisionless boson spheres do not have this limitation.
Finally, we look at isolated dwarf galaxies of visible radii of about 10 Kpc. They show rotational velocities ranging between 50 and 150 km s^{1} (see, for example, Kent 1987). Fermions with masses between 15 g^{−1/4} and 30 g^{−1/4} eV could build more extensive degenerate structures that at 10 kpc would show masses in the range (0.6−5.2) × 10^{10} M_{⊙}. These would lead to rotational velocity contributions, v_{h}, similar to those shown by the visible dwarf galaxies. From the calculation of M(r) for degenerate fermion clusters (by integrating Eq. (75) with the boundary conditions Eqs. (83) and (85)), it is deduced that M(r) ∝ r^{3} for approximately r < 0.4R_{c}, so v_{h}(r) ∝ r. Incidentally, this behavior of v_{h}(r) is what is needed to describe the rotation velocity in some of those galaxies (see, for example, NGC 247, NGC 3109, NGC 2259 from Kent 1987).
6. Dependence on m of the degeneracy of a cosmological background of fermions
In this section we consider that the degenerate fermions decouple from radiation at an expansion parameter a_{D}. For simplicity, we assume that the decoupling occurs instantaneously. Before the decoupling, the fermions are in equilibrium with the radiation, so their temperatures are equal, i.e., T(a ≤ a_{D}) = T_{γ}(a ≤ a_{D}). We assume that at a_{D}, such fermions have a kinetic energy E_{D}, a chemical potential μ_{D} and a temperature T_{D} = T(a_{D}) = T_{γ}(a_{D}). Because we are considering degenerate fermions, μ_{D} ≫ K_{B}T_{D}. We also assume that during a_{NR} < a ≤ a_{0} degenerate fermions are nonrelativistic, and that their temperature at a_{NR} is T_{NR} = T(a_{NR}).
We distinguish between fermions that are ultrarelativistic at a_{D} (hot and warm dark matter) and nonrelativistic fermions at a_{D} (cold dark matter). The degeneracy condition μ_{D}/K_{B}T_{D} ≫ 1 together with the condition a_{NR} > a_{D} for the former and a_{NR} < a_{D} for the latter must lead to restrictions on the mass of fermions. Obtaining these restrictions is the aim of this section.
Fermions that decouple when they are ultrarelativistic could be e, μ and τ neutrinos. Direct observational bounds on neutrino masses, found kinematically, are: m_{νe} < 10 eV (Weinheimer et al. 1999; Lobashev et al. 1999; Particle Data Group 2000); m_{ντ} < 170 keV (Assmagan et al. 1996); and m_{ντ} < 18 MeV (Buskulic et al. 1995; Passalacqua 1997). In the absence of degeneracy, eneutrinos decouple at K_{B}T_{D}(ν_{e}) ≈ 2.0 MeV, and μ and τ neutrinos at K_{B}T_{D}(ν_{μ}) ≈ K_{B}T_{D}(ν_{τ}) ≈ 3.5 MeV (see, for example, Freese et al. 1983). For semidegenerate neutrinos, T_{D} is greater and depends on their degree of degeneration (see, for example, Kang & Steigman, 1992). Heavy neutrinos could be an example of fermions that are nonrelativistic when they decouple. For a review of neutrinos in cosmology, see, for example, Dolgov (2002a).
6.1. Degenerate fermions that decouple when they are ultrarelativistic
At a_{D}, before the decoupling, the degenerate fermions, which are in equilibrium with the radiation, follow the Fermi distribution function given by Eq. (8) (96)where E_{D} and μ_{D} are related with the momentum and the Fermi momentum at a_{D}, p_{D} and P_{F}_{D}, by If fermions are ultrarelativistic at a_{D}, p_{F}_{D}c/mc^{2} ≫ 1, only fermions with p_{D}c ≫ mc^{2} contribute appreciably to physical quantities. Therefore, from Eqs. (97) and (98) and then Eq. (96) expressed as a function of the momentum, p_{D}, is given by (101)At a_{D}, after the decoupling, the fermions travel along geodesics in spacetime. During this free propagation, the distribution function of the fermions is conserved (the expansion of the universe causes the proper volume occupied by any number of these fermions to increase as a^{3} while their volume in the momentum space is redshifted by a^{3}). This allows us to obtain their distribution function after the decoupling, . Thus, taking into account that as the universe expands, the fermion momentum, p, decreases as (102)the conservation of the distribution function implies , so (103)Finally, as the fermion temperature, T, cools as (104)even though the fermions are no longer in thermodynamical equilibrium with radiation (see, for example, Padmanabhan 1995), Eq. (103) reads as (105)Equation (105) tells us that the Fermi momentum, p_{F}, is (106)Hence, from Eq. (106), the fermions will be nonrelativistic when p_{F}c/mc^{2} = (μ_{D}/mc^{2})(T/T_{D}) ≪ 1, and T_{NR} can then be estimated taking p_{F}c/mc^{2} = 1; thus, (107)Dividing now Eq. (107) by the fermion temperature at present, T_{0} = T(a_{0}), and taking into account that T/T_{0} = a_{0}/a, (108)For a > a_{NR}, degenerate fermions are nonrelativistic. Then, their number density and energy density, n_{F} and ρ_{F}, are given by where the distribution function given by Eq. (105) has been used. Assuming a value for Ω_{F0} = ρ_{F0}/ρ_{c0}, the energy density at present will be (111)with (112)so, equating Eqs. (110) at a_{0} and (111), (113)Finally, using Eq. (113) in (107), we have (114)and, as T/T_{0} = a_{0}/a, (115)i.e., (116)As an example, at z = 1000 (a relatively short epoch after recombination when the growth of perturbations is fully given by linear theory and late enough for the decaying modes of perturbations to have decayed away), i.e., at a/a_{0} = 1/(1 + z) = 9.99 × 10^{4}, degenerate fermions are nonrelativistic if (117)At present, the fermion temperature, T_{0}, differs from that of radiation, T_{γ0}, due to different particle annihilation and decoupling occurring after a_{D}. If only e^{+}e^{−} annihilate, T_{0} = (4/11)^{1/3}T_{γ0} (see, for example, Padmanabhan 1995). Thus, expressing the temperature of fermions at present, T_{0}, as (118)with (119)the expansion factor of the universe at the fermion decoupling, a_{D}, can be expressed as (120)Thus, as degenerate fermions are relativistic when decoupling if a_{D}/a_{0} < a_{NR}/a_{0}, we deduce from Eqs. (116) and (120) that m must fulfill (121)If the fermions were nondegenerate, we could estimate a_{NR}/a_{0} assuming K_{B}T_{NR} ≈ mc^{2}/3. Thus, (122)Therefore, nondegenerate fermions would be relativistic when decoupling (a_{NR}/a_{0} > a_{D}/a_{0}) if (123)The degeneracy parameter ξ_{D} = μ_{D}/(K_{B}T_{D}) can be estimated from Eqs. (113) and (118). Thus, (124)Therefore, only fermions with mass m fulfilling (125)will be degenerate (ξ_{D} ≫ 1). Evidently, the degeneracy condition (125) is more restrictive for m than the relativistic condition at a_{D} given in (121).
Cosmological constraints on the neutrino degeneracy parameter began to be analyzed at the end of the past century (see, for example: Lesgourgue & Pastor 1999; Pastor & Lesgourgue 2000; Lesgourgue & Liddle 2001; Orito et al. 2002; Dolgov et al. 2002b). A relic neutrino degeneracy leads to a different amount of neutrinos and antineutrinos. Such a lepton asymmetry could have been imprinted on the cosmological data (Lesgourgue & Pastor 1999, 2006). Studies about the effects of neutrino mass and neutrino degeneracy on the CMB, on the LSS, and on the BBN have been made in recent years, and their results have been compared with cosmological observations (such as data from the Wilkinson Microwave Anisotropy Probe (WMAP) on anisotropies of the CMB, from the BBN, and from the Sloan Digital Sky Survey (SDSS) on the dark matter power spectrum). Nevertheless, current constraints are still weak. WMAP could allow values of ξ_{ν} up to 4 (Shiraishi et al. 2009); however the BBN imposed bounds on  ξ_{ν}  smaller than 1.1 (see, for example: Serpico & Raffet 2005; Lattanci et al. 2005; Hamann et al. 2008; Popa & Vasile 2008; Shiraishi et al. 2009; Castorina et al. 2012).
We did not specify the type of degenerate fermion or its interactions; therefore, no comparison with cosmological data is made. Indeed, the aim of this subsection is simply to derive an upper limit for the mass of degenerate fermions that decouple when they are ultrarelativistic (as a function of their decoupling temperature and their temperature and energy density at present). We have also derived an expression for the degeneracy parameter (as a function of m, T_{0} and ρ_{F0}).
6.2. Degenerate fermions that decouple when they are nonrelativistic
If degenerate fermions are nonrelativistic at a_{D}, p_{F}_{D}c/mc^{2} ≪ 1, and hence, p_{D}c/mc^{2} ≪ 1. Therefore, from Eq. (97), their kinetic energy, E_{D}, is related to their momentum, p_{D}, by (126)Therefore, using Eq. (126) in (96), the distribution function followed by these degenerate fermions at a_{D}, before the decoupling, is (127)Thus, taking into account the conservation of the distribution function together with Eqs. (102) and (104), the distribution function after the decoupling, , will be (128)Hence, the Fermi momentum, p_{F}, is now (129)According to Eq. (129), fermions are nonrelativistic when p_{F}c/mc^{2} = (2μ_{D}/mc^{2})^{1/2}(T/T_{D}) ≪ 1; so T_{NR}, estimated taking p_{F}c/mc^{2} = 1, is (130)and therefore, (131)Using Eq. (128) for nonrelativistic degenerate fermions, their number density and energy density, n_{F} and ρ_{F}, are given by Thus, equating (133) at a_{0} with (111), (134)Now, using Eq. (134) in (130), we have (135)and, therefore, (136)i.e., (137)Equations (135)−(137) are the same equations as derived for degenerate fermions that are ultrarelativistic at a_{D} (see Eqs. (114)−(116)). Finally, taking into account that degenerate fermions are nonrelativistic at a_{D} when a_{D}/a_{0} > a_{NR}/a_{0}, we have, from Eqs. (137) and (120), that m must fulfill (138)From Eqs. (134) and (118), the degenerate parameter, ξ_{D} is now, (139)So using (138) in (139), (140)Therefore, according to Eqs. (140), fermions coming from a cosmological background that decouple when they are nonrelativistic cannot be completely degenerate.
7. Jeans mass for Newtonian degenerate fermions
The Jeans length for a background of Newtonian degenerate fermions can be derived by studying a perturbation solution from the unperturbed state, which is taken as spatially uniform in the neighborhood of the perturbation. The simplest procedure is to decompose the perturbation of a magnitude χ into a system of plane waves with wave vector K and instability growth ω (see, for example, Binney & Tremaine 1987; in this paper we follow the work by Membrado & Aguerri 1996); i.e. χ = χ_{u} + χ_{p}exp [ωt + iK·r], χ being, in this section, the energy density, ρ (related with the density number, n, by ρ = mc^{2}n), velocity, v, and gravitational potential, φ. Magnitudes with subindex u describe the unperturbed solution, while those with subindex p describe small amplitudes of the perturbation. Inserting the perturbation and the unperturbed solution in the hydrodynamical equations and in the Poisson equation (143)and neglecting secondorder terms, we have
Equations (144)−(146) are therefore the equation of continuity, the equation of motion, and the Poisson equation. In Eq. (145), P is the pressure that fulfills an equation of state P = P(ρ). To simplify, we assume that the unperturbed state is at rest (v_{u} = 0).
We calculate the Jeans mass of a cosmological background of degenerate fermions at an expansion factor fulfilling a_{NR} < a < a_{0}, therefore we consider that the fermions are nonrelativistic.
For degenerate fermions that are ultrarelativistic at a_{D}, their pressure, P_{UR}, as a function of their energy density, ρ_{UR}, can be calculated from (147)where is given by (105) and v = p/m is the fermion velocity. Thus, taking into account that ρ_{UR} is given by Eq. (110), the equation of state is (148)The pressure, P_{NR}, of degenerate fermions that are nonrelativistic when decoupling is (149)where is given by (128) and v = p/m. So using ρ_{NR} from Eq. (133) leads to the same equation of state as that shown in (148), i.e., (150)Thus from Eq. (148) or (150) (151)Equations (144)−(146) together with (151) have a nontrivial solution if (152)Thus, there exists a critical value of K_{J} determined by the condition ω = 0. The Jeans length is then given by (153)and the Jeans mass, M_{J} = (4πρ_{u}/3c^{2})(λ_{J}/2)^{3}, is (154)If λ > λ_{J}3, i.e. if M > M_{J}, then ω^{2} > 0 and a cluster can be created.
Equations (153) and (154) can be used to estimate the Jeans length and mass at different redshifts. For this purpose, we assume that ρ_{u} is the background energy density of degenerate fermions at some z. Taking into account that we are dealing with nonrelativistic matter, ρ(z) = ρ(z = 0)(1 + z)^{3} (see, for example, Padmanabhan 1995), and using (111), we have (155)Therefore, We now consider redshift z = 50, when the nonlinear collapse could have started (see for example, Padmanabhan 1995). At this redshift, the Jeans mass and radius of a cosmological background of mc^{2} ≈ 20−30 g^{−1/4} eV degenerate fermions is and (see Eqs. (157) and (156)). These results could be consistent with a typical galaxy dark halo. In addition, these results based on cosmological background densities at z = 50 are coincident with those shown in Sect. 5, where no cosmological information was used (only the relation among fermion mass and bound degenerate cluster mass and radius). For those fermion masses, a Jeans mass of 10^{15} M_{⊙} is reached at , when R_{J} ≈ (2−0.7) kpc; therefore, dark halos of clusters of galaxies cannot be composed of such degenerate fermions. For M_{J} = 10^{10} M_{⊙}, and R_{J} ≈ 58−31 kpc.
At redshift z = 50, a background of degenerate fermions of mc^{2} ≈ 2.8 g^{−1/4} eV shows (according to Eq. (88), 2.8 g^{−1/4} eV fermions cannot build bound degenerate clusters of ), but a size of , smaller than the observed radii of typical clusters of galaxies. At z = 50, the fermion mass of ≈50 g^{−1/4} eV would lead to and ; this result could be consistent with the dark halo of small spiral galaxies; however, for such a fermion mass, the redshift and the Jeans radius for 10^{11} M_{⊙} is and R_{J} ≈ 3.8 kpc.
Finally, we can insert Eq. (157) in (156); thus, we find (158)Comparing the Jeans radius, R_{J} = λ_{J}/2, with the radius of a cluster of Newtonian degenerate fermions, R_{c} = 5.1 b, with b, given by Eq. (79), it can be seen that R_{J} = 0.9 R_{c}. Therefore, the Jeans radius for Newtonian degenerate fermions is approximately the radius of a bound system of degenerate fermions.
8. Conclusions
We have considered structures composed of fully degenerate fermions in a background of dark energy assumed to be the cosmological constant. This has led to the determination of masses allowed for those clusters as a function of fermion mass; i.e., for each fermion mass, m, a minimum allowed cluster mass, M_{m}, and a maximum allowed cluster mass, M_{M} have been deduced.
The minimum fermion mass that is able to build stable clusters of degenerate fermion is 0.024 g^{−1/4} eV, for which M_{m} = M_{M}.
Clusters composed of degenerate fermions with masses smaller than 0.1 g^{−1/4} eV are relativistic and simultaneously suffer the repulsive effects of the cosmological constant (for these masses, clusters with M_{m} are relativistic and those with M_{M} are affected by the cosmological constant).
We have shown that fermions of mc^{2} ≈ (15−30) g^{−1/4} eV could build bound degenerate dark haloes with contributions to the rotational velocities as those needed to reproduce the rotational curves of galaxies. However, these fermions could not compose a degenerate dark halo of 10^{15} M_{⊙} with a size on the order of or greater than the visual radius of clusters of galaxies. Nevertheless, the dark halo in galaxy clusters could be composed of dark galaxies.
If degenerate dark clusters of masses much smaller than 10^{11} M_{⊙} were the dark haloes of dwarf galaxies, their radius should be much greater than those of degenerate clusters of 10^{11} M_{⊙}. However, the tidal effects of larger neighboring galaxies would cause their sizes to be truncated and their masses to be reduced.
We saw that the dark halo of dwarf spheroidal galaxies cannot be modeled by degenerate fermions like those proposed here to describe the dark halo of typical galaxies. In fact, more compact selfgravitating spheres than those provided by degenerate fermions are necessary. Selfgravitating boson spheres might be an answer. In these systems, all bosons occupy the minimum Hartree orbital and the radius–mass relation of the cluster is of the type R ~ M^{1}.
We showed that the values of rotational velocities of isolated small galaxies of ~10 kpc can be reproduced with dark clusters of masses similar to those proposed for typical galaxies, i.e., on the order of ~10^{11} M_{⊙}.
We found degeneracy conditions for fermions that are relativistic when they decouple from the radiation and for fermions that are nonrelativistic at the decoupling. We saw that cold dark matter fermions cannot be fully degenerate.
We also derived the Jeans mass and radius of a cosmological background of Newtonian degenerate fermions as a function of the redshift. We saw that the Jeans masses and radii for degenerate fermions of mc^{2} ≈ (20−30) g^{−1/4} eV at z ≈ 50 are similar to the masses and radii derived (without considering any cosmological information) for the bound degenerate clusters that would be able to reproduce the rotation curves of galaxies. However, these degenerate fermions could not give rise to the dark halo of clusters of galaxies (a Jeans mass of about 10^{15} M_{⊙} would appear at z ≈ 10^{4}, with a Jeans radius on the order of 1 kpc).
References
 Assamagan, K., Brönnimann, Ch., Daum, M., et al. 1996, Phys. Rev. D, 53, 6065 [NASA ADS] [CrossRef] [Google Scholar]
 Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton University Press) [Google Scholar]
 Bond, J. R. 1988, in The Early Universe, eds. W. G. Unruh, & G. W. Semenoff (Dordrecht: Reidel), Proc. NATO Adv. Study Inst., 283 [Google Scholar]
 Bothun, G. D., & Thompson, I. B. 1988, AJ, 96, 877 [NASA ADS] [CrossRef] [Google Scholar]
 Buskulic, D., Casper, D., de Bonis, I., et al. (ALEPH Collaboration) 1995, Phys. Lett. B, 349, 585 [NASA ADS] [CrossRef] [Google Scholar]
 Caldwell, D. O. 1995, Nucl. Phys. B, 38, 394 [CrossRef] [Google Scholar]
 Caldwell, R. R., Dave, R., & Steinhardt, P. J. 1998, Phys. Rev. Lett., 80, 1582 [NASA ADS] [CrossRef] [Google Scholar]
 Carr, B. J. 1994, ARA&A, 32, 531 [NASA ADS] [CrossRef] [Google Scholar]
 Castorina, E., Urbano, F., Massimiliano, L., et al. 2012, Phys. Rev. D, 86, 023517 [NASA ADS] [CrossRef] [Google Scholar]
 Coble, K., Dodelson, S., & Frieman, J. A. 1997, Phys. Rev. D, 55, 1851 [NASA ADS] [CrossRef] [Google Scholar]
 Dekel, A., & Silk, J. 1986, ApJ, 303, 39 [NASA ADS] [CrossRef] [Google Scholar]
 de Vaucouleurs, G., & Freeman, K. C. 1972, Vistas Astron., 14, 163 [NASA ADS] [CrossRef] [Google Scholar]
 Diemand, J., Kuhlen, M., & Madau, P. 2007, ApJ, 667, 859 [NASA ADS] [CrossRef] [Google Scholar]
 Dolgov, A. D. 2002a, Phys. Rep., 370, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Dolgov, A. D., Hansen, S. H., Pastor, S., et al. 2002b, Nucl. Phys. B, 632, 363 [NASA ADS] [CrossRef] [Google Scholar]
 Drimmel, R., & Spergel, D. 2001, ApJ, 556, 181 [NASA ADS] [CrossRef] [Google Scholar]
 Efstathiou, G. 1992, MNRAS, 256, 43 [NASA ADS] [CrossRef] [Google Scholar]
 Efstathiou, G., Bridle, S. L., Lasenby, A. N., et al. 1999, MNRAS, 303, L47 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrara, A., & Tolstoy, E. 2000, MNRAS, 313, 291 [NASA ADS] [CrossRef] [Google Scholar]
 Freese, K., Kolb, E., & Turner, S. 1983, Phys. Rev. D, 27, 1689 [NASA ADS] [CrossRef] [Google Scholar]
 Glendenning, N. K. 1997, Compact Stars. Nuclear Physics, Particle Physics, and General Relativity (New York: SpringerVerlag) [Google Scholar]
 Governato, F., Brook, C., Mayer, L., et al. 2010, Nature, 463, 203 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Guo, Q., White, S., Li, C., & BoylanKolchin, M. 2010, MNRAS, 404, 1111 [NASA ADS] [Google Scholar]
 Hamann, J., Lesgourgues, J., & Mangano, G. 2008, JCAP, 0803, 004 [NASA ADS] [CrossRef] [Google Scholar]
 Hoeft, M., Yepes, G., Gottlöber, S., & Springel, V. 2006, MNRAS, 371, 401 [NASA ADS] [CrossRef] [Google Scholar]
 Honma, M., & Sofue, Y. 1996, PASJ, 45, L103 [NASA ADS] [CrossRef] [Google Scholar]
 Jeans, J. H. 1919, Problems of Cosmology and Stellar Dynamics (Cambridge) [Google Scholar]
 Kang, H., & Steigman, G. 1992, Nucl. Phys. B, 372, 494 [NASA ADS] [CrossRef] [Google Scholar]
 Kent, S. M. 1987, AJ, 93, 816 [NASA ADS] [CrossRef] [Google Scholar]
 Kent, S. M., Dame, T. M., & Fazio, G. 1991, ApJ, 378, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Kerins, E. J., & Carr, B. J. 1994, MNRAS, 266, 775 [NASA ADS] [Google Scholar]
 Kim, S., StaveleySmith, L., Dopita, M. A., et al. 1998, ApJ, 503, 674 [NASA ADS] [CrossRef] [Google Scholar]
 Kippenhahn, R., & Weigert, A. 1994, Stellar Structure and Evolution (New York: SpringerVerlag) [Google Scholar]
 Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, 522, 82 [NASA ADS] [CrossRef] [Google Scholar]
 Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 [NASA ADS] [CrossRef] [Google Scholar]
 Koposov, S., Belokurov, V., Evans, N. W., et al. 2008, ApJ, 686, 279 [NASA ADS] [CrossRef] [Google Scholar]
 Landau, L. D., & Lifshitz, E. M. 2002, Statistical Physics, 5 Course of Theoretical Physics (Oxford: ButterworthHeinemann) [Google Scholar]
 Larson, R. B. 1974, MNRAS, 169, 229 [NASA ADS] [CrossRef] [Google Scholar]
 Lattanzi, M., Ruffini, R., & Vereshchagin, G. V. 2005, Phys. Rev. D, 72, 063003 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, J., & Koh, I. 1996, Phys. Rev. D, 53, 2236 [NASA ADS] [CrossRef] [Google Scholar]
 Lesgourgues, J., & Liddle, A. R. 2001, MNRAS, 327, 13072 [CrossRef] [Google Scholar]
 Lesgourgues, J., & Pastor, S. 1999, Phys. Rev. D, 60, 103521 [NASA ADS] [CrossRef] [Google Scholar]
 Lesgourgues, J., & Pastor, S. 2006, Phys. Rept., 429, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Li, C., & White, S. D. M. 2009, MNRAS, 398, 2177 [NASA ADS] [CrossRef] [Google Scholar]
 Lobashev, V. M., Aseev, V. N., Beselev, A. I., et al. 1999, Phys. Lett. B, 460, 227 [NASA ADS] [CrossRef] [Google Scholar]
 Martinez, G. D., Minor, Q. E., Bullock, J., et al. 2011, ApJ, 738, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Matos, T., & UreñaLópez, L. A. 2000, Class. Quant. Grav., 17, L75 [NASA ADS] [CrossRef] [Google Scholar]
 Matos, T., & UreñaLópez, L. A. 2001, Phys. Rev. D, 63, 063506 [NASA ADS] [CrossRef] [Google Scholar]
 Membrado, M., & Aguerri, J. A. L. 1996, Int. J. Mod. Phys., 5, 257 [NASA ADS] [Google Scholar]
 Membrado, M., & Pacheco, A. F. 2012, Europhys. Lett., 100, 39004 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Membrado, M., Pacheco, A. F., & Sañudo, J. 1989a, Phys. Rev. A, 39, 4207 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Membrado, M., Abad, J., Pacheco, A. F., & Sañudo, J. 1989b, Phys. Rev. D, 40, 2736 [NASA ADS] [CrossRef] [Google Scholar]
 Moore, B., Ghigna, S., Governato, F., et al. 1999, ApJ, 524, L19 [NASA ADS] [CrossRef] [Google Scholar]
 Moscoso, I. 1993, Nucl. Phys. B, 32, 365 [Google Scholar]
 Orito, M., Kajino, T., Mathews, G. J., et al. 2002, Phys. Rev. D, 65, 123504 [NASA ADS] [CrossRef] [Google Scholar]
 Pacheco, A. F., Sañudo, J., & Seguí, A. 1986, Phys. Lett B, 168, 223 [NASA ADS] [CrossRef] [Google Scholar]
 Padmanabhan, T. 1995, Structure Formation in the Universe (Cambridge: Cambridge University Press) [Google Scholar]
 Particle Data Group 2000, EPJ C, 15, 1 [Google Scholar]
 Passalacqua, L. 1997, Nucl. Phys. B, 55, 435 [CrossRef] [Google Scholar]
 Pastor, S., & Lesgourgues, J. 2000, Nucl. Phys. B, 81, 47 [CrossRef] [Google Scholar]
 Peacock, S. 1999, Cosmological Physics (Cambridge: Cambridge University Press) [Google Scholar]
 Pelupessy, F. I., van der Werf, P. P., & Icke, V. 2004, A&A, 422, 55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Perlmutter, S., Aldering, G., Deustua, S., et al. 1997, 191st AAS Meeting #85.04, BAAS, 29, 1351 [NASA ADS] [Google Scholar]
 Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Popa, L. A., & Vasile, A. 2008, JCAP, 0806, 028 [NASA ADS] [CrossRef] [Google Scholar]
 Raychaudhuri, A. K., Banerji, S., & Banerjee, A. 1992, General Relativity, Astrophysics and Cosmology (New York: SpringerVerlag) [Google Scholar]
 Sahni, V., & Wang, L. 2000, Phys. Rev. D, 62, 103517 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Sawala, T., Scannapieco, C., Maio, U., & White, S. 2010, MNRAS, 402, 1599 [NASA ADS] [CrossRef] [Google Scholar]
 Sawala, T., Guo, Q., Scannapieco, C., Jekins, A., & White, S. 2011, MNRAS, 413, 659 [NASA ADS] [CrossRef] [Google Scholar]
 Schmidt, B., Sunzeff, N. B., Phillips, M. M., et al. 1998, ApJ, 507, 46 [NASA ADS] [CrossRef] [Google Scholar]
 Serpico, P. D., & Raffelt, G. G. 2005, Phys. Rev. D, 71, 127301 [NASA ADS] [CrossRef] [Google Scholar]
 Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs, and Neutron Stars (New York: John Wiley & Sons) [Google Scholar]
 Shiraishi, M., Ichikawa, K., Ichiki, K., Sugiyama, N., & Yamaguchi, M. 2009, JCAP, 0907, 005 [NASA ADS] [CrossRef] [Google Scholar]
 Simon, J. D., & Geha, M. 2007, ApJ, 670, 313 [NASA ADS] [CrossRef] [Google Scholar]
 Simon, J. D., Geha, M., Minor, Q. E., et al. 2011, ApJ, 733, 46 [NASA ADS] [CrossRef] [Google Scholar]
 Somerville, R. S. 2002, ApJ, 572, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. N., Verde, L., Peris, H., et al. 2003, ApJS, 148, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Springel, V., Wang, J., Vogelsberger, M., et al. 2008, MNRAS, 391, 1685 [NASA ADS] [CrossRef] [Google Scholar]
 Steinhardt, P. J. 1996, Nature, 382, 768 [NASA ADS] [CrossRef] [Google Scholar]
 Stinson, G. S., Dalcanton, J. J., Quinn, T., Kaufmann, T., & Wadsley, J. 2007, ApJ, 667, 170 [NASA ADS] [CrossRef] [Google Scholar]
 Strukov, I. A., Skulachev, D. P., & Klypin, A. A. 1987, Large Scale of the Universe, IAU Symp. 130, eds. J. Audouze, M.C. Pelletar, & H. Szalay (Dordrecht: Kluwer), 93 [Google Scholar]
 Tollerud, E. J., Bullock, J. S., Strigari, L. E., & Willman, B. 2008, ApJ, 688, 277 [NASA ADS] [CrossRef] [Google Scholar]
 Tremaine, S., & Gunn, J. E. 1979, Phys. Rev. Lett., 42, 407 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Trimble, V. 1987, ARA&A, 25, 425 [NASA ADS] [CrossRef] [Google Scholar]
 Turner, M. S., & White, M. 1997, Phys. Rev. D, 56, 4439 [NASA ADS] [CrossRef] [Google Scholar]
 van der Marel, R. P., Alves, D. R., Hardy, E., & Suntzeff, N. B. 2002, ApJ, 124, 2639 [Google Scholar]
 Wasserman, I., & Salpeter, E. E. 1994, ApJ, 433, 670 [NASA ADS] [CrossRef] [Google Scholar]
 Weinberg, S. 1972, Gravitation and Cosmology (New York: Wiley) [Google Scholar]
 Weinheimer, Ch., Degen, B., Bleile, A., et al. 1999, Phys. Lett. B, 460, 219 [NASA ADS] [CrossRef] [Google Scholar]
 Zaritsky, D. 1999, AJ, 118, 2824 [NASA ADS] [CrossRef] [Google Scholar]
 Zavala, J., Jing, Y. P., Faltenbacher, A., et al. 2009, ApJ, 700, 1779 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Mass limits, M_{M} (upper continuous line) and M_{m} (lower continuous lines), versus fermion mass, m, for dark clusters of degenerate fermions. is assumed. For details see the text. 

Open with DEXTER  
In the text 
Fig. 2 Cluster radii, R, versus fermion mass, m, for several cluster masses, M. is assumed. Horizontal lines correspond to 1, 10, 30 kpc and 1 Mpc. For details see the text. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.