Abstract

When the Large Hadron Collider resumes operations in 2021, several experiments will directly measure the motion of antihydrogen in free fall for the first time. Our current understanding of the universe is not yet fully prepared for the possibility that antimatter has negative gravitational mass. This paper proposes a model of cosmology, where the state of high energy density of the big bang is created by the collapse of an antineutrino star that has exceeded its Chandrasekhar limit. To allow the first neutrino stars and antineutrino stars to form naturally from an initial quantum vacuum state, it helps to assume that antimatter has negative gravitational mass. This assumption may also be helpful to identify dark energy. The degenerate remnant of an antineutrino star can today have an average mass density that is similar to the dark energy density of the ΛCDM model. When in hydrostatic equilibrium, this antineutrino star remnant can emit isothermal cosmic microwave background radiation and accelerate matter radially. This model and the ΛCDM model are in similar quantitative agreement with supernova distance measurements. Therefore, this model is useful as a purely academic exercise and as preparation for possible future discoveries.

1. Introduction

Type-Ia supernovae (SNe Ia) occur when an accreting white dwarf star exceeds the Chandrasekhar limit and collapses until the released potential energy detonates the star with carbon fusion. Due to their similar initial conditions, SNe Ia can serve as standard candles. Twenty years ago, this allowed observations of SNe Ia to show an accelerating expansion of matter on cosmological scales [13]. The big bang model successfully describes this expansion with an expansion of the underlying metric. It assumes that the universe is homogeneous and isotropic on large scales, which is known as the cosmological principle. Mathematically, the Friedman–Lemaître–Robertson–Walker (FLRW) metric upholds the cosmological principle by uniformly changing the metric of space with a scale factor that varies in time. Today, a concordance from various observations defines the ΛCDM model, which divides the energy content of the universe into 31% matter and 69% dark energy, where dark energy is commonly thought to be the constant energy density of the quantum vacuum [4, 5]. This translates to a dark energy density of ρΛ≈6 × 10−30 g/cm3, which is 10120 times smaller than expected [6]. Therefore the identity of dark energy is a major unsolved puzzle. Alternative models posit a negative gravitational mass [710], negative inertial mass [11] or relic neutrino condensation [12, 13]. Despite these alternatives, the ΛCDM model is so far our most successful description of the late universe.

While the ΛCDM model is the most widely accepted model of the late universe, it does not address several key puzzles associated with the early universe. First, when we trace the expansion of all observable matter (∼1055 g) backwards in time, we encounter a collective state of high energy density that is associated with the big bang. Current models treat the big bang either as the beginning of space and time [14], or as part of a continuous bounce [15, 16]. While it serves as initial condition for the ΛCDM model, the origin of the big bang is a major unsolved puzzle. A second problem arises when we assume that this big bang state initially contained equal amounts of matter and antimatter, while unequal amounts are observed in the universe today. To account for this matter-antimatter asymmetry, it is necessary to find a mechanism for baryogenesis [17, 18]. A third puzzle is associated with the cosmic microwave background radiation (CMB). The ΛCDM model assumes that the CMB was emitted by initially hot and dense matter that has cooled sufficiently to become transparent to radiation ∼105 years after the big bang [19]. However, the CMB is more isotropic than expected, which is known as the horizon problem. The theory of cosmological inflation addresses this by introducing a period of exponentially accelerating expansion up to 10−32 s after the big bang [2022]. This could allow any two regions of the CMB in the early universe to have been thermalized. This also addresses the question of why our expanding metric appears to be spatially flat, known as the flatness problem. However, inflation suffers from problems such as the entropy problem or the multiverse problem [23, 24]. Moreover, the standard model interpretation of the CMB appears to be in tension with recent cosmology-independent measurements of the expansion rate at both low redshifts (z < 0.15) and high redshifts (1.4 < z < 5.1) at the ∼4σ level [2527]. Lastly, a dipolar anisotropy in the expansion rate of nearby galaxies (z∼0.1) violates the assumptions of homogeneity and isotropy of the ΛCDM model at the 4σ level [28]. These puzzles motivate a search for new models, with the ΛCDM model as the benchmark.

In the present work, a model of cosmology is proposed that attempts to address the origin of the big bang and dark energy. A degenerate self-gravitating gas of antineutrinos, which we will call an antineutrino star, collapses when its mass exceeds the Chandrasekhar limit, [29, 30]. The small neutrino mass motivates a model of cosmology, where the collapse of an antineutrino star creates the state of high energy density of the big bang with a minimum of new physics. As initial condition of the proposed model, we choose a quantum vacuum state due to its minimal entropy. This state is gravitationally unstable and organically forms spatially separated neutrino stars and antineutrino stars when we assume that antimatter has negative gravitational mass. After collapse in an energetic event as described below, a fraction of the antineutrino gas eventually returns to effective hydrostatic equilibrium. If viewed from the core, a degenerate antineutrino star remnant could today emit isothermal cosmic microwave background radiation and radially accelerate matter.

The main finding is that both the new model and the ΛCDM model describe redshift-distance measurements with comparable quantitative accuracy. The best-fit parameters of the new model give a density of the antineutrino star that is similar to the dark energy density of the ΛCDM model, and constrain the electron neutrino mass to high statistical precision. The new model is qualitatively consistent with CMB anisotropies [31] and large-scale structures [32, 33], which presently challenge the ΛCDM model’s assumptions of homogeneity and isotropy. These results encourage future work to further develop and test the presented model.

2. Antineutrino Star Model of the Early Universe

Similar to a white dwarf star, a degenerate gas of antineutrinos collapses when its mass exceeds the Chandrasekhar limit. In this section we will show that this collapse can transform energy at the scale of the mass-energy content of the known universe. In order to create the conditions for this event to occur, it becomes helpful to assume that antimatter has negative gravitational mass. This assumption will be tested at CERN and, if true, may also help identify dark energy.

2.1. Chandrasekhar’s Equation of State

To characterize a degenerate antineutrino gas in effective hydrostatic equilibrium, we first make three assumptions. First, we ignore thermal or radiation pressure of the antineutrino gas by assuming it is highly degenerate with temperature T/TF ≪ 1, where TF is the Fermi temperature. Second, we assume that the neutrino is a Dirac fermion, which means that it is not its own antiparticle. This standard model assumption is being tested by the search for neutrinoless double-beta decay [34, 35]. Third, we assume that all neutrinos in the star are electron flavored neutrinos with effective inertial mass . Note that due to neutrino oscillations [36, 37], free electron neutrinos have an effective mass , where Uei are the Pontecorvo–Maki–Nakagawa–Sakata leptonic mixing matrix elements and mi are eigenstates of definite mass (i = 1, 2, 3, respectively). We will show below that this third assumption is reasonable in the context of the Schwinger mechanism (see equation (16)) if the electron neutrino mass is much less than the muon or tau neutrino masses.

With these assumptions we can use Chandrasekhar’s equation of state for degenerate matter, derived from hydrostatic equilibrium of gravitational and degeneracy pressures [29, 30, 38]. We thus get the equations of density,and pressure,where x is dimensionless and proportional to the Fermi momentum pF,and K has dimensions of energy density,and the other constants take their usual meaning. We want to solve the following equation for hydrostatic equilibrium:where m(r) is the gravitational mass enclosed inside a radius r and is given by

The spherically symmetric gravitational potential φ(r) is given by

For regions outside an antineutrino star (r > R) the potential simplifies towhere M < 0 is the gravitational mass of the antineutrino star. The radius R of the star can be found where P(R) = 0, while the total mass of the star is given by M = m(R),

Chandrasekhar’s equation of state above does not take into account general relativistic effects, which we can ignore for simplicity if the radius of the star is much larger than its Schwarzschild radius, RS/R ≪ 1. Note that we assume that inertial masses (mi = |m|) of matter and antimatter are equal and positive, and that the gravitational mass ( = m) of matter is positive ( = 1) and of antimatter is negative ( = −1); for generality, we do not use overbar notation. Note that the gauge is fixed to zero at the star’s center, φ(0) = 0, so that φ(r) < 0 everywhere else.

The above equations can be solved numerically for a given central density ρ0 and effective neutrino mass . To facilitate this, we recast equations (5)–(7) above using dimensionless units for radial position η = r/a, mass μ = m/b, and potential ϕ = φ/c2, where

This gives for hydrostatic equilibriumand for the mass enclosedwhile the potential becomes

Using the boundary conditions μ(0) = 0, central density x(0) = x0, and fixing the gauge ϕ(0) = 0 allows us to solve the above three equations. For illustrative purposes, the mass-radius relationship for an antineutrino star is plotted in Figure 1. The constants a and b depend on the neutrino mass as follows:These large-scale factors suggest that neutrino stars and antineutrino stars may be relevant for cosmology [39].

2.2. Chandrasekhar Limit

The estimated mass-energy content of the universe at the time of the big bang is ∼1055 g [5]. In the early universe, this mass-energy was concentrated in a state of high energy density. Possible origins of the big bang have been discussed by several authors [1416]. One possibility, which seems to have been overlooked in earlier discussions, is that the big bang was created by the gravitational collapse of an antineutrino star. When a white dwarf star exceeds the Chandrasekhar limit, its own gravitational pressure overwhelms the degeneracy pressure of electrons and it collapses in a supernova [29]. Similarly, the limiting mass of a degenerate gas of electron antineutrinos occurs at the dimensionless values of η = 1.06 and μ = 0.356, which corresponds respectively to

The electron neutrino mass was recently constrained to  eV/c2 by the KATRIN experiment [40, 41]. This gives a lower limit of  g, which is only four orders of magnitude below the estimated mass-energy content of the universe at the time of the big bang (∼1055 g). Moreover, an effective electron neutrino mass of  meV/c2 would correspond to a Chandrasekhar limit of  g. Due to these cosmological length and mass scales, the study of neutrino stars and antineutrino stars falls within the realm of cosmology. A small neutrino mass motivates the ansatz of this paper that the collapse of an antineutrino star created the state of high energy density of the big bang with a minimum of new physics.

2.3. Quantum Vacuum Instability

The hypothesis that the collapse of an antineutrino star created the state of high energy density of the big bang immediately raises the question on the origin of the first neutrino stars and antineutrino stars. We attempt to address this qualitatively in the following. Since the nature of the quantum vacuum is still poorly understood, we warn the reader that this subsection is more speculative than the others. Before the first neutrino and antineutrino stars are formed, it is reasonable to assume that the initial condition of the universe was an infinite volume of quantum vacuum due to its low number of degrees of freedom. The quantum vacuum contains a sea of virtual particle-antiparticle pairs going into and out of existence. The universe today is no longer in this low-entropy quantum vacuum state. To explain this, we assume that matter and antimatter gravitationally repel, which causes the quantum vacuum to be gravitationally unstable by the following mechanism. Short-lived perturbations in the particle-antiparticle density create a weak and fluctuating gravitational field on small scales (this field is established by the particles themselves). By the Schwinger mechanism, this field has a nonzero probability of creating real particles by separating virtual ones before they can annihilate [9, 42]. For example, the pair creation rate per unit volume and time in a constant local gravitational field gradient, , is

According to the above equation, even the weak and fluctuating gravitational field of particles in the quantum vacuum can create real particles with nonzero probability. The exponential dependence on the effective mass (m) strongly favors creation of neutrino-antineutrino pairs compared to more massive particles of the Standard Model [9, 43]. It also favors cold neutrinos over hot neutrinos, which allows them to bind gravitationally. Thus, the assumption of negative gravitational mass of antimatter enables the gradual formation of mutually repulsive neutrino stars and antineutrino stars.

2.4. Qualitative Description of the Collapse of an Antineutrino Star

Whenever an accreting antineutrino star exceeds its mass limit, it collapses in a “neutrinonova.” At sufficiently high temperatures and densities, antineutrinos could transform most of their kinetic energy via high-energy collisions into equal quantities of baryonic matter and antimatter. This is qualitatively similar to the conversion of large amounts of kinetic energy into equal quantities of matter and antimatter in experiments such as the Large Hadron Collider at CERN. The collapse could subsequently be reversed in a big bounce. While the details of the bounce mechanism are unknown at this point, they would likely involve nuclear fusion and other physical processes similar to the bounce event of supernovae.

During the subsequent expansion of material, baryonic matter and antimatter start to annihilate faster than they are created. We assume that a small remnant of baryonic matter survives due to baryogenesis [17, 18]. While the details of baryogenesis remain an active area of research, negative gravitational mass has been proposed as a possible origin of CP violation [44], which is one of the requirements for baryogenesis. Similar to a supernova bounce, a fraction of the original antineutrino star is accelerated to escape velocities from the shock wave associated with the bounce. The remaining antineutrinos and matter expand as two adiabatic ideal gases initially in thermal equilibrium. During adiabatic expansion, the temperature (T) decreases with an increase in volume (V) as T ∝ V−1/3 for a relativistic gas and T ∝ V−2/3 for a nonrelativistic gas. Thus, baryonic matter undergoes nucleosynthesis until density and temperature decrease sufficiently to “freeze out” certain reactions. This process could produce light elements comparable to big bang nucleosynthesis (BBN) [45]. Note that a neutrinonova is not associated with an expansion of the underlying metric, but an expansion of material (as in a supernova). Since baryonic particles are much more massive than antineutrinos, they become nonrelativistic at a much higher temperature than antineutrinos, and thus begin structure formation in a much smaller volume. The antineutrinos that have not been accelerated to escape velocities eventually form a degenerate self-gravitating gas with sub-Chandrasekhar mass , and reestablish thermal equilibrium a sufficient time afterwards (t ≫ R/c). Newly formed galaxies are subsequently radially accelerated from initial proximity to the center of the antineutrino star remnant. Observers in the rest frame and inside of this antineutrino star in hydrostatic equilibrium would detect isotropic black body radiation, which we identify as the CMB. We call this the ATLAS (AnTineutrino Lepton gAS) model (in Greek mythology Atlas is a titan who is able to lift up the celestial spheres), and refer to the antineutrino star as ATLAS-1 (see Figure 2 for a summary).

To summarize this section, the collapse of an antineutrino star can explain the energetic event commonly known as the big bang. It also suggests that an antineutrino star remnant in hydrostatic equilibrium may exist in the universe today. This degenerate antineutrino star is detectable by isothermal background radiation and the large-scale motion of galaxies. In the following section, we quantitatively show that this antineutrino star remnant could account for the large-scale motion of galaxies with similar precision as the ΛCDM model.

3. Antineutrino Star Model of the Late Universe

The ΛCDM model can mathematically describe the large-scale motion of galaxies with only two parameters, the fractional matter density (Ωm) and dark energy density (ΩΛ). As we will show below, an antineutrino star model can competitively describe the large-scale motion of galaxies, while also accounting for the physical origin of the big bang and dark energy. The equation of state of an antineutrino star is defined by two unknown parameters, namely, the central density (ρ0) and the effective electron neutrino mass . An initial average expansion velocity of baryonic matter () as a third parameter improves the physical basis of our model at low redshifts (z < 0.04) and could qualitatively explain the low redshift contribution to the 4.4σ Hubble tension [25, 26].

3.1. Definition of Our Observational Frame

An antineutrino star establishes a background potential that is the dominant contribution to the apparent velocities of galaxies. Therefore, we can infer our approximate position empirically from our velocities relative to other galaxies and relative to the rest frame of the antineutrino star, which is the rest frame of the CMB. The velocities of galaxies relative to us scale approximately isotropically with distance at a rate of H0 ≈ 70 km·s−1Mpc−1 [25, 26]. The velocity of the Local Group relative to the CMB is comparably small at  km·s−1 [46]. These velocities suggest we are approximately at rest and near (but not at) the center of the star. This can be physically explained by the small potential gradient near the center, which causes matter initially close to the center to remain effectively at rest, while matter far from the center accelerates down the gravitational potential hill. This empirical inference of our observational frame guides our derivation of the distance-redshift relationship below.

3.2. Derivation of a Distance-Redshift Relationship

We will use the Schwarzschild metric to determine the redshift of light emitted by a galaxy in free fall from initial proximity to the center of an antineutrino star. We may identify ourselves as Schwarzschild observers, who are by definition at rest where the gauge is fixed to zero, namely at the star’s center. When allowing for negative gravitational mass, we need to define two different metrics for matter and for photons.

By symmetry, we assume that photons in a gravitational potential undergo blueshift or redshift independently of the matter or antimatter nature of the gravitational source. We will thus use a Schwarzschild metric for photons that is agnostic to the type of matter,where M/|M| is the sign of the gravitational mass of the source, t is coordinate time, dΩ2 = sin θ2dϕ2 + dθ2 is the angular path element in spherical coordinates, and φ(r) is the gravitational potential of the antineutrino star given by equations (7) and (8). Note that the above equation is only an approximation of the metric inside the star (r < R), since we used Chandrasekhar’s Newtonian model for simplicity to find φ(r). It is worth repeating here that this approximation is accurate when changes in gravitational potential inside the star satisfy 2|φ(r)|/c2 ≪ 1. Outside of the star (r > R), the above Schwarzschild metric accurately captures all further changes in gravitational potential relative to the central observer up to 2|φ(r)|/c2 ≲ 1.

Since proper time is zero for photons, the velocity of distant photons moving radially (dΩγ = 0) isTo summarize the above result, photons are assumed to behave in a way that is agnostic to the matter or antimatter nature of the gravitational source.

We assume that matter experiences a Schwarzschild metric that has been minimally extended to capture the interaction of negative and positive gravitational masseswhere m/|m| is the sign of the gravitational mass of a test particle. Here we have introduced a sign-dependence that transforms the metric’s dependence on gravitational potential (φ) into a dependence on potential energy, V(r) = (r). Note that the above expression allows matter and antimatter to experience two different space-time metrics that are equal where the gauge is fixed to zero.

We can now use energy conservation [47] to find the time dilation factor for galaxies with average initial velocity undergoing free fall from the center of the star,

The coordinate velocity  = dr/dt is found similarly,

The above equations allow us to calculate the redshift seen by a Schwarzschild observer [47],

Therefore, redshift is caused by a combination of apparent radial velocity and time dilation of free-falling sources in the gravitational potential φ. Specifically, contributions of to redshift dominate at low-z and contributions from φ(r) via time dilation dominate at high-z (see Figure 3). The distance modulus of distant SNe Ia in receding galaxies iswhere r(z) is the distance from the center of the antineutrino star in Megaparsec [1].

3.3. Cosmological Parameters

We can compare the theoretical distance modulus to observed distance moduli for a given set of cosmological parameters, . Goodness of fit is determined with a χ2 statistic,where is the measurement uncertainty in the observed distance modulus μ0,i. As observational data, the “Union2.1” catalog (http://supernova.lbl.gov/union/; date accessed: 08/09/2017) of 580 SNe Ia is used [48]. The best-fit parameter values are  meV/c2, ρ0 = 1.60 × 10−29 g/cm3, and  = 6.18 × 10−3c. These give , where dof stands for a model’s degrees of freedom. For reference, the ΛCDM model has a comparable fit of [5]. In Figure 4 the theoretical distance modulus of the ATLAS model with the above best-fit parameters is plotted together with the distance-redshift data in a Hubble diagram.

For comparison, the theoretical distance modulus of the concordance ΛCDM model [5] is shown with H0 = 67.74 km·s−1Mpc−1. Note that the central density of the antineutrino star ρ0 is comparable to the dark energy density of the ΛCDM model. Therefore the ATLAS model can account for the overall expansion of matter, with the density profile of an antineutrino star in equilibrium effectively acting as a dark energy density (see Figure 4).

Following Riess et al. [1], the probability density function (PDF) for a given cosmological parameter is quantified with Bayes’ theorem, which gives a PDF for the electron neutrino mass ofwhere μ0 represents all measured distance moduli, which are assumed to be independent and normally distributed. This gives a mass to one standard deviation of  meV/c2 (see Figure 5).

Note that this uncertainty is purely due to statistical error and does not include possible systematic errors, which can be caused by model assumptions or instrument calibration [48]. This mass is consistent with the present experimental upper limit of  eV/c2 [41], and may be tested by future or ongoing experiments. Measurements of neutrino mixing angles can give the muon and tau neutrino masses, for example, the most recent (NuFIT 3.2 (2018), http://www.nu-fit.org/?q=node/166; date accessed: 06/07/2018) best-fit values give  meV/c2 [49].

4. Discussion

The antineutrino star parameters are approximately consistent with the assumptions of the ATLAS model. The best-fit muon and tau neutrino masses are larger than the electron neutrino mass by a factor of >3. This is consistent with the simplifying assumption that the antineutrino star consists only of electron neutrinos due to the Schwinger mechanism (see section IIA). The best-fit electron neutrino mass gives a Chandrasekhar limit ofwith a limiting radius ofwhich is consistent with the model of the early universe (see Section 2.2). The re-formed antineutrino star has mass and radius RS/R = 0.410 or

Note that Chandrasekhar’s equation of state for the above best-fit parameters naturally gives rise to these cosmological scales. Due to its large size, a degenerate antineutrino star can be detected by the large-scale motion of galaxies and an isothermal microwave background [19].

The best-fit value for the initial velocity,  = 6.18 ± 0.55 × 10−3, can be identified as the root-mean-square speed of thermalized protons, since protons are likely to be the most abundant matter particle by mass in the early universe. This would correspond to a proton temperature of T0 = 11.9 ± 2.1 keV, which is consistent with temperatures at which nucleosynthesis ends due to a decrease in fusion reaction rate with a decrease in density [50, 51]. Moreover, this initial average expansion velocity is responsible for the difference between the ΛCDM model and the ATLAS model at low redshifts (z < 0.04). This difference corresponds to a faster average expansion rate at low redshifts than predicted by the ΛCDM model and could qualitatively explain the 4.4σ Hubble tension reported by Riess et al. [25, 26] (see Figure 4).

The best-fit parameters are also consistent with an antineutrino star that is degenerate, giving a central Fermi temperature TF,0 = 34.1 K > TCMB,0 = 2.73 K [52]. However, the above values for RS/R and TCMB/TF suggest that the systematic error due to model assumptions is ∼20% and dominates over the statistical error. Future work could account for general relativistic effects and the temperature ratio T/TF to improve the accuracy of the antineutrino star model. Since these effects are small and dominate at high-z, where supernova data is still sparse, they are ignored in the present model for simplicity.

The ATLAS model relies on testable assumptions, in particular on the assumption of negative gravitational mass. While there are strong theoretical arguments against negative gravitational mass (for a review, see [53]), there have not yet been any conclusive direct experimental tests. The assumption of negative gravitational mass is in principle testable on cosmological scales, galactic scales, and laboratory scales. On cosmological scales, this assumption is consistent with observations as described above: (i) it allows for the formation of the first neutrino and antineutrino stars from a quantum vacuum state, and (ii) it accounts for the overall expansion of matter in the late universe. On galactic scales, it predicts detectable antimatter emission from hot accreting matter close to the event horizon of black holes [54, 55]. For example, there exists a positron excess in the galactic center [56] and in cosmic rays [5760], which could be emitted by accreting compact objects. On laboratory scales, several experiments are currently undergoing upgrades at the GBAR, ALPHA-g, and AEGIS laboratories at CERN [6163]. Using antiprotons from Large Hadron Collider (LHC) operations resuming in 2021, these experiments will produce antihydrogen and directly observe its motion in free fall. Therefore, laboratory experiments are soon able to test a key assumption of the ATLAS model. In the event that experiments at CERN discover negative gravitational mass, the ATLAS model provides a helpful alternative to the ΛCDM model, which is not optimized for this possibility.

The ATLAS model can be empirically distinguished from the ΛCDM model with a variety of other empirical tests, such as observations that challenge the ΛCDM model’s assumption of isotropy and homogeneity on cosmological scales. For instance, observations of the largely isotropic Hubble expansion and small CMB dipole suggest that we are relatively close to, but not at, the exact center of the star. This off-center location could qualitatively explain anisotropies, such as hemispherical asymmetries in the power spectrum, detected at the ∼3σ level in the CMB [31]. These observations challenge the isotropy assumption of the ΛCDM model. In an effort to explain these observations, many CMB anisotropies have been found to become less significant with the introduction of rotation (specifically, by using a Bianchi VIIh model, albeit one that is decoupled from the standard ΛCDM model) [31, 64]. Future work could investigate whether this rotational axis in the CMB corresponds to the rotational axis of an antineutrino star.

The ATLAS model predicts that any hemispherical asymmetries in the CMB [31] and any anomalies due to rotation of the antineutrino star disappear when viewed from rest at the exact center of the star. Therefore, it may be possible to locate this central position with the study of CMB anisotropies alone. For example, the intersection of the axis of rotation found in the CMB with the plane on which hemispherical asymmetries are minimized could help locate the central region of the antineutrino star. Independently, this central region could coincide with an underdensity of matter in our cosmic neighborhood such as the so-called “dipole repeller” at z≈0.05 [65]. The dipole repeller is a void that appears to have contributed to roughly half of our velocity () relative to the CMB [65]. If the center of the antineutrino star were indeed in our cosmic neighborhood (z ≲ 0.05) and contributing significantly to , then one would also expect to see an asymmetry in the expansion rate of nearby galaxies. Specifically, one would expect an asymmetry that is aligned with . This is consistent with observations at the ∼2σ level of a dipole anisotropy in the distribution of nearby (z ≲ 0.1) radio galaxies that is aligned with [66, 67]. This is also consistent with recent observations of a similarly aligned dipole anisotropy in the expansion rate of nearby supernovae (z∼0.1) [28].

Moreover, one would expect the gravitational interaction between galaxy clusters and the antineutrino gas to generate correlations between the distribution of large-scale structures and temperature fluctuations in the CMB. This could qualitatively explain the observed correlations between anisotropies in the CMB and anisotropies in the matter distribution, known as baryon acoustic oscillations [4, 68].

Furthermore, recall that structure formation in the presented model of the early universe begins in a smaller volume for baryons than for antineutrinos due to their large mass difference (see Figure 2). The initial baryonic structures will be relatively weakly bound and on a supercluster mass scale due to their initially high temperature. As they continue their formation, these structures are dispersed by their acceleration in the gravitational potential of the antineutrino star (see Figure 6). This places evolved structures, whose old age is difficult to reconcile with the ΛCDM model, at high redshifts. This is also qualitatively consistent with the existence of large supercluster-size structures in the distant universe that are no longer gravitationally bound. For example, the Hercules–Corona Borealis Great Wall has been detected at the 1 − 3σ level with a size of 7–10 Gly [32, 33]. This challenges the homogeneity assumption of the ΛCDM model, which predicts structures not larger than ∼1.21 Gly [69, 70]. Therefore, the ATLAS model appears qualitatively consistent with observations that challenge the fundamental assumption of homogeneity and isotropy of the ΛCDM model.

Lastly, the ATLAS model proposes a possible origin of the big bang and dark energy, which are two major puzzles of the ΛCDM model. However, this work relies primarily on SNe Ia data to constrain its respective cosmological parameters (, ρ0, ). This work is still at an early stage, and comparable to early work on the ΛCDM model [13]. During the twenty years since the discovery of dark energy, the ΛCDM model interpretation of SNe Ia data has been in general concordance with additional fits to observations such as baryon acoustic oscillations [4, 68]. We do not attempt these fits yet; they are beyond the scope of the present work. This work provides a viable cosmological model for the possibility that experiments at CERN discover negative gravitational mass. The present work encourages future work to further develop and test the ATLAS model while relying on the ΛCDM model as a benchmark.

5. Conclusions

The entropy of the universe increases relative to an initial state of low entropy. A natural candidate for this initial state is the quantum vacuum. With the assumption that antimatter has negative gravitational mass, which will be tested at CERN, this vacuum gradually decays into neutrinos and antineutrinos. The collapse of an antineutrino star is a possible explanation for the energetic event commonly known as the big bang. After collapse and subsequent net creation of matter via baryogenesis, an antineutrino gas adiabatically expands and partially re-forms into an antineutrino star remnant in hydrostatic equilibrium. If viewed from its core, this star can today emit the isothermal CMB radiation and radially accelerate matter. This addresses the nature of dark energy, and removes the horizon and flatness problems without invoking cosmological inflation. The above ATLAS model is in good quantitative agreement with distance-redshift measurements. The electron neutrino mass is constrained as a cosmological parameter to  meV/c2. The model is qualitatively consistent with existing observations of large structures in the distant universe and anisotropies in the CMB. The presented alternative cosmological model is helpful for the possibility that experiments at CERN discover negative gravitational mass. Moreover, the apparent consistency of the presented model with observational data motivates future work to further develop and test this nascent model while relying on the ΛCDM model as a benchmark.

Data Availability

The data used in this work are publicly available at http://www.nu-fit.org/?q=node/166 and http://supernova.lbl.gov/union/.

Conflicts of Interest

The author has no conflicts of interest to declare.

Acknowledgments

The author gratefully thanks Mark Morris, Viktor Linders, Elizabeth Mills, Rainer Sachs, Ariella Machness, Paul Neiser, and Troy Carter for helpful questions and comments on the manuscript. This work used computational resources associated with the Hoffman2 Shared Cluster provided by the UCLA Institute for Digital Research and Education.