Next Article in Journal
Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers
Next Article in Special Issue
The Variation of G and Λ in Cosmology
Previous Article in Journal
Dark Atoms of Nuclear Interacting Dark Matter
Previous Article in Special Issue
Non-Canonical Dark Energy Parameter Evolution in a Canonical Quintessence Cosmology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Primordial Axion Stars and Galaxy Halo Formation

by
Alexander I. Nesterov
Department of Physics, CUCEI, University of Guadalajara, Guadalajara 44430, Jalisco, Mexico
Universe 2024, 10(9), 369; https://doi.org/10.3390/universe10090369
Submission received: 1 August 2024 / Revised: 9 September 2024 / Accepted: 9 September 2024 / Published: 12 September 2024
(This article belongs to the Special Issue Dark Energy and Dark Matter)

Abstract

:
Primordial axion stars, hypothetical stars formed from axions, could play an essential role in forming galaxy halos. These stars could have originated in the early universe shortly after the Big Bang. We show that the ultralight axions forming primordial stars can act as the initial seeds for galaxy halos.

1. Introduction

The axions are one of the favorite cold dark matter (CDM) candidates [1,2,3,4,5,6,7,8,9]. Initially, they were invented to solve the strong CP problem: Why does quantum chromodynamics (QCD) not break the CP -symmetry? The Peccei–Quinn (PQ) theory, solving this problem, leads to the existence of very light axions [10,11,12,13]. Since the axions have a tiny mass and extraordinarily weak coupling to matter and radiation, it would be challenging to detect them. Efforts toward detecting axions include a wide range of technologies [14,15,16,17].
The massless axion arises due to a spontaneous breaking of the U ( 1 ) PQ symmetry. The axion acquires a mass during the QCD phase-transition, following the PQ phase-transition. In terms of the axion decay constant, f a , the axion mass induced by QCD instantons can be written as m a 6 μ eV 10 12 GeV / f a [12,13,18,19,20]. More precise estimation of the axion mass yields [21]
m a 5.7 μ eV 10 12 GeV f a .
Cosmological and astrophysical constraints restrict the value of the decay constant and leave the allowed QCD axion mass window as 6 μ eV m 2000 μ eV [8,22,23]. Recent calculations of the axion mass, based on high-temperature lattice QCD, lead to the more narrow window: 50 μ eV m a 1.5 × 10 3 μ eV [24].
There are two sources of cold light axions: vacuum misalignment and string decay and domain-wall decay [19,25,26,27]. In the first case, inflation occurs after the PQ phase transition. The axion field is homogenized over enormous distances, and the vacuum misalignment yields axions with the zero-momentum mode and zero velocity dispersion. In the second scenario, the axion field will have random values in different causally connected universe volumes when the PQ phase transition occurs after inflation. This results in the appearance of a network of cosmic strings. The strings disappear after the QCD phase transition, radiating massive axions with momentum of order the Hubble expansion rate. The cold LAs, produced by string decay, have a wide distribution of momentums in their vacuum state and an extremely small velocity dispersion of their excitations [28]. The numerical lattice simulations show that the mass of QCD axions produced in the post-inflationary scenario is restricted to 0.5 meV m a 20 meV  [27]. Recent studies refine this result, leading to axion’s mass prediction in the range 40 μ eV m a 180 μ eV [26].
Primordial axion stars formed by QCD axions with mass m a 10 6 eV, also known as axion clumps, are formed in the early universe shortly after the Big Bang when axions condensed into a dense, gravitationally bound Bose–Einstein condensate. These clumps can range in size and mass, potentially forming compact objects with significant gravitational influence [29,30,31,32,33,34,35].
The axions with the mass range 10 33 eV m a 10 18 eV , known as ultralight axions (ULA’s), are considered promising candidates for dark matter [36,37,38,39,40,41]. The mass of ULA’s can be estimated as follows [40,41]:
m a = 5.3 × 10 19 eV 10 16 GeV f a 4 .
Current analyses of the Lyman-alpha forest data have placed the following lower bounds on the mass of scalar DM: m 10 21 eV . Masses below this threshold would result in a universe with insufficient small-scale structure compared to what is observed [42,43,44].
In the modern understanding of galaxy formation, every galaxy forms within a dark matter halo. The formation of galaxy halos is a critical aspect of galaxy formation, providing the gravitational scaffold around which baryonic matter (ordinary matter) accumulates to form stars and galaxies. Primordial axion stars, formed in the early universe, could act as seeds for the formation of larger structures, leading to the formation of galaxies surrounded by halos. Their gravitational influence would attract surrounding dark matter and baryonic matter, facilitating the growth of density perturbations that eventually led to galaxy formation [45].
In this paper, we study the formation of the DM stars by ULAs and the relationship between primordial axion stars and galaxy halos. The essential features of our approach, which distinguish it from others, are the following: (1) We assume that primordial axion stars were formed in the early universe considerably before the galaxies’ formation. (2) The axion density does not vanish at infinity (a usual boundary condition in previous studies). We assume that the star shell’s density equals the average density of the DM.
We show that ULA stars, with a size and mass comparable with observable data for typical galaxies of our universe, are formed in the radiation-dominated era. Our findings reveal that primordial axion stars could act as the initial seeds for galaxy halos. The suggested model predicts the formation of solitonic cores surrounded by extended halos. Another feature is the formation of clusters, with a significant density, that have not merged at the center. It can explain why some dwarf spheroidal galaxies contain globular clusters.
The paper is organized as follows. Section 1 discusses the cosmological production of axions and the formation of axion stars. In Section 2, we introduce a non-relativistic model describing axion stars. Section 3 discusses the obtained results. In Conclusion, we summarize the findings of our model and future perspectives. We use the natural unit convention, = c = 1 , throughout the paper.

2. Cosmological Production of Axions

We assume the universe is described by the Λ CDM model [46,47]. In this model, the universe comprises photons, neutrinos, ordinary (baryon) matter, cold dark matter, and dark energy. The latter is defined by the cosmological constant, Λ , and is responsible for the observable acceleration in the expansion of our universe. The axion DM can be described by a model with the action given by [19,38,40]
S = g d 4 x 1 2 g μ ν ϕ , μ ϕ , ν V ( ϕ ) .
The corresponding equation of motion is the Klein–Gordon equation:
ϕ = V ( ϕ ) ϕ .
Qualitatively, the potential V ( ϕ ) generated by QCD instantons may be written as [11,19]
V ( ϕ ) = ( m a f a ) 2 1 cos ( ϕ / f a ) ,
where m a is the axion mass. The potential V ( ϕ ) is the periodic function with a period defined by a shift-invariance, ϕ ϕ + 2 π f a . There are N degenerate vacua, and the axion field has a range 0 < ϕ < 2 π f a [19].

2.1. Non-Relativistic Limit of the Klein–Gordon Equation

In the non-relativistic limit, the Klein–Gordon equation has the form of the Schrödinger equation with the effective potential describing the interacting forces in the system in the non-relativistic limit [29,30,48,49,50]. By writing
ϕ = 1 2 m a ψ e i m a t + ψ * e i m a t ,
one can show that the non-relativistic theory is described by the effective action for the complex field ψ :
S eff = d 4 x i 2 ψ * ψ ˙ ψ ˙ * ψ 1 2 m a ψ * · ψ V eff ( ψ * ψ ) 1 8 π G Φ · Φ m a ψ * ψ Φ ,
where Φ ( r ) is the Newtonian gravitational potential, and V e f f ( | ψ | ) is the effective potential obtained as the non-relativistic limit of the instanton potential Equation (5). The equations of motion are given by
i ψ t = 1 2 m 2 ψ + m a Φ ψ + V eff ψ * ψ , 2 Φ = 4 π G m a | ψ | 2 .
For the axion field, the effective instanton non-relativistic potential potential was derived in Refs. [48,51] and takes the form
V eff ψ * ψ = 1 2 m a ψ * ψ + m a 2 f a 2 1 J 0 n 1 / 2 ,
where J 0 ( z ) is a Bessel function, n = ψ * ψ / n 0 and n 0 = m a f a 2 / 2 .
Figure 1 represents the normalized effective potential, V 0 ( x ) = 1 J 0 ( x ) , with the first term subtracted:
V 0 = 1 2 m n 0 V eff 1 2 m a ψ * ψ
The stable solutions for the axion stars with n 1 (dense stars) correspond to the minima of the effective potential [29].

2.2. Formation of Axion Stars

Axion stars are characterized by their high density and stability, maintained by a balance between gravitational attraction and quantum pressure due to the axions’ wave-like nature. They can be categorized into two types:
  • Dilute Axion Stars ( n c o r e 1 ). These stars are less dense and more extended, and their size is determined by the axion mass and self-interaction strength.
  • Dense Axion Stars ( n c o r e 1 ). These are more compact and can undergo complex dynamical processes, including potential collapse into black holes and phase transition to the dilute stars under certain conditions.
  • The stability and structure of axion stars are determined by the balance between gravitational attraction and the quantum pressure arising from the Heisenberg uncertainty principle.

3. Model

The non-relativistic axions can form a Bose–Einstein condensate (BEC). The equations that describe a static BEC of axions can be obtained by considering an axion field ψ ( r , t ) whose time dependence is ψ ( r ) exp ( i μ t ) , where μ is the chemical potential [29]. They can be expressed conveniently as coupled equations for ψ ( r ) and the gravitational potential Φ ( r ) :
2 ψ ^ = 2 m a 2 Φ μ m a 1 2 + J 1 ( | ψ ^ | ) | ψ ^ | ψ ^ , 2 Φ = 4 π G m a n 0 | ψ ^ | 2 ,
where ψ ^ = ψ / n 0 , G = 1 / m p 2 , and m p denotes the Planck mass ( m p = 1.2 × 10 19 GeV ). To further simplify matters, we rescale coordinates and functions: r r * = r / r 0 , ψ ψ ( r * ) and Φ γ Φ ( r * ) , where r 0 = γ / m a and γ = m p / f a .
Hereafter, we drop the ∗’s in rescaled functions and coordinates. Then, the dimensionless equations can be recast as follows:
2 ψ ^ = 2 Φ γ ε ψ ^ , 2 Φ = 4 π | ψ ^ | 2 ,
where
ε = μ ¯ + 1 2 J 1 ( | ψ ^ | ) | ψ ^ | ,
and we set μ ¯ = μ / m a .
In Table 1, we present various length and mass scales to facilitate recalculation in physical units. In particular, the mass of the axion star can be obtained using the following formula:
M = 3 M 0 | ψ ^ | 2 d 3 V * ,
where M 0 = 4 π m a n 0 r 0 3 / 3 , and the integration is performed over dimensionless variables.
Our model is based on the previous studies of the rotating axion stars made in [29,31,32,49,52]. In general ψ may consist of a superposition of spherical harmonics Y l m ,
ψ ^ ( r , θ , φ ) = = 0 m = l f m r Y m ( θ , φ ) .
For simplicity, we assume the following ansatz:
ψ ^ ( r , θ , ϕ ) = R ( r ) Y l m ( θ , φ ) .
After substitution of ψ ^ ( r , θ , ϕ ) , one can rewrite Equation (12) as
Δ r R l ( l + 1 ) r 2 R = 2 Φ γ μ ¯ γ 2 + γ J 1 ( R | Y l m | ) R | Y l m | R , 2 Φ = 4 π R 2 | Y l m | 2 ,
where Δ r denotes the radial part of the Laplacian,
Δ r = 1 r 2 r r 2 r .
Note that (17) still depends on angular variables θ and φ . Following the procedure presented in [49], we take the angular average of these equations to obtain the following equations of motion:
Δ r R = 2 Φ ˜ + l ( l + 1 ) 2 r 2 γ μ ¯ γ 2 + γ J l m ) R ,
Δ r Φ ˜ = R 2 ,
where
Φ ˜ = 1 4 π Φ ( r , θ , φ ) d Ω , J l m = 1 4 π J 1 ( R | Y l m | ) R | Y l m | d Ω .
We used the normalization of the spherical harmonics, Y m Y m * d Ω = δ δ m m , to perform averaging. For l = 0 , 1 , we have obtained the analytical expressions for J l m :
J 1 0 = J 0 ( a 0 R ) J 1 ( a 0 R ) a 0 R + π 2 J 1 ( a 0 R ) H 0 ( a 0 R ) J 0 ( a 0 R ) H 1 ( a 0 R ) ,
J 1 ± 1 = 1 cos ( a 1 R ) a 1 2 R 2 , J 0 0 = J 1 ( R ) R ,
where a 0 = 3 / 4 π , a 1 = 3 / 8 π and H n ( z ) denotes the Struve function.
For a given function R ( r ) , the gravitational potential can be obtained by integrating Equation (20). The solution can be written as follows:
Φ ˜ ( r ) = 1 r 0 r x 2 R 2 ( x ) d x r r h x R 2 ( x ) d x ,
where r h is the radius of the star. From here it follows that Φ ˜ ( 0 ) = Φ ˜ 0 and Φ ( 0 ) = 0 , where
Φ ˜ 0 = 0 r h x R 2 ( x ) d x .
We impose the following boundary conditions to solve numerically the coupling system of the differential Equations (19) and (20):
R ( 0 ) = R 0 , R ( 0 ) = R 1 , Φ ˜ ( 0 ) = Φ 0 , Φ ˜ ( 0 ) = 0 .
The radius of the axion star is obtained as the solution of the equation
Φ ˜ ( r h ) = M * ( r h ) r h ,
where M * ( r h ) = 0 r h x 2 R 2 ( x ) d x is the rescaled mass of the star, M * = M / M 0 , and
Φ ˜ ( r h ) = 1 r h 0 r h x 2 R 2 ( x ) d x .
The chemical potential is subject to the following boundary conditions:
μ ¯ = 1 γ Φ ˜ + l ( l + 1 ) 2 r 2 + J l m 1 2 | r = r h = 0
Redshift and universe chronology. The redshift for different stages of the universe’s evolution can be estimated as follows:
  • Energy-dominated era: z = ( 0 1.46 )
  • Matter-dominated era: z = ( 1.46 4.86 × 10 3 )
  • Radiation-dominated era: z = ( 4.86 × 10 3 5.92 × 10 15 )
We use the Λ CDM model [46,47] to estimate the redshift. For the dark energy-dominated era, the redshift is given by
1 + z = a ( t ) a ( t m ) = sinh 2 / 3 3 2 H 0 Ω Λ t sinh 2 / 3 3 2 H 0 Ω Λ t m ,
where Ω Λ = Λ / ( 3 H 0 2 ) , and H 0 is the Hubble constant. This expression is quite accurate for t > t m = 10 Myr ( 10 7 y ) = 3.15 × 10 14 s . For a radiation-dominated era, we have a ( t ) t , and for a matter-dominated universe, a ( t ) t 2 / 3 . Combining all results, we obtain
z t r t i t m t r 2 / 3 sinh 2 / 3 3 2 H 0 Ω Λ t sinh 2 / 3 3 2 H 0 Ω Λ t m .
For an estimation, we take the following values: t m = 9 billion years (end of the matter-dominated era and beginning of the dark energy-dominated era), t r = 47,000 years (end of the radiation-dominated era), and t i = 10 12 s (beginning of the radiation-dominated era).

4. Results and Discussion

Following the procedure outlined in Section 3, we have solved numerically the equations of motion (19) and (20) for the axion mass, m a = 5.3 × 10 19 eV, and the boundary conditions presented below.
  • Dense stars
    DES ( I ) ) : R ( 0 ) = 6 , R ( 0 ) = 0 , Φ ˜ ( 0 ) = 1.5 , Φ ˜ ( 0 ) = 0 , ( l = 0 , m = 0 ) , DES ( II ) ) : R ( 0 ) = 20 , R ( 0 ) = 0 , Φ ˜ ( 0 ) = 3.15 , Φ ˜ ( 0 ) = 0 , ( l = 0 , m = 0 ) .
  • Dilute stars
    DIS ( I ) ) : R ( 0 ) = 0.1 , R ( 0 ) = 0 , Φ ˜ ( 0 ) = 0.055 , Φ ˜ ( 0 ) = 0 , ( l = 0 , m = 0 ) , DIS ( II ) ) : R ( 0 ) = 0 , R ( 0 ) = 10 , Φ ˜ ( 0 ) = 0.001 , Φ ˜ ( 0 ) = 0 , ( l = 1 , m = 0 ) , DIS ( III ) : R ( 0 ) = 0 , R ( 0 ) = 0.12 , Φ ˜ ( 0 ) = 0.5 , Φ ˜ ( 0 ) = 0 , ( l = 1 , m = ± 1 ) .
The boundary of the star is fixed by conditions imposed by Equation (27).
The results of the numerical solution of Equations (19) and (20) are presented in Figure 2, Figure 3 and Figure 4 and summarized in Table 2. As one can see, the axion stars have an explicitly expressed core with the radius, r c o r e r h . Inside the stellar shell, the density oscillates around some non-zero value. The oscillations in graphs reflect the oscillating character of the Bessel functions.
We assume that the primordial axion stars are formed in the relatively late universe compared with the inflation and the Big Bang nucleosynthesis periods (around hours after the Big Bang). At the star’s boundary, the axion particle density is n b = n D M ( z ) , where n D M ( z ) is the average particle density of DM at the redshift z. To obtain the redshift, we use the relation, z = ( n D M ( z ) / n D M ) 1 / 3 , where n D M = ρ DM / m a , and ρ DM = 1.3 × 10 6 GeV / cm 3 is the observed average density of the DM obtained from the data presented in Ref. [53].
The current size of the primordial axion stars can be estimated using the following relation: d = 2 z r h . We obtain
  • DES (I): d = 693.3 kpc ( l = 0 , m = 0 , M = 1.57 × 10 10 M )
  • DES (II): d = 567.8 kpc ( l = 0 , m = 0 , M = 1.99 × 10 9 M )
  • DIS (I): d = 244.3 kpc ( l = 0 , m = 0 , M = 9.56 × 10 8 M )
  • DIS (II): d = 30.6 kpc ( l = 1 , m = 0 , M = 4 × 10 7 M )
  • DIS (III): d = 100 kpc ( l = 1 , m = ± 1 , M = 3.92 × 10 10 M )
These results are in agreement with an estimation of the halo size for typical galaxies:
  • Spiral Galaxies: 61.3 to 92 kpc in diameter
  • Elliptical Galaxies: 92 to 306.6 kpc in diameter
  • Dwarf Galaxies: up to 9.2 kpc in diameter
The mass of galaxies can vary significantly depending on their type and size:
  • Spiral Galaxies: M ( 10 10 10 12 ) M .
  • Elliptical Galaxies: M ( 10 12 10 13 ) M .
  • Dwarf Galaxies: M ( 10 7 10 9 ) M .
Considering the distribution of different types of galaxies, the average mass of galaxies in our universe can be roughly estimated as follows: M 10 11 M . As one can see, the mass of axion stars presented in Table 2 is inside the range of the observable data. Below, we discuss these results in detail.
Elliptical galaxies. Elliptical galaxies are characterized by their ellipsoidal shape. They range from nearly spherical to highly elongated shapes. The stars in elliptical galaxies typically orbit in random directions, which contributes to their more rounded appearance. The best candidates for forming elliptical galaxies are dilute axion stars with spherical symmetry ( l = 0 , m = 0 ). The current size of the presented primordial axion star is an order of the visible halo of the typical elliptical galaxy. These stars are formed at the beginning of the radiation-dominated era ≈150 years after the Big Bang.
Dwarf galaxies. Dwarf galaxies are small galaxies that contain significantly fewer stars compared to giant galaxies like the Milky Way. They are among the most common types of galaxies in the universe and play a crucial role in studying galaxy formation and evolution. Dwarf galaxies are thought to be the building blocks of more giant galaxies. Over time, they can merge with one another or be absorbed by larger galaxies.
Dwarf galaxies exhibit uniform central masses and shallow density profiles that imply the density decreases slowly with distance from the center. The most crucial feature of dwarf galaxies is that, in many cases, they are often dominated by DM so that the DM mass can be several times greater than the mass of its visible components. Dwarf galaxies can be of various types: elliptical, spheroidal, and irregular. Less common are dwarf spiral galaxies, which exhibit spiral structures.
The essential features of dwarf galaxies, like solitonic cores surrounded by extended halos, in our model are presented by spherically- symmetric primordial dense axion stars (see Figure 3). These stars are formed at the beginning of the radiation-dominated era, at high redshifts, around z 6.76 × 10 6 (≈3.6 h after the Big Bang). In Table 3 and Figure 5, we summarize the results of the numerical solution for DES (I)—dense star with the boundary conditions (I) and DES (II)—dense star with the boundary conditions (II). The first minimum of the axion density in Figure 3 defines the radius of the star core.
Figure 5 presents the dimensionless axion field density of the star core, ρ / ρ , where ρ denotes the mean density of DM. All halos demonstrate a distinct inner core with uniform central distribution of masses. Initially, both stars have close sizes. The ratio of the current sizes of star cores is 13.6 times. This significant difference is explained by different redshifts related to times of star formation: DES (I) (≈3.6 h) and DES(II) (≈5.9 h) after the Big Bang. Figure 6 shows that besides the central core, there are clusters with a significant density that have not merged at the center. It can explain why some dwarf spheroidal galaxies contain globular clusters.
Spiral galaxies. Spiral galaxies, one of the most common galaxies in the universe, have a complex structure consisting of a thin rotating disk and a spherical halo. According to our model, the primordial axion stars that could be seeds of spiral galaxies are the stars with angular momentum l = 1 and m ± 1 (Figure 2c). The current size of the axion star is an order of the visible halo of the spiral galaxy. Formation of the rotating axion stars begins in the middle of the radiation-dominated era ≈ 28,000 years after the Big Bang. The axion stars with l = 1 and m = 0 are possible candidates for the formation of dwarf spiral galaxies. These stars are formed at the beginning of the matter-dominated era 1.22 × 10 5 years after the Big Bang.
Our model predicts the formation of solitonic cores surrounded by extended halos. Another feature is the formation of clusters that do not merge at the center. It can give a potential explanation for why the Fornax dwarf spheroidal galaxy, a satellite of the Milky Way, contains globular clusters [54]. Our findings qualitatively agree with the recent numerical simulations of the ψ D M model presented in Refs. [55,56]. The difference is in the shape of the core profile. In our model, the central distribution of the density is uniform. The smoother density profile in the ψ D M model shows a rather fast decay compared to ours. The possible explanation for this discrepancy is the complicated dynamics of the primordial axion star created in the early universe. Over time, the axion star (a dwarf galaxy in the future) could attract more matter and change its profile to the modern one. One needs a more profound study on this issue.

5. Conclusions

The suggested model assumes that ULA stars, formed in the early universe, could act as seeds for forming larger structures that would become galaxies in the future. Their gravitational field would attract surrounding DM and baryonic stars, eventually creating a galaxy.
Primordial axion stars offer potential insights into the nature of DM and its role in galaxy formation. As ULA stars attract additional axions, they contribute to the overall density of DM halos. This increased density can enhance the halo’s gravitational potential, making it more effective at attracting and retaining baryonic matter. Primordial ULA stars can play a crucial role in forming the observable universe’s shape by acting as seeds for galaxy formation. One needs a more profound study of UlA dynamics to make our model more realistic.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study.

Acknowledgments

We want to express our sincere gratitude to the referees for their valuable comments and suggestions, which significantly improved the quality of our manuscript.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Preskill, J.; Wise, M.B.; Wilczek, F. Cosmology of the invisible axion. Phys. Lett. B 1983, 120, 127–132. [Google Scholar] [CrossRef]
  2. Abbott, L.; Sikivie, P. A cosmological bound on the invisible axion. Phys. Lett. B 1983, 120, 133–136. [Google Scholar] [CrossRef]
  3. Dine, M.; Fischler, W. The not-so-harmless axion. Phys. Lett. B 1983, 120, 137–141. [Google Scholar] [CrossRef]
  4. Kawasaki, M.; Nakayama, K. Axions: Theory and Cosmological Role. Annu. Rev. Nucl. Part. Sci. 2013, 63, 69–95. [Google Scholar] [CrossRef]
  5. Feng, J.L. Dark Matter Candidates from Particle Physics and Methods of Detection. Annu. Rev. Astron. Astrophys. 2010, 48, 495–545. [Google Scholar] [CrossRef]
  6. Bertone, G.; Hooper, D. History of dark matter. Rev. Mod. Phys. 2018, 90, 045002. [Google Scholar] [CrossRef]
  7. O’Hare, C.A.J. Cosmology of axion dark matter. In Proceedings of the 1st General Meeting and 1st Training School of the COST Action COSMIC WISPers—PoS (COSMICWISPers), Bari and Lecce, Italy, 5–14 September 2023; p. 040. [Google Scholar] [CrossRef]
  8. Di Luzio, L.; Giannotti, M.; Nardi, E.; Visinelli, L. The landscape of QCD axion models. Phys. Rep. 2020, 870, 1–117. [Google Scholar] [CrossRef]
  9. Chadha-Day, F.; Ellis, J.; Marsh, D.J.E. Axion dark matter: What is it and why now? Sci. Adv. 2022, 8, eabj3618. [Google Scholar] [CrossRef]
  10. Peccei, R.; Quinn, H. CP Conservation in the Presence of Pseudoparticles. Phys. Rev. Lett. 1977, 38, 1440–1443. [Google Scholar] [CrossRef]
  11. Peccei, R.D.; Quinn, H.R. Constraints imposed by CP conservation in the presence of pseudoparticles. Phys. Rev. D 1977, 16, 1791–1797. [Google Scholar] [CrossRef]
  12. Weinberg, S. A New Light Boson? Phys. Rev. Lett. 1978, 40, 223–226. [Google Scholar] [CrossRef]
  13. Wilczek, F. Problem of Strong P and T Invariance in the Presence of Instantons. Phys. Rev. Lett. 1978, 40, 279–282. [Google Scholar] [CrossRef]
  14. Rosenberg, L.J. Dark-matter QCD-axion searches. Proc. Natl. Acad. Sci. USA 2015, 112, 12278–12281. [Google Scholar] [CrossRef]
  15. Arvanitaki, A.; Huang, J.; Van Tilburg, K. Searching for dilaton dark matter with atomic clocks. Phys. Rev. D 2015, 91, 015015. [Google Scholar] [CrossRef]
  16. Arvanitaki, A.; Dimopoulos, S.; Van Tilburg, K. Sound of Dark Matter: Searching for Light Scalars with Resonant-Mass Detectors. Phys. Rev. Lett. 2016, 116, 031102. [Google Scholar] [CrossRef]
  17. Irastorza, I.G.; Redondo, J. New experimental approaches in the search for axion-like particles. Prog. Part. Nucl. Phys. 2018, 102, 89–159. [Google Scholar] [CrossRef]
  18. Turner, M.S. Windows on the axion. Phys. Rep. 1990, 197, 67–97. [Google Scholar] [CrossRef]
  19. Sikivie, P. Axion Cosmology. In Axions: Theory, Cosmology, and Experimental Searches; Kuster, M., Raffelt, G., Beltrán, B., Eds.; Springer: Berlin/Heidelberg, Germany, 2008; pp. 19–50. [Google Scholar]
  20. Chang, S.; Hagmann, C.; Sikivie, P. Studies of the motion and decay of axion walls bounded by strings. Phys. Rev. D 1998, 59, 023505. [Google Scholar] [CrossRef]
  21. di Cortona, G.G.; Hardy, E.; Vega, J.P.; Villadoro, G. The QCD axion, precisely. J. High Energy Phys. 2016, 2016, 34. [Google Scholar] [CrossRef]
  22. Kim, J.E.; Carosi, G. Axions and the strong CP problem. Rev. Mod. Phys. 2010, 82, 557–601. [Google Scholar] [CrossRef]
  23. Sikivie, P. Invisible axion search methods. Rev. Mod. Phys. 2021, 93, 015004. [Google Scholar] [CrossRef]
  24. Borsanyi, S.; Fodor, Z.; Guenther, J.; Kampert, K.H.; Katz, S.D.; Kawanai, T.; Kovacs, T.G.; Mages, S.W.; Pasztor, A.; Pittler, F.; et al. Calculation of the axion mass based on high-temperature lattice quantum chromodynamics. Nature 2016, 539, 69–71. [Google Scholar] [CrossRef] [PubMed]
  25. Gorghetto, M.; Hardy, E.; Villadoro, G. More axions from strings. SciPost Phys. 2021, 10, 050. [Google Scholar] [CrossRef]
  26. Buschmann, M.; Foster, J.W.; Hook, A.; Peterson, A.; Willcox, D.E.; Zhang, W.; Safdi, B.R. Dark matter from axion strings with adaptive mesh refinement. Nat. Commun. 2022, 13, 1049. [Google Scholar] [CrossRef] [PubMed]
  27. Hindmarsh, M.; Lizarraga, J.; Lopez-Eiguren, A.; Urrestilla, J. Approach to scaling in axion string networks. Phys. Rev. D 2021, 103, 103534. [Google Scholar] [CrossRef]
  28. Yamaguchi, M.; Kawasaki, M.; Yokoyama, J. Evolution of Axionic Strings and Spectrum of Axions Radiated from Them. Phys. Rev. Lett. 1999, 82, 4578–4581. [Google Scholar] [CrossRef]
  29. Braaten, E.; Mohapatra, A.; Zhang, H. Dense Axion Stars. Phys. Rev. Lett. 2016, 117, 121801. [Google Scholar] [CrossRef]
  30. Braaten, E.; Zhang, H. Colloquium: The physics of axion stars. Rev. Mod. Phys. 2019, 91, 041002. [Google Scholar] [CrossRef]
  31. Schiappacasse, E.D.; Hertzberg, M.P. Analysis of dark matter axion clumps with spherical symmetry. J. Cosmol. Astropart. Phys. 2018, 2018, 037. [Google Scholar] [CrossRef]
  32. Visinelli, L.; Baum, S.; Redondo, J.; Freese, K.; Wilczek, F. Dilute and dense axion stars. Phys. Lett. B 2017, 777. [Google Scholar] [CrossRef]
  33. Jetzer, P. Boson stars. Phys. Rep. 1992, 220, 163–227. [Google Scholar] [CrossRef]
  34. Eby, J.; Leembruggen, M.; Street, L.; Suranyi, P.; Wijewardhana, L.C.R. Global view of QCD axion stars. Phys. Rev. D 2019, 100, 063002. [Google Scholar] [CrossRef]
  35. Bautista, B.; Degollado, J.C. Static axion stars revisited. Front. Astron. Space Sci. 2024, 11, 1346820. [Google Scholar] [CrossRef]
  36. Diez-Tejedor, A.; Marsh, D.J. Cosmological production of ultralight dark matter axions. arXiv 2017, arXiv:1702.02116. [Google Scholar] [CrossRef]
  37. Lee, J.-W. Brief History of Ultra-light Scalar Dark Matter Models. EPJ Web Conf. 2018, 168, 06005. [Google Scholar] [CrossRef]
  38. Hui, L.; Ostriker, J.P.; Tremaine, S.; Witten, E. Ultralight scalars as cosmological dark matter. Phys. Rev. D 2017, 95, 043541. [Google Scholar] [CrossRef]
  39. Visinelli, L. Light axion-like dark matter must be present during inflation. Phys. Rev. D 2017, 96, 023013. [Google Scholar] [CrossRef]
  40. Marsh, D.J.E. Axion cosmology. Phys. Rep. 2016, 643, 1–79. [Google Scholar] [CrossRef]
  41. Kim, J.E.; Marsh, D.J.E. An ultralight pseudoscalar boson. Phys. Rev. D 2016, 93, 025027. [Google Scholar] [CrossRef]
  42. Iršič, V.; Viel, M.; Haehnelt, M.G.; Bolton, J.S.; Becker, G.D. First Constraints on Fuzzy Dark Matter from Lyman-α Forest Data and Hydrodynamical Simulations. Phys. Rev. Lett. 2017, 119, 031302. [Google Scholar] [CrossRef]
  43. Armengaud, E.; Palanque-Delabrouille, N.; Yèche, C.; Marsh, D.J.E.; Baur, J. Constraining the mass of light bosonic dark matter using SDSS Lyman-α forest. Mon. Not. R. Astron. Soc. 2017, 471, 4606–4614. [Google Scholar] [CrossRef]
  44. Kobayashi, T.; Murgia, R.; De Simone, A.; Iršič, V.; Viel, M. Lyman-α constraints on ultralight scalar dark matter: Implications for the early and late universe. Phys. Rev. D 2017, 96, 123514. [Google Scholar] [CrossRef]
  45. Wechsler, R.H.; Tinker, J.L. The Connection Between Galaxies and Their Dark Matter Halos. Annu. Rev. Astron. Astrophys. 2018, 56, 435–487. [Google Scholar] [CrossRef]
  46. Frieman, J.A.; Turner, M.S.; Huterer, D. Dark Energy and the Accelerating Universe. Annu. Rev. Astron. Astrophys. 2008, 46, 385–432. [Google Scholar] [CrossRef]
  47. Ryden, B. Introduction to Cosmology; Cambridge University Press: New York, NY, USA, 2017. [Google Scholar]
  48. Braaten, E.; Mohapatra, A.; Zhang, H. Nonrelativistic Effective Field Theory for Axions. Phys. Rev. 2016, D94, 076004. [Google Scholar] [CrossRef]
  49. Davidson, S.; Schwetz, T. Rotating drops of axion dark matter. Phys. Rev. D 2016, 93, 123509. [Google Scholar] [CrossRef]
  50. Eby, J.; Mukaida, K.; Takimoto, M.; Wijewardhana, L.C.R.; Yamada, M. Classical nonrelativistic effective field theory and the role of gravitational interactions. Phys. Rev. D 2019, 99, 123503. [Google Scholar] [CrossRef]
  51. Eby, J.; Suranyi, P.; Vaz, C.; Wijewardhana, L.C.R. Axion Stars in the Infrared Limit. J. High Energy Phys. 2015, 3, 080, Erratum J. High Energy Phys. 2016, 11, 134. [Google Scholar] [CrossRef]
  52. Miniati, F.; Gregori, G.; Reville, B.; Sarkar, S. Axion-Driven Cosmic Magnetogenesis during the QCD Crossover. Phys. Rev. Lett. 2018, 121, 021301. [Google Scholar] [CrossRef]
  53. Planck Collaboration; Aghanim, N.; Akrami, Y.; Arroja, F.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; et al. Planck 2018 results—I. Overview and the cosmological legacy of Planck. Astron. Astrophys. 2020, 641, A1. [Google Scholar] [CrossRef]
  54. Binney, J.; Tremaine, S. Galactic Dynamics; Princeton University Press: Princeton, NJ, USA, 2011. [Google Scholar]
  55. Schive, H.Y.; Chiueh, T.; Broadhurst, T. Cosmic structure as the quantum interference of a coherent dark wave. Nat. Phys. 2014, 10, 496–499. [Google Scholar] [CrossRef]
  56. Schive, H.Y.; Liao, M.H.; Woo, T.P.; Wong, S.K.; Chiueh, T.; Broadhurst, T.; Hwang, W.Y.P. Understanding the Core-Halo Relation of Quantum Wave Dark Matter from 3D Simulations. Phys. Rev. Lett. 2014, 113, 261302. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Effective potential: V 0 = 1 J 0 ( x ) .
Figure 1. Effective potential: V 0 = 1 J 0 ( x ) .
Universe 10 00369 g001
Figure 2. Angular distribution of the axion field density: (a) l = 0 , m = 0 ; (b) l = 1 , m = 0 ; (c) l = 1 , m = ± 1 .
Figure 2. Angular distribution of the axion field density: (a) l = 0 , m = 0 ; (b) l = 1 , m = 0 ; (c) l = 1 , m = ± 1 .
Universe 10 00369 g002
Figure 3. Dense stars. Radial distribution of the axion field density. Inset: log–log plot.
Figure 3. Dense stars. Radial distribution of the axion field density. Inset: log–log plot.
Universe 10 00369 g003
Figure 4. Dilute stars. Radial distribution of the axion field density. Inset: log–log plot.
Figure 4. Dilute stars. Radial distribution of the axion field density. Inset: log–log plot.
Universe 10 00369 g004
Figure 5. Log−log plot. Radial density profiles of halos normalized to the mean DM density. (a) Halo profiles of primordial axion stars at z 6.76 × 10 6 . (b) Current halo profiles. All halos possess a distinct inner core with uniform central distribution of masses.
Figure 5. Log−log plot. Radial density profiles of halos normalized to the mean DM density. (a) Halo profiles of primordial axion stars at z 6.76 × 10 6 . (b) Current halo profiles. All halos possess a distinct inner core with uniform central distribution of masses.
Universe 10 00369 g005
Figure 6. Log plot. Radial density profiles of halos normalized to the mean DM density. Radial density profiles of halos normalized to the mean DM density. Some clusters have not merged at the center.
Figure 6. Log plot. Radial density profiles of halos normalized to the mean DM density. Radial density profiles of halos normalized to the mean DM density. Some clusters have not merged at the center.
Universe 10 00369 g006
Table 1. Acronyms: f a —coupling constant; M —sun mass; n D M —particle DM density.
Table 1. Acronyms: f a —coupling constant; M —sun mass; n D M —particle DM density.
Axion Mass [eV] f a [GeV] r 0 [pc] M 0 [ M ] n 0 [ 1 / cm 3 ] n * = n DM / n 0
10 21 4.8 × 10 16 1.6 6.86 × 10 13 1.5 × 10 44 8.67 × 10 21
5.3 × 10 19 10 16 1.448 × 10 2 3.59 × 10 8 3.45 × 10 45 7.11 × 10 23
Table 2. Acronyms: n h = n b / n 0 —normalized axion density n b at the star’s boundary to n 0 . DES (I)—dense star with the boundary conditions (I); DES (II)—dense star with the boundary conditions (II); DIS (I)—dilute star with the boundary conditions (I); DIS (II)—dilute star with the boundary conditions (II); DIS (III)—dilute star with the boundary conditions (III).
Table 2. Acronyms: n h = n b / n 0 —normalized axion density n b at the star’s boundary to n 0 . DES (I)—dense star with the boundary conditions (I); DES (II)—dense star with the boundary conditions (II); DIS (I)—dilute star with the boundary conditions (I); DIS (II)—dilute star with the boundary conditions (II); DIS (III)—dilute star with the boundary conditions (III).
Star TypeAngular MomentumAxion Mass m a = 5.3 × 10 19 eV
Star Radius, r h Star Mass, M Density, n h Redshift, z
DES (I) l = 0 , m = 0 55 [0.796 pc]43.8 [ 1.57 × 10 10 M ] 5.86 × 10 6 4.35 × 10 5
DES (II) l = 0 , m = 0 2.9 [0.042 pc]5.56 [ 1.99 × 10 9 M ] 0.022 6.76 × 10 6
DIS (I) l = 0 , m = 0 100 [1.448 pc]2.66 [ 9.56 × 10 8 M ] 4.27 × 10 8 8.43 × 10 4
DIS (II) l = 1 , m = 0 410 [5.94 pc]0.11 [ 4 × 10 7 M ] 1.21 × 10 12 2572
DIS (III) l = 1 , m = ± 1 550 [7.96 pc]109.5 [ 3.92 × 10 10 M ] 1.77 × 10 11 6291
Table 3. Acronyms: DES (I)—dense star with the boundary conditions (I); DES (II)—dense star with the boundary conditions (II).
Table 3. Acronyms: DES (I)—dense star with the boundary conditions (I); DES (II)—dense star with the boundary conditions (II).
Star TypeStar Core Radius, r core Star Mass
InitialCurrentCore Mass, M core Total Mass, M total
DES (I) 1.5 × 10 3 pc0.65 kpc 1.4 × 10 6 M 1.57 × 10 10 M
DES (II) 1.3 × 10 3 pc8.81 kpc 1.61 × 10 7 M 1.99 × 10 9 M
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nesterov, A.I. Primordial Axion Stars and Galaxy Halo Formation. Universe 2024, 10, 369. https://doi.org/10.3390/universe10090369

AMA Style

Nesterov AI. Primordial Axion Stars and Galaxy Halo Formation. Universe. 2024; 10(9):369. https://doi.org/10.3390/universe10090369

Chicago/Turabian Style

Nesterov, Alexander I. 2024. "Primordial Axion Stars and Galaxy Halo Formation" Universe 10, no. 9: 369. https://doi.org/10.3390/universe10090369

APA Style

Nesterov, A. I. (2024). Primordial Axion Stars and Galaxy Halo Formation. Universe, 10(9), 369. https://doi.org/10.3390/universe10090369

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop