Cosmological Reflection of Particle Symmetry

The standard model involves particle symmetry and the mechanism of its breaking. Modern cosmology is based on inflationary models with baryosynthesis and dark matter/energy, which involves physics beyond the standard model. Studies of the physical basis of modern cosmology combine direct searches for new physics at accelerators with its indirect non-accelerator probes, in which cosmological consequences of particle models play an important role. The cosmological reflection of particle symmetry and the mechanisms of its breaking are the subject of the present review.


Introduction
The laws of known particle interactions and transformations are based on the gauge symmetry-extension of the gauge principle of quantum electrodynamics to strong and weak interactions.Starting from isotopic invariance of nuclear forces, treating proton and neutron as different states of one particle-nucleon; the development of this approach led to the successful creation of the modern standard model of the elementary particle, involving symmetry between different particles and ascribing their difference to the mechanisms of symmetry breaking.However successful the standard model (SM) is in describing particle properties and interactions, it is not sufficient to provide the basis for modern inflationary cosmology with baryosynthesis and dark matter/energy, and it should be extended to resolve its internal problems, such as divergence of mass of the Higgs boson (which may be resolved by supersymmetric extensions of the SM) or CP violation in Quantum Chromodynamics (QCD) (a popular solution which involves additional Peccei-Quinn symmetry).The aesthetic argument for the extension of the SM comes from the possibility to unify strong and electroweak interactions in the framework of Grand Unified Theories (GUT).The discovery of nonzero mass of neutrino-reflected in the experimentally-detected neutrino oscillations-has already moved physics beyond the SM, since in the SM, neutrinos are massless.
Extensions of the standard model involve new symmetries and new particle states.Noether's theorem relates the exact particle symmetry to conservation of respective charge.If the symmetry is strict, the charge is strictly conserved.The lightest particle bearing this charge is stable.Born in the Universe, they should be present in the form of dark matter, corresponding to ∼25% of the total cosmological density.This form of matter (see, e.g., References [1][2][3][4][5][6] for review and reference) should be stable, explain the measured dark matter density, and decouple from plasma and radiation at least before the beginning of the matter dominated stage.Formation of the large scale structure of the universe from small initial density fluctuations is one of the most important reasons for the nonbaryonic nature of the dark matter that is decoupled from matter and radiation, and provides the effective growth of these fluctuations before recombination.It implies dark matter candidates from the physics beyond the standard model (see References [6][7][8][9][10] for recent review).
On the other hand, the initial density fluctuations coming from the very early universe are also originated from physics beyond the standard model.Mechanisms of symmetry breaking induce new fundamental physical scales in particle theory.If the symmetry is spontaneously broken, it is restored when the temperature exceeds the corresponding scale.Such high temperatures should have naturally arisen at the early stages of cosmological evolution.In the course of cosmological expansion, the temperature decreased, and the transition to the phase with broken symmetry took place, which may be reflected in observable cosmological consequences.
It makes the Big Bang Universe a natural laboratory of particle physics, not only due to possibility of the creation of hypothetical stable particles in the early universe, but also owing to the reflection of the hierarchy of particle symmetry breaking in cosmological phase transitions in their observable effects.
In the old Big Bang scenario, cosmological expansion and its initial conditions were given a priori [11,12].In modern cosmology, expansion of the universe and its initial conditions are related to inflation [13][14][15][16][17], baryosynthesis [18,19], and nonbaryonic dark matter (see review in References [20][21][22][23]).The global properties of the universe, as well as the origin of its large scale structure are considered as the result of the process of inflation.The matter content of the modern universe is also originated from the physical processes: the baryon density is the result of baryosynthesis, and the nonbaryonic dark matter represents the relic species of physics beyond the standard model.
These basic elements of modern cosmology are related with physics beyond the standard model, and in various aspects reflect particle symmetry and mechanisms of its breaking.New symmetries are implied to protect the stability of dark matter candidates.CP violation and baryon charge nonconservation are needed for the generation of baryon excess in a baryon symmetric universe.New fields are needed to drive inflation and play the role of the inflaton.Such fields were associated in the old inflational scenario [14] with the Higgs mechanism of symmetry breaking in Grand Unified Theories (GUT).Inflation finds a relationship with supergravity models in some recent approaches [24].
Whatever the extension of the Standard Model which describes these necessary basic elements of modern cosmology, such extensions inevitably contain some additional model-dependent cosmological consequences, and here we would like to discuss [25] various forms of such additional cosmological reflections of the fundamental particle symmetry.The presented list of nontrivial examples of such reflections-being far from complete-is challenging for the development of astrophysical, astroparticle, cosmological, and collider probes for new physics.
In Section 2 we present examples of the cosmological pattern of particle symmetry and its breaking, from various types of particle dark matter candidates to primordial nonlinear structures.We then consider primordial black holes as universal theoretical probes for cosmological consequences of particle theory (Section 3).We relate the observed broken symmetry of the three known families to various types of dark matter embedded in a unique framework of horizontal unification as well as discuss a possibility for stable charged species of new quarks and leptons to form dark matter, hidden in neutral dark atoms (Section 4).In Section 5 we consider the simplest case of heavy stable −2 charged lepton-like particles surrounded by a helium nuclear shell that form nuclear-interacting O-helium (OHe) dark atoms.The qualitative advantages of using this OHe scenario to explain several puzzles of direct and indirect searches for dark matter challenge its test at the Large Hadron Collider (LHC).The conclusive Section 6 considers cosmological probes of fundamental particle structure in the context of cosmoparticle physics, studying the fundamental relationship of micro-and macro-worlds.

Cosmological Pattern of Particle Physics
In the following, we will specify possible links between fundamental particle structure and its cosmological effects.
Most of the known particles are unstable.For a particle with the mass m, the particle physics time scale is t ∼ 1/m (Henceforth,if it is not otherwise specified, we use the units h = c = k = 1), so in the particle world we refer to particles with lifetime τ 1/m as metastable.To be of cosmological significance in the Big Bang Universe, a metastable particle should survive after t ∼ (m Pl /m 2 ), when the temperature of the universe T fell below T ∼ m, and particles go out of thermal equilibrium.It means that the particle lifetime should exceed t ∼ (m Pl /m) • (1/m), and such a long lifetime should be explained by the existence of an (approximate) symmetry.From this viewpoint, cosmology is sensitive to the conservation laws reflecting strict or nearly strict symmetries of particle theory.
The stability of our ordinary (baryonic) matter gives us an example of this relationship, being protected by the conservation of electric and baryon charges.We assume that an electron is absolutely stable, owing to the conservation of electric charge, while the stability of a proton is conditioned by the conservation of baryon charge.According to the SM, the properties of ordinary matter reflect the fundamental physical scales of electroweak and strong interactions.The mass of an electron is originated from the Higgs mechanism of electroweak symmetry breaking, whereas the mass of a proton reflects the scale of QCD confinement.
New stable particles, corresponding to the new strict symmetry, should be present in the universe and play the role of cosmological dark matter.
However, there is no strict symmetry between various quarks and leptons.Symmetry breaking implies the difference in particle masses.The particle mass spectrum reflects the hierarchy and structure of symmetry breaking.
The mechanism of the spontaneous breaking of particle symmetry also has a cosmological impact.Heating of the condensed matter leads to the restoration of its symmetry.When the heated matter cools down, phase transition to the phase of broken symmetry takes place.In the course of the phase transitions-corresponding to the given type of symmetry breaking-topological defects can form.One can directly observe the formation of such defects in liquid crystals or in superfluid He.In the same manner, the mechanism of the spontaneous breaking of particle symmetry implies the restoration of the underlying symmetry in the early universe at high temperatures.
When temperature decreases in the course of cosmological expansion, transitions to the phase of broken symmetry can lead-depending on the symmetry breaking pattern-to the formation of topological defects in the very early universe.Defects can represent new forms of stable particles (as in the case of magnetic monopoles [26][27][28][29][30][31]), or extended macroscopic structures, such as cosmic strings [32,33] or cosmic walls [34].

Cosmoarcheology of New Physics
Cosmoarcheology considers the results of observational cosmology as the sample of the experimental data on the possible existence and features of hypothetical forms of matter predicted by particle theory.It undertakes a Gedanken Experiment with these forms, assuming some theoretical framework for the origin and evolution of the universe [4].
One can specify new forms of matter by their net contribution to cosmological density and by their possible influence on parameters of matter and radiation.
If new forms of matter lend a dominant contribution to the density of the universe, they determine the dynamics of expansion in that period.The contribution of subdominant new forms of matter to the total density is always small.The simplest examples are dominant and subdominant forms of dark matter.
The effects of new forms of matter can be time dependent, being characterized by their time-scale.Particle decays, evaporation of primordial black holes or the development of gravitational instability are examples of such time dependent effects.
The cosmological structures predicted by particle theory can have inhomogeneous distribution in space.The amplitude of density fluctuations δ ≡ δ / measures the level of inhomogeneity relative to the total density, .The partial amplitude δ i ≡ δ i / i measures the level of fluctuations within a particular component with density i contributing to the total density = ∑ i i .The case δ i ≥ 1 within the considered i-th component corresponds to its strong inhomogeneity.Strong inhomogeneity is compatible with the smallness of total density fluctuations if the contribution of the inhomogeneous component into the total density is small: i , so that δ 1 (see Reference [35] for review).The presence of new forms of matter can influence the properties of matter and radiation either indirectly (i.e., changing the cosmological equation of state) or via direct interaction with matter and radiation.In the first case, only the dominant forms of matter are relevant.In the second case, effects of even subdominant forms of matter are accessible to observational data.The detailed analysis of the sensitivity of cosmological data to various phenomena of new physics are presented in References [1,25].
It should be noted that the parameters of physics beyond the SM can relate the accessible effects to direct experimental research at the accelerators to the cosmological effects of new physics.So, in the model of horizontal unification [36][37][38][39], the top quark (or B-meson) physics fixes the parameters, describing the dark matter, while in supersymmetric models, experimental searches for unstable SUSY (supersymmetric) particles at the LHC put constraints on the parameters of SUSY dark matter candidates [40].

Cosmophenomenology of New Stable Particles
To study the imprints of new physics in astrophysical data, cosmoarcheology implies the forms and means in which new physics leaves such imprints.So, the important tool of cosmoarcheology in linking the cosmological predictions of particle theory to observational data is the Cosmophenomenology of new physics.It studies the possible hypothetical forms of new physics (which may appear as cosmological consequences of particle theory) and their properties, which can result in observable effects.

Freezing out
The simplest primordial form of new physics is the gas of new stable massive particles, originated from the early universe.For particles with mass m, at high temperature T > m, their number density is n ∼ T 3 , and the equilibrium condition is valid at t ∼ mPl/T 2 if their annihilation cross section σ > 1/(Tm Pl ) is sufficiently large to support the equilibrium.At T < m, such particles go out of equilibrium and their relative concentration freezes out.This is the main idea of the calculation of primordial abundance for Weakly Interacting Massive Particles (WIMPs, see References [1,4] for details).
If ordinary particles are among the products of WIMP annihilation, even their small fraction can annihilate in the galaxy, causing significant effects in cosmic rays and gamma background.This effect-first revealed in Reference [41] and then proven for even subdominant fraction of annihilating dark matter in Reference [42]-is now the basis of indirect dark matter searches in cosmic rays [40].
The process of WIMP annihilation to ordinary particles, considered in tchannel, determines their scattering cross section on ordinary particles, and thus relates the primordial abundance of WIMPs to their scattering rate in ordinary matter.Forming a nonluminous massive halo of our galaxy, WIMPs can penetrate terrestrial matter and scatter on nuclei in underground detectors.The strategy of direct WIMP searches implies the detection of recoil nuclei from this scattering.
The process inverse to the annihilation of WIMPs corresponds to their production in the collisions of ordinary particles.It should lead to effects of missing mass and energy-momentum, being the challenge of the experimental search for the production of dark matter candidates at accelerators-e.g., at the LHC.

Stable Relics: Decoupling
More weakly interacting and/or more light species decouple from plasma and radiation, being relativistic at T m, when n i.e., at After decoupling, these species retain their equilibrium distribution until they become non-relativistic at T < m.Conservation of partial entropy in the cosmological expansion links the modern abundance of these species to the number density of relic photons with the account for the increase of the photon number density due to the contribution of heavier ordinary particles, which were in equilibrium in the period of decoupling.
For example, primordial neutrino decouple in the period when relativistic electron-positron plasma was present in equilibrium.The account for the increase of the number density of relic photons due to electron-positron annihilation at T < m e -where m e is the mass of an electron-results in the well-known prediction of Big Bang cosmology [11,12] where n ν ν is the modern number density of one species of primordial left-handed neutrinos (and the corresponding antineutrinos) and n γ = 400 cm −3 is the number density of Cosmological Microwave Background (CMB) photons at the modern CMB temperature T = 2.7 K. Multiplying the predicted modern concentration of neutrinos by their mass, we obtain their contribution to the total density.This contribution should not exceed the total density, which gave the early cosmological upper limits on neutrino mass.For a long time, it seemed possible that relic neutrinos could be the dominant form of cosmological dark matter, and the corresponding neutrino-dominated universe was considered as the physical ground of the Hot Dark Matter scenario of large scale structure formation.Experimental discovery of neutrino oscillations, together with stringent upper limits on the mass of an electron neutrino, exclude this possibility.Moreover, even neutrino masses in the range of 1 eV lead to features in the spectrum of density fluctuations that are excluded by the observational data of CMB.Right-handed neutrinos and left-handed antineutrinos, involved in the seesaw mechanism of neutrino mass generation, are sterile relative to ordinary weak interaction.If these species were in thermal equilibrium in the early universe, they should decouple much earlier than ordinary neutrinos in the period when there were much more particle species (leptons, quarks, gluons, etc.) in the equilibrium, what leads to the primordial abundance of sterile neutrinos much smaller than the ordinary ones.Therefore, cosmological constraints permit sterile neutrinos with mass in the keV range.We refer to the Reference [43] for a recent review of models of sterile neutrinos and their possible effects.

Stable Relics: SuperWIMPs
The maximal temperature which is reached in inflationary universe is the reheating temperature, T r , after inflation.So, the very weakly interacting particles with the annihilation cross section as well as very heavy particles with mass m T r (6) cannot be in thermal equilibrium, and the detailed mechanism of their production should be considered to calculate their primordial abundance.
In particular, the thermal production of the gravitino in the very early universe is proportional to the reheating temperature T r , which puts an upper limit on this temperature from constraints on primordial gravitino abundance [44][45][46][47][48][49][50].

Self-Interacting Dark Matter
Extensive hidden sector of particle theory can provide the existence of new interactions, which only new particles possess.Historically, one of the first examples of such self-interacting dark matter was presented by the model of mirror matter.Mirror particles-first proposed by T. D. Lee and C. N. Yang [51] to restore the equivalence of left-and right-handed co-ordinate systems in the presence of P-and C-violation in weak interactions-should be strictly symmetric by their properties to their ordinary twins.After the discovery of CP-violation, it was shown by I. Yu.Kobzarev, et al. in Reference [52] that mirror partners cannot be associated with antiparticles and should represent a new set of symmetric partners for ordinary quarks and leptons with their own strong, electromagnetic, and weak mirror interactions.It means that there should exist mirror quarks, bound in mirror nucleons by mirror QCD forces and mirror atoms, in which mirror nuclei are bound with mirror electrons by mirror electromagnetic interaction [53,54].If gravity is the only common interaction for ordinary and mirror particles, mirror matter can be present in the universe in the form of elusive mirror objects, having symmetric properties with ordinary astronomical objects (gas, plasma, stars, planets, etc.), but causing only gravitational effects on ordinary matter [55,56].
Even in the absence of any other common interaction except for gravity, the observational data on primordial helium abundance and the upper limits on local dark matter seem to exclude mirror matter evolving in the universe in a fully symmetric way in parallel with the ordinary baryonic matter [57,58].The symmetry in the cosmological evolution of mirror matter can be broken either by initial conditions [59,60] or by breaking mirror symmetry in the sets of particles and their interactions as it takes place in the shadow world [61,62], arising in the heterotic string model.We refer to References [2,63,64] for a current review of mirror matter and its cosmology.
If new particles possess new y-charge, interacting with massless bosons or intermediate bosons with sufficiently small mass (y-interaction), for slow y-charged particles a Coulomb-like factor of "Gamov-Sommerfeld-Sakharov enhancement" [65][66][67] should be added in the annihilation cross section where v is relative velocity and α y is the running gauge constant of the y-interaction.This factor may not be essential in the period of particle freezing out in the early universe (when v was only a few times smaller than c), but can cause strong enhancement in the effect of annihilation of nonrelativistic dark matter particles in the galaxy.

Subdominant Dark Matter
If charge-symmetric stable particles (and their antiparticles) represent the only subdominant fraction of cosmological dark matter, a more detailed analysis of their distribution in space, of their condensation in galaxies, of their capture by stars, Sun, and Earth, as well as the effects of their interaction with matter and of their annihilation provides more sensitive probes for their existence.
In particular, fourth-generation hypothetical stable neutrinos with mass about 50 GeV should be the subdominant form of modern dark matter, contributing less than 0.1% to the total density [41,42].However, direct experimental search for cosmic fluxes of weakly interacting massive particles (WIMPs) may be sensitive to existence of such components (see [68][69][70][71][72][73][74][75][76] and references therein).It was shown in References [77][78][79][80] that annihilation of fourth generation neutrinos and their antineutrinos in the galaxy is severely constrained by the measurements of gamma-background, cosmic positrons, and antiprotons.Fourth generation neutrino annihilation inside the Earth should lead to the flux of underground monochromatic neutrinos of known types, which can be traced in the analysis of the already existing and future data of underground neutrino detectors [79,[81][82][83].

Charged Stable Relics: Dark Atoms
New particles with electric charge and/or strong interaction can form anomalous atoms and be contained in ordinary matter as anomalous isotopes.For example, if the lightest fourth generation quark is stable, it can form stable charged hadrons, serving as nuclei of anomalous atoms of, for example, anomalous helium [84][85][86][87][88][89].Therefore, stringent upper limits on anomalous isotopes, especially on anomalous hydrogen, put severe constraints on the existence of new stable charged particles.However, as we discuss in Section 4, stable doubly charged particles cannot only exist, but even dominate in cosmological dark matter, being effectively hidden in neutral "dark atoms" [90].

Decaying Dark Matter
Decaying particles with lifetime τ, exceeding the age of the universe (t U , τ > t U ) can be treated as stable.By definition, primordial stable particles survive to the present time and should be present in the modern universe.The net effect of their existence is given by their contribution to the total cosmological density.However, even the small effect of their decay can lead to a significant contribution to cosmic rays and gamma background [91].Leptonic decays of dark matter are considered as a possible explanation of the cosmic positron excess, measured in the range above 10 GeV by PAMELA [92], FERMI/LAT [93], and AMS02 [94] (see Reference [95] for a review of the AMS02 experiment).

Charge Asymmetry of Dark Matter
The fact that particles are not absolutely stable means that the corresponding charge is not strictly conserved and the generation of particle charge asymmetry is possible, as is assumed for ordinary baryonic matter.At sufficiently strong particle annihilation cross section, excessive particles (antiparticles) can dominate in the relic density, leaving an exponentially small admixture of their antiparticles (particles) in the same way that primordial excessive baryons dominate over antibaryons in the baryon asymmetric universe.In this case, asymmetric dark matter does not lead to a significant effect of particle annihilation in the modern universe and can be searched for either directly in underground detectors or indirectly by effects of decay or condensation and structural transformations of, e.g., neutron stars (see Reference [96] for recent review and references).If particle annihilation is not strong enough, primordial pairs of particles and antiparticles dominate over excessive particles (or antiparticles), and this case has no principle difference from the charge symmetric case.In particular, for very heavy charged leptons (with mass above 1 TeV), like "tera electrons" [97], discussed in Section 4.2, their annihilation due to electromagnetic interaction is too weak to provide effective suppression of primordial tera electron-positron pairs relative to primordial asymmetric excess [84].

Unstable Particles
Primordial unstable particles with a lifetime less than the age of the universe (τ < t U ) cannot survive to the present time.However, if their lifetime is sufficiently large to satisfy the condition τ (m Pl /m) • (1/m), their existence in the early universe can lead to direct or indirect traces [98].
Weakly interacting particles decaying to invisible modes can influence large scale structure formation.Such decays prevent the formation of structure if they take place before the structure is formed.Invisible products of decays after the structure is formed should contribute to the cosmological dark energy.The unstable dark matter [39,[99][100][101][102][103][104][105][106] implied weakly interacting particles that form structure in the matter-dominated stage and then decay to invisible modes after the structure is formed.
The cosmological flux of decay products contributing to the cosmic and gamma ray backgrounds represents the direct trace of unstable particles [98,107].If the decay products do not survive to the present time, their interaction with matter and radiation can cause indirect trace in the light element abundance [46][47][48]108] or in the fluctuations of thermal radiation [109].
If the particle lifetime is much less than 1 s, the multi-step indirect traces are possible, provided that particles dominate in the universe before their decay.In the dust-like stage of their dominance, black hole formation takes place, and the spectrum of such primordial black holes traces the particle properties (mass, frozen concentration, lifetime) [110][111][112].The particle decay in the end of the dust-like stage influences the baryon asymmetry of the universe.Cosmophenomenological chains link the predicted properties of even unstable new particles to the effects accessible in astronomical observations.Such effects may be important in the analysis of the observational data.

Phase Transitions
Parameters of new stable and metastable particles are also determined by a pattern of particle symmetry breaking.This pattern is reflected in a succession of phase transitions in the early universe.First order phase transitions proceed through bubble nucleation, which can result in black hole formation (see References [113,114] for review and references).Phase transitions of the second order can lead to the formation of topological defects, such as walls, string, or monopoles.The observational data puts severe constraints on magnetic monopole [28] and cosmic wall production [34], as well as on the parameters of cosmic strings [32,33].The structure of cosmological defects can be changed in a succession of phase transitions.More complicated forms, like walls-surrounded-by-strings, can appear.Such structures can be unstable, but their existence can leave a trace in the nonhomogeneous distribution of dark matter and give rise to large scale structures of nonhomogeneous dark matter like archioles [115][116][117].This effect should be taken into account in the analysis of the cosmological effects of weakly interacting slim particles (WISPs) (see Reference [118] for current review) that can play the role of cold dark matter in spite of their small mass.
A wide class of particle models possesses a symmetry breaking pattern, which can be effectively described by pseudo-Nambu-Goldstone (PNG) field and which corresponds to the formation of unstable topological defect structure in the early universe (see Reference [114] for review and references).The Nambu-Goldstone nature in such an effective description reflects the spontaneous breaking of global U(1) symmetry, resulting in continuous degeneracy of vacua.The explicit symmetry breaking at smaller energy scale changes this continuous degeneracy by discrete vacuum degeneracy.The character of formed structures is different for phase transitions, taking place on post-inflationary and inflationary stages.

Large Scale Correlations of Axion Field
At high temperatures, such a symmetry breaking pattern implies the succession of second order phase transitions.In the first transition, continuous degeneracy of vacua leads-at scales exceeding the correlation length-to the formation of topological defects in the form of a string network; in the second phase transition, continuous transitions in space between degenerated vacua form surfaces: domain walls surrounded by strings.This last structure is unstable, but, as was shown in the example of the invisible axion [115][116][117], it is reflected in the large scale inhomogeneity of the distribution of the energy density of coherent PNG (axion) field oscillations.This energy density is proportional to the initial phase value, which acquires a dynamical meaning of the amplitude of the axion field, when axion mass is switched on as a result of the second phase transition.
The value of phase changes by 2π around string.This strong nonhomogeneity of phase leads to corresponding nonhomogeneity of the energy density of coherent PNG (axion) field oscillations.The usual argument (see Reference [119] and references therein) is essential only at scales corresponding to the mean distance between strings.This distance is small, being on the order of the scale of the cosmological horizon in the period when PNG field oscillations start.However, since the nonhomogeneity of phase follows the pattern of the axion string network, this argument misses large scale correlations in the distribution of oscillations' energy density.
Indeed, numerical analysis of the string network (see the review in [120]) indicates that large string loops are strongly suppressed, and the fraction of about 80% of string length (corresponding to long loops) remains virtually the same in all large scales.This property is the other side of the well known scale invariant character of the string network.Therefore, the correlations of energy density should persist on large scales, as was revealed in References [115][116][117].
The large scale correlations in topological defects and their imprints in primordial inhomogeneities is the indirect effect of inflation, if phase transitions take place after the reheating of the universe.Inflation provides, in this case, the equal conditions of phase transition, taking place in causally disconnected regions.

Primordial Seeds for Active Galactic Nuclei
If the phase transitions take place at the inflational stage, new forms of primordial large scale correlations appear.The example of global U(1) symmetry-broken spontaneously in the period of inflation and successively broken explicitly after reheating-was considered in Reference [121].In this model, spontaneous U(1) symmetry breaking at the inflational stage is induced by the vacuum expectation value ψ = f of a complex scalar field Ψ = ψ exp (iθ), also having an explicit symmetry breaking term in its potential V eb = Λ 4 (1 − cos θ).The latter is negligible in the period of inflation if f Λ, so there appears a valley relative to values of phase in the field potential in this period.Fluctuations of the phase θ along this valley-being of the order of ∆θ ∼ H/(2π f ) (here H is the Hubble parameter at the inflational stage)-change in the course of inflation, its initial value within the regions of smaller size.Owing to such fluctuations, for the fixed value of θ 60 in the period of inflation with e-folding N = 60 corresponding to the part of the universe within the modern cosmological horizon, strong deviations from this value appear at smaller scales, corresponding to later periods of inflation with N < 60.If θ 60 < π, the fluctuations can move the value of θ N to θ N > π in some regions of the universe.After reheating, when the universe cools down to temperature T = Λ, the phase transition to the true vacuum states, corresponding to the minima of V eb , takes place.For θ N < π, the minimum of V eb is reached at θ vac = 0; whereas in the regions with θ N > π, the true vacuum state corresponds to θ vac = 2π.For θ 60 < π in the bulk of the volume within the modern cosmological horizon θ vac = 0.However, within this volume there appear regions with θ vac = 2π.These regions are surrounded by massive domain walls, formed at the border between the two vacua.Since regions with θ vac = 2π are confined, the domain walls are closed.After their size equals the horizon, closed walls can collapse into black holes (BHs).
The mass range of formed BHs is constrained by fundamental parameters of the model, f and Λ.The maximal BH mass is determined by the condition that the wall does not dominate locally before it enters the cosmological horizon.Otherwise, local wall dominance leads to a superluminal a ∝ t 2 expansion for the corresponding region, separating it from the other part of the universe.This condition corresponds to the mass [114] The minimal mass follows from the condition that the gravitational radius of BH exceeds the width of wall, and it is equal to [114,122] Closed wall collapse leads to primordial GW spectrum, peaked at with energy density up to At f ∼ 10 14 GeV this primordial gravitational wave background can reach Ω GW ≈ 10 −9 .For the physically reasonable values of 1 < Λ < 10 8 GeV (12) the maximum of the spectrum corresponds to In the range from tens to thousands of Hz, such background may be a challenge for Laser Interferometer Gravitational-Wave Observatory (LIGO) experiment.Another profound signature of the considered scenario are gravitational wave signals from merging of BHs in the primoridal black hole (PBH) cluster.These effects can provide a test of the considered approach in Evolved Laser Interferometer Space Antenna (eLISA) experiment.
This mechanism can lead to the formation of primordial black holes of a whatever large mass (up to the mass of active galactic nuclei (AGNs) [123,124], see for latest review Reference [35]).Such black holes appear in the form of primordial black hole clusters, exhibiting fractal distribution in space [114,122,125].This can shed new light on the problem of galaxy formation [114,124].

Antimatter in Baryon Asymmetric Universe?
Primordial strong inhomogeneities can also appear in the baryon charge distribution.The appearance of antibaryon domains in the baryon asymmetrical universe (reflecting the inhomogeneity of baryosynthesis) is the profound signature of such strong inhomogeneity [126].In an example of a model of spontaneous baryosynthesis (see Reference [127] for review), the possibility of the existence of antimatter domains surviving to the present time in the inflationary universe with inhomogeneous baryosynthesis was revealed in.
The mechanism of spontaneous baryogenesis [127][128][129] implies the existence of a complex scalar field χ = ( f / √ 2) exp (θ) carrying the baryonic charge.The U(1) symmetry-which corresponds to the baryon charge-is broken spontaneously and explicitly.The explicit breakdown of U(1) symmetry is caused by the phase-dependent term The possible baryon and lepton number violating interaction of the field χ with matter fields can have the following structure [127]: where fields Q and L represent a heavy quark and lepton, coupled to the ordinary matter fields.In the early universe, at a time when the friction term induced by the Hubble constant becomes comparable with the angular mass m θ = Λ 2 f , the phase θ starts to oscillate around the minima of the PNG potential and decays into matter fields, according to (15).The coupling (15) gives rise to the following [127]: as the phase starts to roll down in the clockwise direction (Figure 1), it preferentially creates an excess of baryons over antibaryons, while the opposite is true as it starts to roll down in the opposite direction.
The fate of such antimatter regions depends on their size.If the physical size of some of them is larger than the critical surviving size L c = 8h 2 kpc [130], they survive annihilation with surrounding matter.The evolution of sufficiently dense antimatter domains can lead to the formation of antimatter globular clusters [131].The existence of such clusters in the halo of our galaxy should lead to the pollution of the galactic halo by antiprotons.Their annihilation can reproduce [132] the observed galactic gamma background in the range of tens-hundreds MeV.The prediction of an antihelium component of cosmic rays [133]-accessible to future searches for cosmic ray antinuclei in PAMELA and AMS II experiments-as well as of antimatter meteorites [134] provides the direct experimental test for this hypothesis.The possibility of the formation of dense antistars within an extension of the Affleck-Dine scenario of baryogenesis and the strategies for their search were considered in [135].
So, the primordial strong inhomogeneities in the distribution of total dark matter and baryon density in the universe is the new important phenomenon of cosmological models, based on particle models with a hierarchy of symmetry breaking.containing phase θ N gets divided (after one e-fold) into e 3 causally disconnected domains of radius H −1 in f l .Each new domain contains almost homogeneous phase value θ N−1 = θ N ± δθ e f f .This process repeats in every domain with every successive e-fold.

Primordial Black Holes as a Cosmological Reflection of Particle Structure
Any object of mass M can become a black hole, being put within its gravitational radius r g = 2GM/c 2 .At present time, black holes can be created only by a gravitational collapse of compact objects with mass more than about three solar mass [137,138].It can be a natural end of massive stars or can result from the evolution of dense stellar clusters.However, in the early universe, there were no limits on the mass of BH.
Ya. B. Zeldovich and I. D. Novikov (see Reference [139]) noticed that if cosmological expansion stops in some region, a black hole can be formed in this region within the cosmological horizon.This corresponds to strong deviation from general expansion and reflects strong inhomogeneity in the early universe.There are several mechanisms for such strong inhomogeneity and formation of Primordial Black Holes (PBHs) [140,141].
Here we outline the role of PBHs as a link in cosmoarcheological chain, connecting cosmological reflections of particle symmetry with observational data.We discuss the way in which the spectrum of PBHs reflects the properties of superheavy metastable particles and of phase transitions on inflationary and post-inflationary stages.We illustrate in Section 3.1 some mechanisms of PBH formation on the stage of dominance of superheavy particles and fields (Section 3.1.3)and from second order phase transition on the inflationary stage.An effective mechanism of BH formation during bubble nucleation provides a sensitive tool to probe the existence of cosmological first order phase transitions by PBHs (Section 3.3).The existence of stable remnants of PBH evaporation can strongly increase the sensitivity of such a probe, and we demonstrate this possibility in Section 3.4 on an example of gravitino production in PBH evaporation.Being formed within the cosmological horizon, PBHs seem to have masses much less than the mass of stars, constrained by the small size of the horizon in the very early universe.
However, if phase transition takes place in the inflationary stage, closed walls of practically any size can be formed, and their successive collapse can give rise to clouds of massive black holes, which can play the role of seeds for galaxies as discussed above in Section 2.4.2.

PBHs from Early Dust-Like Stages
The possibility of forming a black hole is highly improbable in a homogeneous expanding universe, since it implies metric fluctuations of order 1.For metric fluctuations distributed according to Gaussian law with dispersion a probability for fluctuation of order 1 is determined by an exponentially small tail of the high amplitude part of this distribution.This probability can be further suppressed in the case of non-Gaussian fluctuations [142].
In the universe with equation of state p = γ (17) with numerical factor γ being in the range the probability of forming a black hole from fluctuations within the cosmological horizon is given by [143] This provides the exponential sensitivity of the PBH spectrum to the softening of the equation of state in the early universe (γ → 0) or to the increase of ultraviolet part of the spectrum of density fluctuations ( δ 2 → 1).These phenomena can appear as a cosmological consequence of particle theory.

Dominance of Superheavy Particles in the Early Universe
Superheavy particles cannot be directly studied at accelerators.If they are stable, their existence can be probed by cosmological tests, but there is no direct link between astrophysical data and the existence of superheavy metastable particles with lifetime τ 1s.It was first noticed in Reference [111] that the dominance of such particles in the universe before their decay at t ≤ τ can result in the formation of PBHs, remaining in the universe after the particles decay and keeping some information on particle properties in their spectrum.This provided (though indirect) a possibility to probe the existence of such particles in astrophysical observations.Even the absence of observational evidence for PBHs is important.It puts restrictions on the allowed properties of superheavy metastable particles, which might form such PBHs on a stage of particle dominance, and thus constrains the parameters of models predicting these particles.
After reheating, at particles with mass m and relative abundance r = n/n r (where n is the frozen out concentration of particles and n r is the concentration of relativistic species) must dominate in the universe before their decay.Dominance of these nonrelativistic particles at t > t 0 , where corresponds to a dust-like stage with equation of state p = 0, at which particle density fluctuations grow as and the development of gravitational instability results in the formation of gravitationally-bound systems which decouple at from general cosmological expansion, when δ(t f ) ∼ 1 for fluctuations entering the horizon at t = t i > t 0 with amplitude δ(t i ).
The formation of these systems can result in black hole formation, either immediately after the system decouples from expansion or as a result of the evolution of the initially formed nonrelativistic gravitationally-bound system.

Direct PBH Formation
If density fluctuation is especially homogeneous and isotropic, it directly collapses to BH as soon as the amplitude of fluctuation grows to 1 and the system decouples from expansion.A probability for direct BH formation in the collapse of such homogeneous and isotropic configurations gives minimal estimation of BH formation in the dust-like stage.
This probability was calculated in Reference [111] with the use of the following arguments.In the period t ∼ t f , when fluctuation decouples from expansion, its configuration is defined by averaged density ρ 1 , size r 1 , deviation from sphericity s, and by inhomogeneity u of internal density distribution within the fluctuation.Having decoupled from expansion, the configuration contracts and the minimal size to which it can contract is r min ∼ sr 1 (24) being determined by a deviation from sphericity where γ 1 , γ 2 , and γ 3 define a deformation of configuration along its three main orthogonal axes.It was first noticed in Reference [111] that in order to form a black hole as a result of such a contraction it is sufficient that the configuration returns to the size which had the initial fluctuation δ(t i ) when it entered the horizon at cosmological time t i .If the configuration is sufficiently isotropic to concentrate its mass in the course of collapse within its gravitational radius, but such a concentration also implies sufficient homogeneity of configuration.Density gradients can result in gradients of pressure, which can prevent collapse to BH.This effect does not take place for the contracting collisionless gas of weakly interacting massive particles, but due to the inhomogeneity of collapse, the particles which have already passed the caustics can free stream beyond the gravitational radius before the whole mass is concentrated within it.Collapse of nearly spherically symmetric dust configuration is described by the Tolman solution.Its analysis [110,112,144,145] has provided a constraint on the inhomogeneity u = δρ 1 /ρ 1 within the configuration.It was shown that for both collisionless and interacting particles, the condition is sufficient for the configuration to contract within its gravitational radius.The probability of direct BH formation is then determined by a product of probability for sufficient initial sphericity W s and homogeneity W u of configuration, which is determined by the phase space for such configurations.In a calculation of W s , one should take into account that the condition (27) implies five conditions for independent components of the tensor of deformation before its diagonalization (two conditions for three diagonal components to be close to each other, and three conditions for nondiagonal components to be small).Therefore, the probability of sufficient sphericity is given by [110][111][112]144,145] W s ∼ δ(t i ) 5 (29) and together with the probability for sufficient homogeneity results in the strong power-law suppression of probability for direct BH formation Though this calculation was originally done in References [110][111][112]144,145] for Gaussian distribution of fluctuations, it does not imply a specific form of the high amplitude tail of this distribution, and thus should not change strongly in a case of non-Gaussian fluctuations [142].
The mechanism [1,2,[110][111][112]144,145] is effective for the formation of PBHs with mass in an interval The minimal mass corresponds to the mass within the cosmological horizon in the period t ∼ t 0 , when particles start to dominate in the universe and it is equal to [1,2,[110][111][112]144,145] The maximal mass is indirectly determined by the condition that fluctuation in the considered scale M bhmax , entering the horizon at t(M bhmax ) with an amplitude δ(M bhmax ), can manage to grow up to nonlinear stage, decouple, and collapse before particles decay at t = τ.For scale-invariant spectrum δ(M) = δ 0 , the maximal mass is given by [114] The probability, given by Equation (31), is also appropriate for the formation of PBHs in the dust-like preheating stage after inflation [1,2,146].The simplest example of such a stage can be given by the use of a model of a homogeneous massive scalar field [1,2].Slow rolling of the field in the period t 1/m (where m is the mass of the field) provides a chaotic inflation scenario, while at t > 1/m, the field oscillates with period 1/m.Coherent oscillations of the field correspond to an average over a period of oscillations in a dust-like equation of state p = 0, at which gravitational instability can develop.The minimal mass in this case corresponds to the Jeans mass of scalar field, while the maximal mass is also determined by the condition that fluctuation grows and collapses before the scalar field decays and reheats the universe.
The probability W PBH (M) determines the fraction of total density corresponding to PBHs with mass M. For δ(M) 1, this fraction (given by Equation 31) is small.This means that the bulk of particles do not collapse directly into black holes, but form gravitationally-bound systems.The evolution of these systems can give a much larger amount of PBHs, but it strongly depends on particle properties.

Evolutional Formation of PBHs
Superweakly interacting particles form gravitationally bound systems of collisionless gas, which resemble modern galaxies with collisionless gas of stars.Such a system can finally collapse to a black hole, but energy dissipation within it-and consequently, its evolution-is a relatively slow process [1,2,147].The evolution of these systems is dominantly determined by the evaporation of particles, which gain velocities exceeding the parabolic velocity of the system.In the case of binary collisions, the evolution timescale can be roughly estimated [1,2,147] as for a gravitationally-bound system of N particles, where the free fall time t f f for a system with density . This time scale can be shorter due to collective effects in collisionless gas [148] and at large N can be on the order of However, since the free fall time scale for gravitationally-bound systems of collisionless gas is on the order of cosmological time t f for the period when these systems are formed, even in the latter case the particles should be very long-living (τ t f ) to form black holes in such a slow evolutional process.The evolutional time scale is much smaller for gravitationally-bound systems of superheavy particles interacting with light relativistic particles and radiation.Such systems have analogy with stars, in which evolution time scale is defined by energy loss by radiation.An example of such particles is superheavy color octet fermions of asymptotically free SU(5) model [149] or magnetic monopoles of GUT models.Having decoupled from expansion, frozen out particles and antiparticles can annihilate in gravitationally-bound systems, but detailed numerical simulation [150] has shown that annihilation cannot prevent the collapse of the majority of mass, and the timescale of collapse does not exceed the cosmological time of the period when the systems are formed.

Spikes from Phase Transitions in the Inflationary Stage
Scale non-invariant spectrum of fluctuations-in which the amplitude of small scale fluctuations is enhanced-can be another factor, increasing the probability of PBH formation.The simplest functional form of such a spectrum is represented by a blue spectrum with a power law dispersion with the amplitude of fluctuations growing at k > 0 to small M. The realistic account for the existence of other scalar fields together with inflaton in the period of inflation can give rise to spectra with distinguished scales, determined by the parameters of the considered fields and their interaction.
In the chaotic inflation scenario, interaction of a Higgs field φ with an inflaton η can give rise to phase transitions in the inflationary stage if this interaction induces positive mass term + ν 2 2 η 2 φ 2 .When in the course of slow rolling, the amplitude of an inflaton decreases below a certain critical value η c = m/ν, the mass term in Higgs potential changes sign, and phase transition takes place.Such phase transitions in the inflationary stage lead to the appearance of a characteristic spike in the spectrum of initial density perturbations.These spike-like perturbations, on scales that cross the horizon (60 ≥ N ≥ 1), e-fold before the end of inflation and reenter the horizon during the radiation or dust-like era and could in principle collapse to form primordial black holes.The possibility of such spikes in the chaotic inflation scenario was first pointed out in Reference [151] and realized in Reference [152] as a mechanism of of PBH formation for the model of horizontal unification [36][37][38][39].
For the vacuum expectation value of a Higgs field and λ ∼ 10 −3 , the amplitude δ of a spike in the spectrum of density fluctuations, generated in phase transition in the inflationary stage is given by [152] δ ≈ 4 9s with where κ ∼ 1.
If phase transition takes place at e-folding N before the end of inflation, the spike re-enters the horizon at the radiation dominance (RD) stage and forms a Black hole of mass where H 0 is the Hubble constant in the period of inflation.
If the spike re-enters the horizon in the matter dominance (MD) stage, it should form black holes of mass

First Order Phase Transitions as a Source of Black Holes in the Early Universe
First order phase transitions go through bubble nucleation, recalling the common example of boiling water.The simplest way to describe first order phase transitions with bubble creation in the early universe is based on a scalar field theory with two non-degenerated vacuum states.Being stable at a classical level, the false vacuum state decays due to quantum effects, leading to the nucleation of bubbles of true vacuum and their subsequent expansion [153].The potential energy of the false vacuum is converted into the kinetic energy of bubble walls, thus making them highly relativistic in a short time.The bubble expands until it collides with another one.As it was shown in References [154,155], a black hole may be created in a collision of several bubbles.The probability of the collision of two bubbles is much higher.The opinion of the absence of black holes (BHs) in such processes was based on strict conservation of the original O(2,1) symmetry.As shown in References [113,156,157], there are ways to break it.Firstly, radiation of scalar waves indicates increasing entropy, and hence the permanent breaking of the symmetry during bubble collision.Secondly, the vacuum decay due to thermal fluctuation does not possess this symmetry from the beginning.The investigations [113,156,157] have shown that BH can also be created with a probability of order unity in collisions of only two bubbles.This initiates an enormous production of BH that leads to essential cosmological consequences discussed below.
Inflation models ended by a first order phase transition hold a dignified position in the modern cosmology of the early universe (see for example [158][159][160][161][162][163][164]).The interest in these models is due to the fact that such models are able to generate the observed large-scale voids as remnants of the primordial bubbles for which the characteristic wavelengths are several tens of Mpc [163,164].A detailed analysis of a first order phase transition in the context of extended inflation can be found in Reference [165].Hereafter, we will be interested only in a final stage of inflation when the phase transition is completed.Remember that a first order phase transition is considered completed immediately after establishing the true vacuum percolation regime.Such a regime is established approximately when at least one bubble per unit Hubble volume is nucleated.Accurate computation [165] shows that first order phase transition is successful if the following condition is valid: Here Γ is the bubble nucleation rate.In the framework of first order inflation models, the filling of all space by true vacuum takes place due to bubble collisions, nucleated at the final moment of exponential expansion.The collisions between such bubbles occur when they have comoving spatial dimension less than or equal to the effective Hubble horizon H −1 end at the transition epoch.If we take H 0 = 100 hKm/s/Mpc in Ω = 1 universe, the comoving size of these bubbles is approximately 10 −21 h −1 Mpc.In the standard approach it is believed that such bubbles are rapidly thermalized without leaving a trace in the distribution of matter and radiation.However, in the previous section it was shown that for any realistic parameters of the theory, the collision between only two bubbles leads to BH creation with a probability close to 100% .The mass of this BH is given by [113,156,157] where γ 1 10 −2 and M bub is the mass that could be contained in the bubble volume at the epoch of collision in the condition of the full thermalization of bubbles.The discovered mechanism leads to a new direct possibility of PBH creation at the epoch of reheating in first order inflation models.In the standard picture, PBHs are formed in the early universe if density perturbations are sufficiently large, and the probability of PBH formation from small post-inflation initial perturbations is suppressed (see Section 3.1).A completely different situation takes place at the final epoch of the first order inflation stage; namely, collision between bubbles of Hubble size in the percolation regime leads to copious PBH formation with masses where M hor end is the mass of the Hubble horizon at the end of inflation.According to (47), the initial mass fraction of this PBH is given by β 0 ≈ γ 1 /e ≈ 6 × 10 −3 .For example, for a typical value of H end ≈ 4 × 10 −6 m pl , the initial mass fraction β is contained in PBHs with mass M 0 ≈ 1 g.
In general, the Hawking evaporation of mini BHs [166] could give rise to a variety of possible end states.It is generally assumed that evaporation proceeds until the PBH vanishes completely [167], but there are various arguments against this proposal (see References [168][169][170][171]).If one supposes that BH evaporation leaves a stable relic, then it is natural to assume that it has a mass of order m rel = km pl , where 1 ≤ k ≤ 10 2 .We can investigate the consequences of PBH formation at the percolation epoch after first order inflation, supposing that the stable relic is a result of its evaporation.As it follows from the above consideration, the PBHs are preferentially formed with a typical mass M 0 at a single time t 1 .Hence, the total density ρ at this time is where β 0 denotes the fraction of the total density corresponding to PBHs in the period of their formation t 1 .The evaporation time scale can be written in the following form where g * is the number of effective massless degrees of freedom.Let us derive the density of PBH relics.There are two distinct possibilities to consider.The universe is still radiation dominated (RD) at τ BH .This situation will hold if the following condition is valid: ρ BH (τ BH ) < ρ γ (τ BH ).It is possible to rewrite this condition in terms of Hubble constant at the end of inflation Taking the present radiation density fraction of the universe to be Ω γ 0 = 2.5 × 10 −5 h −2 (h being the Hubble constant in the units of 100 km • s −1 Mpc −1 ), and using the standard values for the present time and time when the density of matter and radiation become equal, we find the contemporary densities fraction of relics It is easily to see that relics overclose the universe (Ω rel >> 1) for any reasonable k and H end > 10 −6 m pl .
The second case takes place if the universe becomes PBHs dominated at period t 1 < t 2 < τ BH .This situation is realized under the condition ρ BH (t 2 ) < ρ γ (t 2 ), which can be rewritten in the form The present day relics density fraction takes the form Thus, the universe is not overclosed by relics, only if the following condition is valid This condition implies that the masses of PBHs created at the end of inflation have to be larger then On the other hand, there are a number of well-known cosmological and astrophysical limits [172][173][174][175][176][177][178] which prohibit the creation of PBHs in the mass range (56) with initial fraction of mass density close to β 0 ≈ 10 −2 .
So, one has to conclude that the effect of the false vacuum bag mechanism of PBH formation makes the coexistence of stable remnants of PBH evaporation with the first order phase transitions at the end of inflation impossible.

PBH Evaporation as Universal Particle Accelerator
Presently, there is no observational evidence proving the existence of PBHs.However, even the absence of PBHs provides a very sensitive theoretical tool to study the physics of the early universe.PBHs represent a nonrelativistic form of matter, and their density decreases with scale factor a as ∝ a −3 ∝ T 3 , while the total density is ∝ a −4 ∝ T 4 in the period of radiation dominance (RD).Being formed within the horizon,a PBH of mass M can be formed not earlier than at If they are formed in the RD stage, the smaller the masses of PBHs, the larger becomes their relative contribution to the total density in the modern MD stage.Therefore, even the modest constraint for PBHs of mass M on their density in units of critical density ρ c = 3H 2 /(8πG) from the condition that their contribution α(M) into the the total density for ρ tot = ρ c does not exceed the density of dark matter converts into a severe constraint on this contribution in the period t f of their formation.If formed in the RD stage at t f = t(M) (given by 57), which corresponds to the temperature T f = m pl m pl /M, PBHs contribute to the total density in the end of the RD stage at t eq , corresponding to T eq ≈ 1 eV, by a factor of a(t eq )/a(t f ) = T f /T eq = m pl /T eq m pl /M larger than in the period of their formation.The constraint on β(M), following from Equation ( 60) is then given by The possibility of PBH evaporation, revealed by S. Hawking [166], strongly influences the effects of PBHs.In the strong gravitational field near the gravitational radius r g of PBH, a quantum effect of the creation of particles with momentum p ∼ 1/r g is possible.Due to this effect, PBHs turn out to be a black body source of particles with temperature (in the units h = c = k = 1) The BH evaporation timescale is τ BH ∼ M 3 /m 4 pl (see Equation (50) and discussion in previous section), and at M ≤ 10 14 g is less than the age of the universe.Such PBHs cannot survive to the present time, and their magnitude Equation (60) should be re-defined and has the meaning of contribution to the total density in the moment of PBH evaporation.For PBHs formed in the RD stage and evaporated in the RD stage at t < t eq , the relationship Equation (62) between β(M) and α(M) is given by [110,179]: The relationship between β(M) and α(M) has a more complicated form if PBHs are formed in early dust-like stages [1,110,144,180], or if such stages take place after PBH formation [1,180].The relative contribution of PBHs to total density does not grow in the dust-like stage, and the relationship between β(M) and α(M) depends on the details of a considered model.Minimal model independent factor α(M)/β(M) follows from the account for enhancement, taking place only during the RD stage between the first second of expansion and the end of the RD stage at t eq , since radiation dominance in this period is supported by observations of light element abundance and the spectrum of CMB [1,110,144,180].
Effects of PBH evaporation make astrophysical data much more sensitive to the existence of PBHs.Constraining the abundance of primordial black holes can lead to invaluable information on cosmological processes, particularly as they are probably the only viable probe for the power spectrum on very small scales which remain far from the Cosmological Microwave Background (CMB) and Large Scale Structures (LSS) sensitivity ranges.To date, only PBHs with initial masses between ∼10 9 g and ∼10 16 g have led to stringent limits (see References [110,168,181,182]) from consideration of the entropy per baryon, the deuterium destruction, the 4 He destruction, and the cosmic-rays currently emitted by the Hawking process [166].The existence of light PBHs should lead to important observable constraints, either through the direct effects of the evaporated particles (for initial masses between 10 14 g and 10 16 g) or through the indirect effects of their interaction with matter and radiation in the early universe (for PBH masses between 10 9 g and 10 14 g).In these constraints, the effects taken into account are those related to known particles.However, since the evaporation products are created by the gravitational field, any quantum with a mass lower than the black hole temperature should be emitted, independently of the strength of its interaction.This could provide a copious production of superweakly interacting particles that cannot be in equilibrium with the hot plasma of the very early universe.This makes evaporating PBHs a unique source of all the species which can exist in the universe.
Following References [1,2,98,180] and [183,184] (but in a different framework and using more stringent constraints), limits on the mass fraction of black holes at the time of their formation (β ≡ ρ PBH /ρ tot ) were derived in Reference [185] using the production of gravitinos during the evaporation process.Depending on whether gravitinos are expected to be stable or metastable, the limits are obtained using the requirement that they do not overclose the universe and that the formation of light nuclei by the interactions of 4 He nuclei with nonequilibrium flux of D,T, 3 He and 4 He does not contradict the observations.This approach is more constraining than the usual study of photo-dissociation induced by photons-photinos pairs emitted by decaying gravitinos.This opened a new window for the upper limits on β below 10 9 g, and correspondingly on various mechanisms of PBH formation [185].Some other aspects of PBH formation and effects are discussed in References [186][187][188].

Symmetry of Known Families
The existence and observed properties of the three known quark-lepton families appeal to the broken SU(3) H family symmetry [36][37][38], which should be involved in the extension of the standard model.This provides the possibility of horizontal unification in the "bottom-up" approach to the unified theory [39].Even in its minimal implementation, the model of horizontal unification can reproduce the main necessary elements of modern cosmology.It provides the physical mechanisms for inflation and baryosynthesis and offers a unified description of candidates for Cold, Warm, Hot, and Unstable Dark Matter.Methods of cosmoparticle physics [1,2] have provided the complete test of this model.Here we discuss the possibilities of linking the physical basis of modern cosmology to the parameters of broken family symmetry.

Horizontal Hierarchy
The approach of References [36][37][38][39] (and its revival in References [189][190][191]) follows the concept of local gauge symmetry SU(3) H , first proposed by Chkareuli [192].Under the action of this symmetry, the left-handed quarks and leptons transform as SU(3) H triplets and the right-handed quarks and leptons as antitriplets.Their mass term transforms as 3 3 = 6 3, and therefore, can only form as a result of horizontal symmetry breaking.
This approach can be trivially extended to the case of n generations, assuming the proper SU(n) symmetry.For three generations, the choice of horizontal symmetry SU(3) H is the only possible choice because the orthogonal and vector-like gauge groups cannot provide different representations for the left-and right-handed fermion states.
In the considered approach, the hypothesis that the structure of the mass matrix is determined by the structure of horizontal symmetry breaking-i.e., the structure of the vacuum expectation values of horizontal scalars carrying the SU(3) H breaking-is justified.
The mass hierarchy between generations is related to the hypothesis of a hierarchy of such symmetry breaking.This hypothesis is called the hypothesis of horizontal hierarchy (HHH) [193][194][195].
The model is based on the gauge SU(3) H flavor symmetry, which is additional to the symmetry of the standard model.It means that there exist eight heavy horizontal gauge bosons, and there are three multiplets of heavy Higgs fields ξ (n) ij (i,j-family indexes, n = 1, 2, 3) in nontrivial (sextet or triplet) representations of SU(3) H .These heavy Higgs bosons are singlets relative to electroweak symmetry, and don't have Yukawa couplings with ordinary light fermions.They have direct coupling to heavy fermions.The latter are singlets relative to electroweak symmetry.Ordinary Higgs φ of the standard model is singlet relative to SU(3) H .It couples left-handed light fermions f i L to their heavy right-handed partners F i R , which are coupled by heavy Higgs multiplets ξ ij with heavy left-handed states F j L .Heavy left-handed states F j L are coupled to right-handed light states f j R by a singlet scalar Higgs field η, which is singlet both relative to SU(3) H and electroweak group of symmetry.The described succession of transitions realizes the Dirac see-saw mechanism, which reproduces the mass matrix m ij of ordinary light quarks and charged leptons f due to mixing with their heavy partners F. This fixes the ratio of vacuum expectation values of heavy Higgs fields, leaving their absolute value as the only main free parameter, which is determined from the analysis of physical, astrophysical, and cosmological consequences.
The SU(3) H flavor symmetry should be chiral in order to eliminate the flavor symmetric mass term.The condition of absence of anomalies implies heavy partners of light neutrinos, and the latter acquire mass by Majorana see-saw mechanism.The natural absence of triple couplings in the heavy Higgs potentials, which do not appear as radiative effects of any other (gauge or Yukawa) interaction, supports additional global U(1) symmetry, which can be associated with Peccei-Quinn symmetry and whose breaking results in the Nambu-Goldstone scalar filed, which shares the properties of axion, Majoron, and singlet familon.

Horizontal Unification
The model provides a complete test (in which its simplest implementation is already ruled out) in a combination of laboratory tests and analysis of cosmological and astrophysical effects.The latter include the study of the effect of radiation of axions on the processes of stellar evolution, the study of the impact of the effects of primordial axion fields and massive unstable neutrinos on the dynamics of the formation of the large-scale structure of the universe, as well as the analysis of the mechanisms of inflation and baryosynthesis based on the physics of the hidden sector of the model.
The model results in physically self-consistent inflationary scenarios with dark matter in the baryon-asymmetric universe.In these scenarios, all steps of cosmological evolution correspond quantitatively to the parameters of particle theory.The physics of the inflaton correspond to the Dirac see-saw mechanism of the generation of the mass of quarks and charged leptons, and the leptogenesis of baryon asymmetry is based on the physics of Majorana neutrino masses.The parameters of axion cold dark matter (CDM), as well as the masses and lifetimes of neutrinos, correspond to the hierarchy of breaking of the SU(3) H symmetry of families.

Stable Charged Constituents of Dark Atoms
New stable particles may possess new U(1) gauge charges and bind by Coulomb-like forces in composite dark matter species.Such dark atoms would look nonluminous, since they radiate invisible light of U(1) photons.Historically, mirror matter (see Section 2.2.4 and References [1,63] for review and references) seems to be the first example of such a nonluminous atomic dark matter.
However, it turned out that the possibility of new stable charged leptons and quarks is not completely excluded, and Glashow's tera-helium [97] has offered a new solution for dark atoms of dark matter.Tera-U-quarks with electric charge +2/3 formed stable (UUU) +2 charged "clusters" that formed with two −1 charged tera-electrons E neutral [(UUU)EE] tera-helium "atoms" that behaved like Weakly Interacting Massive Particles (WIMPs).The main problem for this solution was to suppress the abundance of positively charged species bound with ordinary electrons, which behave as anomalous isotopes of hydrogen or helium.This problem turned out to be unresolvable [84], since the model [97] predicted stable tera-electrons E − with charge −1.As soon as primordial helium is formed in the Standard Big Bang Nucleosynthesis (SBBN), it captures all the free E − in the positively charged (HeE) + ion, preventing any further suppression of positively charged species.Therefore, in order to avoid the overproduction of anomalous isotopes, stable particles with charge −1 (and corresponding antiparticles) should be absent, so that stable negatively charged particles should have charge −2 only.
Since all these models also predict corresponding +2 charge antiparticles, the cosmological scenario should provide a mechanism for their suppression, which can naturally take place in the asymmetric case, corresponding to an excess of −2 charge species, O −− .Then, their positively charged antiparticles can effectively annihilate in the early universe.
If new stable species belong to non-trivial representations of the electroweak SU(2) group, sphaleron transitions at high temperatures can provide the relationship between baryon asymmetry and an excess of −2 charge stable species, as was demonstrated in the case of WTC in References [199,[207][208][209][210][211].

Problem of Tera-Fermion Composite Dark Matter
Glashow's Tera-helium universe was the first inspiring example of the composite dark matter scenario.SU(3) c × SU(2) × SU(2) × U(1) gauge model [97] aimed to explain the origin of the neutrino mass and to solve the problem of strong CP-violation in QCD.New extra SU(2) symmetry acts on three heavy generations of tera-fermions linked with the light fermions by CP transformation.SU(2) symmetry breaking at the TeV scale makes tera-fermions much heavier than their light partners.Tera-fermion mass spectrum is the same as for light generations, but all the masses are scaled by the same factor of about 10 6 .Thus, the masses of the lightest heavy particles are in the tera-eV (TeV) range, explaining their name.
Glashow's model [97] takes into account that very heavy quarks Q (or antiquarks Q) can form bound states with other heavy quarks (or antiquarks) due to their Coulomb-like QCD attraction, and the binding energy of these states substantially exceeds the binding energy of QCD confinement.Then, stable (QQq) and (QQQ) baryons can exist.
According to Reference [97], primordial heavy quark U and heavy electron E are stable and may form a neutral (UUUEE) "atom" with a (UUU) hadron as nucleus and two E − s as "electrons".The gas of such "tera-helium atoms" was proposed in Reference [97] as a candidate for a WIMP-like dark matter.
The problem of such a scenario is the inevitable presence of "products of incomplete combustion", and the necessity to decrease their abundance.
Unfortunately, as was shown in Reference [84], this picture of a Tera-helium universe cannot be realized.
When ordinary 4 He is formed in Big Bang Nucleosynthesis, it binds all the free E − into positively charged ( 4 HeE − ) + "ions".This creates a Coulomb barrier for any successive E − E + annihilation or any effective EU binding.It removes the possibility of suppressing the abundance of unwanted tera-particle species (like (eE + ), ( 4 HeEe), etc.).For instance, the remaining abundance of (eE + ) and ( 4 HeE − e) exceeds the terrestrial upper limit for anomalous hydrogen by 27 orders of magnitude [84].

Composite Dark Matter from Almost Commutative Geometry
The AC-model is based on the specific mathematical approach of unifying general relativity, quantum mechanics, and gauge symmetry [87,198].This realization naturally embeds the standard model, both reproducing its gauge symmetry and Higgs mechanism with the prediction of a Higgs boson mass.The AC model is in some sense alternative to SUSY, GUT, and superstring extensions of the standard model.The AC-model [87] extends the fermion content of the Standard model by two heavy particles-SU(2) electro-weak singlets, with opposite electromagnetic charges.Each of them has its own antiparticle.Having no other gauge charges of the standard model, these particles (AC-fermions) behave as heavy stable leptons with charges −2e and +2e, called A −− and C ++ , respectively.
Similar to the Tera-helium Universe, AC-lepton relics from intermediate stages of a multi-step process towards a final (AC) atom formation must survive in the present universe.In spite of the assumed excess of particles (A −− and C ++ ), the abundance of relic antiparticles ( Ā++ and C−− ) is not negligible.There may also be a significant fraction of A −− and C ++ which remains unbound after the recombination process of these particles into (AC) atoms took place.As soon as 4 He is formed in Big Bang nucleosynthesis, the primordial component of free anion-like AC-leptons (A −− ) is mostly trapped in the first three minutes into a neutral O-helium atom 4 He ++ A −− .O-helium is able to capture free C ++ , creating (AC) atoms and releasing 4 He back.In the same way, the annihilation of antiparticles speeds up.C ++ -O-helium reactions stop when their timescale exceeds a cosmological time, leaving O-helium and C ++ relics in the universe.The catalytic reaction of O-helium with C ++ in the dense matter bodies provides successive (AC) binding that suppresses terrestrial anomalous isotope abundance below the experimental upper limit.Due to the screened charge of AC-atoms, they have WIMP-like interaction with ordinary matter.Such WIMPs are inevitably accompanied by a tiny component of nuclear interacting O-helium.

Stable Charged Techniparticles in Walking Technicolor
The minimal walking technicolor model [200][201][202][203][204][205] has two techniquarks-i.e., up U and down D-that transform under the adjoint representation of a SU(2) technicolor gauge group.The six Goldstone bosons UU, UD, DD, and their corresponding antiparticles carry a technibaryon number, since they are made of two techniquarks or two anti-techniquarks.This means that if there are no processes violating the technibaryon number, the lightest technibaryon will be stable.
The electric charges of UU, UD, and DD are given in general by q + 1, q, and q − 1, respectively, where q is an arbitrary real number.The model additionally requires the existence of a fourth family of leptons-i.e., a "new neutrino" ν and a "new electron" ζ.Their electric charges are in terms of q, respectively, (1 − 3q)/2 and (−1 − 3q)/2.
There are three possibilities for a scenario of dark atoms of dark matter.The first one is to have an excess of Ū Ū (charge −2).The technibaryon number TB is conserved, and therefore UU (or Ū Ū) is stable.The second possibility is to have an excess of ζ that also has −2 charge and is stable, if ζ is lighter than ν and technilepton number L is conserved.In both cases, stable particles with −2 electric charge have substantial relic densities and can capture 4 He ++ nuclei to form a neutral techni-O-helium atom.Finally there is a possibility to have both L and TB conserved.In this case, the dark matter would be composed of bound atoms ( 4 He ++ ζ −− ) and (ζ −− (UU) ++ ).In the latter case, the excess of ζ −− should be larger than the excess of (UU) ++ ), so that WIMP-like (ζ −− (UU) ++ ) is subdominant at the dominance of nuclear interacting techni-O-helium.
The technicolor and the Standard Model particles are in thermal equilibrium as long as the timescale of the weak (and color) interactions is smaller than the cosmological time.The sphalerons allow violation of TB, of baryon number B, of lepton number L, and L , as long as the temperature of the universe exceeds the electroweak scale.It was shown in [199] that there is a balance between the excess of techni(anti)baryons, ( Ū Ū) −− , technileptons ζ −− , or of the both over the corresponding particles (UU and/or ζ ++ ), and the observed baryon asymmetry of the universe.It was also shown the there are parameters of the model at which this asymmetry has proper sign and value, explaining the dark matter density.

Stable Particles of Fourth Generation Matter
Though precision data on the parameters of the standard model did not exclude [212] the existence of the fourth generation of quarks and leptons, the LHC data on the 125 GeV Higgs boson exclude its full strength coupling to fourth family fermions.However, since these fermions should have much larger mass than quarks and leptons of the three known families, it would be reasonable to assume that the main contribution in their mass generation comes from another, more heavy Higgs field.It leads to the suppression of their coupling to a 125 GeV Higgs boson.Then, the limits on the deviation of its production cross-section and probabilities of decay modes from the prediction of the Standard model convert in the constraints on this suppressed coupling [213].The analysis of the application of the LHC results to the direct searches for stable fourth generation quarks and leptons, as well as of the precision SM data for suppressed Higgs couplings need special study.
The fourth generation follows from heterotic string phenomenology, and its difference from the three known light generations can be explained by a new conserved charge, possessed only by its quarks and leptons [67,80,88,196,214].Strict conservation of this charge makes the lightest particle of the fourth family (neutrino) absolutely stable, but it was shown in References [67,80,214] that this neutrino cannot be the dominant form of dark matter.The same conservation law requires the lightest quark to be long living [88,196].In principle, the lifetime of U can exceed the age of the universe if m U < m D [88,196].Provided that sphaleron transitions establish excess of Ū, antiquarks at the observed baryon asymmetry ( Ū Ū Ū) can be formed and bound with 4 He in atom-like state of O-helium [88].
In the successive discussion of OHe dark matter, we generally do not specify the type of −2 charged particle, denoting it as O −− .However, one should note that the AC model does not provide OHe as the dominant form of dark matter, so that the quantitative features of the OHe-dominated universe are not related to this case.

Dark Atoms with Helium Shell
Here we concentrate on the properties of OHe atoms, their interaction with matter, and the qualitative picture of OHe cosmological evolution [87,88,199,209,[215][216][217] and observable effects.We show from the following References [90,218] that the interaction of OHe with nuclei in underground detectors can explain the positive results of dark matter searches in DAMA/NaI (see for review Reference [69]) and DAMA/LIBRA [70] experiments by annual modulations of the radiative capture of O-helium, resolving the controversy between these results and the results of other experimental groups.
In all the considered forms of O-helium, O −− behaves either as lepton or as a specific "heavy quark cluster" with strongly-suppressed hadronic interaction.Therefore, O-helium interaction with matter is determined by nuclear interaction of He.These neutral primordial nuclear interacting species can play the role of a nontrivial form of strongly interacting dark matter [219][220][221][222][223][224][225][226][227], giving rise to a Warmer than Cold dark matter scenario [207,208,215].

OHe Atoms and Their Interaction with Nuclei
The structure of an OHe atom follows from the general analysis of the bound states of O −− with nuclei.
Consider a simple model [228][229][230] in which the nucleus is regarded as a sphere with uniform charge density and in which the mass of the O −− is assumed to be much larger than that of the nucleus.Spin dependence is also not taken into account, so both the particle and nucleus are considered as scalars.Then, the Hamiltonian is given by for short distances r < R and for long distances r > R, where α is the fine structure constant, R = d n A 1/3 ∼ 1.2A 1/3 /(200 MeV) is the nuclear radius, Z is the electric charge of the nucleus, and Z o = 2 is the electric charge of the negatively charged particle O −− .Since Am p M o , the reduced mass is 1/m = 1/(Am p ) + 1/M o ≈ 1/(Am p ).
For small nuclei, the Coulomb binding energy is like in hydrogen atom and is given by For large nuclei, O −− is inside the nuclear radius and the harmonic oscillator approximation is valid for the estimation of the binding energy For the intermediate regions between these two cases with the use of trial function of the form ψ ∼ e −γr/R , variational treatment of the problem [228][229][230] gives where the function F(a) has limits and where a = ZZ o αAm p R. For 0 < a < 1 the Coulomb model gives a good approximation, while at 2 < a < ∞ the harmonic oscillator approximation is appropriate.
In the case of OHe a = ZZ o αAm p R ≤ 1, which proves its Bohr-atom-like structure, assumed in References [88,199,[209][210][211].The radius of Bohr orbit in these "atoms" [88,215] is r o ∼ 1/(Z o Z He αm He ) ≈ 2 × 10 −13 cm.However, the size of the He nucleus rotating around O −− in this Bohr atom turns out to be on the order of and even a bit larger than the radius r o of its Bohr orbit, and the corresponding correction to the binding energy due to non-point-like charge distribution in He is significant.
The Bohr atom-like structure of OHe seems to provide a possibility to use the results of atomic physics for the description of OHe interaction with matter.However, the situation is much more complicated.The OHe atom is similar to hydrogen, in which the electron is hundreds of times heavier than the proton, so that it is a proton shell that surrounds an "electron nucleus".Nuclei that interact with such "hydrogen" would interact first with the strongly interacting "protonic" shell, and such interaction can hardly be treated in the framework of perturbation theory.Moreover, in the description of OHe interaction, accounting for the finite size of He-which is even larger than the radius of the Bohr orbit-is important.One should consider, therefore, the analysis presented below as only a first step approaching the true nuclear physics of OHe.
The approach of References [207,215] assumes the following picture of OHe interaction with nuclei: OHe is a neutral atom in the ground state, perturbed by the Coulomb and nuclear forces of the approaching nucleus.The sign of OHe polarization changes with the distance: at larger distances, Stark-like effect takes place-nuclear Coulomb force polarizes OHe so that the nucleus is attracted by the induced dipole moment of OHe, while as soon as the perturbation by nuclear force starts to dominate, the nucleus polarizes OHe in the opposite way so that He is situated more closely to the nucleus, resulting in the repulsive effect of the helium shell of OHe.When helium is completely merged with the nucleus, the interaction is reduced to the oscillatory potential of O −− with a homogeneously charged merged nucleus with the charge Z + 2.
Therefore OHe-nucleus potential can have a qualitative feature, presented in Figure 2: the potential well at large distances (regions III-IV) is changed by a potential wall in region II.The existence of this potential barrier is crucial for all the qualitative features of the OHe scenario: it causes suppression of reactions with the transition of the OHe-nucleus system to levels in the potential well of the region I, provides the dominance of elastic scattering while transitions to levels in the shallow well (regions III-IV) should dominate in reactions of OHe-nucleus capture.The proof of this picture implies accurate and detailed quantum-mechanical treatment, which was started in Reference [231].With the use of perturbation theory, it was shown that OHe polarization changes sign as the nucleus approaches OHe (as is given in Figure 3), but the perturbation approach was not valid for the description at smaller distances, while the estimates indicated that this change of polarization may not be sufficient for the creation of the potential, given by Figure 2. If the picture of Figure 2 is not proven, one may need more sophisticated models retaining the ideas of the OHe scenario which involve more elements of new physics, as proposed in Reference [232].
On the other hand, O-helium-being an α-particle with screened electric charge-can catalyze nuclear transformations, which can influence primordial light element abundance and cause primordial heavy element formation.It is especially important for the quantitative estimation of the role of OHe in Big Bang Nucleosynthesis and in stellar evolution.These effects need a special detailed and complicated study of OHe nuclear physics, and this work is under way.
The qualitative picture of OHe cosmological evolution is presented below, following References [87,88,90,199,207,209,215,216], and is based on the idea of the dominant role of elastic collisions in OHe interaction with baryonic matter.[231] in the framework of perturbation theory.Note that here R denotes the distance between OHe and the nucleus and not radius of the nucleus, as it was in Equations ( 65), (66), (68), and (69).

Large Scale Structure Formation by OHe Dark Matter
Due to elastic nuclear interactions of its helium constituent with nuclei in the cosmic plasma, the O-helium gas is in thermal equilibrium with plasma and radiation in the Radiation Dominance (RD) stage, while the energy and momentum transfer from plasma is effective.The radiation pressure acting on the plasma is then transferred to density fluctuations of the O-helium gas, and transforms them in acoustic waves at scales up to the size of the horizon.
At temperature T < T od ≈ 1S 2/3 3 eV, the energy and momentum transfer from baryons to O-helium is not effective [88,199] because where m o is the mass of the OHe atom, and S 3 = m o /(1 TeV).Here and v = 2T/m p is the baryon thermal velocity.Then, O-helium gas decouples from plasma.It starts to dominate in the universe after t ∼ 10 12 s at T ≤ T RM ≈ 1 eV, and O-helium "atoms" play the main dynamical role in the development of gravitational instability, triggering large scale structure formation.The composite nature of O-helium determines the specifics of the corresponding dark matter scenario.
At T > T RM , the total mass of the OHe gas with density ρ d = (T RM /T)ρ tot is equal to within the cosmological horizon l h = t.In the period of decoupling T = T od , this mass depends strongly on the O-helium mass S 3 , and is given by [199] M where M is the solar mass.O-helium is formed only at T o and its total mass within the cosmological horizon in the period of its creation is M o = M od (T od /T o ) 3 = 10 37 g.In the RD stage before decoupling, the Jeans length λ J of the OHe gas was restricted from below by the propagation of sound waves in plasma with a relativistic equation of state p = /3, being on the order of the cosmological horizon and equal to λ J = l h / √ 3 = t/ √ 3.After decoupling at T = T od , it falls down to λ J ∼ v o t, where v o = √ 2T od /m o .Though, after decoupling, the Jeans mass in the OHe gas correspondingly falls down one should expect a strong suppression of fluctuations on scales M < M o , as well as adiabatic damping of sound waves in the RD plasma for scales M o < M < M od .It can provide some suppression of small scale structure in the considered model for all reasonable masses of O-helium.The significance of this suppression and its effect on structure formation needs a special study in detailed numerical simulations.In any case, it cannot be as strong as the free streaming suppression in ordinary Warm Dark Matter (WDM) scenarios, but one can expect that, qualitatively, we deal with the Warmer Than Cold Dark Matter model.At temperature T < T od ≈ 1S 2/3 3 keV, the energy and momentum transfer from baryons to O-helium is not effective [88,207,215], and O-helium gas decouples from plasma.It starts to dominate in the universe after t ∼ 10 12 s at T ≤ T RM ≈ 1 eV, and O-helium "atoms" play the main dynamical role in the development of gravitational instability, triggering large scale structure formation.The composite nature of O-helium determines the specifics of the corresponding warmer than cold dark matter scenario.
Being decoupled from baryonic matter, the OHe gas does not follow the formation of baryonic astrophysical objects (stars, planets, molecular clouds, etc.) and forms dark matter halos of galaxies.It can be easily seen that O-helium gas is collisionless for its number density, saturating galactic dark matter.Taking the average density of baryonic matter, one can also find that the galaxy as a whole is transparent for O-helium in spite of its nuclear interaction.Only individual baryonic objects like stars and planets are opaque for it.

Anomalous Component of Cosmic Rays
O-helium atoms can be destroyed in astrophysical processes, giving rise to the acceleration of free O −− in the galaxy.
O-helium can be ionized due to nuclear interaction with cosmic rays [88,211].Estimations [88,233] show that for the number density of cosmic rays n CR = 10 −9 cm −3 during the age of a galaxy, a fraction of about 10 −6 of the total amount of OHe is disrupted irreversibly, since the inverse effect of recombination of free O −− is negligible.Near the solar system, it leads to the concentration of free , and due to effect of Solar modulation these particles initially can hardly reach Earth [208,233].Their acceleration by Fermi mechanism or by collective acceleration forms the power spectrum of the O −− component at the level of O/p ∼ n O /n g = 3 × 10 −10 S −1 3 , where n g ∼ 1 cm −3 is the density of baryonic matter gas.At the stage of red supergiant, stars have the size ∼ 10 15 cm, and during the period of this stage∼ 3 × 10 15 s, up to ∼ 10 −9 S −1 3 of O-helium atoms per nucleon can be captured [208,233].In the Supernova explosion, these OHe atoms are disrupted in collisions with particles in the front of the shock wave and the acceleration of free O −− by the regular mechanism gives the corresponding fraction in cosmic rays.However, this picture needs detailed analysis, based on the development of OHe nuclear physics and numerical studies of OHe evolution in the stellar matter.
If these mechanisms of O −− acceleration are effective, the anomalous low Z/A component of −2 charged O −− can be present in cosmic rays at the level O/p ∼ n O /n g ∼ 10 −9 S −1 3 , and be within the reach for PAMELA and AMS02 cosmic ray experiments.
In the framework of the Walking Technicolor model, the excess of both stable ζ −− and (UU) ++ is possible [208]; the latter being two to three orders of magnitude smaller than the former.This leads to the two-component composite dark matter scenario, with the dominant OHe accompanied by a subdominant WIMP-like component of (ζ −− (UU) ++ ) bound systems.Technibaryons can be metastable, and decays of (UU) ++ can provide an explanation for anomalies observed in the high energy cosmic positron spectrum by PAMELA, FERMI-LAT, and AMS02.

Positron Annihilation and Gamma Lines in Galactic Bulge
Inelastic interaction of O-helium with matter in the interstellar space and its de-excitation can give rise to radiation in the range from a few keV to a few MeV.In the galactic bulge with radius r b ∼ 1 kpc, the number density of O-helium can reach the value n o ≈ 3 × 10 −3 /S 3 cm −3 , and the collision rate of O-helium in this central region was estimated in [211]: At the velocity of v h ∼ 3 × 10 7 cm/ s, energy transfer in such collisions is ∆E ∼ 1 MeVS 3 .These collisions can lead to the excitation of O-helium.If the 2S level is excited, pair production dominates over two-photon channel in the de-excitation by E0 transition, and positron production with the rate 3 × 10 42 S −2 3 s −1 is not accompanied by strong gamma signal.According to Reference [234], this rate of positron production for S 3 ∼ 1 is sufficient to explain the excess in positron annihilation line from bulge, measured by INTEGRAL (see Reference [235] for review and references).If OHe levels with nonzero orbital momentum are excited, gamma lines should be observed from transitions (n > m) E nm = 1.598MeV(1/m 2 − 1/n 2 ) (or from the similar transitions corresponding to the case I o = 1.287MeV) at the level 3 × 10 −4 S −2 3 ( cm 2 s MeVster) −1 .

O-Helium in the Terrestrial Matter
The evident consequence of the O-helium dark matter is its inevitable presence in the terrestrial matter, which appears opaque to O-helium and stores all its in-falling flux.
After they fall down to the terrestrial surface, the in-falling OHe particles are effectively slowed down due to elastic collisions with matter.They then drift, sinking down towards the center of the Earth with velocity Here A med ∼ 30 is the average atomic weight in terrestrial surface matter, n = 2.4 × 10 24 /A is the number of terrestrial atomic nuclei, σv is the rate of nuclear collisions, and g = 980 cm/ s 2 .
Near the Earth's surface, the O-helium abundance is determined by the equilibrium between the in-falling and down-drifting fluxes.
At a depth L below the Earth's surface, the drift timescale is t dr ∼ L/V, where V ∼ 400 S 3 cm/ s is the drift velocity and m o = S 3 TeV is the mass of O-helium.This means that the change of the incoming flux, caused by the motion of the Earth along its orbit, should lead at the depth L ∼ 10 5 cm to the corresponding change in the equilibrium underground concentration of OHe on the timescale t dr ≈ 2.5 × 10 2 S −1 3 s.The equilibrium concentration-which is established in the matter of underground detectors at this timescale-is given by with ω = 2π/T, T = 1yr and t 0 the phase.So, there is an averaged concentration given by n and the annual modulation of concentration characterized by the amplitude Here V h is the speed of Solar System (220 km/s), V E is the speed of the Earth (29.5 km/s), and n 0 = 3 × 10 −4 S −1 3 cm −3 is the local density of O-helium dark matter.

OHe in the Underground Detectors
The explanation [90,215,218] of the results of the DAMA/NaI [69] and DAMA/LIBRA [70] (see Reference [71] for the latest review of these results) experiments is based on the idea that OHe-slowed down in the matter of detector-can form a few keV bound state with a nucleus, in which OHe is situated beyond the nucleus.Therefore, the positive result of these experiments is explained by annual modulation in reaction of radiative capture of OHe A + ( 4 He ++ O −− ) → [A( 4 He ++ O −− )] + γ (81) by nuclei in DAMA detector.
To simplify the solution of the Schrödinger equation, the potential was approximated in References [207,215] by a rectangular potential, presented in Figure 2. The solution of the Schrödinger equation determines the condition under which a low-energy OHe-nucleus bound state appears in the shallow well of region III and the range of nuclear parameters was found at which OHe-sodium binding energy is in the interval 2-4 keV.
The rate of radiative capture of OHe by nuclei can be calculated [215,218] with the use of the analogy with the radiative capture of a neutron by a proton, accounting for: (i) absence of M1 transition that follows from the conservation of orbital momentum; and (ii) suppression of E1 transition in the case of OHe.Since OHe is isoscalar, the isovector E1 transition can take place in the OHe-nucleus system only due to the effect of isospin nonconservation, which can be measured by the factor f = (m n − m p )/m N ≈ 1.4 × 10 −3 , corresponding to the difference of mass of neutron, m n , and proton, m p , relative to the mass of nucleon, m N .In the result, the rate of OHe radiative capture by a nucleus with atomic number A and charge Z to the energy level E in a medium with temperature T is given by Formation of an OHe-nucleus bound system leads to energy release of its binding energy, detected as ionization signal.In the context of our approach, the existence of annual modulations of this signal in the range 2-6 keV and the absence of such an effect at energies above 6 keV means that binding energy E Na of the Na-OHe system in the DAMA experiment should not exceed 6 keV, being in the range 2-4 keV.The amplitude of annual modulation of the ionization signal can reproduce the result of DAMA/NaI and DAMA/LIBRA experiments for E Na = 3 keV.Accounting for the energy resolution in DAMA experiments [239] can explain the observed energy distribution of the signal from a monochromatic photon (with E Na = 3 keV) emitted in OHe radiative capture.
At the corresponding nuclear parameters, there is no binding of OHe with iodine and thallium [215].
It should be noted that the results of the DAMA experiment also exhibit the absence of annual modulations of the energy above 1 MeV.Energy release in this range should take place if the OHe-nucleus system comes to the deep level inside the nucleus.This transition implies tunneling through dipole Coulomb barrier, and is suppressed below the experimental limits.
For the chosen range of nuclear parameters, reproducing the results of DAMA/NaI and DAMA/LIBRA, the results of Reference [215,240] indicate that there are no levels in the OHe-nucleus systems for heavy nuclei.In particular, there are no such levels in Xe, which seems to prevent direct comparison with DAMA results in the XENON100 experiment [75].The existence of such levels in Ge and the comparison with the results of CDMS [72][73][74] and CoGeNT [76] experiments need special study.According to Reference [215], OHe should bind with O and Ca, which is of interest for interpretation of the signal observed in the CRESST-II experiment [241].
In thermal equilibrium, OHe capture rate is proportional to the temperature.Therefore, it looks like it is suppressed in cryogenic detectors by a factor on the order of 10 −4 .However, for the size of cryogenic devices (less than a few tens of meters), OHe gas in them has the thermal velocity of the surrounding matter, and this velocity dominates in the relative velocity of the OHe-nucleus system.This gives the suppression relative to room temperature only ∼ m A /m o .Then, the rate of OHe radiative capture in cryogenic detectors is given by Equation (82), in which room temperature T is multiplied by factor m A /m o .Note that in the case of T = 70 K in the CoGeNT experiment, relative velocity is determined by the thermal velocity of germanium nuclei, which leads to enhancement relative to cryogenic germanium detectors.

Discussion
The existence of heavy stable particles is one of the popular solutions for the dark matter problem.Usually they are considered to be electrically neutral.However, dark matter can potentially be formed by stable heavy charged particles bound in neutral atom-like states by Coulomb attraction.Analysis of the cosmological data and atomic composition of the universe gives constraints on the particle charge, showing that only −2 charged constituents-being trapped by primordial helium in neutral O-helium states-can avoid the problem of overproduction of the anomalous isotopes of chemical elements, which are severely constrained by observations.A cosmological model of O-helium dark matter can even explain puzzles of direct dark matter searches.
The proposed explanation is based on the mechanism of low-energy binding of OHe with nuclei.Within the uncertainty of nuclear physics parameters, there exists a range at which OHe binding energy with sodium is in the interval 2-4 keV.Annual modulation in radiative capture of OHe to this bound state leads to the corresponding energy release observed as an ionization signal in DAMA/NaI and DAMA/LIBRA experiments.
Accounting for the of the numerical results to the values of nuclear parameters and for the approximations made in the calculations, the presented results can be considered only as an illustration of the possibility of explaining puzzles of the dark matter search in the framework of the composite dark matter scenario.An interesting feature of this explanation is a conclusion that the ionization signal may be absent in detectors containing light (e.g., 3 He) or heavy (e.g., Xe) elements.Therefore, tests of results of the DAMA/NaI and DAMA/LIBRA experiments by other experimental groups can become a very nontrivial task.Recent indications of a positive result in the matter of the CRESST detector [241], in which OHe binding is expected together with absence of signal in xenon detector [75], may qualitatively favor the presented approach.For the same chemical content, an order of magnitude suppression in cryogenic detectors can explain why a possible positive effect in the CoGeNT experiment [76] can be compatible with the constraints of the CDMS experiment.
The present explanation contains distinct features, by which it can be distinguished from other recent approaches to this problem (see [25] for review and references).
An inevitable consequence of the proposed explanation is the appearance of anomalous superheavy isotopes in the matter of underground detectors, having mass roughly m o larger than ordinary isotopes of the corresponding elements.
It is interesting to note that in the framework of the presented approach, positive results of the experimental search for WIMPs by effect of their nuclear recoil would be a signature for a multicomponent nature of dark matter.Such OHe+WIMPs multicomponent dark matter scenarios naturally follow from the AC model [87] and can be realized in models of Walking technicolor [208].
Stable −2 charge states (O −− ) can be elementary, like AC-leptons or technileptons, or look like technibaryons.The latter, composed of techniquarks, reveal their structure at much higher energy scale and should be produced at the LHC as elementary species.The signature for AC leptons and techniparticles is unique and distinctive, which allows their separation from other hypothetical exotic particles.
Since simultaneous production of three U Ū pairs and their conversion in two doubly charged quark clusters UUU is suppressed, the only possibility to test the models of composite dark matter from fourth generation in the collider experiments is a search for the production of stable hadrons containing single U or Ū-like Uud and Ūu/ Ūd.
The presented approach sheds new light on the physical nature of dark matter.Specific properties of dark atoms and their constituents are challenging for the experimental search.The development of the quantitative description of OHe interaction with matter confronted with the experimental data will provide the complete test of the composite dark matter model.It challenges the search for stable doubly charged particles at accelerators and cosmic rays as direct experimental probe for charged constituents of the dark atoms of dark matter.

Conclusions
The mystical Ouroboros (self-eating-snake) illustrates the main problem of modern fundamental physics: The theory of the universe is based on the predictions of particle theory, which in turn need cosmology for their test.Indeed, our modern understanding of the structure and evolution of the universe involves phenomena of inflation, baryosynthesis, and dark matter and dark energy.Cosmoparticle physics [1][2][3]242,243] offers the way out of this wrong circle.It studies the fundamental basis and mutual relationship between micro-and macro-worlds in the proper combination of physical, astrophysical, and cosmological signatures.Some aspects of this relationship, which arise in the extension of nuclear symmetry to physics beyond the standard model, is the subject of this review.

Figure 1 .
Figure1.The inflational evolution of the phase (taken from Reference[136]).The phase θ 60 sits in the range [π, 0] at the beginning of inflation and makes Brownian step δθ e f f = H in f l /(2π f e f f ) at each e-fold.The typical wavelength of the fluctuation δθ is equal to H −1 in f l .The whole domain H −1 in f l

Figure 2 .Figure 3 .
Figure 2. The potential of the O-Helium (OHe)-nucleus system and its rectangular well approximation.